boxsteel 发表于 2003-3-11 08:51:51

【求助】关于“单元生死”的疑问

ANSYS用“焓”来计算相变,但ANSYS如何处理焊接的熔化过程呢?
: 焊缝的熔化过程:固态-熔化-熔点-液态-熔点-凝固。
: 在加热和冷却过程中,分别在熔点处,由于相变潜热的存在,熔池此时
: 应该温度不变,那么使用ANSYS如何解决这一过程的模拟呢?
: 我做了个简单的例子,即将达到一定温度的单元杀死,可是却发现
: 在后面的载荷作用下,前面已经被杀死的单元的温度还是变化的?
: 看帮助,是说被杀死的单元的载荷矩阵不是近似为0么,那么它的温度是否
: 应该不变化才对呢?
BTW:精华区里讲单元生死的例子,似乎没涉及分析相变过程中单元生死的应用。所有才发这个帖子。因为我需要在相变开始的时候杀死单元,相变结束的时候激活单元,相变开始可以温度来判断,那么相变结束用什么来判断呢?

那位知道,请多多指教!!!

boxsteel 发表于 2003-3-12 11:26:06

回复: 【求助】关于“单元生死”的疑问

大家多多帮助啊!!!
做焊接的没有人考虑这个问题么?ANSYS虽然通过H来考虑相变,但是相变过程中的温度不变,只考虑H是不能做到这一点的吧?
多谢!!!!

edodo 发表于 2003-3-12 11:38:20

回复: 【求助】关于“单元生死”的疑问

你要将你的log文件贴出来,不然理论上大家都知道,可是你操作中的错误就不清楚了!!

boxsteel 发表于 2003-3-12 12:31:05

回复: 【求助】关于“单元生死”的疑问

/prep7
/UNITS,si

et,1,solid70
mp,c,1,544.28
mp,kxx,1,15.06
mp,dens,1,4500

block,0,30,0,15,0,3

esize,1
vmesh,all
FINISH

/SOLU
!PARAMETER TABLE TO INPUT HEAT FLUX AS A MOVING SOURCE
!SPEED OF HEAT SOURCE IN X DIRECTION v=(200/60)m/s
!GAUSS VARIATION ASSUMED IN X & Y DIRECTION AS q=qmax/EXP(K*R**2)
!TOTAL HEAT INPUT Q=3.14*qmax/K
!DIMENSIONING TABLE NAMED "FLUX2"
v=1
L=30
NUM_TIME=L/v+1
NUM_X=L/v+1

*DIM,FLUX2,TABLE,NUM_TIME,NUM_X,7,TIME,X,Y

*DO,K,1,7,1

*DO,I,1,NUM_TIME,1
FLUX2(I,0,K)=(I-1)
*ENDDO

*DO,J,1,NUM_X,1
FLUX2(0,J,K)=(J*v-v)
*ENDDO

Q=1000
KAPA=1/4
qmax=((Q*KAPA)/3.14159)
PERCENT=0.98
RADIUS=ABS(SQRT(1/KAPA*LOG(1/(1-PERCENT))))

!N=

*DO,ANTIME,1,NUM_TIME,1
*DO,N,0,15,1
XDIST=(N*v)
*IF,(ANTIME-N),GT,0,THEN
*IF,SQRT((XDIST)**2+((RADIUS/6)*(K-1))**2),LE,RADIUS,THEN
!
GAUSS=qmax/EXP(KAPA*((XDIST)**2+((RADIUS/6)*(K-1))**2))
FLUX2(ANTIME,ANTIME-N,K)=GAUSS*10e5
!
*ELSE
*ENDIF
*ENDIF
!
*IF,(ANTIME+N),LE,NUM_TIME,THEN
*IF,SQRT((XDIST)**2+((RADIUS/6)*(K-1))**2),LE,RADIUS,THEN
!
GAUSS=qmax/EXP(KAPA*((XDIST)**2+((RADIUS/6)*(K-1))**2))
FLUX2(ANTIME,ANTIME+N,K)=GAUSS*10e5
!
*ELSE
*ENDIF
*ENDIF
*ENDDO
*ENDDO
!
!INPUT Y DIRECTION DIVISIONS
!---------------------------
FLUX2(0,0,K)=(K-1)*(RADIUS/6)
*ENDDO

