卒論日誌

日付時間帯作業時間内容立会
4/221:00~4:303.5hvi練習後藤さん、野田さん、河原さん、近藤さん
4月合計3.5h
5/234:00~9:305.5hsalome
5/2713:00~18:005h現場見学
5/281:00~3:002hsalome
5/3012:30~13:301hsalome
5/3016:00~18:002hsalome
5月合計20.5h
6/213:00~14:301.5h英語後藤さん、野田さん、河原さん、近藤さん
6/64:30~6:302hgnuplot
6/615:30~17:001.5hgnuplot
6/816:00~18:302.5hgnuplot
6/913:20~14:201hgnuplot
6/1315:30~17:001.5hgnuplot
6/1616:00~18:002.0hsalome
6/195:00~9:304.5hsalome
6/2016:00~18:002hsalome
6/2313:30~14:301hsalome
6/2315:30~18:303hparavis
6/274:30~7:303hsalome
6/2712:30~13:301hsalome
6/2715:30~16:301hsalome
6/3013:30~14:301h英語
6/3016:00~18:302.5hsalome
6月合計31h
7/35:00~8:003hsalome
7/413:00~14:301.5hsalome
7/712:30-14:302hsalome
7/1112:30-14:302hsalome
7/1413:00-14:001hsalome
7/1416:30-18:001.5hsalome
7/2112:30-14:302hsalome
7/2512:30-14:302hオープンキャンパス準備
7/289:30-14:004.5hsalome
7/2913:00-15:002hオープンキャンパス準備
7/309:00-17:008hオープンキャンパス
7月合計29.5h
8/17:30-8:301hsalome
8/113:30-14:301hsalome
8/413:30-17:003.5hsalome
8/913:30-salome
8/275:00-7:002hsalome-摩擦/LaTeX
8月合計9h
9/143:30-6:002.5hsalome-摩擦/接触
9/151:30-6:004.5hsalome-摩擦/接触
9/161:00-6:005hsalome-摩擦/接触後藤さん、野田さん、河原さん、近藤さん
9/171:30-5:304hsalome-摩擦/接触
9/192:30-7:004.5hsalome-摩擦/接触
9/2111:30-15:003.5hsalome-摩擦/接触
9/229:00-15:306.5hsalome-摩擦/接触
9/2412:30-21:006.5hsalome-摩擦/接触
9/2513:00-24:009hsalome-摩擦/接触
9/2614:00-23:009hsalome/TeX
9/2713:00-21:308.5hsalome/TeX
9/2813:00-4:301.5salome/TeX
9月合計65h
10/413:30-19:306h荷重で接触解析
10/511:00-15:304.5h荷重で接触解析
10/614:00-19:305.5h荷重で接触解析
10/711:00-17:004.5h接触解析
10/1014:00-18:004h接触解析
10/1115:00-20:005h接触解析
10/1312:00-18:006h接触解析
10/1410:30-18:306h接触解析
10/179:00-18:009h実験室の片付け/接触解析
10/1815:00-20:005h接触解析
10/1910:30-18:007.5h接触解析
10/2013:00-22:309.5h接触解析
10/2111:00-13:002h接触解析
10/2415:00-19:304.5h接触解析/調べ物
10/2515:00-17:302.5h接触解析/調べ物
10/2613:00-17:004h接触解析
10/279:00-17:008h現場見学/接触解析
10/2810:00-14:304.5h接触解析
10月合計96.5h
11/114:00-20:006h接触解析
11/215:30-19:003.5h接触解析
11/315:00-20:005h接触解析
11/410:00-13:003h接触解析
11/910:00-17:307.5h接触解析
11/1012:00-19:007h接触解析
11/1110:00-18:006.5h接触解析後藤さん、野田さん、河原さん、近藤さん
11/1514:00-18:004h接触解析/調べ物
11/1715:00-21:006h接触解析
11/1811:00-13:002h接触解析
11/2215:00-22:007h接触解析
11/2417:00-24:007h接触解析
11/2512:00-13:001h接触解析
11/2913:30-20:006.5h接触解析
11/3013:00-20:307.5h接触解析
11月合計73.5h
11月までの合計328.5h
12/115:00-18:303.5h接触解析
12/211:00-13:002h接触解析
12/614:30-22:007.5h接触解析
12/712:30-18:005.5h接触解析
12/817:00-22:305.5h接触解析
12/912:00-13:00/17:00-19:303.5h接触解析
12/1216:00-23:007h接触解析/調べ物
12/1315:00-21:306.5h接触解析
12/1412:30-18:305.5h接触解析/Tex
12/1515:30-23:007.5hTex/調べ物
12/1612:00-13:00/17:00-25:009hTex/調べ物
12/1815:00-21:006hTex
12/1914:00-21:007h中間発表関係
12/207:30-8:301h中間発表関係
12/2618:30-19:301h調べ物
12月合計78h
12月までの合計406.5h
1/1213:30-22:008.5h接触解析
1/1310:00-13:003h接触解析
1/171:00-4:003h接触解析
1/1816:00-24:308.5h接触解析
1/1914:00-18:004h接触解析
1/2010:00-13:00/16:00-19:306.5h接触解析
1/2316:00-24:308.5h概要作成
1/2410:00-17:307.5h概要作成
1/2512:00-18:306.5h概要作成/調べ物後藤さん、野田さん、河原さん、近藤さん
1/2717:00-20:003h接触解析/調べ物
1/281:00-6:005h接触解析/調べ物
1/292:00-6:00/15:30-18:307h接触解析/調べ物
1/3015:00-21:306.5hデータ整理
1/3113:00-20:307.5hデータ整理
1月合計77.5h
1月までの合計484h
2/114:30-18:304hデータ整理
2/214:00-22:008hデータ整理/接触解析
2/88:00-18:0010hデータ整理/スライド作成
2/912:00-18:006hデータ整理/スライド作成
2/1010:00-19:009h発表練習/スライド作成
2/136:30-18:0011.5hスライド作成
2/149:30-17:007.5hスライド作成
2月合計56h
合計540h

