跳轉至

非線性靜力分析 (黏彈性)

本分析使用 tutorial/07_viscoelastic_cylinder 的資料。

分析對象

分析對象與非線性靜力分析(超彈性其1)相同,為圓棒的1/8模型。

項目 內容 備註 參照
分析類型 非線性靜力分析(彈塑性) !SOLUTION,TYPE=NLSTATIC
節點數 629
元素數 432
元素類型 8節點六面體一次元素 !ELEMENT,TYPE=361
材料性質名稱 MAT1 !MATERICAL,NAME=MAT1
材料性質 ELASTIC,VISCOELASTIC !ELASTIC !VISCOELASTIC
邊界條件 約束, 強制位移
矩陣解法 CG/SSOR !SOLVER,METHOD=CG,PRECOND=1

分析內容

對圓棒施加軸向拉伸位移,執行應力鬆弛分析。分析控制資料如下所示。

分析控制資料 cylinder.cnt

#  Control File for FISTR
## Analysis Control
!VERSION
 3
!SOLUTION, TYPE=NLSTATIC
!WRITE,VISUAL
!WRITE,RESULT
## Solver Control
### Boundary Conditon
!BOUNDARY, GRPID=1
 LOADS, 3, 3, -7.0
 FIX,   3, 3, 0.0
 XSYMM, 1, 1, 0.0
 YSYMM, 2, 2, 0.0
### STEP
!STEP, TYPE=VISCO, CONVERG=1.0e-5
 0.2, 2.0
 BOUNDARY, 1
### Material
!MATERIAL, NAME=MAT1
!ELASTIC
 206900.0, 0.29
!VISCOELASTIC
 0.5, 1.0
### Solver Setting
!SOLVER,METHOD=CG,PRECOND=1,ITERLOG=YES,TIMELOG=YES
 10000, 1
 1.0e-8, 1.0, 0.0
## Post Control
!VISUAL,metod=PSR
!surface_num=1
!surface 1
!output_type=VTK
!END

分析步驟

執行 FrontISTR 的執行指令 fistr1

cd FrontISTR/tutorial/07_viscoelastic_cylinder
fistr1 -t 4

(以4個執行緒執行)

分析結果

使用REVOCAP_PrePost製作加入 von Mises 應力等值線的變形圖,並示於圖4.7.1。 這是2秒後(第10步)的分析結果。此外,以下列出分析結果日誌檔案的一部分,作為分析結果的數值資料。

變形與 von Mises 應力的分析結果

圖4.7.1 變形與 von Mises 應力的分析結果

分析結果日誌 0.log

 fstr_setup: OK
#### Result step=     0
 ##### Local Summary @Node    :Max/IdMax/Min/IdMin####
 //U1    0.0000E+00         1  0.0000E+00         1
 //U2    0.0000E+00         1  0.0000E+00         1
 //U3    0.0000E+00         1  0.0000E+00         1
 //E11   0.0000E+00         1  0.0000E+00         1
 //E22   0.0000E+00         1  0.0000E+00         1
 //E33   0.0000E+00         1  0.0000E+00         1
 //E12   0.0000E+00         1  0.0000E+00         1
 //E23   0.0000E+00         1  0.0000E+00         1
 //E31   0.0000E+00         1  0.0000E+00         1
 //S11   0.0000E+00         1  0.0000E+00         1
 //S22   0.0000E+00         1  0.0000E+00         1
 //S33   0.0000E+00         1  0.0000E+00         1
 //S12   0.0000E+00         1  0.0000E+00         1
 //S23   0.0000E+00         1  0.0000E+00         1
 //S31   0.0000E+00         1  0.0000E+00         1
 //SMS   0.0000E+00         1  0.0000E+00         1
 ##### Local Summary @Element :Max/IdMax/Min/IdMin####
 //E11   0.0000E+00         1  0.0000E+00         1
 //E22   0.0000E+00         1  0.0000E+00         1
 //E33   0.0000E+00         1  0.0000E+00         1
 //E12   0.0000E+00         1  0.0000E+00         1
 //E23   0.0000E+00         1  0.0000E+00         1
 //E31   0.0000E+00         1  0.0000E+00         1
 //S11   0.0000E+00         1  0.0000E+00         1
 //S22   0.0000E+00         1  0.0000E+00         1
 //S33   0.0000E+00         1  0.0000E+00         1
 //S12   0.0000E+00         1  0.0000E+00         1
 //S23   0.0000E+00         1  0.0000E+00         1
 //S31   0.0000E+00         1  0.0000E+00         1
 //SMS   0.0000E+00         1  0.0000E+00         1
 ##### Global Summary @Node    :Max/IdMax/Min/IdMin####
 //U1    0.0000E+00         1  0.0000E+00         1
 //U2    0.0000E+00         1  0.0000E+00         1
 //U3    0.0000E+00         1  0.0000E+00         1
 //E11   0.0000E+00         1  0.0000E+00         1
 //E22   0.0000E+00         1  0.0000E+00         1
 //E33   0.0000E+00         1  0.0000E+00         1
 //E12   0.0000E+00         1  0.0000E+00         1
 //E23   0.0000E+00         1  0.0000E+00         1
 //E31   0.0000E+00         1  0.0000E+00         1
 //S11   0.0000E+00         1  0.0000E+00         1
 //S22   0.0000E+00         1  0.0000E+00         1
 //S33   0.0000E+00         1  0.0000E+00         1
 //S12   0.0000E+00         1  0.0000E+00         1
 //S23   0.0000E+00         1  0.0000E+00         1
 //S31   0.0000E+00         1  0.0000E+00         1
 //SMS   0.0000E+00         1  0.0000E+00         1
 ##### Global Summary @Element :Max/IdMax/Min/IdMin####
 //E11   0.0000E+00         1  0.0000E+00         1
 //E22   0.0000E+00         1  0.0000E+00         1
 //E33   0.0000E+00         1  0.0000E+00         1
 //E12   0.0000E+00         1  0.0000E+00         1
 //E23   0.0000E+00         1  0.0000E+00         1
 //E31   0.0000E+00         1  0.0000E+00         1
 //S11   0.0000E+00         1  0.0000E+00         1
 //S22   0.0000E+00         1  0.0000E+00         1
 //S33   0.0000E+00         1  0.0000E+00         1
 //S12   0.0000E+00         1  0.0000E+00         1
 //S23   0.0000E+00         1  0.0000E+00         1
 //S31   0.0000E+00         1  0.0000E+00         1
 //SMS   0.0000E+00         1  0.0000E+00         1

 ...

