24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ymh0427

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 47.3
散金: 10
红花: 1
帖子: 230
在线: 91.9小时
虫号: 3431778
注册: 2014-09-21
性别: GG
专业: 金属材料的凝固与结晶学

[求助] ansys 焊接应力场问题已有2人参与

下面是我做的焊接应力场，温度场是对的，但是算应力场时老出现下面的问题，求大神解决！！！
附上命令流。
finish
/filname
/filname,al_fe_weld
/title,al_fe_weld
/units,si
btol,1e-7
/prep7
len=0.1
wid=0.05
thic=0.0015
lnum=20
i=100
u=10
v=0.5/60
pi=acos(-1)
eta=0.8
k=1.2e5
them0=25
qmax=k*eta*u*i/pi
radius=0.06
et,1,solid70
mp,dens,1,8030
mp,dens,2,4850
mp,alpx,1,1.78e-5
mp,alpx,2,9.36e-6
mp,kxx,1,16.3
mp,kxx,2,7.44
mp,c,1,502
mp,c,2,544
k,1,0,0,0
k,2,0,0,0.1
l,1,2
rectng,0,0.0005,0,0.001
rectng,0,0.05,0,0.001
rectng,0,-0.05,-0.0005,0
rectng,0,-0.05,0,0.001
rectng,0,-0.05,0,0.0015
pcirc,0,0.0005,270,360
wpoffs,0,0.001,0
pcirc,0,0.0005,0,90
wpcsys
vdrag,all,,,,,,1
vptn,all
vglue,all
nummrg,all
numcmp,volume
vplot
/view,,2,4,3
/replot
allsel
lsel,s,,,88,91
lesize,all,,,50,10
allsel
lsel,a
lsel,s,,,10
lsel,a,,,44
lesize,all,,,50,0.1
allsel
lsel,s,,,92
lsel,a,,,95
lsel,a,,,20
lsel,a,,,65
lsel,a,,,93
lsel,a,,,96
lesize,all,,,50,10
allsel
lsel,s,length,,0.1
lesize,all,,,100
allsel
esize,0.0005
vsel,s,,,1
vsel,a,,,4
mat,1
vsel,all
vsel,s,,,2,3
vsel,a,,,5,7
mat,2
vsel,all
vsweep,all,,,0
allsel
finish
/solu
esel,all
*get,nemax,elem,,num,max
*get,nemin,elem,,num,min
alls
antype,trans
trnopt,full
nropt,full,,on
pred,on
timint,0,struct
timint,1,thern
tintp,0.005,,,1,0.5,0.2
tref,25
t=0
dt0=1e-6
dt1=0.5
tinc=len/(lnum*v)
time,dt0
timint,off
kbc,1
nsubst,1
esel,all
nsle,s
ic,all,temp,25
alls
outres,all,all
solve
esel,all
eplot
*do,im,0,lnum,0.1
vsel,s,,,4
vsel,a,,,3
vsel,a,,,7
eslv,r
nsle
asel,s,area,,30
asel,a,area,,47
asel,a,area,,52
nsla,s,1
esln,r,0
c=v*tinc*im
b=thic
a=-0.0005
tm=dt1+im*tinc
time,tm
antype,4,vtrest
timint,on
autots,on
kbc,1
deltim,0.1,0.1,0.2
*do,i,nemin,nemax,1
*if,esel(i),eq,1,then
xsy=centrx(i)
ysy=centry(i)
zsy=centrz(i)
rr=abs(sqrt((xsy-a)*(xsy-a)+(ysy-b)*(ysy-b)+(zsy-c)*(zsy-c)))
*if,rr,le,radius,then
qr=qmax*exp(-k*rr*rr)
sfe,i,nmface(i),hflux,,qr
*endif
*endif
*enddo
asel,s,mat,,1
asel,a,mat,,2
asel,u,,,30
asel,u,,,47
asel,u,,,52
nsla,r
sf,all,conv,200,293
allsel
outres,all,all
solve
esel,all
vsel,s,,,4
vsel,a,,,3
vsel,a,,,7
eslv,r
nsle
asel,s,area,,30
asel,a,area,,47
asel,a,area,,52
nsla,r,1
esln,r,0
*do,i,1,6
sfedele,all,i,hflux
*enddo
allsel
*enddo
esel,all
eplot
save
finish
/filname
/filname,heat_stress
/prep7
allsel
esel,s,,type,2
etchg,tts
mptemp,1,25,500,1000,1500,2000
mpdata,ex,1,1,1.93e11,1.5e11,0.7e11,0.1e11,0.01e11
mpdata,ex,2,1,1.02e11,0.5e11,0.08e11,0.001e11,0.0001e11
tb,bkin,1,5
tbtemp,25,1
tbdata,1,1200e6,0.193e11
tbtemp,500,2
tbdata,1,933e6,0.15e11
tbtemp,1000,3
tbdata,1,435e6,0.07e11
tbtemp,1500,4
tbdata,1,70e6,0.01e11
tbtemp,2000,5
tbdata,1,7e6,0.001e11
tb,bkin,2,5
tbtemp,25,1
tbdata,1,800e6,0.102e11
tbtemp,500,2
tbdata,1,400e6,0.05e11
tbtemp,1000,3
tbdata,1,70e6,0.008e11
tbtemp,1500,4
tbdata,1,1e6,0.0001e11
tbtemp,2000,5
tbdata,1,0.1e6,0.00001e11
mp,nuxy,1,0.29
mp,nuxy,2,0.3
finish
/solu
antype,trans
solcontrol,off
neqit,100
lumpm,0
allsel
tref,25
autots,on
kbc,1
lnsrch,1
trnopt,full
nropt,full,,off
allsel
*do,im,0,12,1
tm=dt1+im
time,tm
ldread,temp,,,tm,0,al_fe_weld,rth
solve
*enddo
finish

wenti.JPG

回复此楼

倚楼听风雨，淡看江湖路

1楼 2015-01-15 16:02:26

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

维度lavery

新虫 (初入文坛)

应助: 0 (幼儿园)
帖子: 1
在线: 15分钟
虫号: 4529384
注册: 2016-03-22

引用回帖:

5楼: Originally posted by ymh0427 at 2015-01-19 10:37:07
多谢，这个问题解决了。但是又出问题了
Temperatures for material properties must be in ascending order. The
MPTEMP command is ignored.
我把命令流改为这样了
mptemp,,,,,,,,
mptemp,1,25
mpt ...

楼主有这个完整的命令流马，可运行的，小弟刚入门，想学习一下