521majie

鐵蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 79
帖子: 52
在線: 26.6小時
蟲號: 3046058

[交流] 一個單層多道的例子，但是總是在進行一道之后，熱源就不再移動了。

找了一個單層多道的例子，但是總是在進行一道之后，熱源就不再移動了。
/UNITS,SI
/CONFIG,NRES,10000
/PREP7
ET,1,SOLID70    !單元類型，1表示此單元編號為1
MP,DENS,1,7930 !定義材料密度
MP,C,1,502    !定義材料比熱
mptemp,1,20,100,500
mpdata,kxx,1,1,12.1,16.3,21.4 !kxx熱導系數(shù)
/VIEW,1,1,1,1
!*************定義基板尺寸***********
lx=0.1
lz=0.1
ly=0.006       !基板大小
!*************定義焊接參數(shù)***********
layer=2       !總的焊接道數(shù)
lap=1/3       !兩道之間搭接量
weld=0.01       !每道焊道寬度
wweld=layer*weld-(layer-1)*lap*weld    !焊道的總寬度
wheight=0.004 !焊道的高度
remain=0.03

LSIZE=0.002
V=0.008       !速度
pi=3.1415926
U=25
I=180
Q=U*I          !電源有效功率（在此，可以對Q通過UI來調(diào)整）
R=0.006       !電源有效熱半徑
effect=0.8
Qmax=effect*Q/(pi*R*R)    !中心處最大熱流
tm=(lx-2*remain)/V       !一道的總時間
tinc=LSIZE/V             !每小段的時間間隔

!*****************建立模型，生成網(wǎng)格**************
block,0,lx,0,ly,0,lz       !創(chuàng)建一個塊體積，可以根據(jù)基板尺寸調(diào)整坐標
wpoff,remain,ly,(lz-wweld)/2  !隨著坐標系軸的變化平移工作平面
block,0,lx-2*remain,0,wheight,0,wweld
vsbw,all,,delete             !對block布爾運算
wpoff,,,wweld
vsbw,all,,delete
wprota,90,90                !轉(zhuǎn)動工作面
vsbw,all,,delete
wpoff,,,lx-2*remain
vsbw,all,,delete
save

vglue,all                   !粘合所有的幾何體

lsel,s,,,9,12,1
lsel,a,,,25,26,1
lsel,a,,,30
lsel,a,,,32
lsel,a,,,38
lsel,a,,,45
lsel,a,,,59
lsel,a,,,73
lesize,all,0.005,,,0.2,,,,1

allsel,all
lsel,s,,,13,16,1
lsel,a,,,19,21,1
lsel,a,,,23,24,1
lsel,a,,,34
lsel,a,,,41,44,1
lsel,a,,,63,64,1
lsel,a,,,77,78,1
lsel,a,,,81,82,1
lesize,all,0.001,,,1,,,,1

allsel,all
lsel,s,,,1
lsel,a,,,3,6,1
lsel,a,,,8
lsel,a,,,27
lsel,a,,,29
lsel,a,,,31
lsel,a,,,33
lsel,a,,,35
lsel,a,,,37
lsel,a,,,40
lsel,a,,,50
lsel,a,,,53
lsel,a,,,55
lsel,a,,,62
lsel,a,,,67
lsel,a,,,69
lsel,a,,,85
lesize,all,0.002,,,1,,,,1

allsel,all
lsel,s,,,46
lsel,a,,,86
lsel,a,,,88
lesize,all,0.001,,,5,,,,1

allsel,all
vsweep,all

/PNUM,DEFA
EPLOT
fini
!*******************進入后處理器*********************
/solu
antype,trans,new
!nlgeom,on
timint,0,struct
timint,1,therm
timint,0,mag
timint,0,elect
tref,25
nropt,auto
autos,on
kbc,0
pred,on
lnsrch,on
!************殺死焊縫去單元***********
nsel,s,loc,y,ly,ly+wheight
esln,s,1
ekill,all
allsel,all

esel,s,live
eplot

!**********施加對流載荷*************
esel,s,live
nsel,s,loc,z,0
nsel,a,loc,z,lz
nsel,a,loc,x,0
nsel,a,loc,x,lx
sf,all,conv,10,25

