采用ANSYS进行疲劳分析的apdl.docx
《采用ANSYS进行疲劳分析的apdl.docx》由会员分享,可在线阅读,更多相关《采用ANSYS进行疲劳分析的apdl.docx(3页珍藏版)》请在冰豆网上搜索。
ANSYS疲劳分析的apdl
/units,si
/title,Fatigueanalysisofcylinderwithflathead
!
***************参数设定***************Di=1000!
筒体内径t=20!
筒体厚度
hc=nint(4*sqrt(Di/2*t)/10)*10!
模型中筒体长度
tp=60!
平板封头厚度
r1=10!
平板封头外测过渡圆弧半径
r2=10!
平板封头内侧应力释放槽圆弧半径exx=2e5!
材料弹性模量
mu=0.3!
材料泊松比
p1=2!
最高工作压力p3=2.88!
水压试验压力
n1=2e4!
最高/最低压力循环次数n2=5!
水压试验次数
!
***************前处理***************
/prep7
et,1,82!
设定单元类型keyopt,1,3,1!
设定周对称选项mp,ex,1,exx!
定义材料弹性模量mp,nuxy,1,mu!
定义材料泊松比
!
*******建立模型*******
k,1,0,0!
定义关键点k,2,Di/2+t,,
k,3,Di/2+t,-(tp+hc)k,4,Di/2,-(tp+hc)k,5,Di/2,-tp
k,6,Di/2-r2,-tp!
定义应力释放槽圆弧中心关键点k,7,0,-tp
l,1,2!
生成线l,2,3
l,3,4
l,4,5
l,5,7
l,7,1
LFILLT,1,2,r1!
生成外测过渡圆弧al,all!
生成子午面
CYL4,kx(6),ky(6),r2,180!
生成应力释放槽面域ASBA,1,2!
面相减
wprot,,,90!
旋转工作平面wpoff,,,kx(6)-3*r2!
移动工作平面asbw,all!
用工作平面切割子午面wprot,,90!
旋转工作平面wpoff,,,tp+r2!
移动工作平面asbw,all!
用工作平面切割子午面
esize,5!
设定单元尺寸MSHKEY,1!
设定映射剖分amesh,1!
映射剖分面1amesh,3!
映射剖分面3esize,2!
设定单元尺寸MSHKEY,0!
设定自由剖分amesh,4!
自由剖分面4fini!
退出前处理
!
***************求解***************
/solu!
筒体端部施加轴向约束dl,3,,uy!
筒体端部施加轴向约束dl,6,,symm!
平板封头对称面施加对称约束
time,1!
载荷步1lsel,s,,,8!
选择内表面各线段lsel,a,,,11,13
lsel,a,,,15
cm,lcom1,line!
生成内表面线组件SFL,all,PRES,p1,!
内表面施加内压alls!
全选
solve!
求解
fini!
退出求解器
!
***************后处理***************
/post1!
进入后处理
FTSIZE,1,2,2,!
设定疲劳评定的位置数、事件数及载荷数
FP,1,1e1,2e1,5e1,1e2,2e2,5e2!
根据疲劳曲线输入S-N数据
FP,7,1e3,2e3,5e3,1e4,2e4,5e4FP,13,1e5,2e5,5e5,1e6,,FP,19,,
FP,21,4000,2828,1897,1414,1069,724FP,27,572,441,331,262,214,159FP,33,138,114,93.1,86.2,,
FP,39,,
!
******水压试验循环******
fs,4760,1,1,1,0,0,0,0,0,0!
储存节点4760对应其第一载荷的应力set,1,last!
读入第一载荷步数据FSNODE,4760,1,2!
储存节点4760对应其第二载荷的应力fe,1,n2,p3/p1!
设定事件循环次数及载荷比例系数
!
******最高/最低压力循环******
fs,4760,2,1,1,0,0,0,0,0,0!
储存节点4760对应其第一载荷的应力set,1,last!
读入第一载荷步数据FSNODE,4760,2,2!
储存节点4760对应其第二载荷的应力FE,2,n1,1,!
设定事件循环次数及载荷比例系数FTCALC,1!
进行疲劳计算(并记录使用系数)fini