热应力间接耦合与直接耦合

happygirl168505

为什么这个分析的热应力直接耦合与间接耦合的应力场结果不一样呢？命令流分别如下，哪点出错了呢？求解释
直接耦合
/prep7
et,1,plane13
KEYOPT,1,1,4       
KEYOPT,1,3,1       
mptemp,1,34
mptemp,2,114
mptemp,3,794
mptemp,4,1034
mptemp,5,1109
mptemp,6,1113
mptemp,7,1121
mptemp,8,1213
mptemp,9,1337
mptemp,10,1409
mptemp,11,1489               
mpdata,kxx,1,,18.6
mpdata,kxx,1,,18.98
mpdata,kxx,1,,23
mpdata,kxx,1,,24.86
mpdata,kxx,1,,25.45
mpdata,kxx,1,,25.53
mpdata,kxx,1,,25.7
mpdata,kxx,1,,26.06
mpdata,kxx,1,,25.29
mpdata,kxx,1,,26.3
mpdata,kxx,1,,27.4
mpdata,enth,1,,123.8*7998e3
mpdata,enth,1,,160.1*7961e3
mpdata,enth,1,,521.3*7639e3
mpdata,enth,1,,666.9*7526e3
mpdata,enth,1,,723.1*7462e3
mpdata,enth,1,,729.8*7445.7e3
mpdata,enth,1,,744.9*7410.3e3
mpdata,enth,1,,891.96*7416e3
mpdata,enth,1,,1084*6927.3e3
mpdata,enth,1,,1137.6*6872.7e3
mpdata,enth,1,,1199.5*6811.7e3
mpdata,dens,1,,7998
mpdata,dens,1,,7961
mpdata,dens,1,,7639
mpdata,dens,1,,7526
mpdata,dens,1,,7462
mpdata,dens,1,,7445.7
mpdata,dens,1,,7410.3
mpdata,dens,1,,7416
mpdata,dens,1,,6927.3
mpdata,dens,1,,6872.7
mpdata,dens,1,,6811.7
mpplot,kxx,1
mpplot,enth,1
MPTEMP,,,,,,,,
mptemp,1,290
mptemp,2,610
mptemp,3,810
mptemp,4,910
mptemp,5,1010
mptemp,6,1110
mptemp,7,1150
mptemp,8,1185
mptemp,9,1191
mptemp,10,1205
mptemp,11,1251
mptemp,12,1271
mpdata,ex,1,,1.725e11
mpdata,ex,1,,1.283e11
mpdata,ex,1,,0.9161e11
mpdata,ex,1,,0.7224e11
mpdata,ex,1,,0.5322e11
mpdata,ex,1,,0.3537e11
mpdata,ex,1,,0.2843e11
mpdata,ex,1,,0.1769e11
mpdata,ex,1,,0.1237e11
mpdata,ex,1,,0.05638e11
mpdata,ex,1,,1.15e9
mpdata,ex,1,,4.031e8
mpdata,prxy,1,,0.2745
mpdata,prxy,1,,0.2745
mpdata,prxy,1,,0.2745
mpdata,prxy,1,,0.2745
mpdata,prxy,1,,0.2745
mpdata,prxy,1,,0.2745
mpdata,prxy,1,,0.2745
mpdata,prxy,1,,0.2877
mpdata,prxy,1,,0.3094
mpdata,prxy,1,,0.3534
mpdata,prxy,1,,0.4029
mpdata,prxy,1,,0.4277
MPTEMP,,,,,,,,
mptemp,1,320
mptemp,2,440
mptemp,3,530
mptemp,4,680
mptemp,5,710
mptemp,6,830
mptemp,7,1130
mptemp,8,1199
mptemp,9,1259
mptemp,10,1304
mpdata,alpx,1,,23.9e-6
mpdata,alpx,1,,27.1e-6
mpdata,alpx,1,,29.06e-6
mpdata,alpx,1,,23.77e-6
mpdata,alpx,1,,23.78e-6
mpdata,alpx,1,,24.01e-6
mpdata,alpx,1,,23.9e-6
mpdata,alpx,1,,20.63e-6
mpdata,alpx,1,,35.08e-6
mpdata,alpx,1,,32.61e-6
save

rectng,0.175,0.255,0,0.1
aplot
save


mat,1
esize,0.001
amesh,1



/SOLU
antype,trans                                       
esel,s,mat,,1                                       
nsle,s
/replot
ic,all,temp,1000
save
lplot
sfl,2,CONV,88,,25                               
sfl,4,CONV,135,,35
mp,reft,1,1000
!施加位移约束
asel,s,mat,,1
esel,s,mat,,1
nsle,s
d,all,uy,0
!施加惯性力约束
omega,0,120,0
/psf,conv,2
allsel
time,300
autots,on
deltim,10,10,30       
outres,all,all
solve

