| 查看: 935 | 回复: 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  热应力分析 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
国自然面上复盘~欢迎讨论
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
基金不中，共勉
已经有5人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有3人回复
-
我不理解！
已经有23人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有10人回复
-
系统查不到
已经有10人回复
-
怎么查啊
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
caemechanics 
至尊木虫 (正式写手)
- 仿真EPI: 1
- 应助: 38 (小学生)
- 贵宾: 0.218
- 金币: 12469.4
- 红花: 70
- 帖子: 734
- 在线: 258小时
- 虫号: 1361907
- 注册: 2011-08-05
- 专业: 固体力学
2楼2011-12-29 11:11:13
3楼2011-12-29 22:43:42
投票:
