Salome-Meca2019による片持ちはりの曲げ解析（ビーム要素）

# 解析対象
+ 矩形断面はり

全長 *L* = 300 mm
幅 *b* = 100 mm
高さ *h* = 20 mm

+ 境界条件

 先端面に集中荷重 *P* = 1000 N
 荷重点の反対側を固定（3方向変位と3軸回りの回転の合計6自由度）


#はりの弾性理論解

+ 材料特性
ヤング率 *E* = 200 GPa
ポアソン比 *&nu;* = 0.3

+ モデル化と理論解

下図のようにモデル化し，z軸方向の直線が，曲げモーメントを与えた後にも直線を保つと仮定する（Bernoulli-Euler theory）と，自由端変位 *v* は以下のようになる．

断面二次モーメント *I* = *bh<sup>3</sup>*/12 = 6.6667×10<sup>4</sup> mm<sup>4</sup>
自由端たわみ *v* = *PL*<sup>3</sup>/(3 *EI*) = 0.67500 mm
固定端反モーメント *M*<sub>max</sub> = PL = 3.0×10<sup>5</sup> N・mm
曲げ応力 *&sigma;*<sub>*b*</sub> = (*M*<sub>max</sub>/*I*) ×(*h*/2) = 45.000 MPa

※ 曲げ応力（bending stress）とは，はりの断面での応力の分布が線形と仮定し，曲げモーメントとつり合うようにしたものであり，実際の表面応力とは限らない．

![modeling.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/608825/1e4e9f9f-b3e1-65d3-3615-645a7aa71123.jpeg)

#メッシュ
+ チモシェンコはり要素（POU_D_T）
+ AFFE_CARA_ELEM で断面形状 RECTANGLE（矩形）と幅 HY = 100，高さ HZ = 20 の寸法を与える
+ グループ
はりの全体：Edgeグループ（Beam）  
左端（固定端）：Nodeグループ（fix）  
右端（自由端，荷重点）：Nodeグループ（load）
+ 長さ方向均等15分割

# コマンドファイル
コマンドファイルは以下の通り．

```DEBUT(LANG='EN')

mesh = LIRE_MAILLAGE(FORMAT='MED',
                     UNITE=20)

model = AFFE_MODELE(AFFE=_F(MODELISATION=('POU_D_T', ),
                            PHENOMENE='MECANIQUE',
                            TOUT='OUI'),
                    MAILLAGE=mesh)

elemprop = AFFE_CARA_ELEM(MODELE=model,
                          POUTRE=_F(CARA=('HY', 'HZ'),
                                    GROUP_MA=('Beam', ),
                                    SECTION='RECTANGLE',
                                    VALE=(100.0, 20.0)))

mater = DEFI_MATERIAU(ELAS=_F(E=200000.0,
                              NU=0.3))

fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=(mater, ),
                                 TOUT='OUI'),
                         MODELE=model)

fix = AFFE_CHAR_MECA(DDL_IMPO=_F(DRX=0.0,
                                 DRY=0.0,
                                 DRZ=0.0,
                                 DX=0.0,
                                 DY=0.0,
                                 DZ=0.0,
                                 GROUP_NO=('fix', )),
                     MODELE=model)

load = AFFE_CHAR_MECA(FORCE_NODALE=_F(FZ=-1000.0,
                                      GROUP_NO=('load', )),
                      MODELE=model)

reslin = MECA_STATIQUE(CARA_ELEM=elemprop,
                       CHAM_MATER=fieldmat,
                       EXCIT=(_F(CHARGE=fix),
                              _F(CHARGE=load)),
                       MODELE=model)

reslin = CALC_CHAMP(reuse=reslin,
                    CONTRAINTE=('SIEF_ELNO', 'SIPO_ELNO', 'EFGE_NOEU', 'SIPO_NOEU'),
                    FORCE=('REAC_NODA', ),
                    RESULTAT=reslin)

IMPR_RESU(FORMAT='MED',
          RESU=_F(MAILLAGE=mesh,
                  NOM_CHAM=('DEPL', 'EFGE_NOEU', 'SIPO_NOEU', 'REAC_NODA'),
                  RESULTAT=reslin),
          UNITE=80)

IMPR_RESU(FORMAT='RESULTAT',
          RESU=(_F(GROUP_NO=('load', ),
                   IMPR_COOR='OUI',
                   NOM_CHAM=('DEPL', ),
                   RESULTAT=reslin),
                _F(GROUP_NO=('fix', ),
                   IMPR_COOR='OUI',
                   NOM_CHAM=('EFGE_NOEU', 'SIPO_NOEU', 'REAC_NODA'),
                   RESULTAT=reslin)),
          UNITE=2)

FIN()
```

# 結果

![BeamElement.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/608825/c6c60cb8-2fde-44fe-a1e5-e762c8d920e7.jpeg)

+ 固定端（fix）のテキスト出力（IMPR_RESU FORMAT='RESULTAT'）から  
*z*方向変位DEPL DZ = 0.67734 mm（理論解 0.6750 mm）
曲げ応力 SIPO_NOEU SMFY = 45.0000 MPa（理論解 45.0 MPa）


![Text.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/608825/f3e4a763-8c4a-3633-113b-10d205bfb632.jpeg)

本研究の一部は科学研究費補助金（18K02963，代表：藤岡照高）の助成を受けて実施した．