複合材料特性計算 戻る
材料工学・複合材料

複合材料(CFRP)特性計算シミュレーター

炭素繊維/ガラス繊維強化プラスチックの微視力学モデルで特性を計算。ハルピン-ツァイモデルによるE1・E2・G12の繊維体積率依存性を可視化。古典積層板理論も対応。

材料選択
繊維の種類
マトリックスの種類
繊維体積率
Vf(繊維体積率) 0.60
単方向層の弾性定数
E1 (GPa)
E2 (GPa)
G12 (GPa)
ν12

計算モデル

E1 = Ef·Vf + Em·(1-Vf) 【則混合則】
E2 = Em·(1+ξ·η·Vf)/(1-η·Vf) 【H-T】
η = (Ef/Em-1)/(Ef/Em+ξ)
ξ=2 (E2), ξ=1 (G12)
縦の破線 = 現在のVf設定値。グラフはVf=0〜0.75の全範囲を表示。

複合材料(CFRP)特性計算シミュレーターとは

🧑‍🎓
このシミュレーターで計算できる「E1」と「E2」って何が違うんですか?同じ材料なのにヤング率が二つあるのが不思議です。
🎓
ざっくり言うと、繊維方向かどうかの違いだね。CFRPのような繊維強化プラスチックは、繊維の並んだ方向(1方向)とそれに垂直な方向(2方向)で硬さが全然違うんだ。E1は繊維方向のヤング率で、E2は横方向のヤング率。シミュレーターの「繊維の種類」を「炭素繊維」から「ガラス繊維」に変えてみると、E1はガラクッと下がるけど、E2はあまり変わらないのがわかるよ。
🧑‍🎓
え、そうなんですか!じゃあ、繊維体積率「Vf」のスライダーを動かすと、E1とE2は同じように変化するんですか?
🎓
そこが面白いところで、計算モデルが全く違うんだ。E1は単純な足し算(則混合則)で、Vfを上げると直線的に増える。でもE2は「ハルピン-ツァイ」という少し複雑な式で計算するから、Vfを上げてもE1ほどは増えない。実際にVfを0から1まで動かしてグラフを見比べてみて。自動車の軽量化設計では、この非線形な挙動を理解して最適なVfを決めるんだ。
🧑‍🎓
なるほど!「ハルピン-ツァイ」って聞きましたが、なんでE2とG12(せん断剛性)で計算式が違うんですか?パラメータ「ξ」が2とか1とかありますよね。
🎓
良いところに気が付いたね。ξ(グザイ)は「強化係数」で、繊維の形状や荷重の掛かり方で効き方が変わる経験的なパラメータなんだ。引っ張り(E2)では繊維が長いほど効くのでξ=2、せん断(G12)ではξ=1がよく使われる。現場では、この値を実験データにフィッティングして微調整することも多いよ。シミュレーターで「マトリックスの種類」をエポキシからポリエステルに変えると、ξの影響の大きさが実感できるはずだ。

物理モデルと主要な数式

繊維方向のヤング率 E1 は、繊維とマトリックスの体積分率に応じた単純な加重平均(則混合則)で計算されます。これは繊維とマトリックスが並列に荷重を分担すると仮定したモデルです。

$$E_1 = E_f V_f + E_m (1 - V_f)$$

$E_f$: 繊維のヤング率, $E_m$: マトリックス(母材)のヤング率, $V_f$: 繊維体積率 (0 ≤ V_f ≤ 1)

横方向ヤング率 E2 や面内せん断剛性 G12 は、ハルピン-ツァイの半経験式で計算されます。繊維とマトリックスの相互作用をパラメータξで表現し、則混合則よりも実験値との一致が良いモデルです。

$$E_2 = E_m \frac{1 + \xi \eta V_f}{1 - \eta V_f}, \quad \eta = \frac{(E_f / E_m) - 1}{(E_f / E_m) + \xi}$$

$E_m$: マトリックスのヤング率, $\xi$: 強化係数(E2の計算では2、G12の計算では1が典型的), $V_f$: 繊維体積率。この式は、繊維の存在が横方向の剛性をどの程度「強化」するかを表しています。

実世界での応用

航空機・宇宙構造物の軽量化設計:CFRPは高い比強度・比剛性が要求される機体の主翼や胴体に広く採用されています。シミュレーターで計算したE1, E2, G12等の基本物性を入力し、積層板理論で多層構造の剛性や強度を予測することで、最適な積層構成を決定します。

自動車部品(ドライブシャフト、ボンネット)の開発:燃費向上のため金属部品をCFRPやGFRPで置き換える際、繊維体積率Vfと繊維配向をどうするかが鍵です。ツァイ-ウー破壊基準と組み合わせて、ねじりや曲げ荷重下での破壊を予測し、安全性を確保します。

