Salome-Meca 2018を用いた部品を連結したAssyの応力解析（SSLV）

複数のSolidを連結（結合）してAssyを作り、Assyの線形解析をする問題を考えます。
複数の部品を配置した後、連結する箇所を定義して連結し解析を実行します。
連結部で荷重や変位を相手側に伝え構造解析を行います。

この資料は[オープンCAE勉強会＠岐阜](http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?OpenCAE%20Users%20Wiki)で公開されているFS氏ご提供の「[SALOME-Mecaの使用法解説](http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?SALOME-Meca%A4%CE%BB%C8%CD%D1%CB%A1%B2%F2%C0%E2)：4.0 部品の連結」をSalome-Meca 2018（Code_Aster 13.6）のAsterStudyモジュールを用いて実行しました。

# ジオメトリ

直方体が２個

# 材料

２種類の材料を定義します。柔らかい材料（Aluminum）と硬い材料（Steel）を定義します。

|材料名|ヤング率|ポアソン比|
|:--|:--|:--|
|Aluminum|70,600MPa|0.345|
|Steel|212,000MPa|0.293|

# 拘束条件と荷重条件


|境界条件|トポロジ|境界条件の種類|条件|
|:--|:--|:--|:--|
|fix|面|変位|DX=DY=DZ=0|
|load|面|圧力荷重|0.1MPa|

![0400BC.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/201403/85e5bbe4-635e-5722-ae8f-b6f9149ec181.png)

# 解析結果

相当応力図を示します。

![0400Assy.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/201403/5d99f12f-1191-4750-7023-fd905447b329.png)

# コマンドファイル

```python
DEBUT(LANG='EN')


mesh = LIRE_MAILLAGE(identifier=u'0:1',
                     FORMAT='MED',
                     UNITE=20)

mesh = MODI_MAILLAGE(identifier=u'1:1',
                     reuse=mesh,
                     MAILLAGE=mesh,
                     ORIE_PEAU_3D=_F(GROUP_MA=('press', )))

model = AFFE_MODELE(identifier=u'2:1',
                    AFFE=_F(MODELISATION=('3D', ),
                            PHENOMENE='MECANIQUE',
                            TOUT='OUI'),
                    MAILLAGE=mesh)

Aluminum = DEFI_MATERIAU(identifier=u'3:1',
                         ELAS=_F(E=70600.0,
                                 NU=0.345))

Steel = DEFI_MATERIAU(identifier=u'4:1',
                      ELAS=_F(E=212000.0,
                              NU=0.293))

fieldmat = AFFE_MATERIAU(identifier=u'5:1',
                         AFFE=(_F(GROUP_MA=('Base', ),
                                  MATER=(Aluminum, )),
                               _F(GROUP_MA=('Bar', ),
                                  MATER=(Steel, ))),
                         MAILLAGE=mesh)

load = AFFE_CHAR_MECA(identifier=u'6:1',
                      DDL_IMPO=_F(DX=0.0,
                                  DY=0.0,
                                  DZ=0.0,
                                  GROUP_MA=('fix', )),
                      LIAISON_MAIL=_F(GROUP_MA_ESCL=('contBase', ),
                                      GROUP_MA_MAIT=('Bar', ),
                                      TYPE_RACCORD='MASSIF'),
                      MODELE=model,
                      PRES_REP=_F(GROUP_MA=('press', ),
                                  PRES=0.1))

reslin = MECA_STATIQUE(identifier=u'7:1',
                       CHAM_MATER=fieldmat,
                       EXCIT=_F(CHARGE=load),
                       MODELE=model)

reslin = CALC_CHAMP(identifier=u'8:1',
                    reuse=reslin,
                    CONTRAINTE=('SIGM_ELNO', 'SIGM_NOEU'),
                    CRITERES=('SIEQ_ELNO', 'SIEQ_NOEU'),
                    RESULTAT=reslin)

IMPR_RESU(identifier=u'9:1',
          FORMAT='MED',
          RESU=_F(NOM_CHAM=('SIGM_NOEU', 'SIEQ_NOEU', 'DEPL'),
                  RESULTAT=reslin),
          UNITE=80)

FIN()
```