Accurate simulation of delamination with a resin-rich layer-dependent penalty stiffness based on structural cohesive elements

あなたは、何千もの小さなカラフルなレゴブロックを接着して非常に強力な壁を作る、という想像をしてみてください。現実の世界では、エンジニアはこれと同様の構造物を、炭素繊維とプラスチック樹脂の層を使用して構築し、航空機の翼やロケットの部品を作り上げています。これらの材料は非常に強力ですが、ある秘密の弱点があります。それは、層同士が、まるで本からページを引き裂くように、時として剥がれてしまう可能性があることです。この剥離は「デラミネーション（層間剥離）」と呼ばれます。もしこれが予期せず発生すると、構造全体が崩壊してしまう可能性があります。これを防ぐために、科学者たちはコンピュータプログラムを使用して、層がどのように分離するかをシミュレーションします。しかし、これらのプログラムには厄介な問題があります。計算を正しく行うためには、非常に詳細な（数百万もの小さなデジタルブロックを使用する）設定が必要であり、そのせいでコンピュータのシミュレーション実行に膨大な時間がかかってしまうのです。もしデジタルブロックが大きすぎると、コンピュータは混乱し、層が不可能な方法で押しつぶされる（圧縮される）と予測してしまい、誤った答えを導き出してしまいます。

この論文は、この混乱を解決するために、層が剥がれる際に層がどのように振る舞うべきかをコンピュータに伝える、よりスマートな方法を考案することで取り組んでいます。著者である Xiaopeng Ai、Christos Kassapoglou、そして Boyang Chen は、層の間の「接着剤」に対する新しい数学的な「ルールブック」を提案しています。彼らは、接着剤の硬さを壁全体の厚さに基づいて推測するのではなく、繊維層の間に自然に存在する、あの微細で乱れた樹脂の層に着目しました。梁（ビーム）が曲がる仕組みを利用した巧妙なトリックを用いることで、彼らはその薄い接着層を通じてストレス（応力）がどのように広がるかを正確に解明しました。彼らの新しい手法により、コンピュータは数学的な間違いを犯すことなく、より大きな、粗いデジタルブロックを使用できるようになりました。これにより、エンジニアは複雑な剥離問題をより速く、かつ高い精度でシミュレーションできるようになり、これらの材料で作られた航空機やロケットの安全性と信頼性を確保できるのです。

問題点： 「グニャグニャ」な接着剤のミス

複合材料を、多くのパンの層（繊維）と、各スライスの間にある少量のバター（樹脂）で作られた巨大なサンドイッチだと考えてみてください。サンドイッチが破れ始めるとき、バターの層が最初に伸びて壊れます。過去には、コンピュータシミュレーションにおいて、このバター層を単純で均一なバネとして扱っていました。問題は、この「バネ」が、コンピュータがサンドイッチをデジタル的なピースにどのように切り分けたかに応じて、硬すぎたり柔らかすぎたりしたことです。

コンピュータが大きなデジタルスライスを使用したとき、特定のミスが発生しました。それは、破断の直前にある層が、まるで自動車事故のように巨大な力で押しつぶされていると予測してしまうことでした。実際には、層は穏やかに剥がれていくはずなのにです。この「グニャグニャ」としたエラーのせいで、正しい答えを得るためには、エンジニアは極めて小さなデジタルスライスを使用しなければならず、シミュレーションの実行に数日かかることもありました。それは、滑らかな曲線を描こうとして、非常に小さくギザギザしたステップ（階段状の線）のみを使用しているようなものでした。正しく見せるためには、何百万ものステップが必要だったのです。

新しいアイデア： ストレスの声を聞く

著者たちは、接着剤の「剛性（硬さ）」を計算する従来の方法が単純すぎると気づきました。彼らは、ストレス（内部の力）がサンドイッチの厚さ全体に均一に広がるわけではなく、上層から下層へ移動するにつれて変化することを知っていました。これを解決するために、彼らは単に接着剤を見るだけでなく、梁が曲がる物理学に着目しました。

