| 查看: 1211 | 回复: 5 | ||
| 【悬赏金币】回答本帖问题，作者面对太阳丶将赠送您 10 个金币 | ||
[求助]
温度场耦合后的应力场结果无法实现不同载荷步的焊缝行走问题，求大神解答 已有1人参与
|
||
|
温度场模拟结果可以实现焊缝的行走，显示不同载荷步下的已走过的活焊缝，显示正常，过程图见图1、2、3。 但经过耦合后，应力场结果无法实现焊缝的行走，只显示所有的焊缝，过程图见4、5、6。 我检查了好几遍命令流，但没找到问题的所在，请大神帮忙解答一下。 命令流如下： !******************************第二部分 应力场分析****************************** /prep7 alls eplot etchg,tts !热模型转化为结构计算模型 lsclear,solid !清除相关热计算载荷 !***********施加位移约束*************** asel,s,area,,106 !选择两侧节点 nsla,s,1 nplot d,all,all,0 asel,a,area,,112 nsla,s,1 nplot d,all,all,0 !施加约束 alls eplot !***********非线性求解选项************** /solu antype,4 trnopt,full !完全瞬态积分法 lumpm,0 !使用与单元相关的质量矩阵公式 nlgeom,on !打开大变形选项 pred,on !打开预测校正 sstif,1 !打开应力刚化功能 solcontrol,on !自动控制 nropt,full,,on !完全牛顿-拉普森方法 cnvtol,u,,0.5 !设置收敛值，位移 !cnvtol,f,,0.05 !设置收敛值，力 cutcontrol,plslimit,10,!设置在一个非线性求解中时间步的减少 eqslv,sparse,,2 !求解器类型 lnsrch,on !线性搜索 pstres,1 !打开预应力选项 autots,on !自动步长 tintp,,,,1,, !向后积分 timint,on !瞬态效果 tref,25 !参考温度 kbc,0 !连续性载荷 neqit,50 !规定每个子步中最大迭代次数,默认为25 !***********死亡************** /solu /auto,1 /replot alls eplot vsel,s,volume,,1,10 eplot eslv,s,all ekill,all alls esel,s,live eplot !**************稳态分析得到的初始应力值************** ldread,temp,1,,,,,rth !读取第一步的温度场温度 /pbf,temp,,1 /replot time,1e-5 deltim,1e-5,1e-5,1e-5 solve !**************焊缝瞬态分析************** alls eplot zmin=-1E-5 zmax=0.005+1E-5 zd=0.005 k=1 !(选择焊缝体号) j=1.93 !(时间) h=1 timint,on *do,m,1,9 !(9道焊缝循环) *do,n,1,30 !一道焊缝30步 /solu vsel,s,volu,,k nslv,s,1 nsel,r,loc,z,zmin,zmax esln,s,all ealive,all alls eplot ldread,temp,h+1,,,,,rth! /pbf,temp,,1 /replot !****************求解********************** time,j deltim,1.93,1.93,1.93 alls solve j=j+1.93 h=h+1 zmin=zmin+zd zmax=zmax+zd *enddo !单道焊缝走完 alls eplot j=j-1.93 !(8分钟冷却) j=j+240 *do,n,151,152 /solu ldread,temp,h+1,,,,,rth time,j deltim,240,240,240 alls solve j=j+240 h=h+1 *enddo j=j-240 j=j+1.93 zmin=zmin-30*zd !减去多加的步数回到最初的起点（不减会走出模型） zmax=zmax-30*zd k=k+1 *enddo !**************************最后一道******************************** *do,n,1,30 /solu vsel,s,volu,,k nslv,s,1 nsel,r,loc,z,zmin,zmax esln,s,all ealive,all eplot ldread,temp,h+1,,,,,rth time,j deltim,1.93,1.93,1.93 alls solve j=j+1.93 h=h+1 zmin=zmin+zd zmax=zmax+zd *enddo j=j-1.93 !(8分钟冷却) j=j+240 *do,n,151,152 /solu ldread,temp,h+1,,,,,rth time,j deltim,240,240,240 alls solve j=j+240 h=h+1 *enddo j=j-240  图1 第61步温度场结果.jpg  图2 第150步温度场结果.jpg  图3 第300步温度场结果.jpg  图4 第20步应力场结果.png  图5 第150步应力场结果.png  图6 第320步应力场结果.png |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有11人回复
2楼2022-04-28 14:08:22
3楼2022-04-28 15:00:02
4楼2022-04-29 15:21:05
有可以参考的命令流么，求发，我依旧没找到这个应力场焊缝不能移动的原因5楼2022-04-29 16:13:43
6楼2022-04-30 13:43:51
40