| 查看: 1564 | 回复: 0 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 ZJ970617 的 5 个金币 | |||
ZJ970617新虫 (初入文坛)
|
[交流]
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 |
回复此楼» 猜你喜欢
基金中了
已经有12人回复
-
关于代码变化问题，想知道的进来
已经有9人回复
-
关于filecode
已经有6人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有7人回复
-
关于豆爷回答的JTJC与%2F数量
已经有6人回复
100