!************定義數(shù)組維數(shù)**********
MAX_TIME=1+tm/tinc
!MAX_X=1+(lx-2*remain)/LSIZE
MAX_X=1+lx/LSIZE
MAX_Z=1+lz/LSIZE

!*************定義table數(shù)組

!******************************************
*DIM,heatflux,TABLE,MAX_X,MAX_Z,MAX_TIME,X,Z,TIME
*do,n,1,layer,1
m=mod(n,2)       !返回n/2的余數(shù)部分
*if,m,eq,1,then
!n=1
      *DO,K,1,MAX_TIME,1
      *DO,I,1,MAX_X,1
      heatflux(I,0,K)=(I-1)*LSIZE
      *ENDDO
      *DO,J,1,MAX_Z,1
      heatflux(0,J,K)=(J-1)*LSIZE
      *ENDDO
      *DO,I,1,MAX_X,1
      *DO,J,1,MAX_Z,1
      XCENTER=remain+V*(K-1)*tinc
      ZCENTER=2*(n-1)*weld/3+(lz-wweld)/2+weld/2
      DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-ZCENTER)**2)
      *IF,DISTANCE,LE,R,THEN
      heatflux(I,J,K)=Qmax/EXP(2*DISTANCE**2/(R**2))
      *ELSE
      heatflux(I,J,K)=0
      *ENDIF
      *ENDDO
      *ENDDO
      heatflux(0,0,K)=(K-1)*tinc
!**********激活單元*********
      esel,s,live
      nsel,s,loc,x,xcenter-0.002,xcenter+0.001
      nsel,r,loc,z,zcenter-weld/2,zcenter+weld/2
      nsel,r,loc,y,ly,ly+wheight
      esln,s,0
      ealive,all
      allsel,all

!**********施加熱流載荷**********
      esel,s,live
      nsel,s,loc,Y,ly
      sf,all,hflux,%heatflux%
      allsel,all

!************求解*****************
      time,(n-1)*(tm+tinc)+k*lsize/v
      nsubst,2
      SOLVE
!*************刪除熱流載荷***********
      nsel,s,loc,y,ly
      sfdele,all,hflux
      allsel,all

      esel,s,live
      eplot

      OUTRES,ALL,ALL,
      /PSF,HFLUX,,2
      /REPLOT
      OUTRES,ALL,ALL,
      *ENDDO

*else
!n=2
      *DO,K,1,MAX_TIME,1
      *DO,I,1,MAX_X,1
      heatflux(I,0,K)=(I-1)*LSIZE
      *ENDDO
      *DO,J,1,MAX_Z,1
      heatflux(0,J,K)=(J-1)*LSIZE
      *ENDDO
      *DO,I,1,MAX_X,1
      *DO,J,1,MAX_Z,1
      XCENTER=lx-remain-V*(K-1)*tinc
      ZCENTER=2*(n-1)*weld/3+(lz-wweld)/2+weld/2
      DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-ZCENTER)**2)
      *IF,DISTANCE,LE,R,THEN
      heatflux(I,J,K)=Qmax/EXP(2*DISTANCE**2/(R**2))
      *ELSE
      heatflux(I,J,K)=0
      *ENDIF
      *ENDDO
      *ENDDO
      heatflux(0,0,K)=(K-1)*tinc
!**********激活單元*********
      esel,s,live
      nsel,s,loc,x,xcenter-0.001,xcenter+0.002
      nsel,r,loc,z,zcenter-weld/2,zcenter+weld/2
      nsel,r,loc,y,ly,ly+wheight
      esln,s,0
      ealive,all
      allsel,all

!**********施加熱流載荷**********
      esel,s,live
      nsel,s,loc,Y,ly
      sf,all,hflux,%heatflux%
      allsel,all

!************求解*****************
      time,(n-1)*(tm+tinc)+k*lsize/v
      nsubst,2
      SOLVE
!*************刪除熱流載荷***********
      nsel,s,loc,y,ly
      sfdele,all,hflux
      allsel,all

      esel,s,live
      eplot

      OUTRES,ALL,ALL,
      /PSF,HFLUX,,2
      /REPLOT
      OUTRES,ALL,ALL,
      *ENDDO
