固有値解析¶

本解析の実施には、tutorial/15_eigen_spring のデータを用います。

解析対象¶

解析対象は、前節の非線形静解析(超弾性その2)と同一のスプリングです。

項目	内容	備考
解析の種別	固有値解析	!SOLUTION, TYPE=EIGEN
節点数	78,771
要素数	46,454
要素タイプ	10節点四面体二次要素	!ELEMENT,TYPE=342
材料物性名	MAT1	!MATERIAL,NAME=MAT1
境界条件	拘束
行列解法	直接法	!SOLVER,METHOD=MUMPS

解析内容¶

図4.4.1に示す拘束面の変位を拘束し、5次までの固有値解析を実施します。解析制御データを以下に示します。

スプリングの形状

図4.4.1　スプリングの形状

解析制御データ `spring.cnt`¶

#  Control File for FISTR
## Analysis Control
!VERSION
 3
!SOLUTION, TYPE=EIGEN
!EIGEN
 5, 1.0E-8, 60
!WRITE,RESULT
!WRITE,VISUAL
## Solver Control
### Boundary Condition
!BOUNDARY
 XFIX, 1, 1, 0.0
 YFIX, 2, 2, 0.0
 ZFIX, 3, 3, 0.0
### Material
# define in mesh file
### Solver Setting
!SOLVER,METHOD=MUMPS
## Post Control
!VISUAL,method=PSR
!surface_num=1
!surface 1
!output_type=VTK
!END

解析手順¶

FrontISTRの実行コマンド fistr1 を実行します。

cd FrontISTR/tutorial/15_eigen_spring
fistr1 -t 4

(4スレッドで実行)

解析結果¶

解析結果データファイルspring.res.0.3を用いて、3次の振動モード（スプリングのy方向圧縮伸長）をREVOCAP_PrePostで作成して図4.15.1に示します。変形倍率を10としています。また、解析結果の数値データとして、解析結果ログファイルに出力された固有振動数リストを以下に示します。

スプリングの3次振動モード

図4.15.1　スプリングの3次振動モード

解析結果ログ `0.log`¶

********************************
*RESULT OF EIGEN VALUE ANALYSIS*
********************************

### NUMBER OF ITERATIONS =       20
### TOTAL MASS =   3.4184E-06

       PERIOD     FREQUENCY  PARTICIPATION FACTOR             EFFECTIVE MASS
  NO.  [Sec]      [HZ]       X          Y          Z          X          Y          Z
  ---  ---------  ---------  ---------  ---------  ---------  ---------  ---------  ---------
    1  2.245E-03  4.454E+02  1.029E+00 -8.894E-02 -7.052E-01  1.301E-06  9.718E-09  6.109E-07
    2  2.239E-03  4.465E+02  6.969E-01  1.075E-01  1.011E+00  6.129E-07  1.460E-08  1.289E-06
    3  1.100E-03  9.087E+02  4.862E-03  1.236E+00 -7.917E-02  4.007E-11  2.591E-06  1.062E-08
    4  1.014E-03  9.858E+02 -2.965E-02  3.385E-01 -6.782E-03  9.823E-10  1.280E-07  5.138E-11
    5  5.525E-04  1.810E+03 -5.867E-01 -5.386E-02 -7.659E-01  2.838E-07  2.391E-09  4.837E-07