| 查看: 5694 | 回复: 9 | |||
[求助]
ANSYS 焊接残余应力模拟 已有1人参与
|
|
热分析时热源荷载是取得生热率和生死单元的方法，这是其中一个荷载步的温度场。主要命令流如下： /CLEAR /FILNAME, exercise10 /title,weldingstress10 /prep7 !**************定义参数及热源******************** L=0.2 LSIZE=0.002 MAX_X=L/LSIZE A= 0.2827E-04 V=0.004 U=25 I=180 K=0.75 HGENUM=(K*U*I)/(A*V) !**************定义单元及材料参数******************** et,1,mesh200 KEYOPT,1,1,6 ET,2,solid70 dof,temp mp,dens,1,7850 mptemp,1,20,100,200,400,600,800 mptemp,7,1000,1200,1500,1700 mpdata,kxx,1,1,48.2,50.7,46.5,38.1,30,26 mpdata,kxx,1,7,26.8,29.2,35,14 mpdata,c,1,1,470,490,515,590,725,835 mpdata,c,1,7,650,654,660,780 mptemp,1,20,200,400,700,750,880 mptemp,7,1000,1200,1500,20000 mpdata,ex,1,1,20.6e10,19e10,17e10,9e10,2e10,1.5e10 mpdata,ex,1,7,1.3e10,1.0e10,0.5e10,0.1e10 mpdata,prxy,1,1,0.28,0.29,0.31,0.35,0.40,0.45 mpdata,prxy,1,7,0.46,0.47,0.48,0.48 tb,biso,1,6,2 tbtemp,20,1 tbdata,,345e6, tbtemp,200,2, tbdata,,250e6, tbtemp,400,3, tbdata,,100e6, tbtemp,700,4, tbdata,,50e6, tbtemp,880,5, tbdata,,30e6, tbtemp,2500,6, tbdata,,20e6, mp,alpx,1,1.06e-5 mpcopy,,1,2 tb,biso,2,6,2 tbtemp,20,1 tbdata,,345e6, tbtemp,200,2, tbdata,,250e6, tbtemp,400,3, tbdata,,100e6, tbtemp,700,4, tbdata,,50e6, tbtemp,880,5, tbdata,,30e6, tbtemp,2500,6, tbdata,,20e6, !**************建立模型******************** rectng,-0.002,0.002,-0.016,0 rectng,-0.0104853,-0.002,-0.016,0 rectng,0.002,0.0104853,-0.016,0 rectng,-0.030,-0.0104853,-0.016,0 rectng,0.0104853,0.030,-0.016,0 rectng,-0.120,-0.030,-0.016,0 rectng,0.030,0.120,-0.016,0 cyl4,-0.002,0,0.0084853,135,,180 cyl4,0.002,0,0.0084853,0,,45 aglue,all k,34,-0.008,0.016 k,35,0.008,0.016 a,29,4,3,33,35,34,29 rectng,-0.008,0.008,0.016,0.036 rectng,-0.008,0.008,0.036,0.180 aglue,all numcmp,all type,1 esize,0.002 amesh,1,2 amesh,5,8 esize,0.005 amesh,4 amesh,11,12 esize,0.002 amesh,3 amesh,9,10 EXTOPT,aclear,1, EXTOPT,attr,0, ESIZE,,100 type,2 mat,1 Vext,5,6,,,,0.2 EXTOPT,aclear,1, EXTOPT,attr,0, ESIZE,,100 type,2 mat,2 Vext,1,4,1,,,0.2 Vext,7,12,1,,,0.2 Nummrg,elem Nummrg,node Nummrg,kp !**************非线性分析求解选项设置******************** /solu antype,trans trnopt,full nropt,full,,on AUTOTS,on lnsrch,on !**************杀死焊缝单元******************** allsel /replot esel,s,mat,,1 ekill,all allsel esel,s,live eplot !**************施加对流荷载******************** allsel asel,s,,,1,4 asel,a,,,7,12 asel,a,,,21,22 asel,a,,,26 asel,a,,,29 asel,a,,,31,33 asel,a,,,35,41 asel,a,,,45,46 asel,a,,,50,52 asel,a,,,54,55 asel,a,,,57,65 esla,s nsle,s sf,all,conv,30,20 !**************施加初始温度条件******************** nsel,all tunif,20 !**************开始热分析******************** timint,on,therm *do,i,1,40 esel,s,mat,,1 nsle,s nsel,r,loc,z,0.002*(i-1),0.002*i esln,s,all ealive,all eplot bfe,all,HGEN,,HGENUM time,0.5*i kbc,1 DELTIM,0.25,0.25,0.5 allsel,all OUTRES,ALL,last tsres,erase SOLVE bfeDELE,all,HGEN,,HGENUM allsel,all *ENDDO !**************冷却过程******************** time,60 AUTOTS,1 NSUBST,20,20,20 KBC,0 TSRES,ERASE solve time,100 AUTOTS,1 NSUBST,8,8,8 KBC,0 TSRES,ERASE solve time,1100 AUTOTS,1 NSUBST,20,20,20 KBC,0 TSRES,ERASE solve time,2100 AUTOTS,1 NSUBST,10,10,10 KBC,0 TSRES,ERASE solve !**************进入应力分析******************** /prep7 et,3,solid45 dof,ux,uy,uz allsel eplot etchg,tts lsclear,solid !**************施加位移边界条件******************** asel,s,,,38 nsla,s,1 dsym,symm,y asel,s,,,64 nsla,s,1 d,all,ux,0 asel,s,,,61 nsla,s,1 d,all,ux,0,,,,uy allsel eplot !**************非线性分析求解选项设置******************** /solu antype,trans trnopt,full nropt,full,,on nlgeom,on AUTOTS,on lnsrch,on neqit,50 !**************杀死焊缝单元******************** allsel /replot esel,s,mat,,1 ekill,all allsel esel,s,live eplot !**************应力分析******************** timint,on *do,i,1,100 esel,s,mat,,1 nsle,s nsel,r,loc,z,0.002*(i-1),0.002*i esln,s,all ealive,all eplot time,0.5*i kbc,1 DELTIM,0.25,0.25,0.5 ldread,temp,i,,,,,rth allsel,all OUTRES,ALL,last tsres,erase SOLVE *ENDDO allsel,all eplot time,60 AUTOTS,1 NSUBST,20,20,20 KBC,0 TSRES,ERASE solve time,100 AUTOTS,1 NSUBST,8,8,8 KBC,0 TSRES,ERASE solve time,1100 AUTOTS,1 NSUBST,20,20,20 KBC,0 TSRES,ERASE solve time,2100 AUTOTS,1 NSUBST,10,10,10 KBC,0 TSRES,ERASE solve 可产生的温度为什么没有扩散，是模型网格还是单元的问题，希望能帮忙修改下，然后就是转应力场分析时有问题，收敛也存在问题。  360截图20160623094633177.jpg |
回复此楼» 收录本帖的淘帖专辑推荐
 ANSYS命令流 |
