【Abaqus】使い方講座その5「Hashin損傷」


# はじめに
使い方講座その5です。
他の回は以下のサイトにまとめています。
https://qiita.com/yuki_2020/items/d613eedb149cbbd16ce8#%E4%BD%BF%E3%81%84%E6%96%B9%E8%AC%9B%E5%BA%A7%E9%80%A3%E8%BC%89

今回は以下の動画を紹介したいと思います。
内容としては、ガラス繊維複合材を重ねた穴あき試験片を引っ張る解析です。

https://www.youtube.com/watch?v=uPVQy7JtZEA


# 学べる事
Hashin損傷
連続体シェル
スタッキング





# パート

## 押し出し

最終的に連続体シェルを使用しますが、パートでは
「3次元」「変形体」「ソリッド」「押し出し」を選択します。
![image.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/feb21e57-d49d-c8fd-4e9a-bfbefcef7c63.jpeg)


そして↓の寸法のパートを作成します。
奥行きは2です。

![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/1c510bc3-5fa1-ee03-67c5-f3d1f5176b67.png)





その際、長方形の左下がスケッチの原点に来るように書き、「拘束の追加」→「固定」で原点から動かないようにしましょう。
![image (4).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/70297ef3-0f6e-f0d6-4c4c-9418f3203115.jpeg)

間違ってもこうしないように
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/07a7f4f9-0d0b-84dd-75e1-382175568b53.png)

::: note warn
ここでスケッチ平面の中心を無視して書くと、のちの平面のオフセットでうまく作れなくなります。
ちゃんと長方形の左下が中心に来るようにしましょう。
:::

## パーティションの作成

次に、パーティションを作成します。

::: note point
パーティションの作成することでメッシュをきれいに切ることができます。
:::

### データム平面
4:20~
「データム平面の作成：主平面からのオフセット」を選択します。
今回は「YZ平面」を選択し、「オフセット」90とします。
もう一つ「YZ平面」で「オフセット」110も作りましょう。
するとこのように黄色い点線で面が作成されます。
![image (5).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/b8a4b9d1-f159-7a6a-959c-0cf026193f43.jpeg)



「セルのパーティション：データム平面の使用」を選択
![image (6).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/21516263-5464-0c4a-c612-1feb33a66a0b.jpeg)

先ほど作成したデータム平面をクリックし、「パーティションの作成」を押します。
2つ目を作るときは、「分割するセルを選択してください」と出るので、切りたいデータム平面を含んでいる方を選択して、同じようにパーティションを作ってください。
完成がこちらです。
![image (7).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/b2bdc814-96d6-05f0-49c6-a3f67eaf8843.jpeg)



### 切断面の定義
6:30~
「セルのパーティション：切断面の定義」を選択
![image (8).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/09d8acac-96ef-f886-cd1d-966f7d6b1416.jpeg)

真ん中のセルを選択し、「どのように平面を選択しますか？」とせるので、「3ポイント」を選択
↓の3点を選択し「パーティションの作成」を押す
![名称未設定 1.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/3e7469b7-c67a-767a-e69c-c0ad93eecc64.jpeg)


するとこうなる。
![image.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/2cf4ad8f-145e-5f9b-bf90-120aec122b88.jpeg)




同様にして↓のように分割する。
![image (1).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/e5b1e80a-847a-a2db-03ca-0edd02bd0dad.jpeg)





# 特性
## 材料特性の作成
7:50~


特性では「機械的」→「弾性」→「弾性」を選択し、
![image.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/4bfc38e0-0a93-0126-0b73-c1ef1f1cdd32.jpeg)


そして複合材の場合、タイプを「工学的定数による異方性」として各種剛性やポワソン比を入力します。
以下をコピペしてください。（タブ区切り）
体積分率45%のガラス繊維複合材だそうです。

```
34600	4300	4300	0.38	0.38	0.4	1700	1700	1450
```




なお、なお、今回は寸法をmmで入力しているため、剛性の単位はMPaとして入力しています。

https://www.lancemore.jp/dyna_db/unit_01/index.html


また、今回の座標系は以下のようになっていると予想されます。
![test1_2.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/f1dfaae7-ea51-5503-be7d-ffaccdf39761.png)


次に、「機械的」→「繊維補強複合材での損傷」→「Hashin損傷」を選択

各強度は以下をコピペ
```
700	400	20	60	35	35
```
![image (3).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/1165c09b-76f3-4678-79d6-320b2e6e5a40.jpeg)

次に「サブオプション」→「損傷発展」から、
![image (5).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/d2efcfc7-982c-3714-0a0d-50f5731891f3.jpeg)


以下をコピペ
```
10	6	0.2	0.4
```
![image (4).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/fb1cb7b8-63d7-ce9f-f973-4679c7204a36.jpeg)



次に「サブオプション」→「損傷安定化」から、粘性係数を決める。
![image (6).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/6b96b685-e26b-7463-76bf-8f01a5b693e7.jpeg)

