24小时热门版块排行榜    

查看: 2518  |  回复: 7

founderyuan

新虫 (初入文坛)

[交流] ansys模拟残余应力的问题,不会写命令流 已有4人参与

/FILNAME,EXERCISE1                              
/TITLE,TRANSIENT THERMAL ANALYSYS TO A SPHERE     
KEYW,PR_THERM,1                                 

/PreP7   


MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,1780

MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,700
MPDATA,C,1,,1050
MPDATA,C,1,,1130
MPDATA,C,1,,1170
MPDATA,C,1,,1210
MPDATA,C,1,,1260
MPDATA,C,1,,1300

MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,700
MPDATA,KXX,1,,96.4  
MPDATA,KXX,1,,101   
MPDATA,KXX,1,,105   
MPDATA,KXX,1,,109   
MPDATA,KXX,1,,113   
MPDATA,KXX,1,,150   




/prep7        
/units,si      
rb=0.001*0.1     
length=0.001*5   
width=0.001*5     
height=0.001*1


et,1,solid70   
et,2,solid70   


block,0,length,-0.5*width,0.5*width,0,height
Block,0,length,-3*Rb,3*Rb,0,height
Block,0,length,-12*Rb,12*Rb,0,height
VOVLAP,ALL   



!!!!!!!!!!!!!!!!!!!!!!!!!!!!!主要焊接部位
VSEL,S, , ,       2
ASLV,S  
LSLA,S  
lesize,all,rb
ALLSEL,ALL  

type,1
vsweep,2



!!!!!!!!!!!!!!!!!!!!!!!!!!!!!上下次要区域
VSEL,S, , ,       4
ASLV,S  
LSLA,S  
lesize,all,3*rb
ALLSEL,ALL  

VSEL,S, , ,       5
ASLV,S  
LSLA,S  
lesize,all,3*rb
ALLSEL,ALL

type,1
vsweep,4
vsweep,5



!!!!!!!!!!!!!!!!!!!!!!!!!!!!!焊接部位与次要区域的链接区域
MSHAPE,1,3D
MSHKEY,0

type,2
vmesh,6
vmesh,7






/Solu
antype,trans,new   
nlgeom,on         
timint,1,struct     
Timint,1,therm      
timint,0,mag     
Timint,0,elect   
Tintp,0.005,-1,0.5,0.2   

nropt,auto,,      
autos,on           
outres,,all         


tref,25                      !!!!!初始条件,室温,25°
SFA,ALL,1,CONV,62.3,25       !!!!!对流换热系数62.3



p=600                       !!激光功率
ef=1                      !!激光吸收率


                     
r=0.001*0.1                      !!激光热源半径      
qmax=(p*ef)/(3.14*r*r)
v=0.004                      !!焊接速度为250mm/min         
tinc=rb/v         



VSEL,S, , ,       2                 !选取焊接部分的节点
ESLV,S  
*get,nemax,elem,,num,max                    
*get,nemin,elem,,num,min
ALLSEL,ALL

max_time=length/rb+1   
      

*do,k,1,max_time,1  
a=v*(k-1)*tinc
c=height

*do,i,nemin,nemax,1
   xsy=centrx(i)               !读取该单元的中心坐标
   ysy=centry(i)
   zsy=centrz(i)
   rr=abs(sqrt((xsy-a)*(xsy-a)+ysy*ysy+(zsy-c)*(zsy-c)))    !该单元中心距离热源中心
   *if,rr,le,r,then                     !该单元中心在加热半径范围内
    qr=qmax*exp(-(rr*rr)/(r*r))         !该单元中心处的热流大小
    sfe,i,,hflux,,qr,qr,qr,qr           !在每个单元所指定的面上施加热源载荷
   *endif
*enddo

time,k*rb/v
nsubst,5                !设定每个时间步长内有5个子步  
solve                    !计算

SFEDELE,ALL,1,HFLUX
*enddo

save                     





这是找别人写的命令流,但是那哥们儿临时有事儿写不了了,我对这个不太懂模拟了一下 好像是有温度场,找不到应力的结果。有没有哪位大神能帮忙给写下 应力场的命令流。跪谢了。十万火急啊~!!  285466480这是小弟的QQ号 有什么问题可以联系。必有重谢。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

founderyuan

新虫 (初入文坛)

别沉啊
2楼2014-05-31 13:24:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

founderyuan

新虫 (初入文坛)

3楼2014-05-31 16:58:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

founderyuan

新虫 (初入文坛)

大神都去哪了
4楼2014-06-02 17:25:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

285561706

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
单元首先是可以热应力分析单元solid45之类   再参数补充弹性模量  泊松比等物性参数 可以参考下《基于温度和应力场的焊接残余应力数值分析》 或者其他“热应力”相关的例子
、简单就好、
5楼2014-06-04 09:38:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

一缕轻烟01

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
你好,你现在的这个搞得怎么样啦,我现在也在学着做残余应力相关的东西,可否请教一下啊,谢谢!
6楼2014-08-12 21:11:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

346507037

铁虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
热应力分析很简单的,只要温度场计算准确,自己就可以多琢磨一下
7楼2019-06-11 18:57:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

曹枣林

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我看楼主是模拟激光焊的吧,我目前在做这方面的。这个是温度场的命令流,没有涉及应力场的命令流。
8楼2019-06-11 19:51:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 founderyuan 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +15 低垂的野花 2026-07-31 20/1000 2026-08-01 14:10 by gltch
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 18/900 2026-08-01 13:56 by gltch
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jzg3xk1po8 2026-07-31 8/400 2026-08-01 13:53 by y4XuShGyllvT
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 13:51 by y4XuShGyllvT
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 13:51 by y4XuShGyllvT
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jzg3xk1po8 2026-07-31 7/350 2026-08-01 13:49 by y4XuShGyllvT
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 13:49 by y4XuShGyllvT
[基金申请] 小木虫看见有人已经知道结果了 +16 1234567wang 2026-07-25 20/1000 2026-08-01 08:57 by type2type
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 8/400 2026-07-31 18:13 by beefly
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[基金申请] 今年的WR进展到哪一步了? +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[有机交流] 这个自发加氧反应的机理是什么? 50+4 YUAN2273 2026-07-25 7/350 2026-07-29 12:27 by xiaobenpu
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
请填处理意见