作者harry750265 (O咖稱)
看板Cad_Cae
標題[問題] ANSYS 關於壓力容器的分析
時間Wed Oct 21 22:45:35 2009
最近在分析壓力容器，照道理來說這CASE因該是線性分析的。
可是在求解時卻會跑去非線性分析。而且還會求解超過十二小時 = =a
想請問這種案例真的是非線性求解麼???
我看網路上有些人說如果要做一個非線性分析時，
一般都會先用線性求解看看，看答案對不對，會不會差太遠，
才去決定要不要做非線性求解。
那這樣說不就代表可以強迫我的case去做線性求解???
還是根本就是我設定有錯誤，才會跑去非線性???
跪求板上高手幫我看看。
為此事稿的頭昏腦脹的新手 敬上
以下為APDL:
finish
/clear
/PREP7
ET,1,solid186
!=== 低碳鋼 AISI1017===========
MP,EX,1,205E3
MP,NUXY,1,0.29
!===============================
K,1,-1200,306,0,
K,2,-1200,300,0,
K,3,-1375,0,0,
K,4,-1369,0,0,
wpoff,-1200,0,0
CSWPLA,11,1,1.748571429,0
L, 1, 3
L, 2, 4
LSTR,1,2
LSTR,3,4
!==========圓頂面================
FLST,2,4,4
FITEM,2,1
FITEM,2,4
FITEM,2,2
FITEM,2,3
AL,P51X
!==========圓柱面================
RECTNG,0,1200,300,306,
!==========環狀複製==============
NUMMRG,KP, , , ,LOW
FLST,2,2,5,ORDE,2
FITEM,2,1
FITEM,2,-2
FLST,8,2,3
FITEM,8,3
FITEM,8,4
VROTAT,P51X, , , , , ,P51X, ,-180,1,
!==========重合座標===============
WPCSYS,-1,0
!=================================
wprot,0,-90,0
CYL4,0,0,325/2,90,325/2+6,180,906
!===========工作平面歸正======================
wprot,0,90,0
!=================================
CSYS,0
k,,0,906+106,0
k,,0,906+100,0
wpoff,0,906,0
CSWPLA,11,1,106/168.5,0
L, 16, 19
L, 17, 20
LSTR,19,20
LSTR,16,17
!=======================
FLST,2,4,4
FITEM,2,23
FITEM,2,30
FITEM,2,31
FITEM,2,32
AL,P51X
NUMMRG,KP, , , ,LOW
!============================
FLST,2,1,5,ORDE,1
FITEM,2,17
FLST,8,2,3
FITEM,8,19
FITEM,8,20
VROTAT,P51X, , , , , ,P51X, ,90, ,
VPLOT
!==============整體重合===============
/VIEW,1,1,2,3
/ANG,1
/REP,FAST
FLST,2,4,6,ORDE,2
FITEM,2,1
FITEM,2,-4
VOVLAP,P51X
!===========刪除多餘體=====================
FLST,2,2,6,ORDE,2
FITEM,2,6
FITEM,2,-7
VDELE,P51X, , ,1
!============映射============================
CSYS,0
FLST,3,5,6,ORDE,4
FITEM,3,1
FITEM,3,5
FITEM,3,8
FITEM,3,-10
VSYMM,Z,P51X, , , ,0,0
FLST,3,10,6,ORDE,2
FITEM,3,1
FITEM,3,-10
VSYMM,X,P51X, , , ,0,0
!==========================================
NUMMRG,KP, , , ,LOW
!===================內壓 0.5 ~ 7.4===========
FLST,2,24,5,ORDE,24
FITEM,2,4
FITEM,2,7
FITEM,2,14
FITEM,2,20
FITEM,2,23
FITEM,2,27
FITEM,2,29
FITEM,2,-30
FITEM,2,33
FITEM,2,-34
FITEM,2,37
FITEM,2,46
FITEM,2,52
FITEM,2,57
FITEM,2,62
FITEM,2,67
FITEM,2,72
FITEM,2,76
FITEM,2,78
FITEM,2,81
FITEM,2,87
FITEM,2,91
FITEM,2,93
FITEM,2,96
SFA,P51X,1,PRES,0.5 !初始壓力0.5，爆破壓力7.4
!=========================================
ALLSEL
!==========mesh============
SMRT,8
MSHAPE,1,3D
MSHKEY,0
FLST,5,20,6,ORDE,2
FITEM,5,1
FITEM,5,-20
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!===============================
!========求解收斂設定=================
ALLSET
FINISH
/SOL
ANTYPE,0
NLGEOM,1 !打開大變形
求解
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.133.13.134
1F：推 shin:有塑性就是非線性 如果不考慮塑性變形那就不用開大變形 10/22 00:40
2F：→ harry750265:可是不開大變形，求解時會說有NODE位移xx超過極限xx 10/22 01:15
3F：→ harry750265:這該怎麼辦?? 10/22 01:15
4F：推 quaintness:幾點建議...debug階段完全不必要建個35萬node的模型 10/22 02:13
5F：→ quaintness:Constraint應該是有問題,材料線性 NLGEOM,0就線性求解 10/22 02:14
6F：→ SCHLAGER:有大的變形或是大的旋轉量也都是非線性 10/23 13:24
7F：→ mises:有確定單位的一致性嗎?你給的壓力是否真的會導致如此大變形? 10/23 22:44
8F：推 lalilaca:你在CAE論壇的發問,我有幫你跑過了..小變形情況下 10/24 01:24
9F：→ lalilaca:並沒有出現你所說過大位移的問題.. 10/24 01:25
10F：→ harry750265:D感謝樓上大大，請問你跑的時候有貼上 NLGEOM,1 這句? 10/24 05:32
11F：→ harry750265:我今天在MODEL上加上FIX兩端，就很順利的跑完了= =a 10/24 05:33
12F：→ harry750265:像是類似這種case是一定要FIX嚜??? 10/24 05:34
13F：→ harry750265:因為我看書上 他們是有FIX某端 10/24 05:34
14F：→ harry750265:在問一下la大你跑的時候大概花了多久呢?是非線性求解? 10/24 05:36
15F：→ harry750265:WOW~LA大你居然還附圖!!真是太貼心了!! 10/24 05:42
16F：→ harry750265:你PO的結果跟我FIX兩端圓頭與圓筒的線的結果是一樣的 10/24 05:44
17F：推 lalilaca:沒阿..我跑的時候,把大變形關掉.. 10/24 18:52
18F：→ lalilaca:跑小變形的分析,因為你用SOLID186且元素不少..大概花了.. 10/24 18:54
19F：→ lalilaca:5-10分鐘吧.. 10/24 18:54
20F：→ harry750265:感謝黑~!我也順利解決了! 10/24 20:58