| 查看: 1611 | 回复: 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 |
回复此楼» 猜你喜欢
越来越担心硕士毕业后没有竞争力了
已经有14人回复
-
论文的挂名怎么评判贡献？基金申请中，排名第一的挂名文章
已经有3人回复
-
申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请
已经有3人回复
-
闲聊
已经有6人回复
-
申博求助
已经有4人回复