| 查看: 841 | 回复: 2 | ||
[求助]
平板加热的热应力分析的结果中，位移与理论不太一致，急求帮助！！！
|
|
各位大神，我现在在做移动热源加热平板的热应力分析，采用间接热应力分析的方法，如图，在板的左下角网格最密的地方施加高斯移动热源（板是对称的，取了一半做分析），然后在结构分析的时候读入热分析的温度场，按照理论，结构分析的结果中所有节点的X方向位移应该都是负值，但计算得到的结果却有很多节点的位移是正值，Y向和Z向也有同样的问题，哪位大神能帮我看看是那方面出了问题么,多谢多谢啊！！！ 急啊！！！急啊！！！ 我的结构分析的命令流如下： /prep7 !删除面8上的表面效应单元 asel,s,area,,8 aclear,all etdel,3 asel,all !将热单元转换成结构单元 etchg,tts !材料热物性 mpdata,alpx,,1,12e-6,12.5e-6,14e-6,14.25e-6,14.65e-6,15.35e-6 !定义线膨胀系数 mpdata,alpx,,7,15.7e-6,15.85e-6,16.4e-6,16.6e-6,17.7e-6,16.9e-6 MPDATA,PRXY,1,1,0.3,0.319,0.338,0.357,0.376,0.395 !定义泊松比 MPDATA,PRXY,1,7,0.414,0.433,0.452,0.471,0.49,0.49 fy=275E+6 !常温下屈服应力 exx=2.1E+11 !常温下杨氏模量 MPTEMP !清除旧的温度场 MPTEMP,,20,100,200,300,400 !定义随温度变化的杨氏模量 MPDATA,EX,1,,exx,exx,0.9*exx,0.8*exx,0.7*exx MPTEMP,,500,600,700,800,900 MPDATA,EX,1,,0.6*exx,0.31*exx,0.13*exx,0.09*exx,0.0675*exx !MP,NUXY,1,0.3 !定义泊松比 !MP,ALPX,1,1.4E-5 !定义热膨胀系数 !定义随温度变化的应力-应变关系 TB,MISO,1,10,3 !共10个温度，每个温度时的应力-应变由3个点描述 TBTEMP,20 !20度时的应力-应变关系 TBPT,,fy/exx,fy TBPT,,0.02,fy TBPT,,0.15,fy TBTEMP,100 !100度时的应力-应变关系 TBPT,,fy/exx,fy TBPT,,0.02,fy TBPT,,0.15,fy TBTEMP,200 !200度时的应力-应变关系 TBPT,,0.807*fy/(0.9*exx),0.807*fy TBPT,,0.02,fy TBPT,,0.15,fy TBTEMP,300 !300度时的应力-应变关系 TBPT,,0.613*fy/(0.8*exx),0.613*fy TBPT,,0.02,fy TBPT,,0.15,fy TBTEMP,400 !400度时的应力-应变关系 TBPT,,0.420*fy/(0.7*exx),0.420*fy TBPT,,0.02,fy TBPT,,0.15,fy TBTEMP,500 !500度时的应力-应变关系 TBPT,,0.360*fy/(0.6*exx),0.360*fy TBPT,,0.02,0.780*fy TBPT,,0.15,0.780*fy TBTEMP,600 !600度时的应力-应变关系 TBPT,,0.180*fy/(0.310*exx),0.180*fy TBPT,,0.02,0.470*fy TBPT,,0.15,0.470*fy TBTEMP,700 !700度时的应力-应变关系 TBPT,,0.075*fy/(0.130*exx),0.075*fy TBPT,,0.02,0.230*fy TBPT,,0.15,0.230*fy TBTEMP,800 !800度时的应力-应变关系 TBPT,,0.050*fy/(0.090*exx),0.050*fy TBPT,,0.02,0.110*fy TBPT,,0.15,0.110*fy TBTEMP,900 !900度时的应力-应变关系 TBPT,,0.0375*fy/(0.0675*exx),0.0375*fy TBPT,,0.02,0.060*fy TBPT,,0.15,0.060*fy fini /solu !antype,trans lsclear,all !从数据库文件中删除载荷步选项 !定义刚性位移约束（是不是这个地方出错？？） d,node(0,0,-thick),uy d,node(0,0,-thick),uz d,node(0,width,-thick),uz asel,s,area,,6,9,3 asel,a,area,,24 da,all,symm asel,all tref,23 !参考温度 nlgeom,on kbc,0 *do,t,1e-6,tm,tinc*10 TIME,T lDREAD,TEMP,,,T,,,RTH, nsubst,2,10,1 alls solve *enddo *do,nc,1,300,10 time,tm+nc*tinc lDREAD,TEMP,,,tm+nc,,,RTH, NSUBST,2,5,1 solve *enddo *do,nc2,1,600,10 time,tm+(nc+nc2)*tinc lDREAD,TEMP,,,tm+nc+nc2,,,RTH, NSUBST,1,5,1 solve *enddo *do,nc3,1,5000,100 time,tm+(nc+nc2+nc3)*tinc lDREAD,TEMP,,,tm+nc+nc2+nc3,,,RTH, NSUBST,1,5,1 solve *enddo OUTRES,ALL,last fini save  热应力分析 |
回复此楼» 猜你喜欢
基金申报
已经有5人回复
-
基金委咋了？2026年的指南还没有出来？
已经有7人回复
-
国自然申请面上模板最新2026版出了吗？
已经有17人回复
-
纳米粒子粒径的测量
已经有8人回复
-
疑惑？
已经有5人回复
-
计算机、0854电子信息（085401-058412）调剂
已经有5人回复
-
Materials Today Chemistry审稿周期
已经有5人回复
-
溴的反应液脱色
已经有7人回复
-
推荐一本书
已经有12人回复
-
常年博士招收(双一流，工科)
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
caemechanics 
至尊木虫 (正式写手)
- 仿真EPI: 1
- 应助: 38 (小学生)
- 贵宾: 0.218
- 金币: 12467.8
- 红花: 70
- 帖子: 734
- 在线: 258小时
- 虫号: 1361907
- 注册: 2011-08-05
- 专业: 固体力学
2楼2011-12-29 11:11:13
3楼2011-12-29 22:43:42