test

modal.hdf

https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca2017

卒論日誌

日付時間帯作業時間内容立会
3年10/1116:00-17:001.0創造工房後藤さん
10/1814:30-16:302.0創造工房後藤さん
10/2112:00-15:203.5創造工房高橋さん
10/2813:00-16:303.5創造工房高橋さん
11/113:00-16:303.5創造工房後藤さん 他
11/612:00-13:001.0創造工房
11/814:30-16:001.5創造工房後藤さん
11/2214:30-16:001.5創造工房後藤さん高橋さん
12/614:30-16:001.5創造工房後藤さん高橋さん
12/2014:30-16:001.5創造工房加藤さん、4年生のみなさん
1/1014:30-16:001.5創造工房後藤さん
1/1514:30-16:001.5創造工房後藤さん
1/3114:30-19:004.5創造工房後藤さん
2/219:00-24:005.0創造工房課題後藤さん
2/714:30-14:500.3創造工房後藤さん
4年7/1716:00-18:002.0卒論課題
7/2014:00-18:004.0ゼミ・卒論課題
7/2024:00-25:451.75卒論課題
7/2111:30-13:001.5卒論課題
7/2216:00-18:002.0卒論課題
7/2414:00-17:003.0卒論課題
7/2714:00-17:003.0ゼミ
8/314:00-17:003.0ゼミ
8/2414:00-15:301.5ゼミ
9/3014:30-18:003.5卒論課題
10/215:00-18:003.0卒論課題
10/514:00-18:004.0ゼミ・卒論課題
10/615:00-16:301.5卒論課題
10/2213:00-18:005.0ゼミ・卒論課題
10/2611:00-18:007.0卒論課題
10/2712:00-20:008.0卒論課題
10/3013:30-16:303.0ゼミ・卒論課題
11/216:00-18:302.5卒論課題
11/420:00-21:301.5卒論課題
11/914:00-18:004.0卒論課題
11/1113:00-18:005.0ゼミ・卒論課題
11/1613:00-14:301.5卒論課題
11/1814:30-18:304.0卒論課題
11/1915:00-16:001.0卒論課題
11/2013:00-17:004.0ゼミ・卒論課題
11/2013:00-19:006.0ゼミ・卒論課題・チューリップ
11/3016:00-17:301.5卒論課題
12/113:00-16:003.0卒論課題・掃除
12/316:00-17:301.5卒論課題
12/413:00-18:005.0ゼミ・卒論課題
12/715:30-18:002.5卒論課題
12/918:00-21:003.0卒論課題
12/1017:00-18:001.0卒論課題
12/1113:00-21:308.5現場見学・ゼミ・卒論課題
12/1519:30-23:003.5卒論課題
12/1617:00-23:306.5卒論課題
12/179:00-16:307.5卒論課題
12/1813:00-18:305.5ゼミ・卒論課題
12/1914:00-17:303.5中間発表準備
12/2014:30-15:301.0中間発表準備
12/2216:00-18:002.0中間発表準備
12/2314:00-25:3011.5中間発表準備
12/2419:00-24:005.0中間発表準備
12/259:00-13:004.0中間発表
1/813:30-18:305.0ゼミ
1/1013:00-16:303.5卒論課題
1/1121:00-26:305.5卒論課題
1/124:00-8:304.5卒論課題
1/1220:00-26:306.5卒論課題
1/1313:00-18:005.0卒論課題
1/1418:00-19:301.5卒論課題
1/1513:00-16:003.0ゼミ・卒論課題
1/1622:00-30:008.0卒論課題
1/182:00-3:301.5卒論課題
1/1815:00-18:303.5卒論課題
1/1821:30-28:006.5卒論課題
1/1916:00-17:301.5卒論課題
1/2014:00-20:006.0卒論課題
1/2213:00-23:0010.0ゼミ・卒論課題
1/2313:00-18:005.0卒論課題
1/2416:00-24:008.0卒論課題
1/2510:30-18:308.0卒論課題

当面の課題

20/5/11(後藤)

5/31

指示された課題が何を指しているかそもそも理解できず、今までなにもせずにいてしまいました。申し訳ありません。 ようやく何をやればいいのかがなんとなくわかったため、及川さんが日誌に載せてくださったものを見たりしながらSalome解析を取り組み始めます。

6/21

