24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者xinxin998899将赠送您 30 个金币

xinxin998899

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3
帖子: 2
在线: 44分钟
虫号: 35412709
注册: 2024-05-06
性别: MM
专业: 金属材料的合金相、相变及

[求助] abqus的二道激光淬火代码

求助大佬帮忙。软件：abqus，激光淬火数值模拟，撰写激光淬火代码，一道激光淬火代码已经写好了，现在差按照原路返回的二道激光淬火代码，只需要增加按照原路径返回的代码，现有代码：

SUBROUTINE DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS,JLTYP,
   1                TEMP,PRESS,SNAME)
C
   INCLUDE 'ABA_PARAM.INC'
   parameter(one=1.d0)
   DIMENSION COORDS(3),FLUX(2),TIME(2)
   CHARACTER*80 SNAME
C    Q,激光有效热功率W
C    v,激光移动速度m/s
C    Rh,加热斑点半径，95%的热量落在以Rh为半径的面积内
C d,当前时刻焊接斑点中心跟焊接初始位置的距离
Q=1000
v=0.01
Rh=0.003
d=v*TIME(2)

x=COORDS(1)
y=COORDS(2)
z=COORDS(3)
C    焊接板厚度为0.01m,焊枪移动从坐标-0.1，0，0.01开始，沿着x方向移动
x0=-0.1
y0=0
            z0=0.003
C
PI=3.1415
C
R=sqrt((x-x0-d)**2+(y-y0)**2)
C    JLTYP＝0，表示为面热源
JLTYP=0
IF (KSTEP. eq. one) THEN
      FLUX(1)=3*Q/(PI*Rh**2)*exp(-3*R**2/Rh**2)
   ENDIF
   RETURN
   END

回复此楼

» 猜你喜欢

1楼2025-01-06 13:12:53

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

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

返回列表