後藤メモ(16/12/9)

紙を挟んでいる写真は、2013/2/18のものだと思います。 この時の実験の写真とデータをk2のzikkenの中に入れておきました。

卒論テーマ

  • プレストレスを受ける鋼板と木材の接触解析

接触解析

非一体化孔無しモデル(片持ち梁/スパン6m)のcommファイル(載荷断面の全節点に一様載荷)

  • ( DEBUT();

kou=DEFI_MATERIAU(ELAS=_F(E=206e+9,

                         NU=0.3,),);

moku=DEFI_MATERIAU(ELAS_ORTH=_F(E_L=6.929e9,

                               E_T=0.27716e9,
                               E_N=0.27716e9,
                               NU_LT=0.016,
                               NU_LN=0.016,
                               NU_TN=0.4,
                               G_LT=0.46193e9,
                               G_LN=0.46193e9,
                               G_TN=0.46193e9,),);

MAIL=LIRE_MAILLAGE(FORMAT='MED',);

MAIL=MODI_MAILLAGE(reuse =MAIL,

                  MAILLAGE=MAIL,
                  ORIE_PEAU_3D=_F(GROUP_MA=('m1r','m2l','m6l','m3l','m2r','m3r','m4l','m','m4r','m5r','m7l','m6r','m7r','m8l',
                                             'm8r',),),);

MODE=AFFE_MODELE(MAILLAGE=MAIL,

                AFFE=_F(TOUT='OUI',
                        PHENOMENE='MECANIQUE',
                        MODELISATION='3D',),);

MATE=AFFE_MATERIAU(MAILLAGE=MAIL,

                  AFFE=(_F(GROUP_MA='ita',
                           MATER=kou,),
                        _F(GROUP_MA=('moku1','moku2','moku3','moku4','moku5','moku6','moku7','moku8',),
                           MATER=moku,),),);

contact=DEFI_CONTACT(MODELE=MODE,

                    FORMULATION='DISCRETE',
                    FROTTEMENT='COULOMB',
                    ITER_GEOM_MAXI=10000,
                    ZONE=(
                    _F(GROUP_MA_MAIT='m1r',
                       GROUP_MA_ESCL='m2l',
                       ALGO_CONT='PENALISATION',
                       E_N=6.929e9,
                       COULOMB=0.4,
                       ALGO_FROT='PENALISATION',
                       E_T=6.929e8,),
                    _F(GROUP_MA_MAIT='m3l',
                       GROUP_MA_ESCL='m2r',
                       ALGO_CONT='PENALISATION',
                       E_N=6.929e9,
                       COULOMB=0.4,
                       ALGO_FROT='PENALISATION',
                       E_T=6.929e8,),
                    _F(GROUP_MA_MAIT='m3r',
                       GROUP_MA_ESCL='m4l',
                       ALGO_CONT='PENALISATION',
                       E_N=6.929e9,
                       COULOMB=0.4,
                       ALGO_FROT='PENALISATION',
                       E_T=6.929e8,),
                    _F(GROUP_MA_MAIT='m',
                       GROUP_MA_ESCL='m4r',
                       ALGO_CONT='PENALISATION',
                       E_N=6.929e9,
                       COULOMB=0.4,
                       ALGO_FROT='PENALISATION',
                       E_T=6.929e8,),
                    _F(GROUP_MA_MAIT='m5r',
                       GROUP_MA_ESCL='m6l',
                       ALGO_CONT='PENALISATION',
                       E_N=6.929e9,
                       COULOMB=0.4,
                       ALGO_FROT='PENALISATION',
                       E_T=6.929e8,),
                    _F(GROUP_MA_MAIT='m7l',
                       GROUP_MA_ESCL='m6r',
                       ALGO_CONT='PENALISATION',
                       E_N=6.929e9,
                       COULOMB=0.4,
                       ALGO_FROT='PENALISATION',
                       E_T=6.929e8,),
                    _F(GROUP_MA_MAIT='m7r',
                       GROUP_MA_ESCL='m8l',
                       ALGO_CONT='PENALISATION',
                       E_N=6.929e9,
                       COULOMB=0.4,
                       ALGO_FROT='PENALISATION',
                       E_T=6.929e8,),
                    _F(GROUP_MA_MAIT='m',
                       GROUP_MA_ESCL='m8r',
                       ALGO_CONT='PENALISATION',
                       E_N=6.929e9,
                       COULOMB=0.4,
                       ALGO_FROT='PENALISATION',
                       E_T=6.929e8,),),);

kotei1=AFFE_CHAR_MECA(MODELE=MODE,

                     DDL_IMPO=(_F(GROUP_MA=('kotei1','kotei2','kotei3','kotei4','kotei5','kotei6','kotei7','kotei8','kotei9',),
                                  DX=0.0,
                                  DY=0,
                                  DZ=0,),
                               _F(GROUP_MA=('kotei19','kotei20',),
                                  DY=0,),),);

kajuu=AFFE_CHAR_MECA(MODELE=MODE,

                    FORCE_FACE=_F(GROUP_MA='ps_ue',
                                  FY=-0.972e6,),);

kajuu2=AFFE_CHAR_MECA(MODELE=MODE,

                     FORCE_FACE=_F(GROUP_MA='ps_sita',
                                   FY=-0.972e6,),);

zurasi=AFFE_CHAR_MECA(MODELE=MODE,

                     FORCE_NODALE=_F(GROUP_NO=('ten_m5','ten_m6','ten_m7','ten_m8','ten_k','ten_m1','ten_m2','ten_m3','ten_m4',),
                                     FZ=-141,),);                #25kN/節点数177

fonc_fy=DEFI_FONCTION(

                     NOM_PARA='INST',
                     VALE=(0,0,
                           10,1,
                           20,1,
                           ),);

fonc_fz=DEFI_FONCTION(

                     NOM_PARA='INST',
                     VALE=(0,0,
                           10,0,
                           20,1,
                           ),);

kankaku1=DEFI_LIST_REEL(DEBUT=0,

                       INTERVALLE=_F(JUSQU_A=10,
                                     NOMBRE=20,),);

kankaku2=DEFI_LIST_REEL(DEBUT=10,

                       INTERVALLE=_F(JUSQU_A=20,
                                     NOMBRE=20,),);

kankaku3=DEFI_LIST_REEL(DEBUT=0,

                       INTERVALLE=_F(JUSQU_A=20,
                                     NOMBRE=40,),);

resu=STAT_NON_LINE(MODELE=MODE,

                  CHAM_MATER=MATE,
                  EXCIT=(_F(CHARGE=kotei1,),
                         _F(CHARGE=kajuu,
                            FONC_MULT=fonc_fy,),
                         _F(CHARGE=kajuu2,
                            FONC_MULT=fonc_fy,),),
                  CONTACT=contact,
                  COMP_ELAS=_F(RELATION='ELAS',),
                  INCREMENT=_F(LIST_INST=kankaku1,),
                  NEWTON=_F(REAC_ITER=1,),
                  CONVERGENCE=_F(RESI_GLOB_MAXI=100,
                                 ITER_GLOB_MAXI=1000,),
                  SOLVEUR=_F(METHODE='MUMPS',
                             SYME='OUI',),
                  ARCHIVAGE=_F(LIST_INST=kankaku3,),);     #結果の出力数

resu=STAT_NON_LINE(reuse =resu,

                  MODELE=MODE,
                  CHAM_MATER=MATE,
                  EXCIT=(_F(CHARGE=kotei1,),
                         _F(CHARGE=zurasi,
                            FONC_MULT=fonc_fz,),
                         _F(CHARGE=kajuu,
                            FONC_MULT=fonc_fy,),
                         _F(CHARGE=kajuu2,
                            FONC_MULT=fonc_fy,),),
                  CONTACT=contact,
                  COMP_ELAS=_F(RELATION='ELAS',),
                  ETAT_INIT=_F(EVOL_NOLI=resu,),
                  INCREMENT=_F(LIST_INST=kankaku2,),
                  NEWTON=_F(REAC_ITER=1,),
                  CONVERGENCE=_F(RESI_GLOB_MAXI=200,
                                 ITER_GLOB_MAXI=10000,),
                  SOLVEUR=_F(METHODE='MUMPS',
                             SYME='OUI',),
                  ARCHIVAGE=_F(LIST_INST=kankaku3,),);     #結果の出力数

resu=CALC_ELEM(reuse =resu,

              MODELE=MODE,
              CHAM_MATER=MATE,
              RESULTAT=resu,
              OPTION=('SIEF_ELNO','SIEQ_ELNO',),);

resu=CALC_NO(reuse =resu,

            RESULTAT=resu,
            OPTION=('SIGM_NOEU','SIEQ_NOEU',),);

IMPR_RESU(FORMAT='MED',

         UNITE=80,
         RESU=_F(MAILLAGE=MAIL,
                 RESULTAT=resu,
                 NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','DEPL',),),);

FIN();

  • ) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/danmen2.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/bara_128_ten2.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/hyou128.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/128_ten_tawami2.png
  • スパンは6m
  • meshの分割数 軸方向128分割、角材は1本当たり幅2分割×高さ5分割、鋼板は1枚当たり幅2分割×高さ14分割とした(1次要素)。

非一体化孔無しモデル(単純梁/スパン6m)のcommファイル(木材先端の上部0.1*0.25に面載荷)

  • ( DEBUT();

kou=DEFI_MATERIAU(ELAS=_F(E=206e+9,

                         NU=0.3,),);

moku=DEFI_MATERIAU(ELAS_ORTH=_F(E_L=6.929e9,

                               E_T=0.27716e9,
                               E_N=0.27716e9,
                               NU_LT=0.016,
                               NU_LN=0.016,
                               NU_TN=0.4,
                               G_LT=0.46193e9,
                               G_LN=0.46193e9,
                               G_TN=0.46193e9,),);

MAIL=LIRE_MAILLAGE(FORMAT='MED',);

MAIL=MODI_MAILLAGE(reuse =MAIL,

                  MAILLAGE=MAIL,
                  ORIE_PEAU_3D=_F(GROUP_MA=('m1r','m2l','m6l','m3l','m2r','m3r','m4l','m','m4r','m5r','m7l','m6r','m7r','m8l',
                                             'm8r',),),);

newmesh=CREA_MAILLAGE(MAILLAGE=MAIL,

                     CREA_POI1=_F(NOM_GROUP_MA='spElmt',
                                  GROUP_NO=('addsp1','addsp2','addsp3','addsp4','addsp5','addsp6','addsp7','addsp8','addsp9',),),);

MODE=AFFE_MODELE(MAILLAGE=newmesh,

                AFFE=(_F(TOUT='OUI',
                         PHENOMENE='MECANIQUE',
                         MODELISATION='3D',),
                      _F(GROUP_MA='spElmt',
                         PHENOMENE='MECANIQUE',
                         MODELISATION='DIS_T',),),);

MATE=AFFE_MATERIAU(MAILLAGE=newmesh,

                  AFFE=(_F(GROUP_MA='ita',
                           MATER=kou,),
                        _F(GROUP_MA=('moku1','moku2','moku3','moku4','moku5','moku6','moku7','moku8',),
                           MATER=moku,),),);

contact=DEFI_CONTACT(MODELE=MODE,

                    FORMULATION='DISCRETE',
                    FROTTEMENT='COULOMB',
                    ITER_GEOM_MAXI=10000,
                    ZONE=(
                    _F(GROUP_MA_MAIT='m1r',
                       GROUP_MA_ESCL='m2l',
                       ALGO_CONT='PENALISATION',
                       E_N=6.929e9,
                       COULOMB=0.4,
                       ALGO_FROT='PENALISATION',
                       E_T=6.929e8,),
                    _F(GROUP_MA_MAIT='m3l',
                       GROUP_MA_ESCL='m2r',
                       ALGO_CONT='PENALISATION',
                       E_N=6.929e9,
                       COULOMB=0.4,
                       ALGO_FROT='PENALISATION',
                       E_T=6.929e8,),
                    _F(GROUP_MA_MAIT='m3r',
                       GROUP_MA_ESCL='m4l',
                       ALGO_CONT='PENALISATION',
                       E_N=6.929e9,
                       COULOMB=0.4,
                       ALGO_FROT='PENALISATION',
                       E_T=6.929e8,),
                    _F(GROUP_MA_MAIT='m',
                       GROUP_MA_ESCL='m4r',
                       ALGO_CONT='PENALISATION',
                       E_N=6.929e9,
                       COULOMB=0.4,
                       ALGO_FROT='PENALISATION',
                       E_T=6.929e8,),
                    _F(GROUP_MA_MAIT='m5r',
                       GROUP_MA_ESCL='m6l',
                       ALGO_CONT='PENALISATION',
                       E_N=6.929e9,
                       COULOMB=0.4,
                       ALGO_FROT='PENALISATION',
                       E_T=6.929e8,),
                    _F(GROUP_MA_MAIT='m7l',
                       GROUP_MA_ESCL='m6r',
                       ALGO_CONT='PENALISATION',
                       E_N=6.929e9,
                       COULOMB=0.4,
                       ALGO_FROT='PENALISATION',
                       E_T=6.929e8,),
                    _F(GROUP_MA_MAIT='m7r',
                       GROUP_MA_ESCL='m8l',
                       ALGO_CONT='PENALISATION',
                       E_N=6.929e9,
                       COULOMB=0.4,
                       ALGO_FROT='PENALISATION',
                       E_T=6.929e8,),
                    _F(GROUP_MA_MAIT='m',
                       GROUP_MA_ESCL='m8r',
                       ALGO_CONT='PENALISATION',
                       E_N=6.929e9,
                       COULOMB=0.4,
                       ALGO_FROT='PENALISATION',
                       E_T=6.929e8,),),);

softSp=AFFE_CARA_ELEM(MODELE=MODE,

                     DISCRET=_F(CARA='K_T_D_N',
                                GROUP_MA='spElmt',
                                VALE=(1e5,1e5,1e5,),),);

kotei1=AFFE_CHAR_MECA(MODELE=MODE,

                     DDL_IMPO=(_F(GROUP_MA=('kotei1','kotei2','kotei3','kotei4','kotei5','kotei6','kotei7','kotei8','kotei9',),
                                  DX=0.0,),
                               _F(GROUP_MA=('y_u','y_s',),
                                  DY=0,),
                               _F(GROUP_MA=('sen1','sen2','sen3','sen4','sen5',),
                                  DZ=0,),),);

kajuu=AFFE_CHAR_MECA(MODELE=MODE,

                    FORCE_FACE=_F(GROUP_MA='ps',
                                  FY=-0.9259259259e6,),);

zurasi=AFFE_CHAR_MECA(MODELE=MODE,

                     FORCE_FACE=_F(GROUP_MA=('saika1','saika2','saika3',),
                                   FZ=-1000000,),);

fonc_fy=DEFI_FONCTION(

                     NOM_PARA='INST',
                     VALE=(0,0,
                           10,1,
                           20,1,
                           ),);

fonc_fz=DEFI_FONCTION(

                     NOM_PARA='INST',
                     VALE=(0,0,
                           10,0,
                           20,1,
                           ),);

kankaku1=DEFI_LIST_REEL(DEBUT=0,

                       INTERVALLE=_F(JUSQU_A=10,
                                     NOMBRE=1,),);

kankaku2=DEFI_LIST_REEL(DEBUT=10,

                       INTERVALLE=_F(JUSQU_A=20,
                                     NOMBRE=10,),);

kankaku3=DEFI_LIST_REEL(DEBUT=0, #ここの3行は必要ない

                       INTERVALLE=_F(JUSQU_A=20,
                                     NOMBRE=200,),);

resu=STAT_NON_LINE(MODELE=MODE,

                  CHAM_MATER=MATE,
                  CARA_ELEM=softSp,
                  EXCIT=(_F(CHARGE=kotei1,),
                         _F(CHARGE=kajuu,
                            FONC_MULT=fonc_fy,),),
                  CONTACT=contact,
                  COMP_ELAS=_F(RELATION='ELAS',
                               DEFORMATION='PETIT',),
                  INCREMENT=_F(LIST_INST=kankaku1,),
                  NEWTON=_F(REAC_ITER=1,),
                  CONVERGENCE=_F(RESI_GLOB_MAXI=1000,
                                 ITER_GLOB_MAXI=1000,),
                  SOLVEUR=_F(METHODE='MUMPS',
                             SYME='OUI',),);

resu=STAT_NON_LINE(reuse =resu,

                  MODELE=MODE,
                  CHAM_MATER=MATE,
                  CARA_ELEM=softSp,
                  EXCIT=(_F(CHARGE=kotei1,),
                         _F(CHARGE=zurasi,
                            FONC_MULT=fonc_fz,),
                         _F(CHARGE=kajuu,
                            FONC_MULT=fonc_fy,),),
                  CONTACT=contact,
                  COMP_ELAS=_F(RELATION='ELAS',
                               DEFORMATION='PETIT',),
                  ETAT_INIT=_F(EVOL_NOLI=resu,),
                  INCREMENT=_F(LIST_INST=kankaku2,),
                  NEWTON=_F(REAC_ITER=1,),
                  CONVERGENCE=_F(RESI_GLOB_MAXI=1000,
                                 ITER_GLOB_MAXI=10000,),
                  SOLVEUR=_F(METHODE='MUMPS',
                             SYME='OUI',),);

resu=CALC_ELEM(reuse =resu,

              MODELE=MODE,
              CHAM_MATER=MATE,
              RESULTAT=resu,
              OPTION=('SIEF_ELNO','SIEQ_ELNO',),);

resu=CALC_NO(reuse =resu,

            RESULTAT=resu,
            OPTION=('SIGM_NOEU','SIEQ_NOEU',),);

IMPR_RESU(FORMAT='MED',

         UNITE=80,
         RESU=_F(MAILLAGE=MAIL,
                 RESULTAT=resu,
                 NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','DEPL',),),);

FIN();

  • ) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/hankai_mk2_dan2.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/hankai_mk2_2.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/hankai_tawa3.png
  • スパンは6m
  • meshの分割数 軸方向63分割、角材は1本当たり幅12分割×高さ5分割、鋼板は1枚当たり幅2分割×高さ14分割とした(1次要素)。
  • 中央付近のたわみは12mm程度になった(グラフはまだ作ってない)。滝田さんのプログラムによると6.5mmくらいで、上図の値も平均するとそのくらいになりそう。

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/hankai_126_mk2_dan2.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/hankai_126_mk2_2.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/hankai_126_mk2_tawa.png

  • スパンは6m
  • meshの分割数 軸方向126分割、角材は1本当たり幅12分割×高さ8分割、鋼板は1枚当たり幅4分割×高さ24分割とした(1次要素)。

非一体化孔有りモデル(単純梁/スパン7m)

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/anaaki73_dan_2.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/anaaki73_yoko_2.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/anaaki73_2.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/anaaki73_tawa.png

  • スパンは7m
    • meshの分割数 軸方向73分割、角材は1本当たり幅12分割×高さ8分割、鋼板は角材と接触している部分は幅4分割×高さ8分割,三角孔の部分は上図の通り四面体要素とした(1次要素)。グラフは下。

非一体化孔無しモデル(単純梁/スパン7m)

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/ananasi73_dan2.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/ananasi73_2.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/73_gura2.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/ananasi73_tawa.png

  • meshの分割数 軸方向73分割、角材は1本当たり幅12分割×高さ8分割、鋼板は1枚当たり幅4分割×高さ24分割(1次要素)。
  • 実験値 実験値は2013/2/18-19のもの(データは"tugite130218"フォルダの中)。FEMの方は補剛材とか対傾構は入れてない。
    • 実験1→J-TEST-1 1回目、補剛材なし(はらみが出たので)
    • 実験2→J-TEST-2 2回目、補剛材なし(クラップゲージではらみ測定)
    • 実験3→J-TEST-hs 補剛材あり
  • 摩擦係数が0.3, 0.4, 0.5, の時の孔ありモデルと孔なしモデルのたわみの比較。 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/0305ana.png

幅員方向に1/2したモデル(中央に載荷)

  • 対傾構なし http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/m.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/m_t.png
  • 対傾構あり http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/mt.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/mt_t.png
  • 中央に載荷した時の各点でのたわみの比較
    • 摩擦係数0.4(左), 摩擦係数0.3(右) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gura_m.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/03_m.png

幅員方向に1/2したモデル(中央から1/8の位置に載荷)

  • 対傾構なし http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/t1.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/t1_t.png
  • 対傾構あり http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/18t.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/18t_t.png
  • 中央から1/8の位置に載荷した時の各点でのたわみの比較
    • 摩擦係数0.4(左), 摩擦係数0.3(右) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gura_18.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/03_18.png

幅員方向に1/2したモデル(中央から1/4の位置に載荷)

  • 対傾構なし http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/t2.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/t2_t.png
  • 対傾構あり http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/14t.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/14t_t.png
  • 中央から1/4の位置に載荷した時の各点でのたわみの比較
    • 摩擦係数0.4(左), 摩擦係数0.3(右) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gura_14_2.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/03_14.png

幅員方向に1/2したモデル(中央から3/8の位置に載荷)

  • 対傾構なし http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/t3.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/t3_t.png
  • 対傾構あり http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/38t.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/38t_t.png
  • 中央から3/8の位置に載荷した時の各点でのたわみの比較
    • 摩擦係数0.4(左), 摩擦係数0.3(右) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gura_38.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/03_38.png

幅員方向に1/2したモデル(木材上部に当分布荷重)

  • 対傾構なし http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/tou.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/tou_tawa.png
  • 対傾構あり http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/tou_t.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/tou_t_tawa.png
  • 当分布荷重を載荷した時の各点でのたわみの比較
    • 摩擦係数0.4(左), 摩擦係数0.3(右) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/toubunpu_gura.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/03_tou.png

1/27-2/3

摩擦係数を変えて計算。データ整理。

1/25-1/27

摩擦係数とかを変えてやってみる。

1/20-1/25

当分布荷重で解く。東北支部に出す概要を作成。

1/13-1/20

対傾構をつける

  • 1/18 接触面の面積を減らしたら計算できた。まだ減らせる場所はあるのでもっと計算を速く出来るようになるかもしれない。 グラフ作成中。
  • 1/17 対傾構を加えたモデルが一応完成した(実物とはちょっと違うけど)。メッシュもある程度綺麗に切ることが出来た。ただ、メモリの関係で計算が出来なかった。要素数は今までとあまり変わっていないので接触面が多いせいかもしれない。現在のモデルは接触面の面積がかなり大きめになっているので次回ここを減らしてやってみる。出来なければどうにかして対傾構部分のみ一体化させる。

12/9-12/16

前やった片持ち梁モデルの結果とティモシェンコ梁の理論値を比較する。単純梁モデルを線荷重でやる。一体化モデルも作る。

  • 斎藤さんとかのデータを確認すると、たわみの値は平均しているようだった。自分の値も平均するとその値に近いのでこれで良いような気もする。

12/2-12/9

meshの分割数を変える。滝田さんのプログラムで理論値を求めFEMと比較する。

  • 12/6 理論値に近い値が出たので、meshをさらに細かくして計算中。

11/25-12/2

実物に近いモデルを作る。

  • 12/1 前日からやっていた計算が終わる気配が無いので一度中断した。meshの分割数を減らし再度計算中。
  • 11/30 計算中。拘束するための線とかがあってmeshが切りづらいので分割の仕方は今までとは若干違う。それと、片持ち梁状のモデルでは拘束面があり物体が安定していたが、今回のモデルでは拘束面がなく不安定なので弱いバネを追加した。

11/18-11/25

一体化モデルで出来たので接触でも出来るか確認する。

  • 11/22 接触のモデルでも出来た。ただ、計算中に収束しにくい部分の分割数を自動で細かくする方法を使いたかったが、上手く機能せず収束せずに終わった。そのため、とりあえず収束条件を緩めに設定して収束させた。

11/11-11/18

完全に一体化したモデルで解析して、エラーの原因が接触させるためのコマンドにあるのかそれ以前の問題なのかを確認する。

11/4-11/11

手計算と比較する。

10/28-11/4

剛性比で荷重を分配すると鋼板と木材のズレは小さくなったが、剛性の値が上村さんの値と合わなかった。また、未だにたわみの値も合わないのでcommファイルの中身をもう一度確認する。上村さんは節点に荷重を載荷しているが、自分は面に載荷しているのでその違いもあるかもしれない。

  • 上村さん\( \frac{E_{s}I_{s}}{{E_{s}I_{s}}+{E_{w}I_{w}}} \)=0.61 \( \frac{E_{w}I_{w}}{{E_{s}I_{s}}+{E_{w}I_{w}}} \)=0.39
  • 自分\( \frac{E_{s}I_{s}}{{E_{s}I_{s}}+{E_{w}I_{w}}} \)=0.426 \( \frac{E_{w}I_{w}}{{E_{s}I_{s}}+{E_{w}I_{w}}} \)=0.574

10/21-10/28

剛性比で荷重を分配して載荷してみる。

10/14-10/21

データをグラフにして比較する。 木材部もバラバラにして接触させて解く。

  • 10/18 木材部も接触させたモデル(軸方向512分割)を計算中。

10/7-10/14

メッシュを細かくして、過去のデータと比較する。

10/3-10/7

変位ではなく荷重を与えた状態で接触解析を行う。

  • 10/4 計算は最後まで行ったが接触できて無いような気がする。
  • 10/6 鋼板と木材の接触面を2箇所に増やし荷重を与えた状態で接触解析ができた。

夏休み

8/11-8/16は帰省します。

  • 8/8 摩擦ありの接触解析(固定した大きめの物体の上に小さめの物体を置き、上から一定の力を加えたまま上の物体を横方向に引っ張る。横方向の力を徐々に大きくし、動き出すときの摩擦力を調べる。)のモデル作成中。
  • 8/16 帰省中のため関連文献などを読んでいた。salomeの方は進展なし。
  • 8/22 引き続きsalome,eficasを編集中。
  • 9/12 commファイル修正中、エラーをなかなか改善できない。
  • 9/15 摩擦有りで接触解析ができた。引き続き別のモデルで接触解析を行う。 LaTeXもやる。

課題

5/16

  • tetrahedoron 線形 線載荷
    • たわみ(初):\( v_{初} = \frac{pl^3}{3EI} = 0.000833333 \mathrm{m} \)
    • たわみ(T):\( v_{T} = \frac{pl^3}{3EI}+ \frac{pl}{KGA} = 0.000859333 \mathrm{m} \)
      length[mm]FEM[m]相対誤差(初)[%]相対誤差(T)[%]
      160.00034306458.860.1
      80.00060983726.829.0
      40.00070794815.017.6
      20.0007873485.528.38
      10.0008364490.3742.66

5/23

  • tetrahedron 線形 面載荷
    length[mm]FEM[m]相対誤差(初)[%]相対誤差(T)[%]
    160.0003425858.960.1
    80.0006092726.929.1
    40.0007070415.217.7
    20.000785945.698.54
    10.000833770.05282.97
  • tetrahedron 二次 面載荷
    length[mm]FEM[m]相対誤差(初)[%]相対誤差(T)[%]
    160.000835800.2962.74
    80.000852232.270.826
    40.000855182.620.483
    20.000856562.790.322
    10.000857162.860.253
  • hexahedoron 線形 面載荷
    length[mm]FEM[m]相対誤差(初)[%]相対誤差(T)[%]
    20.0008342130.1062.92
    10.0008507042.081.00
    0.50.0008554412.650.453
  • hexahedoron 二次 面載荷
    length[mm]FEM[m]相対誤差(初)[%]相対誤差(T)[%]
    20.0008556422.680.43
    10.0008567542.810.30
    0.50.0008554412.650.453

5/30

  • 固定端の応力分布 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gazou.png

6/6

  • 固定端の応力分布 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gazou2.png
  • 固定端から一要素となりの断面の応力分布 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gazou3.png
  • 梁の真ん中の断面の応力分布 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gazou4.png
  • 3箇所の比較(data2:固定端,data2mm:固定端の一要素となり,data25mm:梁の真ん中) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gazou5.png
  • 手計算との比較
    • 固定端 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gazou6.png
    • 固定端の一要素となり http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gazou7.png
    • 梁の真ん中 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gazou8.png

6/13

  • 固定端 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gazou12.png
  • 固定端の一要素となり http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gazou10.png
  • 梁の真ん中 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gazou11.png

6/20

  • 固定端 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gazou14.png
  • 固定端の一要素となり http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gazou15.png
  • 梁の真ん中 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2016/ebi/gazou17.png

gnuplot関係

  • 凡例非表示
    • ( unset key
    • )
  • y軸反転
    • ( set yrange [:]
    • )
           [20:0]とかにすればy軸が20-0まで表示。x軸も同じ。
  • xfigに出力
    • ( set term fig color set output 'ファイル名.fig' plot 'データ名'
    • )    文字の大きさを変えるなら
    • ( set term fig color fontsize 20
    • )
           とか。
           gnuplotを終了して
    • ( xfigp&
    • )
           でxfigを立ち上げる。
         File→Openで開く。

トップ   編集 凍結 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2020-01-20 (月) 11:16:08