今回は解析を早めるためにもすべて0.001とする。（0.0001が一般的）
![image (7).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/a7a549ee-8a75-d4b9-bf6c-8aefc3b0a616.jpeg)

:::note point
動画によれば、損傷進展が速すぎるので、粘性係数を決めて解析が収束するようにするとのこと。
基本的にステップで決める　Time Periodは1にするため、それよりも十分に小さい10^-4を指定することが一般的で、そうすることで真の解からのずれが小さくなるとのこと
:::

:::note point
Hashin損傷の強度の考え方については、動画の11:00~から詳しく説明しているので、理解していない人はぜひ見てみてください。
縦圧縮強度は縦引張強度の0.5~0.7倍などの情報も。
:::


## 材料特性の割り当て
18:23~

「複合材レイアップの作成」を選択
![image (8).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/357f31a7-d81c-eaca-9521-21c7a3f13142.jpeg)


今回は[0/+45/-45/90]sの積層を模擬するようです。
初期プライ数を4にし、「連続体シェル」を選択します。

![image (9).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/2c071780-83f4-7756-b17e-252e15745b74.jpeg)


その後
「領域」について作成したパートをすべて選択し、
「材料特性」についてもすべて先ほど選択した特性を指定します。
「要素厚」については、今回2 mmの板を8層にするので、0.25mmを指定します。
「回転角度」は上から0、45、-45、90を入力します。
そして「計算荒れた特性を対称にする」にチェックを入れます。

その他については触りません。

「ツール」→「クエリ」→「プライスタック方向プロット」で見ると、ちゃんと積層されていることがわかります。
![image (10).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/5b897779-7244-d678-6780-c261dcc66f27.jpeg)

この見方について詳しく知りたいならこちらを見てください。https://qiita.com/yuki_2020/items/4031daf7b30f1dedd234#%E7%A9%8D%E5%B1%A4%E3%81%AE%E7%A2%BA%E8%AA%8D

# アセンブリ

作ったパートを配置するだけで特に特別なことはしません。


# ステップ


## ステップの作成
21:20~

ステップの作成から、「Static, General」を選択
![image (11).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/861cb1f0-ad55-ccfa-b0eb-b9467d1f8992.jpeg)

「基本」のタブで、自動安定化は「散逸エネルギ比を指定する」を選択。
![image (12).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/9161c8b2-ca12-86af-590b-259d344751bd.jpeg)



「インクリメント」のタブでは、
最大インクリメント数を100→10000に増やし、

時間増分値は
初期値を1→0.001
最小を1E-05→1E-010
に変更
![image (16).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/19e322ba-2565-8fa7-119e-71a570337221.jpeg)



## 一般解法コントロール
22:00~


「ステップ」の上のタブ内の「その他」→「一般解法コントロール」→「編集」→「step-1」を選択します。
![image (15).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/8e4fd492-9a7a-223b-ca02-b4ba1a8d1c5c.jpeg)

警告が出てくるので、そのまま進みます。
![image (14).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/d4905abb-4aa4-06c5-4ef8-e1d45e0921c6.jpeg)


「指定」を選択し、「時間増分」の「時間増分」のタブでI0を20に、IRを40に変更します。
![image (16).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/9b54a4d9-baac-b45b-bffb-89a1e22a2d78.jpeg)

次に、「時間増分」の「次」のタブでIAを20に変更
![image (17).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/ae38f518-10f0-d9d3-7fe9-46ccb0f853f4.jpeg)


## フィールド出力
22:20~


「フィールド出力マネージャー」を選択
![image (20).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/b0f6b4ab-3ea4-fac8-4c91-15282da75a92.jpeg)

step-1の部分の「作成」をダブルクリックする。
![image (21).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/19d5ddfb-8dbb-6568-25fb-cc58acc0b941.jpeg)


「頻度」を「等分時間間隔」にし、「間隔」を100に設定します。
さらに、「疲労/破壊」から以下にチェックを入れます。

* DAMAGEC, 圧縮損傷
* DAMAGEC, 引張り損傷率
* DAMAGEC, 繊維引張り損傷
* DAMAGEC, 繊維圧縮損傷
* DAMAGEC, 母材引張り損傷
* DAMAGEC, 母材圧縮損傷
* SDEG, 剛性低下率
* HSNFTCRT, Hashin 繊維引張り損傷発生基準
* HSNFTCRT, Hashin 繊維圧縮損傷発生基準
* HSNFTCRT, Hashin 母材引張り損傷発生基準
* HSNFTCRT, Hashin 母材圧縮損傷発生基準
![image (22).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/f1b3939e-a940-b7a8-3ad9-526024843630.jpeg)

さらに、「状態/場/ユーザ/時間」のタブで、「STATUSXFEM, xfem 要素の状態」にチェックを入れる
![image (23).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/42c77979-a19b-e18a-d554-27058528482f.jpeg)



