yezi 发表于 2002-8-6 19:52:40

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

我做的移动激光束扫描金属表面, 求金属表面的温度场以及最终残余应力
其中激光束(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=
!
*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 发表于 2002-8-6 19:58:02

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

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

bigbird 发表于 2002-8-6 20:40:20

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

這是一個已參數化的程式.....
其實他已寫的很清楚.........
____________________________________________________
!
---->敘述句...也就是標注...以下的步驟是在做什麼或是某參數的意義
____________________________________________________
*DO
....
.....
*ENDDO
---->作一個do loop
____________________________________________________
*IF
.....
.....
*ENDIF
---->作一各命令句...當參數等於或大於或小於....就會執行或停止等動作
_____________________________________________________
p=123
q=456
y=789
這就是所說的參數化.....
你可隨意指定一個值.....此程式就會依據你給的參數值下去跑....
用來建模非常方便.....
_____________________________________________________

這個程式是大致寫的還不錯..........唯一...*else是多餘的
還有....do loop最好這樣寫
*do
空格......
空格*do
      空格空格.....
空格*enddo
空格.....
*enddo
也就是把包含兩個 loop區隔出來.....
這樣比較不會眼花....容易出錯...
______________________________________________
這樣應該都懂了吧...^^
懂了幫我測試一下我的問題吧.....^^!!!

yezi 发表于 2002-8-7 09:09:57

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

谢谢!

你说的问题是关于MUTIPRO的吗?
我现在用的是5.7
哪里能DOWNLOAD 6.0? 我装一下,希望能有帮助

(我接触ANSYS刚两个月,复杂的问题可能回答不了,还请原谅哦!!!)

bigbird 发表于 2002-8-7 11:39:10

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

^^...沒關係....6.1以前的版本都可以用.......
我的問題很簡單....
就是幫我測試一下....run run看...
寫二十個以上...跑四次以上....
看mutipro這個指令還能不能讀到......
目前有四台電腦.....已經確定不行...
想多測試幾台....看是不是ANSYS本身的問題...^^!!

(我接触ANSYS刚两个月, 复杂的问题可能回答不了, 还请原谅哦!!!)
......有這程度已經很好囉^^

xbzhang 发表于 2003-11-6 19:54:52

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

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

zlf81 发表于 2004-7-1 18:58:25

vivian1010 发表于 2006-7-31 16:00:29

之前的命令流与移动热源有关么?如果相关也EMAIL一份完整的命令流给我,不胜感激:)
lrain1010@sina.com谢谢哦!

walabc086828 发表于 2013-3-8 16:08:41


也给我发一份吧
谢谢 lz

walabc086828 发表于 2013-3-8 16:09:29

忘记给你 邮箱了20lpf@163.com
页: [1]
查看完整版本: 【讨论】有关移动HEAT SOURCE 加载的COMMAND LINE, 其逻辑看不明白