最近在做回流过程的热应力耦合,发现模拟出来的跟实际有点差别,所以就用了一个铜块的模型来做实验。模型是2D的,我发现用三角形网格和用四边形网格画出来的回流过程铜块的应变结果很不一样,三角形的网格还比较符合实际,四边形网格就很离谱,感觉整个模型在绕中心旋转。并且我试用了几种单元类型,包PLANE13,PLANE2,PLANE223,发现PLANE2的变形云图形状最好。我想请教一下网格形状和单元类型对于结构分析的影响是怎么样的?模型复杂后怎么样选择最合理的单元类型?为什么四边形网格和三边形网格的结果差的这么大?下面是模拟的结果和命令流
先是三角形网格命令流和结果,总位移云图结果比较符合实际情况,铜块在降温过程中是沿着X,Y方向均匀收缩的
FINISH
/CLEAR
/TITLE,2D Model
/TRIAD,OFF
/UNITS,uMKS
/PREP7
ET,1,PLANE2
!1 Cu
MP,EX,1,0.128E+6
MP,ALPX,1,0.165E-4
MP,DENS,1,0.896E-14
MP,KXX,1,3.98E+8
MP,C,1,0.38456E+15
MP,PRXY,1,0.36
RECTNG,0,24890,0,7950 !A1 Cu
NUMCMP,AREA
NUMCMP,NODE
/PNUM,AREA,1
MSHAPE,1,2D
MSHKEY,0
ESIZE,313
MAT,1 !Cu
AMESH,1
/PNUM,AREA,1
/PNUM,VOLU,1
/SOLU
ANTYPE, STATIC
AUTOTS, ON
TREF,156
BFUNIF,TEMP,25
ALLSEL
SOLVE
FINISH
接下来是四边形网格,整个云图好像绕着中心旋转了,和实际降温过程铜块收缩情况不符
FINISH
/CLEAR
/TITLE,2D Model
/TRIAD,OFF
/UNITS,uMKS
/PREP7
ET,1,PLANE2
!1 Cu
MP,EX,1,0.128E+6
MP,ALPX,1,0.165E-4
MP,DENS,1,0.896E-14
MP,KXX,1,3.98E+8
MP,C,1,0.38456E+15
MP,PRXY,1,0.36
RECTNG,0,24890,0,7950 !A1 Cu
NUMCMP,AREA
NUMCMP,NODE
/PNUM,AREA,1
MSHAPE,0,2D
MSHKEY,1
ESIZE,100
MAT,1 !Cu
AMESH,1
/PNUM,AREA,1
/PNUM,VOLU,1
/SOLU
ANTYPE, STATIC
AUTOTS, ON
TREF,156 !初始温度156度,降到20度
BFUNIF,TEMP,25
ALLSEL
SOLVE
FINISH
最后再加一个PLANE223单元,四边形网格的,整个变形就离谱的不行了
FINISH
/CLEAR
/TITLE,2D Model
/TRIAD,OFF
/UNITS,uMKS
/PREP7
ET,1,PLANE223
KEYOPT,1,1,11
KEYOPT,1,3,2
!1 Cu
MP,EX,1,0.128E+6
MP,ALPX,1,0.165E-4
MP,DENS,1,0.896E-14
MP,KXX,1,3.98E+8
MP,C,1,0.38456E+15
MP,PRXY,1,0.36
RECTNG,0,24890,0,7950 !A1 Cu
NUMCMP,AREA
NUMCMP,NODE
/PNUM,AREA,1
SMRTSIZE,1
MAT,1 !Cu
AMESH,1
/PNUM,AREA,1
/PNUM,VOLU,1
/SOLU
ANTYPE, STATIC
AUTOTS, ON
TREF,156 !初始温度156度,降到20度
BFUNIF,TEMP,25
ALLSEL
SOLVE
FINISH
![网格形状对热应力耦合的影响]()
PLANE2 三边形网格
![网格形状对热应力耦合的影响-1]()
PLANE2 四边形网格
![网格形状对热应力耦合的影响-2]()
PLANE223 四边形网格 |