24小时热门版块排行榜    

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

阿七不笨

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
ymh0427: 金币+8, ★★★很有帮助 2015-01-19 10:37:16
试一下用CUTCONTROL改变塑形应变增率
CAE
4楼2015-01-16 13:53:49
已阅   回复此楼   关注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的回帖

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0805 316求调剂 +6 大雪深藏 2026-03-18 6/300 2026-03-24 22:13 by peike
[考研] 311求调剂 +3 冬十三 2026-03-24 3/150 2026-03-24 21:31 by peike
[考研] 求调剂323材料与化工 +3 1124361 2026-03-24 3/150 2026-03-24 20:24 by peike
[考研] 289求调剂 +8 硕星赴 2026-03-23 8/400 2026-03-24 20:17 by peike
[考研] 调剂 +4 13853210211 2026-03-24 4/200 2026-03-24 19:44 by ms629
[考研] 306求0703调剂一志愿华中师范 +10 纸鱼ly 2026-03-21 11/550 2026-03-24 17:22 by qingfeng258
[考研] 291求调剂 +3 HanBeiNingZC 2026-03-24 3/150 2026-03-24 16:34 by barlinike
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +4 曼111 2026-03-24 5/250 2026-03-24 15:44 by 星空星月
[考研] 一志愿华东理工大学081700,初试分数271 +5 kotoko_ik 2026-03-23 6/300 2026-03-24 10:29 by 学术搬砖er
[考研] 求调剂 +7 十三加油 2026-03-21 7/350 2026-03-23 23:48 by 热情沙漠
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +7 vv迷 2026-03-22 7/350 2026-03-23 23:44 by Txy@872106
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
[考研] 293求调剂 +3 涛涛Wjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-20 3/150 2026-03-22 16:00 by ColorlessPI
[考研] 319求调剂 +4 小力气珂珂 2026-03-20 4/200 2026-03-22 15:53 by ColorlessPI
[考研] 材料与化工(0856)304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
信息提示
请填处理意见