<gust1977>版主，你好，帮忙吧，谢谢！
版主，你好，我在做热结构耦合分析，进行瞬态温度场应力场分析，包括机械应力，但结果不正确，您帮我看看怎么回事好吗？程序如下，我用的是 直接耦合的方式（自己认为），你觉得呢？选用solid 5单元，如果采取间接耦合，进行瞬态温度场分析后，再进行应力场分析，您觉得应该怎么编写程序？瞬态温度场分析的程序也是您以前写的！快要交论文了，呵呵，有时间就快点答复，谢谢您了！
别的大虾们也帮忙啊！！！
Re:<gust1977>版主，你好，帮忙吧，谢谢！
程序如下：：：! define drum material parameter
/prep7
MPTEMP,1,297,373,473,573,673,773
MPDATA,DENS,1,,7330
MPDATA,EX,1,1,148e9,145e9,141e9 ,140e9 ,137e9 ,132e9
MPDATA,PRXY,1,1,0.31,0.31 ,0.29 ,0.30 ,0.29 ,0.27
UIMP,1,REFT,,,297
MPDATA,ALPX,1,1,1.010E-005,1.098E-005,1.141E-005,1.207E-005,1.266E-005,1.299E-005
MPDATA,KXX,1,1, 40, 40.0,40.3 ,38.3 ,36.5 ,35.1
MPDATA,C,1,1,480,494 ,536,557,565 ,573
！creat drum model
et,1,solid5
k,1,0,0,-0.038
k,2,0,0,0.200
k,3,0.200,0,0.200
k,4,0.215,0,0.200
k,5,0.215,0,0.185
k,6,0.230,0,0.185
k,7,0.230,0,0.142
k,8,0.215,0,0.142
k,9,0.215,0,0
k,10,0.200,0,0
k,11,0.176,0,-0.038
k,12,0.176,0,-0.024
k,13,0.117,0,-0.038
k,14,0.117,0,-0.024
k,15,0.200,0,0.185
k,16,0.215,0,0.185
LSTR, 3, 4
LSTR, 4, 5
LSTR, 5, 6
LSTR, 6, 7
LSTR, 7, 8
LSTR, 8, 9
LSTR, 9, 11
LSTR, 11, 13
LSTR, 13, 14
LSTR, 14, 12
LSTR, 12, 10
LSTR, 10, 3
LPLOT
LSTR, 12, 11
LSTR, 10, 9
LSTR, 8, 5
LSTR, 15, 16
LSEL,S,,,1
*DO,I,2,12,1
LSEL,A,,,i
*ENDDO
AL,ALL
LSEL,S,,,13
*DO,I,14,16,1
LSEL,A,,,I
*ENDDO
ASBL, 1,ALL
APLOT
FLST,2,5,5,ORDE,2
FITEM,2,2
FITEM,2,-6
FLST,8,2,3
FITEM,8,1
FITEM,8,2
VROTAT,P51X, , , , , ,P51X, ,360, ,
VSEL,ALL
CM,DRUM,VOLU
CMSEL,A,DRUM
/MREP,EPLOT
!mesh drum
cmsel,s,drum
!*
vatt, 1, , 1, 0
ESIZE,0.010,0
vsweep,drum
/solu
anty,trans
outres,basic,all
csys,1
tunif,297
NSEL,S,LOC,Z,-0.038
D,ALL,,,,,,UZ,,,,,
ALLSEL,ALL
NSEL,S,LOC,X,0.2
NSEL,R,LOC,Z,0,0.185
SF,all,PRES,1.4e6
allsel,all
OMEGA,0,0,57,0
DOMEGA,0,0,-12.1,
*SET,tm_start,1e-16
*SET,tm_end,4.7
*SET,tm_inc,0.1
NSEL,S,LOC,X,0.215
NSEL,R,LOC,Z,0,0.142
NSEL,A,LOC,X,0.230
SF,ALL,CONV,62.3,297
NSEL,S,LOC,X,0.215,
NSEL,R,LOC,Z,0.185,0.200
SF,ALL,CONV,62.3,297
NSEL,S,LOC,X,0.215,0.230
NSEL,R,LOC,Z,0.142
NSEL,A,LOC,Z,0.185
SF,ALL,CONV,62.3,297
ASEL,S,,,16
ASEL,A,,,38
ASEL,A,,,60
ASEL,A,,,80
SFA,ALL,1,CONV,62.3,297
*do,tm, tm_start, tm_end, tm_inc
time,tm
deltim,tm_inc
NSEL,S,LOC,X,0.200,
NSEL,R,LOC,Z,0,0.185
sf,all,hflux,1656508+(8589-1656508)/4.7*tm
nsel,all
solve
*enddo
Re:<gust1977>版主，你好，帮忙吧，谢谢！
为了方便大虾们容易的读懂程序，说明如下：采用的单元 solid 5
制动鼓内表面施加随时间变化的热流密度，外表面施加散热边界条件
内表面还受制动压力（常数）
另外 还施加了惯性载荷，即因制动鼓的转动的载荷；
版主，麻烦您仔细看看，期待您的回复。。。
Re:<gust1977>版主，你好，帮忙吧，谢谢！
为了方便大虾们容易的读懂程序，说明如下：采用的单元 solid 5
制动鼓内表面施加随时间变化的热流密度，外表面施加散热边界条件
内表面还受制动压力（常数）
另外 还施加了惯性载荷，即因制动鼓的转动的载荷；
版主，麻烦您仔细看看，期待您的回复。。。
Re:<gust1977>版主，你好，帮忙吧，谢谢！
gus1977,您好，不在吗？？您没有看到我的留言吗？？请回复
谢谢
Re:<gust1977>版主，你好，帮忙吧，谢谢！
呵呵，最近比较忙，你忙着交论文，我也在忙我的出站报告，所以没上来，大致看了一下你的程序，我也看不出哪个地方有问题，呵呵。还是请大家一起来看看把，我这周报告必须要出第一稿，不好意思，如果你能够等到下周的话，我仔细帮你看看，可以吗Re:<gust1977>版主，你好，帮忙吧，谢谢！
好的,等到下周吧,多谢了!一定不要忘了帮我看看,你也可以帮我重新编写一个,因为瞬态温度场的求解程序是你编写的,呵呵,不知道你是否记得,就是施加一个随时间变化的热流密度,要是接着下去求这个随时间变化的热流密度引起的热应力怎么编写呢????Re:<gust1977>版主，你好，帮忙吧，谢谢！
你的程序我大致看了一下，主要原因可能是你没有施加位移约束，因为你在直接耦合时，要计算应力和应变，如果没有边界约束，很可能会导致结果不正确。其实要进行间接耦合的话，比较简单，而且容易收敛。就是首先只计算温度场，将solid5的Ux、UY、Uz自由度关闭，在第一阶段只计算温度场。得到温度场结果后，将solid5的温度自由度关闭，开通位移自由度，将得到的温度场结果作为边界条件加上去，计算应力应变场即可。具体作法，论坛里有不少例子，同时，我也可以给你一个算例，我会给你发到信箱里，你自己去研究吧，我最近很忙，程序不能给你编了，你看着我给你的例子，自己去作，应该不难
页:
[1]