24小时热门版块排行榜

>论坛更新日志 (4201)
>文献求助 (447)
>虫友互识 (361)
>导师招生 (357)
>考博 (195)
>博后之家 (176)
>硕博家园 (160)
>休闲灌水 (148)
>论文投稿 (129)
>基金申请 (92)
>考研 (69)
>招聘信息布告栏 (60)
>教师之家 (56)
>公派出国 (45)
>绿色求助(高悬赏) (43)
>找工作 (38)

返回列表

founderyuan

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1.5
帖子: 7
在线: 56分钟
虫号: 2826623
注册: 2013-11-25

[交流] 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号有什么问题可以联系。必有重谢。

回复此楼