## 出力履歴
23:05~

「出力履歴マネージャ」を選択し、
![image (24).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/82e335a4-51ea-3fd5-6089-a8549f25202b.jpeg)

step-1の作成の部分をダブルクリックをする。
![image (25).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/bfd3e3bd-4e49-36e4-79bb-6fe29425405a.jpeg)

先ほどと同じく、「頻度」を「等分時間間隔」にし、「間隔」を100に設定します。
![image.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/cdec301c-8e02-fb3a-3d6f-b968c008fbed.jpeg)

# 荷重
23:15~

:::note point
stepは先ほど作ったstep-1です。
:::


「境界状態の作成」から「変位/回転」を選択
![image (1).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/0fd775bc-707f-57be-7086-c95f68b1edfe.jpeg)

左端を選択し、
![image (2).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/e7a5b901-390c-6d5e-be45-5914a2aa08d2.jpeg)


すべてにチェックを入れます。（なぜ完全固定しないかはわかりませんが）
![image (3).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/cc37d62b-2ab8-3dca-bec0-513bab97141f.jpeg)


次に反対側を引っ張ります。
同様にして反対側の面を選択し、
![image (4).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/83e5a46a-dfa6-0e71-b4ce-1af5386be386.jpeg)


すべてにチェックを入れて「U1」のみ1にします。
![image (5).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/5b72dc6f-b122-d721-4d15-e7cdad60afb9.jpeg)




# メッシュ


## 要素タイプ
24:25~

「要素タイプの割り当て」を選択し、アッセンブリ全体を選択します。
![image (6).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/1005d5be-87d2-f9c9-33a9-faf42cc3ae30.jpeg)


「ファミリ」から「連続体シェル」を選択
動画では最初間違えて「粘着間隙圧」を選択していますが解析できず、後々29:50~で修正しています。
「2次精度」を「はい」
「アワーグラス制御」を「強化型」
「要素の削除」を「はい」
「最大低下率」を「0.999」
とします。
![image (7).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/7f87349d-06d5-8bce-e037-9c8a33a9495c.jpeg)

:::note point
16:50~あたりから、なぜ0.999を使用するのか損傷について説明をされているので興味のある人は見てください。
:::

## スタック方向の確認
27：50~

さらに、「スタック方向の割り当て」を選択すると、
![image (8).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/60d03524-e07b-cce9-515f-19d4a2a85d65.jpeg)

このようにいろいろな面がブラウンとパープルになっていることがわかります。
ブラウンの部分が上面、パープルの部分が下面を表しているため、間違った割り当てになっています。
![image (9).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/f4e8852d-214b-96b5-4dc1-ccb65e4f99b5.jpeg)

:::note warn
このまま間違ったスタック方向で解析するとおかしな結果が出ます。
動画では28:30~からこの状態での解析を実際にやっており、結果だけ見たい方は31:00~見てください。
:::


## スタック方向の修正
32:40~

「クエリ」→「メッシュスタック方向」を押すと、↓のように方向が合わない要素間の境界を赤くハイライトして見ることができます。

![image (10).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/c2283c6d-3f7f-1808-ce46-f7541ad76d41.jpeg)


「スタック方向の割り当て」を選択し、すべてのセルを選択し、スタックの上面となる部分を適当に1箇所選択します。
![image (11).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/367a2a8f-4fdc-21c6-b6c9-51341afefcaf.jpeg)

「参照方向を割り当てますか」と出るので、「はい」を押すときれいに割り当てられます。
![image (12).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/16fe14b9-e5d0-074d-d5ab-3d53f6140853.jpeg)



# ジョブ

特に変更はしません。


# 結果


最初にやった結果はこうなりました。
![image (14).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/1b92dd1b-b1db-0b29-c956-eb10af4b443b.jpeg)


というのも、結果で1コマ進めるだけでこうなっていました。
![image (15).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/00965fa4-e1ec-f097-9b7e-be6e87696454.jpeg)

原因として、「ステップ」の初期値を動画では0.001のところ、0.01としてしまっていたためと考えて直し、さらにメッシュを真ん中で細かくしたところ、1コマ進めて少しは良くなりましたが、だめでした。
![image (17).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/80e25188-8e61-a9e0-2f51-43ead2e33a3f.jpeg)

おそらく、メッシュが粗かったため、初期で変な変形が起きてしまったと考えられます。


初期の変な変形をなくすために、「メッシュ」の「要素タイプの割り当て」の「最大低下率」を0.999のところ、0.5にして、損傷が発生しても剛性は半分までしか落ちないようにしたところ、うまくいきました。
![image (18).jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/993303/bcbfa181-3491-aca4-0743-c90abda63f0f.jpeg)



# 最後に

今回はHashin基準の使い方を説明しました。
同時に連続体シェルの使い方もわかり、かなり勉強になりました。
これで複合材の解析はあらかたできるようになったのではないでしょうか。

また新しい内容があったら連載したいと思います。