» 猜你喜欢
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有18人回复
-
物理学I论文润色/翻译怎么收费?
已经有210人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有20人回复
-
考博自荐
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- ANSYS 模拟焊接残余应力
已经有5人回复
- ansys模拟焊接金属的残余应力，从哪些方面分析结果呢？？
已经有1人回复
- 中国一重焊接研究所人员招聘
已经有11人回复
- ansys模拟残余应力的问题，不会写命令流
已经有7人回复
- ANSYS中如何实现在焊接残余应力基础之上模拟计算热应力？？？
已经有1人回复
- ansys焊接热应力分析温度场是对的，残余应力分布是错的
已经有11人回复
- Ansys中一个模型的残余应力导入另一模型问题
已经有0人回复
- ansys焊接残余应力模拟
已经有6人回复
- ANSYS中如何导入节点的应力
已经有1人回复
- ANSYS中应力强度因子与J积分的计算[转载]
已经有5人回复
- 求ANSYS用于焊接分析的生死单元的例子
已经有4人回复
yini0302
木虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 1941.9
- 散金: 69
- 红花: 1
- 帖子: 232
- 在线: 244小时
- 虫号: 2178416
- 注册: 2012-12-10
- 性别: GG
- 专业: 机械结构强度学
2楼2016-06-23 12:47:54
xiaobai1707
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 532.1
- 帖子: 121
- 在线: 38小时
- 虫号: 3663776
- 注册: 2015-01-26
- 性别: GG
- 专业: 机械结构强度学
3楼2016-06-23 23:39:28
4楼2016-06-24 14:00:27
5楼2016-06-24 14:04:42
6楼2016-06-24 18:35:35
JLU_大母冲
木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 4383.9
- 散金: 29
- 红花: 3
- 帖子: 1071
- 在线: 72.6小时
- 虫号: 2367514
- 注册: 2013-03-22
- 专业: 资源化工
7楼2016-06-26 06:51:48
8楼2017-11-23 16:48:17
9楼2017-12-08 00:02:40
10楼2018-08-10 10:36:40
5