| 查看: 4518 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
ANSYS中采用生死单元计算的热应力，发现被杀死的单元参与了受力，求指教！ 已有1人参与
|
|||
|
各位ANSYS大神： 小弟初次接触热分析，在将计算出来的温度场，施加到温度单元SOLID70转换成实体单元SOLID45后，发现杀死的单元也参与了受力，还有一个问题，是我采用LDREAD命令进行温度场的提取和施加，计算时每一个时间的温度场，都是从1开始然后到指定的时间步（小弟疑虑是该处的时间步导致了被杀死的单元参与了受力，不知道怀疑的对还是不对），由于网格划分的命令流太长，此处就不附加了，现将温度单元转化为实体单元后的命令流附上，请大神帮小弟解惑，小弟在此感激不尽，献上我的大半金币（500），谢谢！ PS:如果命令流中取消掉*IF.....*ELSE这一段命令流的话，第二层的单元就是被杀死的，没有参与受力，如果采用*IF.....*ELSE语句将第二层激活，第二层单元就会在被杀死的状态下参与受力。 转化单元后的命令流： !!!!!求解应力场 /PREP7 ETCHG,TTS HOUR=168 *DIM,HER,ARRAY,1 *CFOPEN,TAN_MO,txt *DO,I,1,HOUR TE00=1.02*(3.55E10)*(1-EXP(-(I/12)*0.09)) *VWRITE,TE00 (F30.5) *ENDDO *CFCLOS *DO,I,1,HOUR MP,ALPX,I,1E-5 !热膨胀系数 *VREAD,HER(1),TAN_MO,txt,,,,,,I-1 (F30.5) MP,EX,I,HER(1) MP,PRXY,I,0.2 MP,DENS,I,2500 *ENDDO !定义边界 NSEL,S,LOC,Y,-0.01,0.01 D,ALL,UY !竖向约束 ALLSEL WPLAN,,15.350,0,-0.225,19.250,1.300,-0.225,19.250,1.300,0.225 NSEL,S,LOC,Z,-0.01,0.01 WPCSYS,-1,0 D,ALL,UY ALLSEL NSEL,S,LOC,Y,-0.01,0.01 NSEL,R,LOC,X,-0.01,0.01 D,ALL,ALL ALLSEL NSEL,S,LOC,X,-0.1,0.1 !对称约束 DSYM,SYMM,X !完成边界的定义 NROPT,FULL !打开牛顿—拉夫 CMSEL,S,DIERCENG !选择已经定义的第二层单元 ESLV,S EKILL,ALL !杀死已经定义的第二层单元 ESEL,S,LIVE NSLE,S,1 NSEL,INVE D,ALL,ALL !约束第二层单元节点的自由度 *DO,I,1,8 TIME,I TIMINT,ON /PREP7 CMSEL,S,DIYICENG ESLV,S EMODI,ALL,MAT,I !修改材料属性 BFEDELE,ALL,TEMP SFADELE,ALL,TEMP *IF,I,GE,6,THEN !如果时间大于等于6，就激活第二层单元 /PREP7 ESEL,S,LIVE NSLE,S NSEL,INVE DDELE,ALL,ALL ! 删除第二层单元节点的约束 CMSEL,S,DIERCENG ESLV,S EALIVE,ALL !激活第二层单元 ALLSEL CMSEL,S,DIYICENG ESLV,S EMODI,ALL,MAT,I ALLSEL CMSEL,S,DIERCENG ESLV,S EMODI,ALL,MAT,I-5 ALLSEL BFEDELE,ALL,TEMP SFADELE,ALL,TEMP LDREAD,TEMP,,,I,,RES1,RTH !提取时间步为I的温度场，并作为荷载施加在结构上 TREF,30 ALLSEL *ELSE CMSEL,S,DIYICENG ESLV,S LDREAD,TEMP,,,I,,RES1,RTH TREF,30 !指定参考温度 ALLSEL *ENDIF /SOLU NLGEOM,ON DELTIM,1,1,1 AUTOTS,ON OUTRES,ALL,ALL SOLVE *ENDDO allsel SAVE FINISH |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有54人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
关于Filecode分析方法
已经有13人回复
-
时间戳又变了8-15
已经有9人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
6楼2017-01-02 13:52:05
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
2楼2017-01-01 08:00:53
3楼2017-01-01 11:19:23
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
4楼2017-01-02 08:05:35