间接耦合
/prep7
et,1,plane55
KEYOPT,1,3,1               
mptemp,1,34
mptemp,2,114
mptemp,3,794
mptemp,4,1034
mptemp,5,1109
mptemp,6,1113
mptemp,7,1121
mptemp,8,1213
mptemp,9,1337
mptemp,10,1409
mptemp,11,1489               
mpdata,kxx,1,,18.6
mpdata,kxx,1,,18.98
mpdata,kxx,1,,23
mpdata,kxx,1,,24.86
mpdata,kxx,1,,25.45
mpdata,kxx,1,,25.53
mpdata,kxx,1,,25.7
mpdata,kxx,1,,26.06
mpdata,kxx,1,,25.29
mpdata,kxx,1,,26.3
mpdata,kxx,1,,27.4
mpdata,enth,1,,123.8*7998e3
mpdata,enth,1,,160.1*7961e3
mpdata,enth,1,,521.3*7639e3
mpdata,enth,1,,666.9*7526e3
mpdata,enth,1,,723.1*7462e3
mpdata,enth,1,,729.8*7445.7e3
mpdata,enth,1,,744.9*7410.3e3
mpdata,enth,1,,891.96*7416e3
mpdata,enth,1,,1084*6927.3e3
mpdata,enth,1,,1137.6*6872.7e3
mpdata,enth,1,,1199.5*6811.7e3
mpdata,dens,1,,7998
mpdata,dens,1,,7961
mpdata,dens,1,,7639
mpdata,dens,1,,7526
mpdata,dens,1,,7462
mpdata,dens,1,,7445.7
mpdata,dens,1,,7410.3
mpdata,dens,1,,7416
mpdata,dens,1,,6927.3
mpdata,dens,1,,6872.7
mpdata,dens,1,,6811.7
mpplot,kxx,1
mpplot,enth,1
MPTEMP,,,,,,,,
mptemp,1,290
mptemp,2,610
mptemp,3,810
mptemp,4,910
mptemp,5,1010
mptemp,6,1110
mptemp,7,1150
mptemp,8,1185
mptemp,9,1191
mptemp,10,1205
mptemp,11,1251
mptemp,12,1271
mpdata,ex,1,,1.725e11
mpdata,ex,1,,1.283e11
mpdata,ex,1,,0.9161e11
mpdata,ex,1,,0.7224e11
mpdata,ex,1,,0.5322e11
mpdata,ex,1,,0.3537e11
mpdata,ex,1,,0.2843e11
mpdata,ex,1,,0.1769e11
mpdata,ex,1,,0.1237e11
mpdata,ex,1,,0.05638e11
mpdata,ex,1,,1.15e9
mpdata,ex,1,,4.031e8
mpdata,prxy,1,,0.2745
mpdata,prxy,1,,0.2745
mpdata,prxy,1,,0.2745
mpdata,prxy,1,,0.2745
mpdata,prxy,1,,0.2745
mpdata,prxy,1,,0.2745
mpdata,prxy,1,,0.2745
mpdata,prxy,1,,0.2877
mpdata,prxy,1,,0.3094
mpdata,prxy,1,,0.3534
mpdata,prxy,1,,0.4029
mpdata,prxy,1,,0.4277
MPTEMP,,,,,,,,
mptemp,1,320
mptemp,2,440
mptemp,3,530
mptemp,4,680
mptemp,5,710
mptemp,6,830
mptemp,7,1130
mptemp,8,1199
mptemp,9,1259
mptemp,10,1304
mpdata,alpx,1,,23.9e-6
mpdata,alpx,1,,27.1e-6
mpdata,alpx,1,,29.06e-6
mpdata,alpx,1,,23.77e-6
mpdata,alpx,1,,23.78e-6
mpdata,alpx,1,,24.01e-6
mpdata,alpx,1,,23.9e-6
mpdata,alpx,1,,20.63e-6
mpdata,alpx,1,,35.08e-6
mpdata,alpx,1,,32.61e-6
save

rectng,0.175,0.255,0,0.1
aplot
save


mat,1
esize,0.001
amesh,1




/SOLU
antype,trans                                       
esel,s,mat,,1                                       
nsle,s
/replot
ic,all,temp,1000
save
lplot
sfl,2,CONV,88,,25                               
sfl,4,CONV,135,,35

timint,on
time,300
autots,on
deltim,10,10,30       
outres,all,all
solve


/prep7
etchg,tts
keyopt,1,3,1
keyopt,1,6,1
mp,reft,1,1000
!施加位移约束
asel,s,mat,,1
esel,s,mat,,1
nsle,s
d,all,uy,0
!施加惯性力约束
omega,0,120,0


/solu
antype,trans
timint,on
*do,i,10,300,10
time,i
autots,on
deltim,10,10,10
outres,all,all

ldread,temp,,,i,,,rth
solve
*enddo

happygirl168505

自己顶一个，哪位大侠解释一下呢？

leiguilin

直接耦合和间接耦合两者分析方法不一样，结果肯定有误差啊，要看你的结果误差多大了！

happygirl168505

leiguilin 发表于 2012-3-3 21:25
直接耦合和间接耦合两者分析方法不一样，结果肯定有误差啊，要看你的结果误差多大了！ ...

现在的问题是应力场结果差别很大啊？温度场的差别还能接受，应力场趋势都不一样啊，怎么回事呢？

leiguilin

happygirl168505 发表于 2012-3-5 20:27
现在的问题是应力场结果差别很大啊？温度场的差别还能接受，应力场趋势都不一样啊，怎么回事呢？ ...

温度场相差很小吗？温度分布是不是基本一样？

happygirl168505

leiguilin 发表于 2012-3-7 23:57
温度场相差很小吗？温度分布是不是基本一样？

温度场差一点点，应力场差别很大

leiguilin

happygirl168505 发表于 2012-3-8 09:31
温度场差一点点，应力场差别很大

如何约束都一样的话，那就不知道为什么了，建议重新试试看！有空了我翻翻书看看！