High-order isogeometric phase-field fracture modeling: a numerical-experimental study

エンジニアが橋梁や航空機、あるいは日常的な道具を設計する際、材料がどのように破壊されるかを予測しなければなりません。最も劇的な破壊の形態は「破断（フラクチャー）」であり、これは固体が突然分裂してしまう現象です。数十年にわたり、科学者たちはこの研究のために主に2つの方法に頼ってきました。一つは実験室での物理的な試験であり、もう一つはコンピュータ・シミュレーションです。物理的な試験は、材料がどのように振る舞うかを示す否定できない証拠となりますが、コストと時間がかかり、かつ破壊を伴います。コンピュータ・シミュレーションは、より迅速で安価な代替手段を提供し、金属片を一つも切断することなく、数千通りのバリエーションをテストすることを可能にします。しかし、長年の課題は両者の間のギャップでした。多くのコンピュータ・モデルは、ガラスのようにパキッと割れる材料を記述するように作られていますが、現実の世界は、壊れる前に曲がったり伸びたりする金属に満ちています。この不一致により、エンジニアは、数学的には単純だが物理的に不正確なモデルか、あるいは複雑だが実行が困難なモデルかの選択を迫られることがよくあります。

イタリアとドイツの大学の研究チームは、最近、驚くべき中間領域を探索しました。彼らは、脆い（もろい）ガラスのような材料のために特別に設計されたコンピュータ・モデルが、延性のある（伸びやすい）金属の挙動をうまく予測できるかどうかを調査しました。対象としたのは、建設や製造において一般的に使用されており、破壊する前に大幅に変形することで知られている特定のアルミニウム合金です。研究者たちは、「脆い」数学的枠組みを用いて、「延性のある」現実世界の事象、具体的には、あらかじめ切り込み（ノッチ）を入れた金属サンプル内での亀裂の発生と進展をシミュレートできるかどうかを検証することにしました。彼らの目的は、金属の曲がりに関するあらゆる微細な詳細を再現することではなく、モデルが2つの重要な事項、すなわち、材料が破壊する前に耐えられる最大荷重と、亀裂が材料を突き抜けていく正確な経路を正確に予測できるかどうかを確認することでした。

実験は、材料強度をテストするための標準的な形状である、ドッグボーン型の金属標本を用いた物理試験から始まりました。研究者たちは、シリコンとマグネシウムを含むEN AW-6060として知られるアルミニウム合金を使用しました。彼らは、厚さ4ミリメートルのものと6ミリメートルのものの2種類の標本を用意しました。各標本には、中央からわずかにオフセンターの位置に、2つの小さな三角形の切り込み（ノッチ）が施されています。これらのノッチは、亀裂が発生しやすい弱点として機能します。チームはこれらの標本を強力な試験機に設置し、破断するまで引き離しました。彼らは、金属の厚さが結果を劇的に変えることを観察しました。より厚い6ミリメートルの標本では、一方のノッチから他方のノッチへと対角線上に進む亀裂が発生し、それらを一つの鋭い線で結びました。対照的に、より薄い4ミリメートルの標本では、2つの別々の亀裂が発生して互いに離れる方向に湾曲し、決して合流することはありませんでした。この挙動の違いは、材料の厚さという物理的な制約がどのように破壊モードを変えるかを示す典型的な例であり、エンジニアリングではよく知られた現象ですが、単純なコンピュータ・モデルで捉えることは困難なものです。

これらの事象をシミュレートするために、研究者たちは「アイソジオメトリック解析」と呼ばれる高度なコンピュータ手法を用いました。これは、従来の格子ベースの手法よりも物体の形状をより正確に表現するために、高度な数学的曲線を使用するものです。彼らは、「フェーズフィールド破断モデル」として知られる特定のタイプのモデルを採用しました。このアプローチでは、亀裂はコンピュータが手動で追跡しなければならない鋭くギザギザした線としては表現されません。代わりに、亀裂は材料が徐々に強度を失っていく、滑らかでぼやけた領域として扱われます。この「フェーズフィールド」変数は、ディマー・スイッチ（調光器）のように機能し、材料が完全に保持されている状態から完全に破壊された状態へと、材料を保持する能力を下げていきます。研究者たちは、もともと脆い材料のために設計されたバージョンのモデルを選択しました。これは、その材料が大きく曲がったり伸びたりすることなく破壊することを前提としているため、大胆な選択でした。なぜなら、彼らがテストしているアルミニウムは延性があり、最終的に破断する前に塑性変形を通じて伸び、エネルギーを吸収するからです。

