24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3492  |  回复: 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的回帖

fangdian

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
ymh0427: 金币+2, ★★★很有帮助 2015-01-15 20:01:59
应该是如下2个原因:
   1)边界条件或者材料属性设置问题。超过了应有的弹性应变速率
   2)时间步太大,建议减小时间步
2楼2015-01-15 18:25:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ymh0427

木虫 (小有名气)

引用回帖:
2楼: Originally posted by fangdian at 2015-01-15 18:25:50
应该是如下2个原因:
   1)边界条件或者材料属性设置问题。超过了应有的弹性应变速率
   2)时间步太大,建议减小时间步

减小时间步长我试过了,我改为*do,im,0,12,0.001都不行。
边界条件或材料属性,不知道怎么改
倚楼听风雨,淡看江湖路
3楼2015-01-15 20:01:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿七不笨

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
ymh0427: 金币+8, ★★★很有帮助 2015-01-19 10:37:16
试一下用CUTCONTROL改变塑形应变增率
CAE
4楼2015-01-16 13:53:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ymh0427

木虫 (小有名气)

引用回帖:
4楼: Originally posted by 阿七不笨 at 2015-01-16 13:53:49
试一下用CUTCONTROL改变塑形应变增率

多谢,这个问题解决了。但是又出问题了
Temperatures for material properties must be in ascending order.  The   
MPTEMP command is ignored.  
我把命令流改为这样了
mptemp,,,,,,,,
mptemp,1,25
mptemp,2,250
mptemp,3,500
mptemp,4,750
mptemp,5,1000
mptemp,6,1500
mptemp,7,1700
mptemp,8,2500
mpdata,kxx,1,,50
mpdata,kxx,1,,47
mpdata,kxx,1,,40
mpdata,kxx,1,,27
mpdata,kxx,1,,30
mpdata,kxx,1,,35
mpdata,kxx,1,,14
mpdata,kxx,1,,14.2
mpdata,c,1,,460
mpdata,c,1,,480
mpdata,c,1,,530
mpdata,c,1,,675
mpdata,c,1,,670
mpdata,c,1,,660
mpdata,c,1,,780
mpdata,c,1,,820
mp,dens,1,7800
mp,alpx,1,1.78e-5
mptemp,,,,,,,,
mptemp,1,25
mptemp,2,100
mptemp,3,200
mptemp,4,300
mptemp,5,400
mptemp,6,500
mptemp,7,600
mpdata,kxx,2,,176
mpdata,kxx,2,,180
mpdata,kxx,2,,184
mpdata,kxx,2,,188
mpdata,kxx,2,,192
mpdata,kxx,2,,198
mpdata,kxx,2,,200
mpdata,c,2,,728
mpdata,c,2,,795
mpdata,c,2,,879
mpdata,c,2,,963
mpdata,c,2,,1089
mpdata,c,2,,1290
mpdata,c,2,,1560
mp,dens,2,2700
mp,alpx,2,9.36e-6
倚楼听风雨,淡看江湖路
5楼2015-01-19 10:37:07
已阅   回复此楼   关注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的回帖

害玩逗

新虫 (初入文坛)

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

你好 请问你第一个问题是怎么解决的
7楼2018-11-26 16:30:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ymh0427 的主题更新
信息提示
请填处理意见