引用回帖: 7樓 : Originally posted by
htbbzzg at 2023-02-22 08:41:29
你好!
簡單看了一下你的命令流,先說一個問題吧:
你所定義的與材料有關(guān)的命令有一個共同的問題,即:每個定義材料性能隨溫度變化的命令中的溫度參數(shù)最多只能是6 個,而你在每個材料性能隨溫度變化的命令中定 ...
老師,您好,按您的指導,溫度云圖的數(shù)據(jù)好看了一些,我后來也修改了一下網(wǎng)格劃分,但是算下來的兩層焊縫的溫度還是有差距,第一層是2200度,第二層是1450度,因為我劃分的兩層焊縫的截面面積是一樣的,不應該會有這么大的差距,而且論文里的兩層焊縫得溫度都是1880度,應該還是我哪里搞錯了,請老師幫我看看。
下面是修改后的命令流
!******************************第一部分進行溫度場分析******************************
/com, units: SI, m - degree C - W - kg -sec - J
!/title,the drive force of solidification crack
/prep7
et,1,solid70
!******************************定義熱分析物理參數(shù)************************************
mptemp,1,0,250,500,750,1000,1500
mptemp,7,1700,2500 !定義材料溫度區(qū)間(C)
mpdata,dens,1,1,7850,7700,7600,7550,7500,7350
mpdata,dens,1,7,7300,7100 !定義材料密度(kg/m3)
mpdata,kxx,1,1,50,45,40,30,26,35
mpdata,kxx,1,7,140,145 !熱傳導系數(shù)(W/(m*K))(考慮熔池對流)
mpdata,c,1,1,458,482,536,677,671,661
mpdata,c,1,7,781,825 !比熱容(Jp/K)
mpdata,enth,1,1,0,250*7700*482,500*7600*536,750*7550*677,1000*7500*671,1500*7350*661
mpdata,enth,1,7,1700*7300*781,2500*7100*825
*dim,conve,table,8,1,1,temp !考慮輻射和對流的總換熱系數(shù)(W/m^2*K),加載時按對流加載 ,定義8行、1列、1面的數(shù)組
conve(0,1,1)=1
conve(1,0,1)=0
conve(1,1,1)=5
conve(2,0,1)=250
conve(2,1,1)=30
conve(3,0,1)=500
conve(3,1,1)=50
conve(4,0,1)=750
conve(4,1,1)=97
conve(5,0,1)=1000
conve(5,1,1)=150
conve(6,0,1)=1500
conve(6,1,1)=300
conve(7,0,1)=1700
conve(7,1,1)=307
conve(8,0,1)=2500
conve(8,1,1)=351
!******************************定義應力分析參數(shù)***************************************
mpdata,ex,1,1,20E10, 17.8E10 ,13.5E10 ,11E10 ,0.7E10 ,0.5E10
mpdata,ex,1,7,0.5E10,0.1E10 !定義彈性模量
tb,bkin,1,6 !定義屈服強度和切變模量
tbtemp,20,1
tbdata,1,235e6,2.06e10
tbtemp,100,2
tbdata,1,220e6,2.03e10
tbtemp,500,3
tbdata,1,160e6,1.47e10
tbtemp,800,4
tbdata,1,75e6,0.198e10
tbtemp,1000,5
tbdata,1,12e6,0.07e10
tbtemp,1500,6
tbdata,1,1e6,0.07e10
mpdata,alpx,1,1,1.1e-5,1.25e-5,1.4e-5,1.3e-5,1.3e-5,1.3e-5
mpdata,alpx,1,7,1.3e-5,1.3e-5 !定義線膨脹系數(shù)(考慮凝固收縮中的固液區(qū))
mpdata,prxy,1,1,0.29,0.29,0.39,0.39
mpdata,prxy,1,7,0.39,0.39 !定義泊松比
!******************************建模型***********************************************
K, 1,0.0875,0,0
K, 2,0.0875,0.016,0
K, 3,0.0045,0.096,0
K, 4,-0.0045,0.096,0
K, 5,-0.0875,0.016,0
K, 6,-0.0875,0,0
K, 7, 0.0045, 0.024,0
K, 8,-0.0045, 0.016,0
K, 9,-0.0025, 0.016,0
K,10,0.0115, 0.016,0
k,11,0.0075, 0.016,0
k,12,0.0025, 0.152/7,0
K,13,-0.0045, 0,0
K,14,0.0115, 0,0
K,15,-0.0045, 0.024,0
K,16,-0.006, 0.016,0
K,17,-0.006, 0,0
K,18,0.013, 0.016,0
K,19,0.013, 0,0
A,1,2,18,19
A,3,4,15,7
A,5,6,17,16
A,8,13,14,10
A,9,11,12
A,11,10,7,12
A,10,14,19,18
A,8,9,7,15
A,16,17,13,8
VEXT, 1, , ,,,0.2
VEXT, 2, , ,,,0.2
VEXT, 3, , ,,,0.2
VEXT, 4, , ,,,0.2
VEXT, 5, , ,,,0.2
VEXT, 6, , ,,,0.2
VEXT, 7, , ,,,0.2
VEXT, 8, , ,,,0.2
VEXT, 9, , ,,,0.2
/facet,norml
vglue,all !粘連所有體
!*******************************網(wǎng)格劃分********
!******************體10和體5
vsel,s,,,10
vsel,a,,,5
esize,0.0025
mshkey,1
mshape,0
mat,1
type,1
vmesh,all
!******************體1********體2********體3
alls
vsel,s,,,1
vsel,a,,,2
vsel,a,,,3
esize,0.01
mshkey,1
mshape,0
mat,1
type,1
vmesh,all
!******************體14
alls
vsel,s,,,14
esize,0.005
mshkey,0
mshape,1
mat,1
type,1
vmesh,all
!******************體11****體12****體13
alls
vsel,s,,,11
vsel,a,,,12
vsel,a,,,13
esize,0.0025
mshkey,0
mshape,1
mat,1
type,1
vmesh,all
alls
!**************************************合并重合點*************************************
nummrg,all
alls
!************************************非線性求解選項****************************************
/solu
antype,4
trnopt,full
pred,on
solcontrol,on
nropt,full,,on
cutcontrol,plslimit,10
lnsrch,on
autots,on
tintp,,,,1,,
timint,on
tref,30
kbc,0
neqit,50
outres,all
!***************************************死亡************************************************
!/solu
!/auto,1
!/replot
alls
eplot
vsel,s,,,10
vsel,a,,,5
eslv,s,1
eplot
ekill,all
alls
esel,s,live
eplot
!************************施加對流邊界條件******************************
esel,s,live
ASEL,S,ext
NSLA,S,1
SF,ALL,CONV,%CONVE%,30
ALLS
!****************************************定義變量*******************************************
!HGEN=Q/A*V*DT!!
!A為焊縫的橫截面積
!V為焊接速度
!DT為每個載荷步的時間
!Q=K*U*I*DT!!
!K為焊接熱源的熱效率
!U為焊接電壓
!I為焊接電流
!*******************************************穩(wěn)態(tài)分析*********************************************
timint,off
alls
ic,all,temp,30
alls
nsla,u,1
sf,all,conv,%conve%,30
time,1e-5
deltim,1e-5,1e-5,1e-5
kbc,1
allsel,all
solve
eplot
!********************************************第一道焊縫瞬態(tài)分析*******************************************
j=0.5 !
zmin=-1E-5
zmax=0.005+1E-5 !
zd=0.005
!***************************************計算生熱率******************************************
K=0.75
U=31
I=270
A=28e-6
V=0.01
DT=0.5 !
HGENUM=(K*U*I)/(A*V*DT)
timint,on
*do,i,1,40 !
/solu
vsel,s,,,5
eslv,s,1
nsle,s,all
nsel,r,loc,z,zmin,zmax
esln,s,all
ealive,all
eplot
bfe,all,hgen,,HGENUM
time,j
deltim,0.5,0.5,0.5
alls
solve
vsel,s,,,5
eslv,s,1
nsle,s,all
nsel,r,loc,z,zmin,zmax
esln,s,all
eplot
bfedele,all,hgen
j=j+0.5
zmin=zmin+zd
zmax=zmax+zd
*enddo
alls
eplot
!*************************************第一道焊后冷卻過程*************************************************
*do,i,41,50 !
/solu
time,j
deltim,1,1,1
alls
solve
j=j+1
*enddo
!********************************************第二道焊縫瞬態(tài)分析*******************************************
j=j+0.5 !
zmin1=-1E-5
zmax1=0.0025+1E-5 !
zd1=0.0025 !
!***************************************計算生熱率******************************************
A1=28e-6
V1=0.005
DT1=0.5
HGENUM1=(K*U*I)/(A1*V1*DT1)
timint,on
*do,i,51,130 !
/solu
vsel,s,,,10
eslv,s,1
nsle,s,all
nsel,r,loc,z,zmin1,zmax1
esln,s,all
ealive,all
eplot
bfe,all,hgen,,HGENUM1
time,j
deltim,0.5,0.5,0.5 !
alls
solve
vsel,s,,,10
eslv,s,1
nsle,s,all
nsel,r,loc,z,zmin1,zmax1
esln,s,all
eplot
bfedele,all,hgen
j=j+0.5
zmin1=zmin1+zd1
zmax1=zmax1+zd1
*enddo
alls
eplot
!*************************************第二道焊后冷卻過程*************************************************
*do,i,131,230
/solu
time,j
deltim,10,10,10
alls
solve
save
j=j+10
*enddo
發(fā)自小木蟲Android客戶端