antype,trans
tunif,20
kbc,1
    
*SET,delt,1
  timint,on
  autots,off

/VIEW, 1 ,1,1,1
/ANG, 1
/REP,FAST
FLST,2,1,5,ORDE,1
FITEM,2,2
/GO
!*
TIME,0.1
AUTOTS,0
DELTIM, , , ,1
KBC,0
!*
TSRES,ERASE
FLST,2,1,5,ORDE,1
FITEM,2,2
/GO
!*
!*
SFA,P51X,1,HFLUX, %FLUX2%
/STATUS,SOLU
SOLVE
OUTRES,ALL,ALL,
/PSF,hflux,,2

/replot
     

OUTRES,ALL,ALL,

FINISH

/POST1

ETABLE,MELTTEMP,TEMP
ESEL,S,ETAB,MELTTEMP,100
finish
/SOLU

EKILL,ALL

ESEL,INVE

antype,trans,rest,,,continue
tunif,20
kbc,1
    
*SET,delt,1
  timint,on
  autots,off

/VIEW, 1 ,1,1,1
/ANG, 1
/REP,FAST
FLST,2,1,5,ORDE,1
FITEM,2,2
/GO
!*
TIME,1
AUTOTS,0
DELTIM, , , ,1
KBC,0
!*
TSRES,ERASE
FLST,2,1,5,ORDE,1
FITEM,2,2
/GO
!*
!*
SFA,P51X,1,HFLUX, %FLUX2%
/STATUS,SOLU
SOLVE
OUTRES,ALL,ALL,
/PSF,hflux,,2

/replot
     

OUTRES,ALL,ALL,

FINISH

boxsteel 发表于 2003-3-12 12:36:36

回复: 【求助】关于“单元生死”的疑问

思路就是先施加第一步载荷,然后杀死在第一步载荷作用下温度超过100的单元,然后施加第二步载荷。
可是结果发现,在第二步载荷作用下,被杀死的单元的温度仍然变化?这是什么原因?帮助中说被杀死的单元的载荷矩阵不是为0么?那么是否应该在该单元被杀死后,其各项属性就应该保持不变呢?
请多指教!!!

咬紧牙关的蚊子 发表于 2003-3-12 13:52:52

回复: 【求助】关于“单元生死”的疑问

温度应该是变化的,这是因为被杀死的单元并不是不参与计算,而是其贡献由于乘上了一个1E-6的因子而被忽略掉了

boxsteel 发表于 2003-3-12 14:08:02

回复: 【求助】关于“单元生死”的疑问

那么用什么办法可以保证其温度不变化呢?不会是强制限制它的温度自由度吧?

yayasmiling 发表于 2003-4-17 10:35:24

回复: 【求助】关于“单元生死”的疑问

是啊
我现在也很苦恼这个问题
明明已经设为杀死的单元,居然还可以传热,而且传热量还很大
有人知道怎么解决这个问题吗?
我用的是ANSYS5.7版的,是不是版本低的原因

boxsteel 发表于 2003-4-17 22:20:33

回复: 【求助】关于“单元生死”的疑问

呵呵,对啊,我试验过,我选择某个单元,杀死,然后求解,查看结果,发现杀死的单元的温度还是随着其他单元一起变化的。也许就是因为被杀死的单元与其他单元连接在一起吧。共有节点,所以其温度不可能保持不变。

cswbww 发表于 2003-4-19 05:16:05

回复: 回复: 【求助】关于“单元生死”的疑问

boxsteel wrote:
那么用什么办法可以保证其温度不变化呢?不会是强制限制它的温度自由度吧?

可以选择被杀死的节点(不与任何活跃单元相连的),然后d,all,all,0

yayasmiling 发表于 2003-4-21 09:35:04

回复: 【求助】关于“单元生死”的疑问

我已经试过了
好像也不行!
页: [1]
查看完整版本: 【求助】关于“单元生死”的疑问