*endif
*enddo
FINISH

回復此樓

» 猜你喜歡

336求調(diào)劑已經(jīng)有5人回復
277材料科學與工程080500求調(diào)劑已經(jīng)有6人回復
333求調(diào)劑已經(jīng)有4人回復
286求調(diào)劑已經(jīng)有8人回復
298求調(diào)劑已經(jīng)有4人回復
306求0703調(diào)劑一志愿華中師范已經(jīng)有5人回復
299求調(diào)劑已經(jīng)有7人回復
26考研一志愿中國石油大學(華東)305分求調(diào)劑已經(jīng)有6人回復
初試 317 已經(jīng)有5人回復
296求調(diào)劑已經(jīng)有3人回復

» 本主題相關商家推薦: (我也要在這里推廣)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2014-05-26 08:58:46

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

簡單回復

ansonal2樓

2014-05-26 09:19 回復

521majie(金幣+1): 謝謝參與

baoer4103樓

2014-05-26 09:34 回復

521majie(金幣+1): 謝謝參與

neu2344樓

2014-05-26 09:36 回復

521majie(金幣+1): 謝謝參與

kingki5樓

2014-05-26 09:36 回復

521majie(金幣+1): 謝謝參與

kiss200855376樓

2014-05-26 09:39 回復

521majie(金幣+1): 謝謝參與

假大空7樓

2014-05-26 09:39 回復

521majie(金幣+1): 謝謝參與

521majie8樓

2014-05-26 09:45 回復

shenrenren9樓

2014-05-26 09:55 回復

521majie(金幣+1): 謝謝參與

henry83510樓

2014-05-26 10:08 回復

521majie(金幣+1): 謝謝參與

鐵觀音茶葉11樓

2014-05-26 10:11 回復

521majie(金幣+1): 謝謝參與

gaomeisdu12樓

2014-05-26 11:12 回復

521majie(金幣+1): 謝謝參與

521majie13樓

2014-05-26 15:18 回復

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 333求調(diào)劑 +4	87639 2026-03-21	4/200	2026-03-21 17:18 by 學員8dgXkO
[考研] 299求調(diào)劑 +4	某某某某位 2026-03-21	4/200	2026-03-21 16:30 by barlinike
[考研] 280求調(diào)劑 +9	咕嚕曉曉 2026-03-18	10/500	2026-03-21 16:01 by ColorlessPI
[考研] 297求調(diào)劑 +11	戲精丹丹丹 2026-03-17	11/550	2026-03-21 15:43 by ColorlessPI
[考研] 332求調(diào)劑 +3	鳳凰院丁真 2026-03-20	3/150	2026-03-21 10:27 by luoyongfeng
[考研] 306求調(diào)劑 +4	chuanzhu川燭 2026-03-18	4/200	2026-03-21 08:25 by laoshidan
[考研] 083200學碩321分一志愿暨南大學求調(diào)劑 +3	innocenceF 2026-03-17	3/150	2026-03-21 02:35 by JourneyLucky
[考研] 354求調(diào)劑 +5	Tyoumou 2026-03-18	8/400	2026-03-21 00:35 by JourneyLucky
[考研] 一志愿西北大學，070300化學學碩，總分287，雙非一本，求調(diào)劑。 +4	晨昏線與星海 2026-03-19	4/200	2026-03-20 22:15 by JourneyLucky
[考研] 本人考085602 化學工程專碩 +19	不知道叫什么！ 2026-03-15	21/1050	2026-03-20 20:48 by zhukairuo
[考研] 289求調(diào)劑 +6	懷瑾握瑜l 2026-03-20	6/300	2026-03-20 20:30 by 學員8dgXkO
[考研] 一志愿吉林大學材料學碩321求調(diào)劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑 +9	手機用戶 2026-03-17	9/450	2026-03-20 14:24 by 無懈可擊111
[考研] 求調(diào)劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學校友
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 0703化學 305求調(diào)劑 +4	FY_yy 2026-03-14	4/200	2026-03-19 05:54 by anny19840123
[考研] 311求調(diào)劑 +11	冬十三 2026-03-15	12/600	2026-03-18 14:36 by 星空星月
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 0703化學調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011