#### Result step=    10
 ##### Local Summary @Node    :Max/IdMax/Min/IdMin####
 //U1    0.0000E+00         1 -7.4530E-01        91
 //U2    0.0000E+00         1 -7.4530E-01        88
 //U3    0.0000E+00         1 -7.0000E+00        38
 //E11  -1.2136E-01        38 -1.2719E-01         7
 //E22  -1.2136E-01        50 -1.2719E-01        13
 //E33   2.3670E-01         7  2.2951E-01        50
 //E12   1.2005E-03        53 -1.2295E-03        10
 //E23   1.2987E-04        72 -2.5210E-03        83
 //E31   1.2987E-04        64 -2.5210E-03        96
 //S11   1.3771E+02        13 -1.3366E+02        50
 //S22   1.3771E+02         7 -1.3366E+02        38
 //S33   6.3193E+04         7  5.9532E+04        50
 //S12   3.5416E+01        53 -3.5955E+01        10
 //S23   1.0704E+01        35 -2.1669E+02        83
 //S31   1.0704E+01        31 -2.1669E+02        96
 //SMS   6.3093E+04         7  5.9628E+04        50
 ##### Local Summary @Element :Max/IdMax/Min/IdMin####
 //E11  -1.2167E-01       759 -1.2687E-01       354
 //E22  -1.2167E-01       768 -1.2687E-01       363
 //E33   2.3619E-01       354  2.3001E-01       768
 //E12   1.1247E-03       765 -1.1492E-03       360
 //E23  -8.9560E-06       758 -2.4072E-03       552
 //E31  -8.9560E-06       771 -2.4072E-03       543
 //S11   1.0437E+02       363 -1.0216E+02       768
 //S22   1.0437E+02       354 -1.0216E+02       759
 //S33   6.2914E+04       354  5.9795E+04       768
 //S12   3.3111E+01       765 -3.3459E+01       360
 //S23  -7.6469E-01       758 -2.0669E+02       552
 //S31  -7.6469E-01       771 -2.0669E+02       543
 //SMS   6.2843E+04       354  5.9864E+04       768
 ##### Global Summary @Node    :Max/IdMax/Min/IdMin####
 //U1    0.0000E+00         1 -7.4530E-01        91
 //U2    0.0000E+00         1 -7.4530E-01        88
 //U3    0.0000E+00         1 -7.0000E+00        38
 //E11  -1.2136E-01        38 -1.2719E-01         7
 //E22  -1.2136E-01        50 -1.2719E-01        13
 //E33   2.3670E-01         7  2.2951E-01        50
 //E12   1.2005E-03        53 -1.2295E-03        10
 //E23   1.2987E-04        72 -2.5210E-03        83
 //E31   1.2987E-04        64 -2.5210E-03        96
 //S11   1.3771E+02        13 -1.3366E+02        50
 //S22   1.3771E+02         7 -1.3366E+02        38
 //S33   6.3193E+04         7  5.9532E+04        50
 //S12   3.5416E+01        53 -3.5955E+01        10
 //S23   1.0704E+01        35 -2.1669E+02        83
 //S31   1.0704E+01        31 -2.1669E+02        96
 //SMS   6.3093E+04         7  5.9628E+04        50
 ##### Global Summary @Element :Max/IdMax/Min/IdMin####
 //E11  -1.2167E-01       759 -1.2687E-01       354
 //E22  -1.2167E-01       768 -1.2687E-01       363
 //E33   2.3619E-01       354  2.3001E-01       768
 //E12   1.1247E-03       765 -1.1492E-03       360
 //E23  -8.9560E-06       758 -2.4072E-03       552
 //E31  -8.9560E-06       771 -2.4072E-03       543
 //S11   1.0437E+02       363 -1.0216E+02       768
 //S22   1.0437E+02       354 -1.0216E+02       759
 //S33   6.2914E+04       354  5.9795E+04       768
 //S12   3.3111E+01       765 -3.3459E+01       360
 //S23  -7.6469E-01       758 -2.0669E+02       552
 //S31  -7.6469E-01       771 -2.0669E+02       543
 //SMS   6.2843E+04       354  5.9864E+04       768