チームは、物理試験と同じ寸法と材料特性を使用してシミュレーションを実行しました。彼らは、モデル内の「内部長さ」として知られる主要なパラメータを調整しました。これは、亀裂周囲のぼやけた損傷ゾーンの幅を制御するものです。この値を注意深く調整することで、彼らはラボで見られた異なる挙動を再現できることを見出しました。内部長さを大きな値に設定すると、シミュレーションは厚い6ミリメートルの標本と一致するように、2つのノッチを対角線上で結ぶ亀裂パターンを生成しました。内部長さを小さくすると、シミュレーションは、2つの別々の湾曲した亀裂を示し、それらは結合しませんでした。これは、薄い4ミリメートルの標本の挙動を完璧に反映していました。これは、単一の比較的単純な数学的枠組みが、金属の伸びる複雑な物理現象を明示的にモデル化することなく、一つのパラメータを調整するだけで、複雑な物理現象を捉えられることを示唆しており、重要な発見でした。

結果は、視覚的に似ているだけでなく、数値的にも精密でした。コンピュータ・モデルは、金属が破断する前に耐えられる最大荷重を、驚くほど正確に予測しました。厚い標本の場合、予測された荷重は、実際に測定された荷重の1パーセント以内でした。薄い標本の場合、予測は3.5パーセント以内でした。これらの数値は、現実の金属で起こる複雑な塑性変形を意図的に無視したシミュレーションとしては、非常に近いものです。研究者たちは、モデルは金属がピーク強度に達した後のゆっくりとした緩やかな軟化を再現することはできないものの、失敗の瞬間と亀裂の経路を予測することには極めて優れていると指摘しました。これにより、設計が構築される前に安全限界を知る必要があるエンジニアにとって、このモデルは強力なツールとなります。

本研究は、脆性破断モデルを使用して延性材料をシミュレートすることは、エンジニアリング設計において有効かつ保守的な戦略となり得ることを結論付けています。著者たちが「保守的」と言うのは、モデルが材料が実際に耐えられる荷重よりも、わずかに低い荷重で破壊を予測する傾向があることを意味しており、これは安全性に関わる用途において安全なアプローチです。このアプローチの成功は、特定の数学的定式化に依存しています。これにより、エネルギーがプロセスの中で早期に散逸し、延性材料のエネルギー吸収を（明示的に計算することなく）模倣することが可能になります。研究者たちは、これは現象論的な妥協であると認めています。つまり、モデルは金属の内部で降伏した後に何が起こるかという全容を語るものではありませんが、最も重要な部分、すなわち、いつ、どのように破壊されるかを伝えているのです。

この研究は、物理的世界を理解するためにコンピュータ・シミュレーションをどのように活用できるかについて、新しい視点を開いています。それは、正しい答えを得るために、必ずしも最も複雑なモデルを必要とするわけではないことを示唆しています。時には、適切な数学的ツールと、その限界に対する注意深い理解を持って適用された、より単純なモデルが、正確で信頼できる予測を提供できるのです。テストされたアルミニウム合金に対して、研究者たちは、損傷ゾーンのスケールを調整するだけで、脆性モデルが、部品が耐えるために必要な力と、亀裂が辿る複雑な経路の両方を正常に予測できることを実証しました。この知見は、金属を用いた設計を行うエンジニアに対し、塑性変形の微視的な詳細をモデル化するという膨大な計算コストに足を取られることなく、効率的な高次シミュレーションを使用して安全性と性能を確保するための、実用的な道筋を提供しています。本研究は、厳格なラボ試験と高度な数値的手法を組み合わせることで、現実世界のエンジニアリング上の課題を解決できるという、力の証明となっています。