厚い本を曲げる場面を想像してください。外側のページは引き伸ばされ、内側のページは押しつぶされ、真ん中のページは中立の状態を保ちます。著者たちは、この同じ論理を使用して、樹脂層の厚さを通じてストレスがどのように変化するかをマッピングしました。彼らは、複雑で変化するストレスパターンを、コンピュータが使用できる単一の正確な剛性数値へと変換する、「スマートな翻訳機」として機能する新しい公式を導き出しました。

決定的なのは、この新しい公式が、コンピュータが同じ種類の層をいくつも一つの大きな「スーパーレイヤー（等価単一層モデル）」としてグループ化する場合でも機能するという点です。従来の方法では、コンピュータはすべての層を個別にモデル化する必要があり、これは遅くて非効率的でした。新しい手法では、精度を損なうことなく層をグループ化できるため、シミュレーションをはるかに高速に実行できます。

証明： 新しいルールブックのテスト

彼らの新しいアイデアが実際に機能するかどうかを確認するために、チームはエンジニアが数学の検証に使用する一連の標準的な「ベンチマーク」問題を用いてテストを行いました。これらには以下が含まれます：

1. ダブルカンチレバービーム (DCB): 本を開くように、二つの層を引き離すテスト。
2. エンドノッチフレックス (ENF): 層を互いに滑らせるテスト。
3. ミックスモードベンディング (MMB): 引っ張りと滑りの両方を組み合わせたテスト。
4. 強化型DCB (R-DCB): 実世界の強化構造をシミュレートするために、外側に余分な層が追加された、より複雑なテスト。

すべてのテストにおいて、彼らは新しい「スマートな」剛性と、従来の「単純な」剛性を比較しました。結果は明白でした：

• メッシュ収束性: 粗い（大きな）デジタルメッシュを使用したとき、従来の方法は失敗しました。それは、ピースの大きさによって異なる答えを出し、層が強烈に押しつぶされていると予測しました。しかし、新しい手法は、デジタルピースが大きくても小さくても、同じ正確な答えを出しました。それは、ズームインしてもズームアウトしても同じように見える地図を持っているようなものでした。
• 精度: DCBテストにおいて、新しい手法は荷重－変位曲線（層を剥がすのに必要な力）の誤差を大幅に減少させました。7.5 mmのメッシュを使用した場合、従来の方法では50%を超える誤差がありましたが、新しい手法では非常に低く抑えられました。
• 実世界との比較: 実物の物理実験と比較された複雑なR-DCBテストにおいて、新しい手法は、従来のメソッドよりも亀裂の経路と破壊に必要な力をより正確に予測しました。具体的には、破壊力の第2ピークを従来の方法よりも14%低く予測しましたが、これは実世界のデータとより密接に一致していました。

これが意味すること

この論文は、複合材料の世界におけるあらゆる問題を解決したと主張しているわけではありません。これは、剥離をシミュレートするために使用されるデジタルツールである「構造的コヒーシブ要素」の「ペナルティ剛性」を改善することに特化しています。著者たちは、シミュレーションと実験との比較を通じて、彼らの新しい公式が従来のものよりも優れていることを示しました。

この新しい方法を使用することで、エンジニアは複雑な3D構造におけるデラミネーションを、より粗いメッシュでシミュレートできます。これは、以前かかっていた時間のわずかな一部で、正確な結果を得られることを意味します。それは、スーパーコンピュータを必要とせずに、低解像度のビデオゲームから滑らかな高精細なゲームへとアップグレードするようなものです。著者たちは、このアプローチが将来的に、層内部の亀裂のような他のタイプの損傷を予測するためにも拡張できる可能性があると示唆していますが、現時点では、彼らは「グニャグニャな接着剤」の問題を解決し、複合材料のデジタルテストをより速く、より信頼性の高いものにすることに成功しました。