【原创】自己做的间接法热-结构耦合的例子，大家看看！

dragon_2002_sd

我做的间接法热-结构耦合的命令流。温度场分析后，转换单元，进行结构分析。读入温度为体载荷。设为瞬态分析。l两个不同的时间（1秒和3e-5秒）结果截然不同！！！恳请各位大虾运行分析一下！
/prep7
et,1,plane55,1
et,2,surf151,
keyopt,2,4,1
keyopt,2,8,3   
mptemp,1,0,400,600,800,1000,2000
mpdata,kxx,1,1,18.9,28.35,34.65,34.65,29.4,3.15
mpdata,c,1,1,1839,7559,9500,10200,8433,4821.6
mpdata,enth,1,1,0,6.42e9,19.3e9,34.17e9,41.64e9,51.03e9
mpdata,hf,1,1,9.66,59.43,84.315,109.2,180.6,537.6
mpdata,dens,1,1,7550,7550,7550,7550,7550,7550
!input material property
  
rectng,-0.1,0.1,0,0.03
wpave,-0.01
wprota,,,90
asbw,all
wpave,0.01
asbw,all
wpave,0.016
asbw,all
wpave,-0.016
asbw,all
wpave,-0.04
asbw,all
wpave,0.04
asbw,all
wpave,0.06
asbw,all
wpave,-0.06
asbw,all
wprota,,90
wpave,-0.1,0.02
asbw,all
wpave,-0.1,0.01
asbw,all
/pnum,area,1
aglue,all
wpstyle,,,,,,,,0
!creat model
  
lsel,s,,,12
lsel,s,,,41
lsel,s,,,12
lsel,s,,,25
lsel,s,,,12
lsel,s,,,11
lesize,all,0.001
allsel
lsel,s,,,18
lsel,a,,,48
lsel,a,,,60
lsel,a,,,17
lsel,a,,,15
lsel,a,,,45
lsel,a,,,58
lsel,a,,,14
lesize,all,0.0015
allsel
lsel,s,,,21
lsel,a,,,51
lsel,a,,,62
lsel,a,,,20
lsel,a,,,24
lsel,a,,,54
lsel,a,,,64
lsel,a,,,23
lplot
lesize,all,0.002
allsel
lsel,s,,,30
lsel,a,,,57
lsel,a,,,66
lsel,a,,,29
lsel,a,,,27
lsel,a,,,56
lsel,a,,,65
lsel,a,,,26
lplot
lesize,all,0.004
allsel
lsel,s,,,3
lsel,a,,,31
lsel,a,,,7
lsel,a,,,1
lsel,a,,,6
lsel,a,,,36
lsel,a,,,16
lsel,a,,,5
lplot
lesize,all,0.005
allsel
lplot
  
lsel,s,,,32
lsel,a,,,33
lsel,a,,,52
lsel,a,,,49
lsel,a,,,42
lsel,a,,,43
lsel,a,,,46
lsel,a,,,55
lsel,a,,,38
lsel,a,,,37
lplot
lesize,all,0.001
allsel
lplot
lsel,s,,,67,76
lplot
lesize,all,0.002
allsel
lplot
lsel,s,,,2
lsel,a,,,4
lsel,a,,,61
lsel,a,,,59
lsel,a,,,19
lsel,a,,,22
lsel,a,,,28
lsel,a,,,63
lsel,a,,,13
lsel,a,,,10
lplot
lesize,all,0.003
allsel
type,1
mat,1
mshkey,1
amesh,all
!creat the model
  
lsel,s,loc,y,0
lsel,a,loc,y,0.03
lsel,a,loc,x,-0.1
lsel,a,loc,x,0.1
lplot
type,2
mat,1
lmesh,all
allsel
eplot
/pnum,line,off
  
/solu
tunif,20
tref,0
esel,s,type,,2
sfe,all,,conv,,-1,20
allsel
eplot
  
*DIM,hflux,TABLE,5,1,1,x, ,  
HFLUX(1,0,1) = -0.016
HFLUX(1,1,1) = 800000
HFLUX(2,0,1) = -0.01  
HFLUX(2,1,1) = 6500000   
HFLUX(3,0,1) = 0  
HFLUX(3,1,1) = 16000000   
HFLUX(4,0,1) = 0.01   
HFLUX(4,1,1) = 6500000   
HFLUX(5,0,1) = 0.016  
HFLUX(5,1,1) = 800000
lsel,s,,,18
lsel,a,,,12
lsel,a,,,15
nsll,s,1
sf,all,hflux,%hflux%
allsel
!施加热流密度
  
