| 查看: 1581 | 回复: 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 |
回复此楼» 猜你喜欢
29号明天会评吗
已经有4人回复
-
国社科又开始会评了，不知道这次命运如何
已经有10人回复
-
梦想
已经有4人回复
-
中青基了要发朋友圈吗？
已经有6人回复
-
基金不中，共勉
已经有12人回复
-
能否申诉？
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上意见出来了
已经有19人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
投票:
10