80 HEX / 189 nodes
Coreform Cubit 2026.8 → Abaqus形式mesh → CalculiX 2.23 → Python / VTK。最大変位 0.167306 mm。AIのE2E実行とエラー回復を確認。
短い自然言語の指示から、Cubitで形状・HEXメッシュを作成し、CalculiXで構造解析、Python + VTKで結果を可視化した実機ワークフローです。今回はその既存デモを再取得し、反力つり合いとメッシュ収束を追加して公開検証ケースにしました。
公開検証基準を満たしました。
5120 C3D8要素で、片持ち梁のEuler–Bernoulli理論解に対する先端変位誤差は 0.883%。固定端反力は 100.0 N で負荷と一致しました。
2026年8月の実機デモでは、AIがCubit、CalculiX、後処理まで自律的に進め、途中の変換エラーにも別手段へ切り替えて完走しました。ただし、最初の80要素モデルはワークフロー確認用で、理論解との数値検証は十分ではありませんでした。そこで同一の形状・材料・荷重条件を固定し、C3D8要素を細分化して収束を確認しました。
Coreform Cubit 2026.8 → Abaqus形式mesh → CalculiX 2.23 → Python / VTK。最大変位 0.167306 mm。AIのE2E実行とエラー回復を確認。
同じ100 × 10 × 10 mm梁、E=210000 MPa、ν=0.30、先端総荷重100 Nでメッシュを細分化し、CalculiX 2.21でも独立再実行。
先端平均|Uz|とEuler–Bernoulli理論解の差を1%未満、固定端の合力と負荷のつり合いを確認。両方を満たしたためPASS。
単位系は mm・N・MPa。検証では形状、材料、境界条件、荷重を変えず、メッシュ密度だけを変更しました。
長さ L = 100 mm
断面 b × h = 10 × 10 mm
断面二次モーメント I = 833.333 mm⁴
E = 210000 MPa
ν = 0.30
x=0面:Ux, Uy, Uz = 0
x=100 mm面に合計 P = 100 Nを-Z方向へ分配。線形静解析。主検証要素は8節点HEXのC3D8(完全積分)。
Euler–Bernoulli梁理論の先端変位 δ = PL³/(3EI) = 0.190476 mm をTargetとしました。80 → 640 → 5120要素で誤差が単調に減少しています。
| Mesh | 先端平均 |Uz| [mm] | 理論解との差 | 固定端反力 Z [N] |
|---|---|---|---|
| 80 C3D8 | 0.166823 | 12.418% | 100.0 |
| 640 C3D8 | 0.183789 | 3.511% | 100.0 |
| 5120 C3D8 | 0.188794 | 0.883% | 100.0 |
下図は80要素の初期デモで、FRDをPython + VTKで直接読み、変形形状に節点von Mises応力を表示したものです。最大節点von Mises応力は47.243 MPaでした。

入力条件と結果だけでなく、元のCubit export、CalculiX入力、細分化ケース、検証値、後処理コードを公開します。将来のAIや第三者が同じ比較契約を再実行できます。
Cubitから出力したmeshと、初期デモで使用したCalculiX C3D8 input。
80 / 640 / 5120 C3D8要素の入力デッキ。形状・材料・境界条件・総荷重は同一です。
検証値CSV・機械可読JSON・元デモのVTK後処理コード。
短い自然言語指示から、Cubit形状・HEXメッシュ、solver input、CalculiX実行、結果読込、VTK/PNG出力まで一連の処理を実機で完走。
同一条件のC3D8メッシュを細分化すると、先端変位が梁理論へ収束し、5120要素で誤差1%未満。全meshで反力つり合いも確認。
局所最大応力の厳密検証、非線形、接触、塑性、複雑形状、任意モデルへの一般化はこのケースのPASS範囲に含めません。