| 查看: 3680 | 回复: 0 | ||
| 【悬赏金币】回答本帖问题，作者xinxin998899将赠送您 30 个金币 | ||
xinxin998899新虫 (初入文坛)
|
[求助]
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 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