及川さんの日誌やYouTubeを参考に片持ち梁の振動解析を試みていますが、④材料物性値 で詰まっています。正直Aster Studyの最初のadd stage with assistant - modal analysis がSalome2018には無くて、これかな?と勘で進めています。怪しいです。 Aster Study以前にジオメトリから、振動解析で載荷面はいるのか などわかっていません…。

7/15

ソリッド単純梁.100/10/20.ヤング率6000.ポアソン比0.4.密度3.8e-0.7.メッシュサイズ1.1次メッシュ

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/aoyama/tamj_0713_soli_1.png

片持ち梁 梁要素 振動

7/17

片持ち梁の梁要素.振動

Salome実験値と理論値が一致しない。Salomeの解析結果が間違っているのか、理論式への代入が間違っているのか。理論式の扱い方の理解も怪しい。

ちなみに上記は梁要素の片持ち梁だが、ソリッド要素片持ち梁とソリッド要素の単純梁のSalomeでの実験値も得られている(うまくいっているかは置いといて。そして多分うまくいっていない。)

7/20

7/17分表の訂正.

e-0.7は10の-7乗という意味らしい.自然対数ではない.

y軸、z軸の2次(?)振動までをSalomeと理論値で比較することができた! 3次以上は単純な曲げを超える(ねじれの要素も含まれる)ので無視した。

片持ち梁 ソリッド要素 振動

7/20

Salome実験値λ理論値
1(y軸曲げ)-20.81881.87520.2962
2(z軸曲げ)-39.95311.87540.593
3(y軸曲げ)-124.44.694127.204
4(ねじれ)-145.725
5(z軸曲げ)-213.8414.694254.407
6-316.751
7-326.558
8-438.839
9-507.218
10-591.374

二次のz軸曲げがSalomeと理論値で誤差がある?

単純梁 ソリッド要素 振動

7/20

Salome実験値λ理論値
1(z軸曲げ)-128.947π113.979
2-157.345
3(y軸曲げ)-229.934227.915
4-350.212
5(z軸曲げ)-442.855455.830
6-499.825
7-522.297
8-605.117
9-634.469
10-653.725

単純梁 梁要素 振動

7/22

  !<S> 例外ユーザーは上げたけれども、not interceptee 。             !
  !ベースはfermeesである。                                     !
  !例外のタイプ:エラー!
  !                                                            !
  ! GROUP_MA Group_2はメッシュの一部を形成しない:メッシュ

Group_2ってどれだろう.

7/24

Salome実験値λ理論値
1(z軸曲げ?)-129.164π113.958
2-258.329
3-278.293
4-356.047
5-556.655
6-628.307
7-697.994
8-712.094
9-835.154
10-1113.86

一個目のほかソリッド要素と梁要素で結果がずれた.

↓曲げ振動の固有振動数の式にパラメータを代入して得られた値

λbhf
π201056.979
π1020113.958
2010227.915
1020455.830
2010512.809
10201025.52

これが上2の表にうまくはまってくれるといいが...

一本のアーチ 梁要素 振動

かじか橋の座標を使って.

これまでと同様に断面を (HY,HZ)=(10,20) にした.

7/24

エラーが出た.

  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  ! <S> Exception user raised but not interceptee.                               !
  ! The bases are fermees.                                                       !
  ! Type of the exception: error                                                 !
  !                                                                              !
  ! Program error.                                                               !
  !                                                                              !
  ! Condition not met:                                                           !
  !     .false.                                                                  !
  ! File                                                                         !
  ! /home/A21173/smeca/yamm_build/V2018.0.1_public_build/V2018_public/tools/src/ !
  ! Code_aster_stable-v136_smeca/bibfor/debug/contex_param                       !
  ! F90, line 147                                                                !
  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

7/27

→解決.

Aster StudyのPre Analysis > ASSEMBLAGEのMaterial fieldとStructural el...とLoadにチェックを入れていなかった.

あとは境界条件のDRY=0を外した.

断面を(HY,HZ)=(10,20)から(HY,HZ)=(0.01,0.02)に変更した.

Salome実験値λ理論値
1-4.02135
2-11.3476
3-13.6050
4-22.5116
5-30.9090
6-37.4207
7-56.0714
8-58.3772
9-78.4644
10-90.7733