antype,trans
time,25                  !焊接时间25s
autots,-1
deltim,0.1,,,1
kbc,1
TIMINT,1
TINTP,0.005, , ,1,0.5,0.2  
lnsrch,1
outres,all,all
lswrite,1
time,40
SFDELE,all,HFLUX  
AUTOTS,1
DELTIM,0.2,0.1,0.3,1
KBC,1  
outres,all,all
lswrite,2
time,100
AUTOTS,1
DELTIM,0.5,0.1,1,1   
KBC,1
outres,all,all
lswrite,3
lssolve,1,3,1
!焊接过程25s
  
/prep7
etchg,tts
mpdata,ex,1,1,1.72e11,1.31e11,1.1e11,0.364e11,0.15e11,0.05e11
mpdata,prxy,1,1,0.26,0.26,0.26,0.26,0.26,0.26
mpdata,alpx,1,1,0,11.2e-6,13.2e-6,15.2e-6,17.2e-6,27.2e-6
tb,bkin,1,6
tbtemp,0,1
tbdata,1,4.6e8,1.72e10
tbtemp,400,2
tbdata,1,4.07e7,1.31e10
tbtemp,600,3
tbdata,1,3.8e8,1.1e10
tbtemp,800,4
tbdata,1,0.45e8,0.364e10
tbtemp,1000,5
tbdata,1,0.36e8,0.15e10
tbtemp,2000,6
tbdata,1,0.35e8,0.05e10
  
/solu
LDREAD,TEMP,,,100, ,'welding','rth',' '!读入温度100s时的温度  
nsel,s,loc,y,0
nsel,r,loc,x,0
d,all,uy,0
  
antype,trans
time,3e-5 !!!!!!!注意时间
AUTOTS,-1   
DELTIM,1e-6,0.5e-6,1.5e-6,1
KBC,0
OUTRES,ALL,ALL,
LNSRCH,1
NROPT,FULL, ,OFF
solve

dragon_2002_sd

不好意思，上面的命令流倒数第5行“KBC,0”改为“KBC,1“。
这是应力图：

dragon_2002_sd

不能显示？

dragon_2002_sd

最后时间的应力图：

dragon_2002_sd

上表面应力分布图：

dragon_2002_sd

最后时刻上表面应力分布图：

dragon_2002_sd

下面的命令流是瞬态求解应力时间1s:
/prep7
et,1,plane55,1
et,2,surf151,
keyopt,2,4,1
keyopt,2,8,3   
mptemp,1,0,400,600,800,1000,2000
mpdata,kxx,1,1,18.9,28.35,34.65,34.65,29.4,3.15
mpdata,c,1,1,1839,7559,9500,10200,8433,4821.6
mpdata,enth,1,1,0,6.42e9,19.3e9,34.17e9,41.64e9,51.03e9
mpdata,hf,1,1,9.66,59.43,84.315,109.2,180.6,537.6
mpdata,dens,1,1,7550,7550,7550,7550,7550,7550
!input material property
  
rectng,-0.1,0.1,0,0.03
wpave,-0.01
wprota,,,90
asbw,all
wpave,0.01
asbw,all
wpave,0.016
asbw,all
wpave,-0.016
asbw,all
wpave,-0.04
asbw,all
wpave,0.04
asbw,all
wpave,0.06
asbw,all
wpave,-0.06
asbw,all
wprota,,90
wpave,-0.1,0.02
asbw,all
wpave,-0.1,0.01
asbw,all
/pnum,area,1
aglue,all
wpstyle,,,,,,,,0
!creat model
  
lsel,s,,,12
lsel,s,,,41
lsel,s,,,12
lsel,s,,,25
lsel,s,,,12
lsel,s,,,11
lesize,all,0.001
allsel
lsel,s,,,18
lsel,a,,,48
lsel,a,,,60
lsel,a,,,17
lsel,a,,,15
lsel,a,,,45
lsel,a,,,58
lsel,a,,,14
lesize,all,0.0015
allsel
lsel,s,,,21
lsel,a,,,51
lsel,a,,,62
lsel,a,,,20
lsel,a,,,24
lsel,a,,,54
lsel,a,,,64
lsel,a,,,23
lplot
lesize,all,0.002
allsel
lsel,s,,,30
lsel,a,,,57
lsel,a,,,66
lsel,a,,,29
lsel,a,,,27
lsel,a,,,56
lsel,a,,,65
lsel,a,,,26
lplot
lesize,all,0.004
allsel
lsel,s,,,3
lsel,a,,,31
lsel,a,,,7
lsel,a,,,1
lsel,a,,,6
lsel,a,,,36
lsel,a,,,16
lsel,a,,,5
lplot
lesize,all,0.005
allsel
lplot
  
