ex4.15(3).inp

ansys有限元分析程序

! 4.15  联轴器膜片多工况分析
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
! 轴向变形引起的弯曲应力计算
FINISH					! 结束以前的命令
/CLEAR					! 清空数据库
/FILNAME, EX4.15(3)	! 定义工程文件名称
/TITLE, EX4.15(3) Membrance of Axial Connector under Axial Displacement.
! （1）设置几何参数
WAIJ=139					! 圆环的外半径D1/2=139mm
NEIJ=95					! 圆环的内半径D3/2=95mm
XIAOJ=11					! 小圆孔的半径d/2=11mm
ZHONGJ=117					! 小圆孔圆心到圆心的距离，中半径D2/2=117mm
THICK=0.4					! 圆环的厚度Thick=0.4mm
Pi=3.1415926				! 圆周率常数
! XIAOJ1=(ZHONGJ-NEIJ-XIAOJ)/2+XIAOJ
! zhouc=Pi/2*XIAOJ		! 1/4小圆孔的弧长
! power=400000e3			! 功率400kW=400E6 N*mm/s
! NR=2960/60*2*Pi			! 用每分钟转数计算转动角速度
! force=power/NR/ZHONGJ/80	! 20层膜片每个螺栓上的受力
! （2）根据对称性，建立实体模型
/PREP7						! 加载前处理模块
*AFUN, DEG					! 声明在三角函数中变量单位为度
PCIRC, NEIJ, WAIJ, 0, 90		! 用内径NEIJ, 外径WAIJ, 从0度到90度，建立1/4圆环
CYL4, ZHONGJ*COS(45), ZHONGJ*SIN(45), , , XIAOJ	! 建立45度方向半径为XIAOJ小圆
CYL4, 0, ZHONGJ, , -90, XIAOJ, 90		! 在(0, ZHONGJ)位置建立半径为XIAOJ的小圆
CYL4, ZHONGJ, 0, , 0, XIAOJ, 180			! 在(ZHONGJ, 0)位置建立半径为XIAOJ的小圆
AOverlap, ALL			! 利用重叠操作生成多个图元
ADele, 2, 4, 1, 1		! 删除面生成三个小圆孔
! （3）选取单元，设置材料属性和实常数
ET, 1, SHELL63			! 定义第1类单元为弹性壳单元SHELL63
MP, EX, 1, 198e3			! 定义第1类材料的弹性模量EX=198E3 MPa
MP, PRXY, 1, 0.3		! 定义第1类材料的泊松比PRXY=0.3
MP, DENS, 1, 7.9e-9	! 定义第1类材料的密度DENS=1.7e-9
R, 1, THICK				! 定义第1类实常数，膜片厚度
! （4）施加右边和上边小孔的固定条件和右上小孔里的压力荷载
DL, 9, , ALL, 0		! 上面小孔右半圆位移固定
DL, 12, , ALL, 0	! 右面小孔上半圆位移固定
DL, 15, , UY, 0		! 右面小孔右侧线段位移竖向固定
DL, 18, , UY, 0		! 右面小孔左侧线段位移竖向固定
DL, 16, , UX, 0		! 上面小孔上侧线段位移水平固定
DL, 17, , UX, 0		! 上面小孔下侧线段位移水平固定
! （5）给右上方小孔施加轴向位移
LSel, S, Line, , 5, 8, 1	! 选择右上方小孔周边的5号到8号线
DL, ALL, , UZ, THICK*7/2	! 对选择集里的这4条线施加1.4mm轴向位移
LSel, ALL		! 选择所有的线
! （6）按照指定单元的大小，将面剖分为面单元
AESize, ALL, 4		! 定义面单元尺寸4mm
AMESH, ALL			! 剖分所有面
Finish				! 退出前处理模块
! （7）进入求解模块求解
/solu
Solve		! 开始求解
Finish		! 退出求解模块
! （8）进入一般后处理模块显示结果
/post1				! 进入一般后处理模块
PLNSol, S, eqv		! 显示节点的Mises应力分布
PLNSol, S, 1			! 显示节点的最大拉应力分布
PLNSol, S, 3			! 显示节点的最大压应力分布
FINISH