| 查看: 1442 | 回复: 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 |
回复此楼» 猜你喜欢
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
假如你的研究生提出不合理要求
已经有5人回复
-
所感
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有28人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复