当面の課題

20/5/11(後藤)

20/5/18

20/5/21

20/5/28

20/6/5

20/6/12

20/6/19

20/6/26

20/6/30

20/7/10

20/7/15

20/7/21

20/7/24

20/8/7

20/10/4

20/10/7

メッシュサイズ変位
60.000154
70.00015381
80.00015373
90.00015379
100.00015365
110.00015372
120.0001538

20/10/8

20/10/9

20/10/12

メッシュサイズ変位
60.000157
70.0001568
80.0001567
90.0001568
100.0001566
110.0001567
120.0001568

20/10/14

20/10/16

メッシュサイズ要素数変位
51458130.000157025
61440520.0001569962
7382440.0001567935
8370090.0001567133
9368040.0001567753
10361160.000156634
11360350.0001567057
12352930.0001567914
1595700.0001533965
2091480.0001537212
2587470.0001540633
3029720.0001336225
3523640.0001379303

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ko.png

20/10/19

20/10/20

20/10/21

メッシュサイズ要素数変位
145819480.0001126382
154431100.0001150793
203354450.0001133839
252556900.0001076753
271312960.0001036925
30832210.000097703671
40636830.000095167519
50549090.000093780481

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/syouban.png

20/10/28

メッシュサイズ要素数変位
36507320.0001525962
42538430.0001528185
52266150.000152853
10573900.0001527818
20135730.0001504562

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/kouz.png

20/10/29

20/10/30

メッシュサイズ要素数変位
36507320.0001525962
42538430.0001528185
52266150.000152853
10573900.0001527818
20135730.0001504562
3059310.0001270705
4039820.0001377993
5034110.0001174439

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/kouz1.png

20/11/2

20/11/3

メッシュサイズ要素数変位
145819480.0001126382
154431100.0001150793
203354450.0001133839
252556900.0001076753
271312960.0001036925
30832210.000097703671
40636830.000095167519
50549090.000093780481
60214130.000076680678
70186830.000073994244
80160840.000072255257
100140870.000070051691
15050700.000042707938

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/syouban1.png

