24小时热门版块排行榜    

查看: 1378  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

346507037

铁虫 (正式写手)

[交流] ansys温度场转应力场计算 已有2人参与

麻烦大神看一下小弟的计算应力场命令流有什么问题吗?
计算一个载荷步的应力场,没有警告,不知道为什么没有计算出来应力值,麻烦大家看看命令流有什么问题吗?是不是求解设置的原因?
/com,structural
/prep7
alls
eplot
lsclear,all    !删除载荷,从数据库中清除实体模型载荷
etchg,tts       !将热单元转换为结构单元

MPTEMP,1,25,100,200,300,400,500,
MPTEMP,7,600,700,800,900,1000,1100, !定义材料属性温度范围
MPDATA,EX,1,1,1.96E11,1.92E11,1.84E11,1.76E11,1.68E11,1.60E11, !弹性模量(Pa)
MPDATA,EX,1,7,1.60E11,1.52E11,1.44E11,1.35E11,1.27E11,1.19E11,

MPDATA,NUXY,1,1,0.30,0.30,0.31,0.31,0.31,0.31,  !泊松比
MPDATA,NUXY,1,7,0.31,0.32,0.32,0.32,0.33,0.33,

MPDATA,ALPX,1,1,15.4E-5,15.8E-5,16.3E-5,17E-5,17.9E-5,19.1E-5,   !热膨胀系数
MPDATA,ALPX,1,7,19.1E-5,15.4E-5,15.8E-5,16.3E-5,17E-5,17.9E-5,


TB,BISO,1,6,2,  
TBTEMP,25   
TBDATA,,3.2E8,1.96E9,,,,
TBTEMP,200  
TBDATA,,1.67E8,1.84E9,,,,
TBTEMP,600  
TBDATA,,1.23E8,1.52E9,,,,
TBTEMP,900
TBDATA,,1.05E8,1.27E9,,,,
TBTEMP,1200
TBDATA,,0.24E8,0.2E9,,,,
TBTEMP,1300
TBDATA,,0.105E8,0.07E9,,,,  


MPTEMP,1,25,100,200,300,400,500,
MPTEMP,7,600,700,800,900,1000,1100, !定义材料属性温度范围
MPDATA,EX,2,1,2.23E11,2.19E11,2.12E11,2.05E11,1.98E11,1.92E11,  !弹性模量(Pa)
MPDATA,EX,2,7,1.84E11,1.76E11,1.67E11,1.53E11,1.43E11,1.43E11,

MPDATA,NUXY,2,1,0.30,0.31,0.32,0.33,0.34,0.35,  !泊松比
MPDATA,NUXY,2,7,0.36,0.37,0.38,0.39,0.40,0.40,  
MPDATA,ALPX,2,1,1.02E-5,1.12E-5,1.21E-5,1.26E-5,1.28E-5,1.36E-5, !热膨胀系数
MPDATA,ALPX,2,7,1.38E-5,1.46E-5,1.54E-5,1.70E-5,1.81E-5,1.81E-5  

TB,BISO,2,6,2,  
TBTEMP,25   
TBDATA,,6.08E8,2.23E9,,,,
TBTEMP,200  
TBDATA,,5.8E8,2.12E9,,,,
TBTEMP,400
TBDATA,,5.12E8,1.98E9,,,,
TBTEMP,600  
TBDATA,,4.7E8,1.84E9,,,,
TBTEMP,800
TBDATA,,2.68E8,1.67E9,,,,   
TBTEMP,1000
TBDATA,,0.78E8,1.43E9,,,,







!*********************施加位移约束***************

csys,11
VSEL,S,MAT,,1
ESLV,S
NSLE,s,1
d,all,ux,0,,,,uy,uz                    
alls
eplot




!********************非线性求解选项**************
/solu
antype,4
trnopt,full            !完全瞬态积分法
lumpm,0                !使用与单元相关的质量矩阵公式
nlgeom,on              !打开大变形选项
pred,on                !打开预测校正
sstif,1                !打开应力刚化功能
solcontrol,on          !自动控制
nropt,full,,on         !完全牛顿-拉普森方法
cnvtol,u,,0.5          !设置收敛值,位移
!cnvtol,f,,0.05           !设置收敛值,力
cutcontrol,plslimit,10,   !设置在一个非线性求解中时间步的减少
eqslv,sparse,,2        !求解器类型
lnsrch,on              !线性搜索
pstres,1               !打开预应力选项   
autots,on              !自动步长        
tintp,,,,1,,           !向后积分        
timint,on              !瞬态效果  
tref,298                !参考温度
kbc,0                  !连续性载荷
neqit,50               !规定每个子步中最大迭代次数,默认为25
timint,on

!*****************杀死单元**************
VSEL,S,MAT,,2                           !*选择增材制造部分体积
ESLV,S                                  !*选择增材制造部分单元
CM,weld_elem,ELEM
EKIll,ALL
ESEL,s,elem,,weld_elem
NSLE,s                                  !*选择增材制造部分节点
CM,weld_node,node
ALLSEL,ALL
ESEL,S,LIVE
EPLOT


!************建立局部坐标系**********
LOCAL,11,0,X0,Y0,Z0,,,,1,1,
csys,11
NSEL,s,node,,weld_node
NSEL,r,loc,x,-a2,a1
NSEL,r,loc,y,-b,b
NSEL,r,loc,z,-0.04E-3,0        !*选择一个载荷步内的节点
CM,cur_node,node
ESLN,s
Ealive,all                   !*激活一个载荷步内的单元
NSEL,s,node,,weld_node
NSEL,u,node,,cur_node
CM,weld_node,node
NSEL,s,node,,cur_node
NSEL,r,loc,x,-a2,a1   !*选择前半热源内的节点
LDREAD,TEMP,,,time,,,RTH

