本人的课题方向是激光熔覆方向的,最近在做激光熔覆的温度场和应力场的仿真,温度场做出来的结果正常,然后用间接热结构耦合的方法,先转换模型的单元类型,再将温度场的计算结果作为载荷进行应力场的仿真,仿真时施加的位移约束是自由位移约束,但是应力场仿真的结果算出来的应力大小比正常值大了一个数量级,我看其他论文里人家算出来的正常值应该是10的9次方,也就是几百兆帕,但是我算出来的是10的10次方,就是几千兆帕了,这很明显太大了,而且实际情况中基板熔覆后两边应该是往上翘的,可是我的模拟结果却是往下翘的,请问各位大佬这是什么情况,怎么解决,下面附上我应力场部分的命令流。仿真结果见下图。
FINISH
/COM,Structural
/PREP7
toffst,273
ETCHG,TTS
MPTEMP,,,,,,,,
MPTEMP,1,0,100,200,300,400,600
MPDATA,EX,1,1,1.985E11,1.93E11,1.85E11,1.76E11,1.67E11,1.59E11
MPDATA,ALPX,1,1,17E-6,17.4E-6,18E-6,18.6E-6,19.1E-6,19.6E-6
MPDATA,PRXY,1,1,0.294,0.295,0.301,0.310,0.318,0.326
MPTEMP,7,800,1200,1300,1500,2000
MPDATA,EX,1,7,1.51E11,0.6E11,0.2E11,0.1E11,0.1E11
MPDATA,ALPX,1,7,20.2E-6,20.7E-6,21.1E-6,21.6E-6,21.6E-6
MPDATA,PRXY,1,7,0.333,0.339,0.342,0.388,0.388
MPTEMP,1,20,100,200,300,400,500
MPDATA,EX,2,1,1.951E11,1.912E11,1.857E11,1.796E11,1.726E11,1.645E11
MPDATA,ALPX,2,1,15.24E-6,15.8E-6,16.43E-6,16.97E-6,17.44E-6,17.85E-6
MPDATA,PRXY,2,1,0.267,0.273,0.29,0.31,0.322,0.313
MPTEMP,7,600,700,800,900,1000,1100
MPDATA,EX,2,7,1.55E11,1.441E11,1.314E11,1.168E11,1.001E11,0.811E11
MPDATA,ALPX,2,7,18.21E-6,18.54E-6,18.83E-6,19.11E-6,19.38E-6,19.66E-6
MPDATA,PRXY,2,7,0.296,0.282,0.262,0.24,0.229,0.223
MPTEMP,13,1200,1420,1460,2000
MPDATA,EX,2,13,0.595E11,0.02E11,0.02E11,0.02E11
MPDATA,ALPX,2,13,19.95E-6,20.7E-6,20.7E-6,20.7E-6
MPDATA,PRXY,2,13,0.223,0.223,0.223,0.223
FINISH
/SOLU !求解
ANTYPE,4 !分析类型:瞬态
TRNOPT,FULL !求解方法:完全的N-R方法
TIMINT,ON !打开瞬态效应
KBC,0
LUMPM, 0 !不使用集中质量矩阵
NLGEOM,ON !使用大变形效应
pstres,1 !打开预应力选项
neqit,50
TREF,20 !定义热应力计算参考温度
NSEL,S,LOC,X,0
NSEL,R,LOC,Y,0
NSEL,R,LOC,Z,BaseHeight
D,ALL,UX,0
D,ALL,UY,0
D,ALL,UZ,0
ALLSEL,ALL
NSEL,S,LOC,X,baselength
NSEL,R,LOC,Y,0
NSEL,R,LOC,Z,BaseHeight
D,ALL,UY,0
D,ALL,UZ,0
ALLSEL,ALL
NSEL,S,LOC,X,0
NSEL,R,LOC,Y,BaseWidth
NSEL,R,LOC,Z,BaseHeight
D,ALL,UZ,0
ALLSEL,ALL
*DO,K,1,101,1
LDREAD,TEMP,,,0.1*K,,'3','rth'
OUTRES,ALL,ALL,
TIME,0.1*K
DELTIM,,,,1 !指定时间步长大小
SOLVE
*ENDDO
![热应力分析的结果比正常值大了一个数量级]()
微信图片_20200623184051.png |