应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » 一个单层多道的例子，但是总是在进行一道之后，热源就不再移动了。
北京石油化工学院2026年研究生招生接收调剂公告
131/1返回列表
查看: 986  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

521majie

铁虫 (小有名气)

[交流] 一个单层多道的例子，但是总是在进行一道之后，热源就不再移动了。

找了一个单层多道的例子，但是总是在进行一道之后，热源就不再移动了。
/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热导系数
/VIEW,1,1,1,1
!*************定义基板尺寸***********
lx=0.1
lz=0.1
ly=0.006         !基板大小
!*************定义焊接参数***********
layer=2          !总的焊接道数
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来调整）
R=0.006          !电源有效热半径
effect=0.8
Qmax=effect*Q/(pi*R*R)     !中心处最大热流
tm=(lx-2*remain)/V         !一道的总时间
tinc=LSIZE/V               !每小段的时间间隔


!*****************建立模型，生成网格**************
block,0,lx,0,ly,0,lz          !创建一个块体积，可以根据基板尺寸调整坐标
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                  !转动工作面
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

!************定义数组维数**********
MAX_TIME=1+tm/tinc
!MAX_X=1+(lx-2*remain)/LSIZE
MAX_X=1+lx/LSIZE
MAX_Z=1+lz/LSIZE

!*************定义table数组

!******************************************
*DIM,heatflux,TABLE,MAX_X,MAX_Z,MAX_TIME,X,Z,TIME
*do,n,1,layer,1
m=mod(n,2)        !返回n/2的余数部分
*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
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2014-05-26 08:58:46
已阅   回复此楼   关注TA 给TA发消息 送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
假大空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   回复  
相关版块跳转 我要订阅楼主 521majie 的主题更新
131/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 321求调剂 +13 认真求上学 2026-04-02 13/650 2026-04-04 18:23 by macy2011
[考研] 材料调剂 +4 一样YWY 2026-04-04 4/200 2026-04-04 18:11 by liu823948201
[考研] 调剂 +7 熊二想上岸 2026-04-04 7/350 2026-04-04 16:16 by dongzh2009
[考博] 申博 +4 IQwQl 2026-04-04 4/200 2026-04-04 13:57 by 工大学长
[考研] 320分人工智能调剂 +6 振—TZ 2026-04-03 6/300 2026-04-04 01:05 by Dream007008
[考研] 化学调剂求助 +6 LULONG1 2026-04-03 6/300 2026-04-03 23:13 by qzxyhcsy
[考研] 一志愿江南大学085501机械工程专硕326分，本科佳木斯大学 +4 顾若浮生 2026-04-03 8/400 2026-04-03 22:31 by zhq0425
[考研] 学硕288调剂!!! +3 小王xw123 2026-04-03 3/150 2026-04-03 21:20 by 啵啵啵0119
[考研] 372分材料与化工（085600）一志愿湖南大学求调剂 +3 蓝笺片 2026-04-03 4/200 2026-04-03 17:58 by Jimmyandyou
[考研] 考研调剂 +8 不爱喝饮料 2026-04-03 8/400 2026-04-03 16:40 by Mistake-J
[考研] 321求调剂 +17 y-yh 2026-04-01 20/1000 2026-04-03 12:57 by y-yh
[考研] 316求调剂 +14 舟自梗 2026-04-01 18/900 2026-04-03 10:28 by linyelide
[考研] 085801 总分275 本科新能源 求调剂 +19 bradoner 2026-04-01 23/1150 2026-04-03 10:07 by linyelide
[考研] 338求调剂，一志愿能源动力，外语是日语203 +5 zzz，，r 2026-04-02 5/250 2026-04-03 09:45 by 蓝云思雨
[考研] 372分材料与化工（085600）一志愿湖南大学求调剂 +5 蓝笺片 2026-04-02 6/300 2026-04-02 21:37 by dongzh2009
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 09:52 by yulian1987
[考研] 本2一志愿C9-333分，材料科学与工程，求调剂 +9 升升不降 2026-03-31 9/450 2026-03-31 18:01 by 无际的草原
[考研] 南京大学化学调剂 +11 景随风 2026-03-29 16/800 2026-03-31 10:14 by herarysara
[考研] 福建理工大学材料学院先进合金团队招收考研调剂学生 +3 大华金商都 2026-03-30 4/200 2026-03-31 01:04 by 方英俊602
[考研] 297 地理学070500 复试求调剂 +3 小圆圈圈ooo 2026-03-30 3/150 2026-03-30 21:05 by 余震yz
信息提示
关闭
请填处理意见
关闭