| 查看: 752 | 回复: 3 | |||
zhuqw金虫 (小有名气)
|
[求助]
网格形状对热应力耦合的影响 已有1人参与
|
|
最近在做回流过程的热应力耦合，发现模拟出来的跟实际有点差别，所以就用了一个铜块的模型来做实验。模型是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 三边形网格  PLANE2 四边形网格  PLANE223 四边形网格 |
回复此楼» 猜你喜欢
青椒八年已不青，大家都被折磨成啥样了？
已经有11人回复
-
限项规定
已经有9人回复
-
免疫学博士有名额，速联系
已经有4人回复
-
交叉科学部支持青年基金，对三无青椒是个机会吗？
已经有5人回复
-
国家基金申请书模板内插入图片不可调整大小?
已经有6人回复
-
国家级人才课题组招收2026年入学博士
已经有5人回复
-
Fe3O4@SiO2合成
已经有6人回复
-
青年基金C终止
已经有4人回复
-
26申博求博导推荐-遥感图像处理方向
已经有4人回复
-
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 2015年南京航空航天大学国家自然科学基金资助清单
已经有14人回复
- 关于deform6.1做金属热成型的实例及deform快捷键的使用和deform入门教程
已经有7人回复
- 大型光学系统中的射线追踪模拟
已经有0人回复
- LS-DYNA DYTRAN ABAQUS ADINA ANSYS NASTRAN 有限元软件比较
已经有0人回复
- 过程装备特殊零部件应力分析
已经有144人回复
- 【资源】COMSOL Multiphysics 有关于化学工程和生物医学工程方面的论文
已经有22人回复
- 【资源】SIMULIA ABAQUS 6.7
已经有7人回复
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
2楼2015-12-09 08:05:23
zhuqw
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1217.4
- 帖子: 77
- 在线: 74.6小时
- 虫号: 2408038
- 注册: 2013-04-08
- 性别: GG
- 专业: 半导体光电子器件
|
1.我在帮助文件里面确实找不到PLANE2，但是用命令流编写的时候却没有报错。。。我也很纠结这个单元到底在这个版本中有没有。。。我也是从书上看到说它是PLANE35的等效单元的。。。 2.因为我就是想验证一下铜块自由收缩，也就是四周没有任何约束的时候是不是应该是均匀的向中心收缩的，所以没有加任何约束。那么为什么用三边形网格画出来的结果就没有出现刚体转动呢？我现在在做复杂的3D模型的时候也出问题了，2D的形变结果比较理想，3D的结果正好和2D相反。。。两种模型施加的底面约束一样，都是UY=0，只是3D的因为只建了一半的模型所以对称面用了对称约束，不知道对称约束和UX=0约束对结果影响大不大？ |
3楼2015-12-09 14:51:37
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
4楼2015-12-10 09:19:21