24小时热门版块排行榜    

查看: 4173  |  回复: 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的回帖

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的回帖
查看全部 7 个回答

fangdian

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
ymh0427: 金币+2, ★★★很有帮助 2015-01-15 20:01:59
应该是如下2个原因:
   1)边界条件或者材料属性设置问题。超过了应有的弹性应变速率
   2)时间步太大,建议减小时间步
2楼2015-01-15 18:25:50
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 2rKfVQg2jQww 2026-11-06 14/700 2026-11-09 04:20 by l48hFd8Nwae2
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 11/550 2026-11-09 04:09 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +10 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 17:18 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-06 13:47 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:36 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-05 10/500 2026-11-06 12:32 by BkA9Bi5cMhkD
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
信息提示
请填处理意见