かじか橋 床版シェル要素とその他梁要素 振動解析

8/3

及川さんからかじか橋のSalomeデータを頂いた.

部材中央に拘束線をつけて解析してみる.

9/30

シェル要素の床版と梁要素のその他部材を結合させたい.

10/9

9/30 続き

 <S> 例外ユーザーが発生しましたが、傍受者ではありません。 !

   ! ベースはフェルミーです。 !
   ! 例外のタイプ:エラー!
   ! !
   ! タイプPOUTREの特性を、そうでないメッシュM155に割り当てます。
   ! このタイプ。 !
   ! !
   ! 助言:                                                                        !
   ! メッシュM155のコマンドAFFE_MODELEの結果を確認してください。

10/22

10/9続き

  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  ! <S> Exception user raised but not interceptee.                !
  ! The bases are fermees.                                        !
  ! Type of the exception: error                                  !
  !                                                               !
  !  the computation of the option:  CARA_SECT_POUT3              !
  !  is not possible for any the types of elements of the LIGREL. !
  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

10/26

10/22続き

AsterStudyでRunは通過した!

しかしParavisを見るとシェルの床版がビーム部材とくっついていない泣

10/27

10/26続き

geometoryを少し変えてやり直してみたら長文のエラー.

でも及川さんいわく,このエラー内容は俺も出たことあるとのこと.通過点だ.

またあした.

  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  ! <S> Exception user raised but not interceptee.                                   !
  ! The bases are fermees.                                                           !
  ! Type of the exception: error                                                     !
  !                                                                                  !
  !   Le shift=20.6546                                                               !
  !   utilisé pour construire la matrice dynamique coïncide avec une valeur propre ! !
  !   Avec l'option 'CENTRE', ce shift vaut FREQ,                                    !
  !   Avec l'option 'BANDE', c'est le milieu de la bande sélectionnée,               !
  !   Avec l'option 'PLUS_PETITE' ou 'TOUT', il prend la valeur 0.                   !
  !                                                                                  !
  !   Malgré la stratégie de décalage du shift, cette matrice dynamique reste        !
  !   numériquement singulière.                                                      !
  !                                                                                  !
  !   -> Risque :                                                                    !
  !   Cette matrice étant abondamment utilisée pour résoudre des systèmes linéaires  !
  !   à chaque itération du processus modal, cette quasi singularité peut fausser    !
  ! les résultats                                                                    !
  !   (mauvais conditionnement matriciel).                                           !
  !                                                                                  !
  !   -> Conseils :                                                                  !
  !   La structure analysée présente probablement des modes de corps rigide.         !
  !                                                                                  !
  !     * si aucun mode de corps rigide n'était attendu :                            !
  !   Vous pouvez modifier les paramètres du solveur linéaire (par exemple METHODE   !
  ! ou NPREC),                                                                       !
  !   ou ceux de l'algorithme de décalage (PREC_SHIFT, NMAX_ITER_SHIFT et            !
  ! SEUIL_FREQ)                                                                      !
  !   pour vérifier qu'il s'agit bien d'une singularité et non d'un problème         !
  ! numérique ponctuel.                                                              !
  !   Si c'est une singularité, vérifiez la mise en donnée du problème :             !
  !   conditions aux limites, maillage (présence de noeuds / mailles orphelin(e)s),  !
  ! unités, ...                                                                      !
  !                                                                                  !
  !    * si ces modes étaient attendus et que vous ne voulez pas les calculer :      !
  !   Utilisez l'option 'BANDE' avec une borne inférieure suffisamment positive (par !
  ! exemple 1.e-1).                                                                  !
  !    * si ces modes étaient attendus et que vous voulez les calculer :             !
  !   - utilisez l'option 'BANDE' avec une borne inférieure légèrement négative (par !
  ! exemple -1.e-1).                                                                 !
  !   - utilisez la méthode 'TRI_DIAG' avec OPTION='MODE_RIGIDE'.                    !
  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

10/30

10/27続き

一旦簡単なモデル(単純梁)で,シェル床版とビーム縦桁をくっつけてみることにした.

1/20

LIAISON_SOLIDEで試みるもおなじみのエラー.

