24小时热门版块排行榜

返回列表

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 ZJ970617 的 5 个金币

ZJ970617

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 29
散金: 3
帖子: 18
在线: 5.1小时
虫号: 17579372
注册: 2019-09-20
性别: GG
专业: 金属材料的合金相、相变及

[交流] Abaqus高斯热源子程序

小弟最近在研究激光仿真，对旋转高斯体热源子程序方面有些小问题，希望能和大家多多交流。

SUBROUTINE DFLUX(FLUX,SOL,JSTEP,JINC,TIME,NOEL,NPT,COORDS,JLTYP,
   1                TEMP,PRESS,SNAME)
C
c    use GlobalVar
   INCLUDE 'ABA_PARAM.INC'
C
   parameter(one=1.d0)
   DIMENSION COORDS(3),FLUX(2),TIME(2)
   CHARACTER*80 SNAME
C    高斯旋转体热源参数，H为热源高度，R0为热源作用有效半径
   Q=2000
   effi=0.8
   P=Q*effi
   R0=0.004
   H=0.005
   v=0.01
   pi=3.1415

x=COORDS(1)
      y=COORDS(2)
z=COORDS(3)

C    初始热源中心坐标
x0=0
y0=0
   z0=0
C    热源形状及热源公式，d当前时刻焊接斑点中心跟焊接初始位置的距离
   d=v*TIME(2)
   Cs=3/R0**2
   K=3*Cs*P/(pi*H*(1-1/(EXP(3)))*EXP(-3*Cs/LOG(H/Z)*((X-d)**2+Y**2)


C    JLTYP=1，体热源
   JLTYP=1
FLUX(1)=3*Cs*P/(pi*H*(1-1/(EXP(3)))*EXP(-3*Cs/LOG(H/Z)*((X-d)**2+Y**2)
   RETURN
   end subroutine DFLUX

回复此楼

» 猜你喜欢

1楼 2020-12-02 10:08:15

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 ZJ970617 的主题更新

返回列表