hangch 发表于 2004-12-31 09:48:58

[求助]移动热源作用下两层材料的应力分析问题

采用间接法来计算应力场,温度场程序代码如下:
/PREP7
ET,1,SOLID70
ET,2,SOLID70
MP,DENS,1,5250
MP,C,1,650
MP,KXX,1,5

MP,DENS,2,7800
MP,C,2,502
MP,KXX,2,16.3
/VIEW,1,1,1,1
BLOCK,0,0.041,0,0.036,0,0.004
BLOCK,0,0.041,0,0.036,0.004,0.0043
VGLUE,ALL
LESIZE,1,,,36
LESIZE,4,,,41
LESIZE,9,,,4
MSHKEY,1
VATT,1,1,1
VMESH,1
LESIZE,8,,,36
LESIZE,5,,,41
LESIZE,25,,,1
MSHKEY,1
VATT,2,1,2
VMESH,3
/PNUM,DEFA
EPLOT
SAVE
FINISH

/SOLU
LSIZE=0.001
V=0.4 !速度
TINC=LSIZE/V !单步时间长度
KK=1 !能量集中系数
Qmax=1800000
R=0.012
MAX_TIME=41
MAX_X=41
MAX_Y=36
   
*DIM,FLUX2,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME
*DO,K,1,MAX_TIME,1
*DO,I,1,MAX_X,1
FLUX2(I,0,K)=(I-1)*LSIZE
*ENDDO
   
*DO,J,1,MAX_Y,1
FLUX2(0,J,K)=(J-1)*LSIZE
*ENDDO
   
   *DO,I,1,MAX_X,1
    *DO,J,1,MAX_Y,1
   XCENTER=0.02
   YCENTER=0.01
   DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2)
   *IF,DISTANCE,LE,R,THEN
   FLUX2(I,J,K)=Qmax/EXP(KK*DISTANCE**2/(2*0.004**2))
   *ELSE
   FLUX2(I,J,K)=0
   *ENDIF
    *ENDDO
   *ENDDO

FLUX2(0,0,K)=(K-1)*TINC
*ENDDO

!开始计算
ANTYPE,TRANS
TUNIF,150 !初始温度150度
KBC,1
   
*SET,DELT,1
   TIMINT,ON
   AUTOTS,OFF
   
   TT=0.05

*DO,TM,1E-6,TT,TINC
/VIEW,1,1,1,1
/ANG,1
/REP,FAST
   
FLST,2,1,5,ORDE,1
FITEM,2,8
   
/GO
!*
TIME,TM
AUTOTS,0
DELTIM, , , ,1
KBC,0
!*
TSRES,ERASE
FLST,2,1,5,ORDE,1
FITEM,2,8
/GO
!*
!*
SFA,P51X,1,HFLUX,%FLUX2%
/STATUS,SOLU
SOLVE
OUTRES,ALL,ALL,
/PSF,HFLUX,,2
   
/REPLOT
      
*ENDDO

OUTRES,ALL,ALL,

FINISH
应力场代码:
/SOLU
LSCLEAR,SOLID

/PREP7
ETCHG,TTS
KEYOPT,1,2,0
KEYOPT,1,4,0
KEYOPT,1,7,0
KEYOPT,1,8,0

KEYOPT,1,2,0
KEYOPT,1,4,0
KEYOPT,1,7,0
KEYOPT,1,8,0

MP,EX,1,3.66E11
MP,ALPX,1,7.962E-6
MP,NUXY,1,0.25

MP,EX,2,1.2E11
MP,ALPX,2,11.5E-6
MP,PRXY,2,0.3

!NSEL,S,LOC,Z,0.004
!CP,1,UZ,ALL
!CP,2,UX,ALL
!CP,3,UY,ALL
ALLSEL,ALL
FINISH

/SOLU
!ANTYPE,TRANS,REST
DK,1,UX,0,,UY,UZ      !固定底面左上角的顶点
!DK,2,UY,0      !固定左侧面左下角的顶点
!DL,1,5,UZ,0    !固定左侧面底部直线
!DA,5,UX,0      !固定左侧面
*DO,tm,1E-6,0.05,0.0025 !设置时间步长和区间
!ANTYPE,TRANS,REST,41
TIME,tm !当前时间为tm
LDREAD,TEMP,,,TIME,0,'sample','RTH' !读入时间tm时的温度分布
DELTIM,,,,20 !初始步长
SOLVE
*ENDDO
FINISH
计算出来的应力却呈现对称状,而热源的轨迹并不对称,不知为什么,请高手指教

hangch 发表于 2004-12-31 09:50:17

Re:[求助]移动热源作用下两层材料的应力分析问题

温度场图

hangch 发表于 2005-1-10 15:05:58

Re:[求助]移动热源作用下两层材料的应力分析问题

没有高人指点吗,太郁闷了

开始进取 发表于 2005-1-10 21:49:14

Re:[求助]移动热源作用下两层材料的应力分析问题

FLST,2,1,5,ORDE,1
FITEM,2,8
请问上面两条语句是什么意思呢?传热理论方面的东西你最好问问热学方面的老师,我估计大家对很专业的东西都不怎么懂得。

开始进取 发表于 2005-1-10 22:38:48

Re:[求助]移动热源作用下两层材料的应力分析问题

FLST,2,1,5,ORDE,1
FITEM,2,8
请问上面两条语句是什么意思呢?传热理论方面的东西你最好问问热学方面的老师,我估计大家对很专业的东西都不怎么懂得。

开始进取 发表于 2005-1-11 23:26:12

Re:[求助]移动热源作用下两层材料的应力分析问题

我知道了。是ANSYS自动生成的。

hangch 发表于 2005-1-17 15:34:12

Re:[求助]移动热源作用下两层材料的应力分析问题

我好可怜啊,没有人帮忙

liujin2000161 发表于 2006-8-16 21:10:40

热源没有移动阿!

cj_hmily 发表于 2006-8-20 10:54:17

FLST,2,1,5,ORDE,1
FITEM,2,8
是用选择菜单在工作面里拾取某些实体(关键点、线等)时的命令流,这种命令流可读性差,如果采用APDL语言的话,最好不要使用.

kitty827 发表于 2006-11-14 10:36:13

好像热源是定的阿,没有移动吧?
页: [1]
查看完整版本: [求助]移动热源作用下两层材料的应力分析问题