/pbf,temp,,1                 !显示材料的温度
/replot  
TM=tm+tinc
time,tm
KBC,1
AUTOTS,ON
nsubst,10
ALLSEL,ALL

solve
FINISH  
SAVE
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
热分析时定义力的收敛准则是无效的,在应力分析时需要重新定义。你用惊叹号将它注销了,就不会再起作用了。
  你的 LDREAD 命令中的参数 time 与热分析时定义的计算终止时间 time 无关,在这里只是作为一个参数而已。由于你没有给它赋值过,默认为零,造成 LDREAD 命令读取不到温度结果,或者即使读取了温度结果也不是你需要的。
  在应力分析部分,如果你只计算一次 (某个指定时间的应力),是否定义 TIME 无关紧要;如果你要计算多个时间的应力,则需要使用 TIME 命令来定义每次的计算终止时间。不过需要将 LDREAD 命令放到每次计算之前,目前这个位置也是不对的 - 只读取一次指定时间的温度,不会随着 TIME 的修改而改变。
5楼2019-09-21 09:07:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

shgao20

专家顾问 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
1  在你的奇特帖子里,已经告诉你: 如下的 NSEL 命令有问题:
      NSEL,s,node,,weld_node       !  其中两个连续的逗号需要去掉一个,否则可能选择不到任何节点。
    不知为什么你不远修改?
  2  你读取温度结果的命令:
      LDREAD,TEMP,,,time,,,RTH
    其中参数 time 没有赋过值,默认为 0,因此可能没有读取任何温度结果,或者即使读取了温度结果,也是错误的。你需要检查一下,是否正确读取了温度结果。
2楼2019-09-20 08:38:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

346507037

铁虫 (正式写手)

引用回帖:
2楼: Originally posted by shgao20 at 2019-09-20 08:38:45
1  在你的奇特帖子里,已经告诉你: 如下的 NSEL 命令有问题:
      NSEL,s,node,,weld_node       !  其中两个连续的逗号需要去掉一个,否则可能选择不到任何节点。
    不知为什么你不远修改?
  2  你读取温 ...

首先感谢您的回答,这里可能忘记改了。疏忽了。time在计算温度场的时候已经赋值,在计算应力场的时候还的重新赋值吗?
3楼2019-09-20 12:49:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

激光增材

新虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
力的收敛值前边是感叹号

发自小木虫Android客户端
4楼2019-09-20 13:41:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+4 NSFC2026我来了 2026-03-10 11/550 2026-03-12 13:13 by NSFC2026我来了
[考研] 304求调剂(085602一志愿985) +9 化工人999 2026-03-09 9/450 2026-03-12 11:32 by 白云朵朵飞
[考研] 314求调剂 +4 无懈可击的巨人 2026-03-12 4/200 2026-03-12 10:49 by vgtyfty
[考研] 泣血叩求调剂恩,愿以丹心报师恩 +5 Iuruoh 2026-03-11 5/250 2026-03-11 16:19 by 暖阳芬芳
[考研] 288求调剂 +13 王晓阳- 2026-03-09 18/900 2026-03-11 14:01 by Linda Hu
[考研] 293求调剂,一志愿陕师大生物学 +3 ??????.?.??? 2026-03-09 3/150 2026-03-11 10:02 by 学员8dgXkO
[考研] 考研材料与化工,求调剂 +7 戏精丹丹丹 2026-03-09 7/350 2026-03-11 09:15 by @飒飒飒飒
[考研] 312求调剂 +5 陌宸希 2026-03-10 5/250 2026-03-11 07:18 by macy2011
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
[考研] 337一志愿华南理工材料求调剂 +5 mysdl 2026-03-07 5/250 2026-03-10 14:54 by ruiyingmiao
[考研] 一志愿天大化工(085600)调剂总分338 +5 蔡大美女 2026-03-09 5/250 2026-03-10 14:44 by ruiyingmiao
[考研] 材料专硕调剂 +7 慕辰123 2026-03-05 10/500 2026-03-10 14:22 by peike
[考研] 材料调剂 +4 xxxcm 2026-03-08 7/350 2026-03-09 10:11 by EBSD
[教师之家] 交大前校长王树国:现在最先进的科技并不在大学实验室,而是在企业研究院 +4 zju2000 2026-03-08 6/300 2026-03-08 19:15 by zju2000
[考研] 347求调剂 +4 浮云满足 2026-03-07 4/200 2026-03-08 16:46 by 星空星月
[考研] 286求调剂 +12 Faune 2026-03-06 14/700 2026-03-08 14:54 by hmn_wj
[考研] 085701环境工程专业,初试305,均过国家A区线 +7 卡卡来了@ 2026-03-07 8/400 2026-03-08 06:40 by 刘兵
[考研] 求调剂,不管什么专业,我是可塑造的人才一枚,希望遇到知己老师捞捞我 +4 13102137290 2026-03-06 5/250 2026-03-07 21:21 by lissomchan
[考研] 化学290求调剂 +3 两颗 西柚 2026-03-07 4/200 2026-03-07 20:55 by macy2011
[考研] 085602 293分求调剂 +3 SivanNano. 2026-03-05 3/150 2026-03-05 19:41 by caszguilin
信息提示
请填处理意见