| 查看: 1551 | 回复: 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 |
回复此楼» 猜你喜欢
明天E口面上会评
已经有8人回复
-
不要再数国自然申请书的 filecode 的分隔符个数了
已经有17人回复
-
内心匮乏
已经有3人回复
-
轴手性化合物的确定
已经有4人回复
-
祈祷青基必中
已经有16人回复
-
生命口会评
已经有12人回复
-
27届辽宁大学应届毕业生申博
已经有3人回复
-
关于如何从代码看上不上会
已经有23人回复
-
chemdraw
已经有6人回复
-
b口会评
已经有6人回复