lsel,s,,,32
lsel,a,,,33
lsel,a,,,52
lsel,a,,,49
lsel,a,,,42
lsel,a,,,43
lsel,a,,,46
lsel,a,,,55
lsel,a,,,38
lsel,a,,,37
lplot
lesize,all,0.001
allsel
lplot
lsel,s,,,67,76
lplot
lesize,all,0.002
allsel
lplot
lsel,s,,,2
lsel,a,,,4
lsel,a,,,61
lsel,a,,,59
lsel,a,,,19
lsel,a,,,22
lsel,a,,,28
lsel,a,,,63
lsel,a,,,13
lsel,a,,,10
lplot
lesize,all,0.003
allsel
type,1
mat,1
mshkey,1
amesh,all
!creat element
  
lsel,s,loc,y,0
lsel,a,loc,y,0.03
lsel,a,loc,x,-0.1
lsel,a,loc,x,0.1
lplot
type,2
mat,1
lmesh,all
allsel
eplot
/pnum,line,off
  
/solu
tunif,20
tref,0
esel,s,type,,2
sfe,all,,conv,,-1,20
allsel
eplot
  
*DIM,hflux,TABLE,5,1,1,x, ,  
HFLUX(1,0,1) = -0.016
HFLUX(1,1,1) = 800000
HFLUX(2,0,1) = -0.01  
HFLUX(2,1,1) = 6500000   
HFLUX(3,0,1) = 0  
HFLUX(3,1,1) = 16000000   
HFLUX(4,0,1) = 0.01   
HFLUX(4,1,1) = 6500000   
HFLUX(5,0,1) = 0.016  
HFLUX(5,1,1) = 800000
lsel,s,,,18
lsel,a,,,12
lsel,a,,,15
nsll,s,1
sf,all,hflux,%hflux%
allsel
  
antype,trans
time,25                  !焊接时间25s
autots,-1
deltim,0.1,,,1
kbc,1
TIMINT,1
TINTP,0.005, , ,1,0.5,0.2  
lnsrch,1
outres,all,all
lswrite,1
time,40
SFDELE,all,HFLUX  
AUTOTS,1
DELTIM,0.2,0.1,0.3,1
KBC,1  
outres,all,all
lswrite,2
time,100
AUTOTS,1
DELTIM,0.5,0.1,1,1   
KBC,1
outres,all,all
lswrite,3
lssolve,1,3,1
!焊接过程25s
  
/prep7
etchg,tts
mpdata,ex,1,1,1.72e11,1.31e11,1.1e11,0.364e11,0.15e11,0.05e11
mpdata,prxy,1,1,0.26,0.26,0.26,0.26,0.26,0.26
mpdata,alpx,1,1,0,11.2e-6,13.2e-6,15.2e-6,17.2e-6,27.2e-6
tb,bkin,1,6
tbtemp,0,1
tbdata,1,4.6e8,1.72e10
tbtemp,400,2
tbdata,1,4.07e7,1.31e10
tbtemp,600,3
tbdata,1,3.8e8,1.1e10
tbtemp,800,4
tbdata,1,0.45e8,0.364e10
tbtemp,1000,5
tbdata,1,0.36e8,0.15e10
tbtemp,2000,6
tbdata,1,0.35e8,0.05e10
  
/solu
LDREAD,TEMP,,,100, ,'welding','rth',' '
nsel,s,loc,y,0
nsel,r,loc,x,0
d,all,uy,0
  
antype,trans
time,1 !!!!!!!注意时间
AUTOTS,-1   
DELTIM,0.1,0.05,0.2,1
KBC,1
OUTRES,ALL,ALL,
LNSRCH,1
NROPT,FULL, ,OFF
solve  
应力几乎分布图保持不变：

dragon_2002_sd

不好意思贴错了！！如图：

dragon_2002_sd

有贴错了！

dragon_2002_sd

上表面应力分布（始终不变）：

ainite

呵呵，做得不错，虽然没时间跑一下看看，还是要支持一下，请斑竹也给点鼓励呀

dragon_2002_sd

to:ainite大侠
谢谢！ainite 大虾，又时间还是帮忙分析一下吧！
个人认为，瞬态结构分析时间为1秒时，结果比较正确。这里再贴上钢板上表面应力分布图：
  

dragon_2002_sd

    难道大家都没有一点儿热情嘛？辛辛苦苦做的。真让人受打击。
    jinl姐姐，ainite大虾，难道你们也没有一点时间运行分析一下嘛？我这里还有个问题等着你们帮忙呢？？？为什么回答得到不同的结果呢？

jnjl