S> Exception user raised but not interceptee. !

  ! The bases are fermees.                                                           !
  ! Type of the exception: error                                                     !
  !                                                                                  !
  ! Problème : la matrice est singulière ou presque singulière :                     !
  !   Lors de la factorisation de la matrice, on a rencontré un problème             !
  !   (pivot nul ou presque nul) à la ligne 37413 qui correspond au degré de liberté !
  ! donné ci-dessus.                                                                 !
  !   Le nombre de décimales "perdues" pour cette ligne est : 13.                    !
  !   Ce nombre de décimales perdues est lié au degré de singularité de la matrice.  !
  !   Plus il est grand, plus le système est singulier.                              !
  !   Quand on a perdu plus de 8 décimales, on estime que l'on a perdu 50% de la     !
  ! précision                                                                        !
  !   des nombres réels (qui ont 15 décimales environ).                              !
  !                                                                                  !
  ! Risques et conseils :                                                            !
  !    * Si la ligne correspond a un degré de liberté physique, il s'agit            !
  ! probablement d'un mouvement                                                      !
  !      de corps rigide mal bloqué.                                                 !
  !      Vérifiez les conditions aux limites.                                        !
  !      Si vous faites du contact, il ne faut pas que la structure ne "tienne" que  !
  ! par le contact.                                                                  !
  !      Vérifiez également les caractéristiques matériaux (module d'Young, ...).    !
  !                                                                                  !
  !    * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans    !
  ! doute d'une condition                                                            !
  !      limite redondante.                                                          !
  !      En particulier, il se peut que la relation linéaire surabondante provienne  !
  ! des conditions de contact.                                                       !
  !      Peut-être devriez vous exclure certains noeuds des conditions de contact    !
  !      (mots clés SANS_NOEUD et SANS_GROUP_NO).                                    !
  !                                                                                  !
  !    * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le       !
  ! solveur MUMPS                                                                    !
  !      car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas  !
  ! définies positives.                                                              !
  !                                                                                  !
  !    * Parfois, en parallèle, le critère de détection de singularité de MUMPS est  !
  ! trop pessimiste ! Il reste néanmoins souvent                                     !
  !      possible de faire passer le calcul complet en relaxant ce critère           !
  ! (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou                              !
  !      en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul    !
  ! sur moins de processeurs.                                                        !
  !                                                                                  !
  !    * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la  !
  ! fissure passe                                                                    !
  !      très près d'un noeud.                                                       !
  !      Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales),  !
  !      vous pouvez relancer le calcul en augmentant le nombre de décimales perdues !
  ! autorisé :                                                                       !
  !      mot-clé NPREC du mot clé facteur SOLVEUR.                                   !
  !      Sinon, contactez l'équipe de développement.  

20/11/6

20/11/8

メッシュサイズ要素数変位
101107650.00014746275
20274750.00014276444
21274190.000142449056
23194020.000142579938
25177920.0001425346
30100500.00012982425
4069120.000105412889

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/3sityuu.png

20/11/11

<ArretCPUError> <MECANONLINE9_1>                                                             !
  !                                                                                              !
  !    Arrêt par manque de temps CPU pendant les itérations de Newton, au numéro d'instant < 1 > !
  !       - Temps moyen par itération de Newton : 81.619164                                      !
  !       - Temps restant                       : 130.515530                                     !
  !    La base globale est sauvegardée. Elle contient les pas archivés avant l'arrêt.            !
  !                                                                                              !
  !    Conseil :                                                                                 !
  !       - Augmentez le temps CPU.                                                              !
  !-------------------------------------       

20/11/12

20/11/13

メッシュサイズ要素数変位
145819480.0001126387
154431100.0001150798
203354450.000113384
252556900.000107676
271312960.0001036928
30832210.000097704057
40636830.000095167927
50549090.000093780519
60214130.000076680739
70186830.000073994344
80160840.000072255329
100140870.0000700518
15050700.000042708863
20035870.00003988833

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/1syouban.png

20/11/16

メッシュサイズ要素数変位
101107650.000148931
20274750.0001490144
21274190.0001489303
23194020.000148931
25177920.0001489472
30100500.00012982425
4069120.000105412889

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/2zisityuu.png

20/11/17

メッシュサイズ要素数変位
20035870.0001210055
15050700.0001209961
100140870.000121258
80160840.0001212511
70186830.0001212428
60214130.0001212307
50549090.0001212352
40636830.000121233
30832210.0001212595
271312960.000121299
252556900.0001213208

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/2zisyouban.png

20/11/18

メッシュサイズ要素数変位
252556900.0001244172
271312960.0001243948
40636830.0001243271
60214130.0001243248
100140870.0001243528
2035870.0001240943

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/2naosi.png

20/11/19

20/11/20

20/11/25

20/11/26

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/nemoto.png

20/11/27

20/11/30

メッシュサイズ要素数変位
Max100 min20421921.3718466667
M100 m10678801.4224408333
M100 m4.51396611.4488306757
M10 m4.53578441.4893606757
M8 m4.54933231.4849333784

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/10sityuu.png

20/12/1

20/12/2

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/poisson1.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/poisson2.png

20/12/3

20/12/4

20/12/7

全体モデル(弾性)
荷重(t)変位(mm)
0.51.329206
1.02.658416
1.53.98762
2.05.31683
2.56.646036
実験値
荷重(t)変位(mm)
0.52.625
1.05.5
1.58.4
2.011.55
2.514.7

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/model.png

20/12/8

荷重(t)変位(mm)
0.52.5460526316
15.0921063158
1.57.6381589474
210.1842136842
2.512.7302684211

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/ogawa/rubber.png

20/12/9

卒論日誌

日付時間帯作業時間内容立会
10/1116:00-17:001hガイダンス後藤さん
10/1814:30-16:302h創造工房実習後藤さん
10/2113:00-16:003h課題
10/2114:30-16:302h創造工房実習後藤さん
10/2813:30-16:002.5h課題
11/114:30-16:001.5h創造工房実習後藤さん
11/611:30-12:301h課題
11/814:30-16:001.5h創造工房実習後藤さん
11/1114:00-15:001h課題
11/1312:00-13:301h課題
11/1514:30-16:001.5h創造工房実習後藤さん
11/2012:00-12:300.5h課題
11/2214:30-16:001.5h創造工房実習後藤さん
11/2712:00-14:002h課題
12/614:30-16:001.5h創造工房実習後藤さん
12/1111:00-13:302.5h課題
12/1314:30-16:001.5h創造工房実習後藤さん
12/1812:00-14:002h課題
12/1914:30-16:302h課題
12/2014:30-16:001.5h創造工房実習
1/1014:30-16:001.5h創造工房実習後藤さん
1/1514:30-15:301h創造工房実習後藤さん
1/2414:30-16:001.5h創造工房実習後藤さん
1/2918:00-19:301.5h課題
1/2921:00-22:001h課題
1/3114:30-16:001.5h創造工房実習後藤さん
1/3116:00-17:301.5h課題
2/117:00-17:300.5h課題
2/314:00-14:300.5h課題
2/316:30-19:002.5h課題
2/320:00-22:002h課題
2/410:30-11:301h課題
2/413:00-16:003h課題
2/512:00-17:005h課題
2/614:30-15:000.5h課題
2/714:30-15:000.5h創造工房実習後藤さん
2/1313:45-14:451h課題
2/1414:30-15:301h創造工房実習後藤さん
2/189:30-10:000.5h春休み課題
2/2614:00-16:302.5h春休み課題
2/2715:00-16:301.5h春休み課題
3/216:45-17:451h春休み課題
3/516:30-17:301h春休み課題
3/716:00-18:302.5h春休み課題
3/916:00-17:301.5h春休み課題
3/1115:30-17:302h春休み課題
3/1614:00-18:304.5h春休み課題
3/1713:00-16:003h春休み課題
3/1813:30-17:304h春休み課題
3/1914:00-17:303.5h春休み課題
3/2316:00-18:302.5h春休み課題
3/2415:00-17:302.5h春休み課題
3/2515:00-16:301.5h春休み課題
3/2616:00-17:301.5h春休み課題

10/18の課題

メッシュ長さ要素数先端変位相対誤差計算者
0.53941216.57938-1.4%佐藤
1722786.443-3.4%吉田
2118175.7382525-13.9%梅田
428624.9428-25.9%青山
88974.0411725-39.4%小川

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

10/25の課題

メッシュ長さ要素数先端変位相対誤差計算者
0.53941216.57938-1.4%佐藤
0.71309166.4781-2.9%佐藤
0.8721016.437-3.5%青山
0.9717186.4314-3.5%吉田
1722786.443-3.4%吉田
1.2655756.443-3.4%小川
1.4410966.3162-5.3%青山
1.5234176.1206-8.2%梅田
1.8117585.7370-13.98%梅田
2118175.7382525-13.9%梅田
428624.9428-25.9%青山
88974.0411725-39.4%小川

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

考察  メッシュ長さ1の時の要素数が0.9、0.8の時よりも多くなっていた。メッシュ長さが大きくなるほど相対誤差が大きくなった。

11/1の課題

メッシュ長さ要素数変位相対誤差計算者
0.53226870.4287747142.9%小川
0.81999680.41646-0.050%吉田
1623600.416513545-0.037%梅田
1.5503590.411953-1.134%青山
2205260.3706311667-11.05%佐藤

理論値 0.416667

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

考察 メッシュを細かくするほど精度が高くなるわけではない。今回はメッシュ長さが1のとき精度が一番良かった。

11/8の課題

メッシュ長さ要素数変位相対誤差計算者
0.53226870.52149585716.07%小川
0.81999680.51454371434.65%吉田
0.9879530.51027538463.78%梅田
1623600.5057687272732.87%青山
1.1832780.50338136362.38%小川
1.2595150.50532672.78%吉田
1.3574550.50384355562.47%青山
1.5503590.500102251.71%梅田
1.8280530.49266685710.20%佐藤
2205260.4765718333-3.08%佐藤

ティモシェンコ梁 0.4916667 理論値 0.416667

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

11/15の課題

わかったこと ・自動メッシュの3次元を選択までの流れ

わからなかったこと ・スケール係数の欄がどこにあるか ・固定面の節点の選択のしかた ・エンティティの追加の欄がどこにあるか

11/22の課題

拡大係数要素数変位相対誤差計算者
0.3574316.422525-3.7%全員
0.4233496.26413-6.0%佐藤
0.5137416.03154-9.5%小川
0.672575.91187-11.3%梅田
0.742455.5716725-16.4%吉田
0.836575.5734-16.4%青山

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

12/6の課題

拡大係数要素数変位相対誤差計算者
0.81225156.917492.6%吉田
1654906.88772.2%梅田
1.2366936.85192251.6%梅田
1.4233416.81171751.0%吉田
1.6153716.6336-1.6%小川
1.8124806.742580.014%佐藤
285896.66499-0.011%梅田

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

節点数で荷重を割るときに精度を上げていれば変位が下がらずに済んだと思う。 有効数字


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