S> Exception user raised but not interceptee. !

  ! The bases are fermees.                                                           !
  ! Type of the exception: error                                                     !
  !                                                                                  !
  !   Le shift=-0.01                                                                 !
  !   utilisé pour construire la matrice dynamique coïncide avec une valeur propre ! !
  !   Avec l'option 'CENTRE', ce shift vaut FREQ,                                    !
  !   Avec l'option 'BANDE', c'est le milieu de la bande sélectionnée,               !
  !   Avec l'option 'PLUS_PETITE' ou 'TOUT', il prend la valeur 0.                   !
  !                                                                                  !
  !   Malgré la stratégie de décalage du shift, cette matrice dynamique reste        !
  !   numériquement singulière.                                                      !
  !                                                                                  !
  !   -> Risque :                                                                    !
  !   Cette matrice étant abondamment utilisée pour résoudre des systèmes linéaires  !
  !   à chaque itération du processus modal, cette quasi singularité peut fausser    !
  ! les résultats                                                                    !
  !   (mauvais conditionnement matriciel).                                           !
  !                                                                                  !
  !   -> Conseils :                                                                  !
  !   La structure analysée présente probablement des modes de corps rigide.         !
  !                                                                                  !
  !     * si aucun mode de corps rigide n'était attendu :                            !
  !   Vous pouvez modifier les paramètres du solveur linéaire (par exemple METHODE   !
  ! ou NPREC),                                                                       !
  !   ou ceux de l'algorithme de décalage (PREC_SHIFT, NMAX_ITER_SHIFT et            !
  ! SEUIL_FREQ)                                                                      !
  !   pour vérifier qu'il s'agit bien d'une singularité et non d'un problème         !
  ! numérique ponctuel.                                                              !
  !   Si c'est une singularité, vérifiez la mise en donnée du problème :             !
  !   conditions aux limites, maillage (présence de noeuds / mailles orphelin(e)s),  !
  ! unités, ...                                                                      !
  !                                                                                  !
  !    * si ces modes étaient attendus et que vous ne voulez pas les calculer :      !
  !   Utilisez l'option 'BANDE' avec une borne inférieure suffisamment positive (par !
  ! exemple 1.e-1).                                                                  !
  !    * si ces modes étaient attendus et que vous voulez les calculer :             !
  !   - utilisez l'option 'BANDE' avec une borne inférieure légèrement négative (par !
  ! exemple -1.e-1).                                                                 !
  !   - utilisez la méthode 'TRI_DIAG' avec OPTION='MODE_RIGIDE'.      

以前出ていた運営に連絡して下さいエラーは解消されたので

進歩といえば進歩か.

原因はなんでしょう.

Geometryの設定か,接合するnodeが多すぎる/距離が近すぎるか,かなぁ.

1/22

やはりGeometryの設定が違っていた.

それを正したらきちんとくっつけることができた!

しかしいったん接合点を少なくしてやってみたので次は点をフルでやってみる.

(Geometryで梁とシェルをコンパウンドして,梁の点がグループ化されてしまって困っている.たぶん縦桁の部分に改善の余地あり.)

1/23

ついに解析が通り,かじか橋のシェルと梁を接合させることができた!

paravisでの振動モードの見方がよくわからないので月曜日に及川さんに聞いてみよう.

単純梁 床版がシェル要素,縦桁が梁要素 解析

11/2-9

10/30続き

impossibility, the node beam N7 should carry the degree of freedom DX

やはり軸と平行方向に接合はできないのか?

11/13

続き

シェルにパーティションで節点を与えて梁要素との結合を試みる.

11/18

続き

the computation of the option: CARA_SECT_POUT3 is not possible for any the types of elements of the LIGREL.

のエラーで足踏み.

11/27

続き

梁は梁の点(V1-10),シェルはシェルの点(vertex_1-10)に分けてやってみる.

そしたら

! one should give one node in the GROUP_NO: V2

ノードが指定されていないとエラー.

ジオメトリでグループを作りなおしてみるが解消せず.

う〜.

12/4

シェルとビームの結合を諦める...

ソリッド縦桁とビームアーチをくっつけてみることにした.

かじか橋 床版と縦桁がソリッド要素,その他梁要素 振動解析

12/11

モデルを作る.

12/15

床版と縦桁はソリッドなのでgeometryでfuseしてくっつけた.

するとメッシュが切れなかった.

12/16

縦桁に接合の節点を1つずつ作った. 横桁も節点1つずつ作った方がいいかも(まだやってない)