スポーツ用品(テニスラケット、ゴルフシャフト)の高性能化:打球感や飛距離を決めるのは複合材料積層板の曲げ剛性やねじり剛性です。設計者はシミュレーターで基本物性を把握した上で、繊維をどの角度に何層重ねるかを繰り返し検討し、理想的な性能を引き出します。

風力発電ブレードの大型化・長寿命化:数十メートルに及ぶ大型ブレードはGFRPが主流です。自重による曲げと遠心力が複合する過酷な環境下で、長期にわたって疲労破壊を起こさない設計が必須です。微視力学モデルは材料開発の段階から信頼性の高い物性値を提供します。

よくある誤解と注意点

まず、「計算結果がそのまま設計値」と思わないことが大事だよ。シミュレーターは「均質で完璧な単層板」を計算している。実際の材料には繊維の波打ちや気泡、界面強度のバラつきがあるから、安全率を見込むのは必須。例えば、ツールでE1=150GPaと出ても、初期設計では120GPaくらいで考えるのが現場の知恵だ。

次に、「繊維体積率Vfは高ければ高いほど良い」という思い込み。確かにE1はVfに比例して上がるけど、Vfを上げすぎると樹脂が繊維を十分に濡らせず、逆に強度が低下したり、製造コストが跳ね上がったりする。自動車部品ではコストと性能のバランスから、Vf=0.5〜0.6がよく選ばれる範囲だ。シミュレーターでVfを0.7以上に設定してみると、E2やG12の上昇が頭打ちになる様子がグラフで確認できるはず。これが「もうこれ以上上げても意味が薄い」サインだ。

最後に、基本物性(E1, E2など)と「強度」を混同しないこと。このツールで計算できるのは「剛性」(変形のしにくさ)であって、「強度」(壊れる限界)ではない。強度は別の破壊基準(ツァイ-ウーなど)で評価する必要がある。剛性が高くても脆くて割れやすい、ということはよくある。例えば、E1が同じ150GPaの二つの材料でも、使用する炭素繊維の種類(高強度型か高弾性型か)で引張強度は全く異なるんだ。

関連する工学分野

このツールで得られる基本物性は、積層板理論(クラシックラミネート理論)への入り口だ。単層のE1, E2, G12, ν12が分かれば、それらを異なる角度(0度、45度、90度など)で何層も重ねたときの全体の剛性(A行列)や曲げ剛性(D行列)を計算できる。飛行機の翼やF1マシンのモノコックは、まさにこの理論に基づいて何十層もの積層構成が最適化されている。

また、これらの物性値は有限要素法(FEA)シミュレーションの必須入力データになる。CAEソフトでCFRP部品を解析する時は、各積層方向ごとに「直交異方性材料」としてこれらの値を設定する。例えば、ドローン用アームを軽量化設計する場合、このツールで材料定数を求め、FEAで変形と固有振動数を評価する、という流れが一般的だ。

さらに深掘りすると、破壊力学疲労解析にも繋がる。ツール内のツァイ-ウー破壊基準は複合材料特有の破壊モード(繊維破断、マトリックス割れ、界面剥離など)を統一的に評価する第一歩。得られた強度比を元に、き裂が進展する寿命予測や信頼性設計を行う分野だ。

発展的な学習のために

まず次のステップは、「則混合則」と「ハルピン-ツァイ則」の式の導出背景を理解すること。則混合則が「並列モデル」と「直列モデル」のどちらに相当するか、図を描いて考えてみよう。ハルピン-ツァイ則のパラメータξが、実は「繊維のアスペクト比(長さ/直径)」や「応力集中」と関係していることを知ると、式の見え方が変わる。

数学的には、テンソルの概念に慣れることが近道だ。複合材料の応力-ひずみ関係は、単なるフックの法則(σ=Eε)ではなく、行列(正確には4階のテンソル)で表される。$$ \begin{bmatrix} \sigma_1 \\ \sigma_2 \\ \tau_{12} \end{bmatrix} = \begin{bmatrix} Q_{11} & Q_{12} & 0 \\ Q_{12} & Q_{22} & 0 \\ 0 & 0 & Q_{66} \end{bmatrix} \begin{bmatrix} \varepsilon_1 \\ \varepsilon_2 \\ \gamma_{12} \end{bmatrix} $$ のような剛性マトリックス[Q]を、ツールで求めたE1, E2, G12, ν12から計算できるようになれば、積層板理論の核心に迫れる。

実務を目指すなら、「製造プロセスが物性に与える影響」を学ぼう。同じ材料仕様でも、オートクレーブ成形とRTM(樹脂転写成形)では繊維体積率や界面状態が変わり、強度が10〜20%も変わることがある。シミュレーターは「理想状態」を計算するツールだから、その結果をどう現実に落とし込むか、という視点がエンジニアの腕の見せ所だ。