如何给刚体添加旋转运动（有proc文件，请帮忙分析）

想做个例子学习其中刚体运动控制的方法。模型中建立了一个平面，采用壳单元划分网格 。一个圆柱面，定义为刚体，想控制圆柱绕自身的中轴旋转。基本做法是定义刚体运动控制点到圆柱中心，然后添加绕z轴的角度位移。结果做出来跟想象差别很大。下面是proc文件，请各位高手给分析一下哪里出了问题。多谢！
*add_points
0 0 0
100 0 0
100 0 100
0 0 100
*set_surface_type quad
*add_surfaces
1
4
3
2
*convert_surfaces
all_existing
# | End of List
*set_surface_type cylinder
*add_surfaces
50 50 0
50 50 100
50
50
*flip_surfaces
all_existing
*geometry_type mech_three_shell
*geometry_value thick
1
*add_geometry_elements
all_existing
*material_type mechanical:isotropic
*material_value isotropic:youngs_modulus
540
*material_value isotropic:poissons_ratio
0.32
*add_material_elements
all_existing
*new_contact_body
*contact_deformable
*add_contact_body_elements
all_existing
*new_contact_body
*add_contact_body_surfaces
2
# | End of List
*cbody_center_rot
50 50 50
*contact_rigid
*contact_option control:position
*contact_value az
1
*apply_type fixed_displacement
*apply_dof x *apply_dof_value x
*apply_dof y *apply_dof_value y
*apply_dof z *apply_dof_value z
*apply_dof rx *apply_dof_value rx
*apply_dof ry *apply_dof_value ry
*apply_dof rz *apply_dof_value rz
*add_apply_nodes
113 114 115 116 117 118 119 120 121
111 112
# | End of List
*add_apply_nodes
1 2 3 4 5 6 7 8 9 10 11
# | End of List
*new_loadcase
*loadcase_type static
*job_class mechanical
*add_loadcase_loads apply1
*job_option large:on
*job_option elasticity:l_strn_upd
*job_class mechanical
*add_job_loadcases lcase1
*update_job
*submit_job 1 *monitor_job