*GET 得到最大TASI-WU值與輸出TXT檔??

maocp

現有複合材料疊層SHELL99  共24層
想得到每一層的最大破壞節點號碼以及最大破壞節點TASI-WU的數值
最後將12個數值儲存起來 來配合C++基因演算法做比較
請問該如何撰寫
下面是我的寫法
不知該如何直接存取最大值出來 ，所以將每層每個節點都存下來
可這樣子效率並不高
麻煩大家幫我改寫，並看看哪裡有錯
以及存檔格式(E13.6) 定義是如何?  我只直接把網路上的格式複製下來使用




NODENUM=855                         !節點總數目     

*DIM,TWSI,ARRAY,TWSI*NODENUM              !定義TWSI矩陣存放大小
*CFOPEN,DISP,TXT,C:\HE,,                          !開啟文件檔儲存位置     

*DO,N,1,NODENUM,1
layer,1
*GET,max_twsi,NODE,TWSI(N),S,TWSI
layer,2
*GET,max_twsi,NODE,TWSI(N),S,TWSI
layer,3
*GET,max_twsi,NODE,TWSI(N),S,TWSI
layer,4
*GET,max_twsi,NODE,TWSI(N),S,TWSI
layer,5
*GET,max_twsi,NODE,TWSI(N),S,TWSI
layer,6
*GET,max_twsi,NODE,TWSI(N),S,TWSI
layer,7
*GET,max_twsi,NODE,TWSI(N),S,TWSI
layer,8
*GET,max_twsi,NODE,TWSI(N),S,TWSI
layer,9
*GET,max_twsi,NODE,TWSI(N),S,TWSI
layer,10
*GET,max_twsi,NODE,TWSI(N),S,TWSI
layer,11
*GET,max_twsi,NODE,TWSI(N),S,TWSI
layer,12
*GET,max_twsi,NODE,TWSI(N),S,TWSI

*ENDDO  

*VWRITE,TWSI(1)                    !寫入
(E13.6)
FINISH

!*************************************

maocp

*DIM,TWSI,ARRAY,TWSI                       !定義TWSI矩陣存放大小
*CFOPEN,DISP,TXT,C:\HE,,                   !開啟文件檔儲存位置     
layer,1                                                   !觀看第一層
allsel                                                       !選取全部節點
nsort,fail,twsi,0,1,all                                !排列大小的方式,破壞準則,tasi-wu,降序,取絕對值,所有節點
*get,max_twsi,sort,0,max                     !得到max_twsi值(自己定的參數),排列,NETNUM(不瞭解 預設),最大值(若薇imax則為最大值的節點號碼)

layer,2
...
...
...

*VWRITE,TWSI(1)                                !寫入
(E13.6)                                                  !依然不瞭解

然後輸出的值為0.0000000E+00

附檔為我的log file檔 到底我錯在哪裡呢?

wjqxlz

楼上说的有道理~

maocp

此篇沉下去了，我不知道哪裡有錯誤一直無法提取數值，都是空白的
一層一層*get  並在最後做寫入動作這樣子
寫入的到底會有幾筆呢?

還是要改用*VGET ??

煩請明指

lscookie

你可以将1L中的各层数据存入一个数组，用数组函数操作得到最值，然后写出即可
参考*vscfun命令，另外nsort操作需要在/post1后处理器中运行才能生效，2L的错误原因查看下nsort命令是否使用正确

maocp

的確是在nsort出現錯誤，可我在GUI介面下操作排列指令依然產生錯誤，附上DB檔(ANSYS12) 可否請優秀斑竹幫我查閱錯誤
此研究抓擷取就卡住快一個月了，欲哭無淚現在，麻煩了謝謝

lscookie

twsi是自己定义的量的话  nsort就会报错，nsort只对系统自带的变量排序，如果你自己定义的量可以用etable 然后esort命令，如果你只想得到最值 就用vscfun

maocp

妳好 這個是我修改過後的結果，可是擷取的數值皆為同一個
我直接在ANSYS GUI介面下執行
*GET,TWSI,NODE,23,S,TWSI
是有法子直接顯示節點23的TASI-WU值
可是改以下面迴圈方式就出錯了 可否再幫我看看呢?
先謝謝各位高手了



NP=12                                                      
NODENUM=4905                           !自己去看所有節點數目所得到的   

*DIM,TWSIVAL,ARRAY,NODENUM   
*DIM,TWSIMAX,ARRAY,NP
*CFOPEN,DISP,DAT,C:\HE,,         


/POST1
*DO,LAY,1,NP,1            
LAYER,LAY                  

    *DO,TWSINODE,1,NODENUM,1  
    *GET,TWSIVAL(TWSINODE),NODE,TWSINODE,S,TWSI
    *ENDDO

*VSCFUN,TWSIMAX(LAY),MAX,TWSIVAL(1)

*ENDDO
*VWRITE,TWSIMAX(1)         
(E13.6)

FINISH

maocp

maocp 发表于 2013-2-6 21:14
*DIM,TWSI,ARRAY,TWSI                       !定義TWSI矩陣存放大小
*CFOPEN,DISP,TXT,C:\HE,,            ...

我知道我這裡的錯誤了
我*GET A處的數值
卻*VWRITE,寫到B處去

maocp

將*GET,TWSIVAL(TWSINODE),NODE,TWSINODE,S,TWSI
抓取Tsai-Wu strength failure criterion.

改為*GET,TWSIVAL(TWSINODE),NODE,TWSINODE,FAIL,TWSR
抓 取Inverse of Tsai-Wu Strength Ratio Index failure criterion.
就行了，可是

*GET,TWSIVAL(TWSINODE),NODE,TWSINODE,FAIL,TWSI
抓取Tsai-Wu Failure Criterion Strength Index failure criterion.
或者*GET,TWSIVAL(TWSINODE),NODE,TWSINODE,S,TWSI
截取出來的值並非實際的值，這點就很奇怪了，差了E40以上，不知道跑去算甚麼東西了