平板焊接之后极限拉伸载荷的计算问题，求教
本帖最后由 ouyangyu3 于 2011-8-15 20:27 编辑FINISH
/CLEAR
/PREP7
/TITLE,YLJS
ET,1,SOLID70
MPTEMP,1,25,500,1000,1500
MPDATA,EX,1,1,1.96E11,1.4275E11,0.7575E11,0.0425E11
MPDATA,EX,2,1,1.96E11,1.4275E11,0.7575E11,0.0425E11
TB,MISO,1,4
TBTEMP,25
TBPT,,0.004,7.84E+8
TBPT,,0.0045,8.685E+8
TBPT,,0.005,9.0E+8
TBPT,,0.009,9.18E+8
TBPT,,0.1,9.18E+8
TBTEMP,500
TBPT,,0.004,5.71E+8
TBPT,,0.005,5.75E+8
TBPT,,0.009,5.88E+8
TBPT,,0.01,5.88E+8
TBTEMP,1000
TBPT,,0.004,3.03E+8
TBPT,,0.005,3.06E+8
TBPT,,0.009,3.13E+8
TBPT,,0.01,3.13E+8
TBTEMP,1500
TBPT,,0.004,1.7E+007
TBPT,,0.005,1.7E+007
TBPT,,0.009,1.8E+007
TBPT,,0.01,1.8E+007
TB,MISO,2,4,
TBTEMP,25
TBPT,,0.004,7.84E+8
TBPT,,0.0045,8.685E+8
TBPT,,0.005,9.0E+8
TBPT,,0.009,9.18E+8
TBPT,,0.1,9.18E+8
TBTEMP,500
TBPT,,0.004,5.71E+8
TBPT,,0.005,5.75E+8
TBPT,,0.009,5.88E+8
TBPT,,0.01,5.88E+8
TBTEMP,1000
TBPT,,0.004,3.03E+8
TBPT,,0.005,3.06E+8
TBPT,,0.009,3.13E+8
TBPT,,0.01,3.13E+8
TBTEMP,1500
TBPT,,0.004,1.7E+007
TBPT,,0.005,1.7E+007
TBPT,,0.009,1.8E+007
TBPT,,0.01,1.8E+007
MP,PRXY,1,0.3
MP,PRXY,2,0.3
MP,ALPX,1,1.06e-5
MP,ALPX,2,1.06e-5
MP,KXX,1,16.3
MP,KXX,2,16.3
MP,C,1,502
MP,C,2,502
MP,MURX,1,1
MP,MURX,2,1
MP,DENS,1,7800
MP,DENS,2,7800
UIMP,1,REFT,25
UIMP,2,REFT,1300
hflengh=0.003
leftlengh=0.4
rightlengh=0.4
width=0.25
depth=0.008
k,1,0,0,0
k,2,leftlengh,0,0
k,3,leftlengh,width,0
k,4,0,width,0
k,5,0,0,depth
k,6,leftlengh,0,depth
k,7,leftlengh,width,depth
k,8,0,width,depth
k,9,leftlengh+hflengh,0,0
k,10,leftlengh+hflengh+rightlengh,0,0
k,11,leftlengh+hflengh+rightlengh,width,0
k,12,leftlengh+hflengh,width,0
k,13,leftlengh+hflengh,0,depth
k,14,leftlengh+hflengh+rightlengh,0,depth
k,15,leftlengh+hflengh+rightlengh,width,depth
k,16,leftlengh+hflengh,width,depth
V,5,1,4,8,6,2,3,7
V,6,2,3,7,13,9,12,16
V,13,9,12,16,14,10,11,15
VGLUE,ALL
NUMCMP,ALL
LESIZE,4,0.005,,,
LESIZE,11,0.005,,,
LESIZE,3,0.004,,,
TYPE,1
MAT,1
MSHAPE,0,3D
MSHKEY,1
VMESH,1
LESIZE,28,0.005,,,
LESIZE,27,0.005,,,
LESIZE,26,0.004,,,
TYPE,1
MAT,1
MSHAPE,0,3D
MSHKEY,1
VMESH,3
LESIZE,19,0.003,,,
TYPE,1
MAT,2
MSHAPE,0,3D
MSHKEY,1
VMESH,2
NUMCMP,ALL
EPLOT
SFA,1,1,CONV,62.5,25
SFA,2,1,CONV,62.5,25
SFA,3,1,CONV,62.5,25
SFA,4,1,CONV,62.5,25
SFA,5,1,CONV,62.5,25
SFA,14,1,CONV,62.5,25
SFA,16,1,CONV,62.5,25
SFA,12,1,CONV,62.5,25
SFA,13,1,CONV,62.5,25
SFA,15,1,CONV,62.5,25
OUTRES,ALL,ALL
VSEL,S,,,2
ESLV,S
CM,TOTAL,ELEM
CMSEL,S,TOTAL
EKILL,ALL
ESEL,S,LIVE
EPLOT
MXT=1300
DT09=2.399
DT01=0.001
T=0
tt=0
/SOL
ANTYPE,4
TRNOPT,FULL
LUMPM,0
ESEL,S,LIVE
OUTRES,ALL,ALL
IC,ALL,TEMP,25
ALLSEL,ALL
*DO,I,1,50
ESEL,S,,,16000+I*2
EALIVE,ALL
ESEL,S,LIVE
EPLOT
T=T+DT09
TIME,T
NSUBST,1
KBC=1
ESEL,S,,,16000+I*2
NSLE,S,ALL
D,ALL,TEMP,MXT
ALLSEL,ALL
SOLVE
T=T+DT01
TIME,T
NSUBST,1
KBC=1
ESEL,S,,,16000+I*2
NSLE,S,ALL
DDELE,ALL,TEMP,MAXT
ALLSEL,ALL
SOLVE
*ENDDO
*DO,I,1,50
ESEL,S,,,15999+I*2
EALIVE,ALL
ESEL,S,LIVE
EPLOT
T=T+DT09
TIME,T
NSUBST,1
KBC=1
ESEL,S,,,15999+I*2
NSLE,S,ALL
D,ALL,TEMP,MXT
ALLSEL,ALL
SOLVE
T=T+DT01
TIME,T
NSUBST,1
KBC=1
ESEL,S,,,15999+I*2
NSLE,S,ALL
DDELE,ALL,TEMP,MAXT
ALLSEL,ALL
SOLVE
*ENDDO
TIME,1200
DELTIM,20
SOLVE
TIME,3600
DELTIM,40
SOLVE
TIME,7200
DELTIM,90
SOLVE
SAVE,'TEMP','db','F:\ANSYS\REAL'
/PREP7
ETCHG,TTS
EPLOT
ALLSEL,ALL
DL,1, ,ALL
DL,3, ,ALL
DL,22, ,ALL
DL,26, ,ALL
DT24=2.4
DT20=20
DT40=40
DT90=90
TT=0
/SOL
ANTYPE,0
ISWRITE,ON
ALLSEL,ALL
EPLOT
CMSEL,S,TOTAL
EKILL,ALL
ESEL,S,LIVE
EPLOT
*DO,I,1,50
ESEL,S,,,16000+I*2
EALIVE,ALL
ESEL,S,LIVE
EPLOT
TT=TT+DT24
ESEL,ALL
TIME,TT
NSUBST,1
LDREAD,TEMP,,, TT, ,'temp','rth','F:\ANSYS\'
SOLVE
*ENDDO
*DO,I,1,50
ESEL,S,,,15999+I*2
EALIVE,ALL
ESEL,S,LIVE
EPLOT
TT=TT+DT24
ESEL,ALL
TIME,TT
NSUBST,1
LDREAD,TEMP,,, TT, ,'temp','rth','F:\ANSYS\'
SOLVE
*ENDDO
*DO,I,1,60
TT=TT+DT20
TIME,TT
NSUBST,1
LDREAD,TEMP,,, TT, ,'temp','rth','F:\ANSYS\'
SOLVE
*ENDDO
*DO,I,1,90
TT=TT+DT40
TIME,TT
NSUBST,1
LDREAD,TEMP,,, TT, ,'temp','rth','F:\ANSYS\'
SOLVE
*ENDDO
*DO,I,1,80
TT=TT+DT90
TIME,TT
NSUBST,1
LDREAD,TEMP,,, TT, ,'temp','rth','F:\ANSYS\'
SOLVE
*ENDDO
SAVE,'result2','db','F:\ansys\real\'
*************************************************************极限强度计算
/SOLU
ANTYPE,0
ALLSEL,ALL
DDELE,ALL,ALL
alls
nsel,s,loc,X,0
d,all,all
alls
nsel,s,loc,X,0.803
cp,1,ux,all
*get,nodenum,node,0,count
f,all,fX,5000000/nodenum
D,ALL,UY,,,,,UZ
alls
RESCONTROL,DEFINE
NLGEOM,on
arclen,on
NSUBST,100
outres,all,all
SOLVE在平板对接焊接，生死单元模拟之后，运用RESCONTROL,进行重启动分析，计算焊接残余应力时候的边界条件，删除与不删除这两种情况下，计算的极限拉伸载荷。都很大，跟在直接计算极限拉伸载荷（是1840KN），相差了几倍的关系，想请教下高手，这个如何解决，是我的命令错了么？
如果运用初始应力施加，ISWRITE输出应力文件，然后ISFILE读入应力值加入模型中计算的极限载荷值跟1840KN接近，但是运用LIST看了下施加的应力场，跟焊接得到的残余应力场有出入，这是怎么回事？计算时候说的是：由于材料定义的是非线性，但是应力场加在了弹性材料上，这个是怎么回事？是不是计算的值不准确了？帮忙解答下。谢谢各位了！ 自己先顶下，各位大侠，多说说自己的想法啊，小弟叩谢了啊！急着呢
页:
[1]