24小时热门版块排行榜    

查看: 3552  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ymh0427

木虫 (小有名气)

[求助] 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

ansys 焊接应力场问题
wenti.JPG
回复此楼

» 收录本帖的淘帖专辑推荐

焊接有限元分析 ANSYS命令流

» 猜你喜欢

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

倚楼听风雨,淡看江湖路
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

维度lavery

新虫 (初入文坛)

引用回帖:
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 ...

楼主有这个完整的命令流马,可运行的,小弟刚入门,想学习一下
6楼2016-03-29 11:58:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ymh0427 的主题更新
信息提示
请填处理意见