部材ごとにメッシュを切って,メッシュでコンパウンドした.

ソリッドの床版と縦桁もAsterStudyで接合できる.

部材ごとにメッシュを切ることで,まとめて切るよりエラーが起きにくい.

12/17

AsterStudyでRunは通ったが,Paravisでみると四面体要素と梁要素がくっついていなかった泣.

結合の節点で,四面体要素の節点と梁要素の節点同士を結合させるやりかたがある?(LIAISON_SOLIDE)

いったん中間発表準備に入るが,年明けからはそれを試してみたい.

線上に点を打つ

http://www.str.ce.akita-u.ac.jp/~gotouhan/2020/aoyama/dotUCHIKATA.png

こうするとSalomeが,点が線上にあるとみなしてくれる.

創造工房

10/18課題

たわみ理論値 6.66667

メッシュ0.5(佐藤)

たわみ実験値相対誤差要素数
6.57938-1.4%394121

メッシュ1(吉田)

たわみ実験値相対誤差要素数
6.443-3.4%72278

メッシュ2(梅田)

たわみ実験値相対誤差要素数
5.73825-13.9%11817

メッシュ4(青山)

たわみ実験値相対誤差要素数
4.94280-25.858%2862

メッシュ8(小川)

たわみ実験値相対誤差要素数
4.404-39.4%897

10/25課題

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/aoyama/kadai1025.png

メッシュ数要素数変位相対誤差作成者
88974.404-39.4%小川
428624.943-25.9%青山
2118175.738-13.9%梅田
1.8117585.737-13.98梅田
1.5234176.121-8.23%梅田
1.4410966.316-5.30%青山
1.2655756.443-3.35%小川
1.0722786.443-3.4%吉田
0.9717186.431-3.53吉田
0.8721016.437-3.48青山
0.71309166.478-2.87%佐藤
0.53941216.579-1.4%佐藤

考察

メッシュ0.8,0.9,1.0で期待通りの数値が出なかった。原因はSalomeでメッシュを作る際にきれいな値でメッシュを割り切れず、切り上げされたなどが考えられる。他はグラフからも期待通りの結果が得られた。

11/1 課題

理論値 0.41667

メッシュ数要素数変位相対誤差作成者
2205260.37063-11.05佐藤
1.5503590.41195-1.13青山
1.0623600.41651-0.037梅田
0.81999680.41646-0.050吉田
0.53226870.428772.9小川

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/aoyama/tangraf.png

考察

試験体ごとに最適なメッシュ長さ(要素数)があり、今回はメッシュ1.0が最も理論値に近い値になった。 プロット数が5つと少ないため、グラフが角ばった形状になった。より多くのメッシュ長さで試験をすることによってよりなめらかなグラフを得られるだろう。

11/8

ティモシェンコ理論値 0.491668 初等理論値 0.41667

メッシュ数要素数変位相対誤差作成者
2.0205260.47657183-3.07佐藤
1.8280530.492666860.203佐藤
1.5503590.500102251.72梅田
1.3574550.503843562.48青山
1.1832780.503381362.38小川
1.0623600.505768732.87青山
0.9879530.510275383.78小川
0.80.514543710.465吉田
0.53226870.521495866.07小川

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/aoyama/tanjun1113.png

考察

グラフを見ると、実験値と初等・ティモシェンコそれぞれの理論値を比較して、ティモシェンコの値の方が実験値と近い値になった。これは単純梁が異方性であることを示している。メッシュ1.8での相対誤差が最小になったため、この単純梁はメッシュ1.8で試験するのがよい。

11/22 課題

拡大係数要素数変位相対誤差作成者
0.836575.5734-16.399青山
0.742455.5717-16.425吉田
0.672575.9119-11.322梅田
0.5137416.03154-9.527小川
0.4233496.26413-6.058佐藤
0.3574316.42253-3.662全員

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/aoyama/macktn.png

12/6 課題

拡大係数要素数変位相対誤差作成者
0.81225156.917492.6%吉田
1.0654906.88772.2%梅田
1.2366936.85192251.6%梅田
1.4233416.81171751.0%吉田
1.6153716.6336-1.6%小川
1.8124806.742580.014%佐藤
2.085896.66499-0.011%梅田

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/aoyama/1206k.png


トップ   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS