24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1415  |  回复: 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, 给个红包,谢谢回帖
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的回帖

shgao20

专家顾问 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
热分析时定义力的收敛准则是无效的,在应力分析时需要重新定义。你用惊叹号将它注销了,就不会再起作用了。
  你的 LDREAD 命令中的参数 time 与热分析时定义的计算终止时间 time 无关,在这里只是作为一个参数而已。由于你没有给它赋值过,默认为零,造成 LDREAD 命令读取不到温度结果,或者即使读取了温度结果也不是你需要的。
  在应力分析部分,如果你只计算一次 (某个指定时间的应力),是否定义 TIME 无关紧要;如果你要计算多个时间的应力,则需要使用 TIME 命令来定义每次的计算终止时间。不过需要将 LDREAD 命令放到每次计算之前,目前这个位置也是不对的 - 只读取一次指定时间的温度,不会随着 TIME 的修改而改变。
5楼2019-09-21 09:07:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 346507037 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研经历 +7 cccchenso 2026-03-29 7/350 2026-03-29 15:00 by 唐沐儿
[考研] 求调剂 +9 张zz111 2026-03-27 10/500 2026-03-29 14:57 by 唐沐儿
[考研] 22408 359分调剂 +4 Qshers 2026-03-27 6/300 2026-03-29 14:16 by Qshers
[考研] 318一志愿吉林大学生物与医药 求调剂 +4 笃行致远. 2026-03-28 4/200 2026-03-29 11:55 by longlotian
[考研] 调剂求院校招收 +6 鹤鲸鸽 2026-03-28 6/300 2026-03-29 08:15 by fmesaito
[考研] 一志愿211 初试270分 求调剂 +8 谷雨上岸 2026-03-23 9/450 2026-03-29 01:00 by 我是小康
[考研] 320分,材料与化工专业,求调剂 +9 一定上岸aaa 2026-03-27 13/650 2026-03-28 15:00 by 神马都不懂
[考研] 0703本科郑州大学求调剂 +3 nhj_ 2026-03-25 3/150 2026-03-28 13:24 by Iveryant
[考研] 286求调剂 +4 丢掉懒惰 2026-03-27 7/350 2026-03-28 08:07 by baoball
[考研] 求调剂 +4 零八# 2026-03-27 4/200 2026-03-27 18:07 by yu221
[考研] 考研化学308分求调剂 +10 你好明天你好 2026-03-23 12/600 2026-03-27 14:43 by shangxh
[考研] 333求调剂 +3 question挽风 2026-03-23 3/150 2026-03-27 11:29 by 不吃魚的貓
[考研] 312求调剂 +9 上岸吧ZJY 2026-03-22 13/650 2026-03-27 11:24 by sanrepian
[考研] 286求调剂 +4 lim0922 2026-03-26 4/200 2026-03-27 10:28 by guoweigw
[考研] 调剂求收留 +7 果然有我 2026-03-26 7/350 2026-03-27 00:26 by wxiongid
[考研] 289求调剂 +17 硕星赴 2026-03-23 17/850 2026-03-26 16:18 by 不吃魚的貓
[考研] 335求调剂 +4 yuyu宇 2026-03-23 5/250 2026-03-23 23:49 by Txy@872106
[考研] 一志愿国科过程所081700,274求调剂 +3 三水研0水立方 2026-03-23 3/150 2026-03-23 23:11 by MajorWen
[考研] 284求调剂 +3 yanzhixue111 2026-03-23 6/300 2026-03-23 22:58 by pswait
[考研] 求老师收我 +3 zzh16938784 2026-03-23 3/150 2026-03-23 12:56 by ztnimte
信息提示
请填处理意见