lsel,u,loc,z,0 ! 去掉筒体线段
LESIZE,all, , ,nme, , , , ,1 ! 设定椭圆封头经线剖分数
vsweep,all ! 扫略剖分x轴正向的部分封头及部分中心接管
alls ! 全选
lsel,s,loc,y,180 ! 选择180度位置线 lsel,r,loc,x,dn3/2,Di ! 再选中心接管以外线 lsel,u,loc,x,dn3/2 ! 去掉接管线段 lsel,u,loc,z,0,hc ! 去掉筒体线段
LESIZE,all, , ,nme, , , , ,1
vsel,s,loc,x,0,dn3/2 csys,0 vsel,u,loc,x,0,Di vsweep,all
lsel,s,radius,,dn1/2-tn1,dn1/2 LESIZE,all, , ,nh, , , , ,1 lsel,r,radius,,dn1/2-tn1 asll vsla aslv vsweep,all cm,vcon1,volum
lsel,s,radius,,dn2/2-tn2,dn2/2 LESIZE,all, , ,nh, , , , ,1 lsel,r,radius,,dn2/2-tn2 asll vsla aslv vsweep,all cm,vcon2,volum
lsel,s,radius,,dn1/2 asll vsla cmsel,u,vcon1 vsweep,all
lsel,s,radius,,dn2/2 asll vsla cmsel,u,vcon2 ! 设定椭圆封头经线剖分数 ! 选择中心接管及其与封头相关区 激活总体直角坐标系 ! 去掉x轴正向的部分
! 扫略剖分接管及其与封头相关区的x轴负向部分 ! 选择半径为接管1半径的线段 设定接管1环向剖分数
! 选择半径为接管1内半径的线段 ! 按线选面 ! 按面选体 ! 按体选面 ! 剖分接管1
! 定义接管1为体组件vcon1 ! 选择半径为接管2半径的线段 设定接管2环向剖分数
! 选择半径为接管2内半径的线段 ! 按线选面 ! 按面选体 ! 按体选面 ! 剖分接管2
! 定义接管2为体组件vcon2 ! 选择半径为接管1外半径的线段 ! 按线选面 ! 按面选体
! 去掉体组件vcon1
! 剖分接管1相关封头及相贯区 ! 选择半径为接管2外半径的线段 ! 按线选面 ! 按面选体
! 去掉体组件vcon2
! ! ! vsweep,all ! 剖分接管2相关封头及相贯区 alls ! 全选
fini ! 退出前处理 /solu ! 进入求解器
csys,0 ! 激活总体直角坐标系 asel,s,loc,y,-hc ! 选择筒体端面 da,all,uy ! 约束轴向位移 nsla,,1 ! 按面选节点
nsel,r,loc,x,0 ! 再选择x为0的节点 d,all,ux ! 约束x方向位移 nsla,,1 ! 按面选节点
nsel,r,loc,z,0 ! 再选择z为0的节点 d,all,uz ! 约束z方向位移
! ********* 以下选择内表面并设定为面组件 ************* alls
lsel,s,radius,,Di/2 ! 选择半径为筒体内半径的线段 lsel,u,loc,y,-hc ! 去掉筒体端部线段 asll ! 按线选面
asel,u,loc,y,0 ! 去掉厚度方向的面 cm,acon,area ! 定义面组件acon
lsel,s,radius,,dn3/2-tn3 ! 选择半径为接管3内半径的线段 asll ! 按线选面
asel,u,loc,y,ky(9) ! 去掉厚度方向的面 asel,u,loc,y,ky(13) ! 去掉厚度方向的面 cmsel,a,acon ! 添加面组件acon cm,acon,area ! 再定义面组件acon
lsel,s,radius,,dn1/2-tn1 ! 选择半径为接管1内半径的线段 lsel,r,loc,y,ky(49),ky(50) ! 再选择接管1端面线段 asll ! 按线选面
cm,acon0,area ! 定义面组件acon0 ksel,s,,,49 ! 选择关键点 lslk ! 按关键点选线
lsel,r,length,,tn1 ! 再选长度为接管1厚度的线段 asll ! 按线选面
cm,acon1,area ! 定义面组件acon1 cmsel,s,acon0 ! 选择面组件acon0
cmsel,u,acon1 ! 从面组件acon0中去掉面组件acon1 cmsel,a,acon ! 添加面组件acon cm,acon,area ! 再定义面组件acon
lsel,s,radius,,dn2/2-tn2 ! 选择半径为接管2内半径的线段 lsel,r,loc,y,ky(48),ky(45) ! 再选择接管2端面线段
asll ! 按线选面
cm,acon0,area ! 定义面组件acon0 ksel,s,,,46 ! 选择关键点 lslk ! 按关键点选线
lsel,r,length,,tn2 ! 再选长度为接管1厚度的线段 asll ! 按线选面
cm,acon1,area ! 定义面组件acon1 cmsel,s,acon0 ! 选择面组件acon0
cmsel,u,acon1 ! 从面组件acon0中去掉面组件acon1 cmsel,a,acon ! 添加面组件acon cm,acon,area ! 再定义面组件acon
! ********* 以下选择接管端面并设定为面组件 *************
lsel,s,radius,,dn1/2-tn1 ! 选择半径为接管1内半径的线段 asll ! 按线选面
asel,r,loc,y,ky(59),ky(58) ! 再选择接管1端面
cm,apax1,area ! 定义接管1端面面组件apax1
lsel,s,radius,,dn2/2-tn2 ! 选择半径为接管2内半径的线段 asll ! 按线选面
asel,r,loc,y,ky(74),ky(75) ! 再选择接管2端面
cm,apax2,area ! 定义接管2端面面组件apax2
lsel,s,radius,,dn3/2-tn3 ! 选择半径为接管2内半径的线段 asll ! 按线选面
asel,r,loc,y,ky(9) ! 再选择接管3端面
cm,apax3,area ! 定义接管3端面面组件apax3 ! ********* 以下选择外表面并设定为面组件 *************
lsel,s,radius,,Di/2+tc ! 选择半径为筒体外半径的线段 lsel,u,loc,y,-hc ! 去掉筒体端部线段 asll ! 按线选面
asel,u,loc,y,0 ! 去掉厚度方向的面 cm,aouter,area ! 定义面组件aouter
lsel,s,radius,,dn3/2 ! 选择半径为接管3外半径的线段 lsel,r,loc,y,ky(9) ! 再选接管3端部线段 asll ! 按线选面
asel,u,loc,y,ky(9) ! 去掉厚度方向的面 cmsel,a,aouter ! 添加面组件aouter cm,aouter,area ! 再定义面组件aouter
lsel,s,radius,,dn1/2 ! 选择半径为接管1外半径的线段 lsel,r,loc,y,ky(55),ky(53) ! 再选接管1端部线段 asll ! 按线选面
asel,u,loc,y,ky(55),ky(53) ! 去掉厚度方向的面
cmsel,a,aouter ! 添加面组件aouter cm,aouter,area ! 再定义面组件aouter
lsel,s,radius,,dn2/2 ! 选择半径为接管2外半径的线段 lsel,r,loc,y,ky(69),ky(71) ! 再选接管2端部线段 asll ! 按线选面
asel,u,loc,y,ky(69),ky(71) ! 去掉厚度方向的面 cmsel,a,aouter ! 添加面组件aouter cm,aouter,area ! 再定义面组件aouter
time,1 cmsel,s,apax1 sfa,all,1,pres,pax21 cmsel,s,apax2 sfa,all,1,pres,pax22 cmsel,s,apax3 sfa,all,1,pres,pax23 cmsel,s,acon sfa,all,1,pres,p2 alls solve save
time,2 alls SFADELE,all,1,pres cmsel,s,apax1 sfa,all,1,pres,pax11 cmsel,s,apax2 sfa,all,1,pres,pax12 cmsel,s,apax3 sfa,all,1,pres,pax13 cmsel,s,aouter sfa,all,1,pres,p1 alls solve save fini ! ****** 后处理 ***********
/post1 LCDEF,1,1 LCDEF,2,2 LCASE,1 LCOPER,sub,2 LCWRITE,3 第一载荷步对应最大工作压力 ! 选择面组件apax1
! 对接管1施加端部平衡面载荷 ! 选择面组件apax2
! 对接管2施加端部平衡面载荷 ! 选择面组件apax3
! 对接管3施加端部平衡面载荷 ! 选择面组件acon ! 施加内压 ! 全选 ! 求解 ! 保存 第一载荷步对应最低工作压力 ! 全选
! 删除所有面载荷 ! 选择面组件apax1
! 对接管1施加端部平衡面载荷 ! 选择面组件apax2
! 对接管2施加端部平衡面载荷 ! 选择面组件apax3
! 对接管3施加端部平衡面载荷 ! 选择面组件aouter ! 施加外压 ! 全选 ! 求解 ! 保存
! 退出求解器 进入后处理
! 定义第一载荷步为load case 1 ! 定义第二载荷步为load case 2 ! 读入load case 1 ! 减去load case 2
! 将计算结果存为load case 3
! ! ! PLNSOL,S,INT,0,1 ! 显示应力云图，读取最大应力强度范围 PATH,a1,2 ! 设定路径
PPATH,1,48630 ! 设定路径第一个节点 PPATH,2,47980 ! 设定路径第二个节点
PRSECT, ,0 ! 读取PL＋Pb＋Q范围判断疲劳曲线类型 fini /post1
! ***** 进行疲劳分析
FTSIZE, 1, 1, 2 ! 设定疲劳评定的位置数、事件数及载荷数 !*
FP,1,1e6,2e6,5e6,1e7,2e7,5e7 FP,7,1e8,1e9,1e10, , , FP,13, , , , , , FP,19, ,
FP,21,194,157,127,113,105,99 FP,27,97,96,94, , , FP,33, , , , , , FP,39, , !*
FL,1,48630,1.0,1.0,1.0,try SET,1,last FSNODE,48630,1,1 SET,2,last FSNODE,48630,1,2 !FE,1,-1 FE,1,3.88e5 FTCALC,1,48630 fini
! 根据疲劳曲线输入S－N数据 ! 定义疲劳分析参数 ! 读入第一载荷数据
! 计算并存储疲劳分析节点的各应力分量! 读入第二载荷数据
! 计算并存储疲劳分析节点的各应力分量! 清除以前的疲劳参数与数据 ! 设定事件循环次数 ! 进行疲劳评定
ANSYS11.0疲劳分析实例2
