【讨论】有关移动HEAT SOURCE 加载的COMMAND LINE, 其逻辑看不明白

yezi

我做的移动激光束扫描金属表面, 求金属表面的温度场以及最终残余应力
其中激光束(HEAT SOURCE)横截面上能量分布为GAUSS分布,  下面是我从同学那里得到的有关处理移动热源的命令流, 但是不太看得明白, 请大家帮忙分析.
  
移动速度v=200/60
总得扫描距离L=30
激光功率Q=360( 不是功率密度)
  
/SOLU
!PARAMETER TABLE TO INPUT HEAT FLUX AS A MOVING SOURCE
!SPEED OF HEAT SOURCE IN X DIRECTION v=(200/60)m/s
!GAUSS VARIATION ASSUMED IN X & Y DIRECTION AS q=qmax/EXP(K*R**2)
!TOTAL HEAT INPUT Q=3.14*qmax/K
!DIMENSIONING TABLE NAMED "FLUX2"
v=(200/60)
L=30
NUM_TIME=L/v+1
NUM_X=L/v+1
!-------------------------------------------------------
*DIM,FLUX2,TABLE,NUM_TIME,NUM_X,7,TIME,X,Y
!DO LOOP FOR ALL PLANES, ie Y DIRECTION
*DO,K,1,7,1
!INPUT FOR TIME COLUMN
*DO,I,1,NUM_TIME,1
FLUX2(I,0,K)=(I-1)
*ENDDO
!INPUT FOR X-DIRECTION DIVISION
*DO,J,1,NUM_X,1
FLUX2(0,J,K)=(J*v-v)
*ENDDO
!-------------------------------
!INPUT GAUSS VARIATION FOR HFLUX
!-------------------------------
!INPUT VARIABLES:
!---------------  
!Q      =    HEAT INPUT     -   J/ SEC
!KAPA  =  CONCENTRATION FACTOR     - 1/mM**2
!PERCENT=  % OF Q CONSIDERED FOR HEAT INPUT AT "RADIUS"  
!-------------------------------
Q=360
KAPA=1/4
qmax=(Q*KAPA)/3.14159
PERCENT=0.98
RADIUS=ABS(SQRT(1/KAPA*LOG(1/(1-PERCENT))))
!
!N=[RADIUS/v]
!
*DO,ANTIME,1,NUM_TIME,1
*DO,N,0,15,1
XDIST=(N*v)
*IF,(ANTIME-N),GT,0,THEN
*IF,(XDIST+ABS((RADIUS-(RADIUS/3)*(K-1)))),LE,RADIUS,THEN  
!
GAUSS=qmax/EXP(KAPA*(XDIST+ABS((RADIUS-(RADIUS/3)*(K-1))))**2)
FLUX2(ANTIME,ANTIME-N,K)=GAUSS
!
*ELSE
*ENDIF
*ENDIF
!
*IF,(ANTIME+N),LE,NUM_TIME,THEN
*IF,(XDIST+ABS((RADIUS-(radius/3)*(K-1)))),LE,RADIUS,THEN
!
GAUSS=qmax/EXP(KAPA*(XDIST+ABS((RADIUS-(RADIUS/3)*(K-1))))**2)
FLUX2(ANTIME,ANTIME+N,K)=GAUSS
!
*ELSE
*ENDIF
*ENDIF
*ENDDO
*ENDDO
!
!INPUT Y DIRECTION DIVISIONS
!---------------------------
FLUX2(0,0,K)=((25-RADIUS)+(K-1)*RADIUS/3)  
*ENDDO

yezi

最好那位大侠看懂以后, 能在COMMAND LINE 后分别加注释,  
BIGBIRD、SIMWE管理员、还有。。。。大家一定要知无不言啊

bigbird

這是一個已參數化的程式.....
其實他已寫的很清楚.........
____________________________________________________
!
---->敘述句...也就是標注...以下的步驟是在做什麼或是某參數的意義
____________________________________________________
*DO
....
.....
*ENDDO
---->作一個do loop
____________________________________________________
*IF
.....
.....
*ENDIF
---->作一各命令句...當參數等於或大於或小於....就會執行或停止等動作
_____________________________________________________
p=123
q=456
y=789
這就是所說的參數化.....
你可隨意指定一個值.....此程式就會依據你給的參數值下去跑....
用來建模非常方便.....
_____________________________________________________
  
這個程式是大致寫的還不錯..........唯一...*else是多餘的
還有....do loop最好這樣寫
*do
  空格......
  空格*do
        空格空格.....
  空格*enddo
  空格.....  
*enddo
也就是把包含兩個 loop區隔出來.....
這樣比較不會眼花....容易出錯...
______________________________________________
這樣應該都懂了吧...^^
懂了幫我測試一下我的問題吧.....^^!!!

yezi

谢谢!  
  
你说的问题是关于MUTIPRO的吗?
我现在用的是5.7
哪里能DOWNLOAD 6.0? 我装一下,希望能有帮助
  
(我接触ANSYS刚两个月,  复杂的问题可能回答不了,  还请原谅哦!!!)

bigbird

^^...沒關係....6.1以前的版本都可以用.......
我的問題很簡單．．．．
就是幫我測試一下．．．．run run看．．．
寫二十個以上．．．跑四次以上．．．．
看mutipro這個指令還能不能讀到．．．．．．
目前有四台電腦．．．．．已經確定不行．．．
想多測試幾台．．．．看是不是ＡＮＳＹＳ本身的問題．．．^^!!
  
(我接触ANSYS刚两个月, 复杂的问题可能回答不了, 还请原谅哦!!!)
......有這程度已經很好囉^^

xbzhang

楼组能否把/solu前的命令流给我一份好吗？我也正在作类似的分析，刚接手，无从下手。先多谢谢了！我的邮件：jbzhang_cn@sina.com

vivian1010

之前的命令流与移动热源有关么？如果相关也EMAIL一份完整的命令流给我，不胜感激：）
lrain1010@sina.com谢谢哦！

walabc086828

也给我发一份吧
谢谢 lz

walabc086828

忘记给你 邮箱了  20lpf@163.com