既然你说１秒的时候应力分布正确，你的命令流应该没什么问题，为什么要算３Ｅ－５时刻呢？
是不是在那么短的时刻内应力分布就应该是这样的呢？对这方面没研究，瞎说的啊，你可以去找找专业的书籍或者问问你老师，用专业的理论来解释你的结果。

dragon_2002_sd

     因为在3e-5秒内钢板要受到一冲击载荷。我想观察分析受冲击后的应力效果。所以先分析钢板在未受冲击时的应力分布。没想到会得到这个结果？
     我觉得，间接耦合法模拟热-结构问题，假如单纯的分析温度产生的应力分布，如焊接应力，在结构分析时，读入温度体载荷后，无论做瞬态分析还是静态分析，应力应该是不变化的或者变化极小，因为在分析的时间内，载荷没有发生改变。因此我感觉在3e-5秒内钢板的应力也是变化很小的，但是得到的结果却截然不同？
    不知我的想法是否正确？

nwpuyl

最近一直没有上来，很想念大家！
这个问题我也遇到过，讨论讨论吧

dragon_2002_sd wrote:
难道大家都没有一点儿热情嘛？辛辛苦苦做的。真让人受打击。  
    jinl姐姐，ainite大虾，难道你们也没有一点时间运行分析一下嘛？

楼主辛苦！
斑竹和ainite很忙，所以体谅体谅啦！
我觉得你的问题已经很明了了，你运行的结果也比较清楚，别人再运行结果也是一样的！
如果没有找到解决问题的正确办法，盲目运行只会浪费机时，对吧！对症下药才是首选嘛

dragon_2002_sd wrote:
     因为在3e-5秒内钢板要受到一冲击载荷。我想观察分析受冲击后的应力效果。所以先分析钢板在未受冲击时的应力分布。没想到会得到这个结果？
      

既然你只考虑结构力的冲击效应，那么在分析未受冲击时的应力分布时－－热应力，读入温度做稳态分析很合理嘛！而且省事省力
只是为了看热应力，就不用考虑结构力了呀，管它冲不冲击呢，干嘛选用瞬态呢？

dragon_2002_sd wrote:
          我觉得，间接耦合法模拟热-结构问题，假如单纯的分析温度产生的应力分布，如焊接应力，在结构分析时，读入温度体载荷后，无论做瞬态分析还是静态分析，应力应该是不变化的或者变化极小，因为在分析的时间内，载荷没有发生改变。因此我感觉在3e-5秒内钢板的应力也是变化很小的，但是得到的结果却截然不同？
    不知我的想法是否正确？

我觉得读入温度体载荷，再做瞬态的结构分析，得出的热应力也应该是瞬态的。温度体载荷作为一种载荷，在瞬态分析中也是随时间变化的，当然可以阶跃也可以线性变化，也可以变化为零（恒定）。楼主所说“因为在分析的时间内，载荷没有发生改变”根据是什么呢？
个人认为，若是做瞬态分析，在相应时刻读入温度载荷进行耦合，那么热应力也和结构力一样可以产生冲击效应－－热冲击
不知斑竹和ainite怎么认为的？大家积极讨论哦

dragon_2002_sd

to:nwpuyl  
   
    载荷没有发生改变我指的是温度没有变化，所以体载荷
没有变化，当然，约束没有变化。因此，得出，载荷没有发生改变，其实，这里，考虑动力学，应该是有变化的，但是分析时间很短，可以认为载荷变化很小。
  
  另外，我说的冲击载荷，并不是因为温度产生的热应力造成的冲击效应，而是要短时间内额外施加的冲击载荷。
  
nwpuyl 大虾，不知你如何分析我提出的问题？

nwpuyl

dragon_2002_sd你好！
不要叫我大虾，我也是小字辈，一起进步啦！
既然温度没有变化，就是载荷不变，当然是稳态过程。我明白你考虑的冲击载荷是额外施加的，所以才建议用稳态查看热应力。
但是你用瞬态分析，就已经把温度载荷当成瞬态的了，那就和原意不符了，虽然你没有定义温度的加载过程，默认是线性加载，比如你算的1秒的例子：温载在1秒内从绝对温度线性加载到读入的温度值。加载时间缩短了，结果当然不一样，所以我才认为温度载荷如果在极断时间内瞬态施加，也会有冲击－－忽略温度冲击就不能这样进行了吧！

lyonfeng

龙哥进展迅速！

dragon_2002_sd

to:nwpuyl
      昨天看了你的回贴，非常高兴。立刻做了回复。回去思考了一下，
明白了许多，想法和你现在的贴子完全一致。   
                  谢谢！互相交流！互相学习！