パワーサイクル試験 — CAE用語解説

自動車のエンジンや駆動系部品に、実際の走行で発生するトルク変動を実験室で再現する試験だ。例えば、エンジンがアイドリング状態から急加速する際の、

良い着眼点だ。パワーサイクルは主に低サイクル疲労（LCF）を対象とする。1回のサイクルで塑性変形が生じるほどの大きな応力振幅が加わる。目安としては、応力振幅が材料の降伏強さの50%以上になるようなサイクルだ。例えば、S45C鋼なら降伏強さは約350MPaだから、175MPa以上の応力振幅が繰り返される試験を想定する。一方、高サイクル疲労（HCF）はエンジンの回転数による振動など、小さな応力が何百万回も加わる現象で、別の試験規格で評価される。

実車で計測したトルクや回転数の時系列データ（実走行データ）を起点にするが、そのままでは試験時間が膨大になる。そこで、統計処理や雨流計数法を使って、損傷に寄与する大きな負荷変動のみを抽出し、試験用の「代表ブロック」を作成する。例えば、市街地、山道、高速道路の各コースデータから、それぞれ損傷の大きいシーケンスを数秒〜数十秒分切り出して連結する。これが試験1サイクルとなり、それを数千〜数万回繰り返す。

多くの場合、準静解析アプローチを取る。つまり、負荷履歴中の複数の特徴的な「タイムステップ」に対して静解析を行い、得られた応力・ひずみ結果を後処理で繋ぎ合わせる。過渡応答（動的）解析は、ギア噛み合い時の衝撃（ラッタ）など、慣性効果が無視できない現象に限定される。Ansys MechanicalやAbaqus/Standardでは、複数の荷重ケースを定義し、各ケースの結果から疲労解析ツール（nCode DesignLifeやFE-SAFE）で寿命を計算するのが標準的なワークフローだ。

最大・最小だけでは不十分だ。材料の塑性や残留応力の影響を評価するには、負荷経路（ローディングパス）が重要になる。例えば、トルクが正転最大→ゼロ→逆転最大→ゼロと変化する場合、単純な往復運動ではなく、応力-ひずみ曲線がヒステリシスループを描く。これを正確に捉えるには、少なくとも各極値点と、経路の中間点を数点追加して解析する必要がある。そうしないと、ひずみ振幅

核心を突く質問だ。弾性解析の応力結果をそのまま使うと、局部応力が材料の降伏強さを超える領域では、実際には起こり得ない非常に高い応力値が算出されてしまう。これを使うと、寿命が過小評価（安全側だが過剰設計）になる。低サイクル疲労の駆動力は「ひずみ振幅」であり、塑性変形を考慮した「真のひずみ」を評価する必要がある。そのため、Chabocheモデルなどの非線形硬化則を用いた弾塑性解析が必須となる。Abaqusでは、このような材料モデルを*PLASTICオプションで定義する。

メッシュ以前に「境界条件と負荷の適用方法」だ。パワーサイクル試験では、試験機のジグで部品をどのように保持し、どこからトルクを加えるかを忠実に再現する。例えば、フランジ部をボルト締結で固定する場合、ボルト穴周辺の接触条件と締結力の初期値設定を誤ると、全体の剛性や応力分布が大きく変わってしまう。まずは単一の静荷重をかけた時の部品の変形モードや反力を計算し、試験機のロードセル値や変位計測値と比較して、FEMモデルの基本挙動を検証する「モデル相関」作業から始める。

最も重要な入力データの一つだ。自社で材料試験（ひずみ制御疲労試験）を行って取得するのが理想だが、コストがかかる。実務では、公開データベースや材料メーカーのカタログを参照する。例えば、日本鉄鋼連盟の「鉄鋼材料の疲労設計データシート」や、MSC FatigueやnCode DesignLifeに内蔵された材料ライブラリを使う。ここで注意すべきは、データが「平滑材」のものである点だ。実際の部品には切欠きがあるため、解析結果に応力集中係数Ktを考慮し、さらに表面粗さや加工硬化の影響を係数で補正する必要がある。

系統的に切り分ける必要がある。第一に「負荷入力」：CAEモデルに与えた力/トルクの振幅と試験機の設定値は一致しているか？ 第二に「応力・ひずみ結果」：試験中に部品に貼ったひずみゲージの測定値と、対応するFEMノードのひずみ値は一致するか？ 第三に「材料データ」：使用したε-N曲線は、実際の材料ロット、熱処理条件を反映しているか？ 第四に「損傷累積則」：Miner則のような線形累積則を使っている場合、負荷順序効果（大きい負荷の後の小さな負荷は損傷が小さいなど）を考慮できておらず、これが誤差の原因になることもある。

両者とも優れた疲労解析専用ソフトウェアだが、統合環境とアルゴリズムの細部に違いがある。nCode DesignLifeはAnsys Workbench環境と深く統合されており、Mechanicalで得られた応力/ひずみの時系列データをシームレスに受け渡せる。一方、MSC Fatigueは元来MSC Nastranと連携が強く、.op2結果ファイルを直接読み込む形式だ。アルゴリズム面では、nCodeはDang Vanの多軸疲労基準の実装が強く、MSC Fatigueは独自の「ストレイン・ライフ・ポストプロセッサ」による詳細な弾塑性修正オプションが特徴と言われる。しかし、基本的な疲労寿命予測精度は、適切に使えば大差ない。

Abaqus/CAE Fatigueは、基本的な高サイクル疲労（応力ベース）や低サイクル疲労（ひずみベース）の評価をAbaqus環境内で完結させたい場合に便利だ。しかし、機能の幅と深さは専用ソフトに及ばない。例えば、複雑な実走行負荷履歴の処理能力、豊富な材料ライブラリ、表面処理や腐食環境の影響を考慮した高度な修正係数、試験データとの相関分析ツールなどは、nCodeやFE-SAFEの方が充実している。プロジェクトで高度な疲労信頼性設計が要求される場合は、Abaqusで応力/ひずみを計算し、nCodeなどで詳細な疲労評価を行う連携が一般的だ。

構造解析（応力計算）までは可能だが、疲労寿命評価までのフルプロセスは難しい。CalculiXやCode_Asterは強力な非線形静解析・過渡応答解析ソルバーであり、適切なモデルを作れば負荷サイクル中の応力・ひずみ履歴を計算できる。問題はその後だ。これらのソフトには疲労専用のポストプロセッサが標準では付属しておらず、出力された時系列応力データを、自分でスクリプト（Pythonなど）を書いて雨流計数法やMiner則を適用し、寿命を計算する必要がある。研究目的や、非常に限定的な部品評価であれば可能だが、実務での効率性と信頼性を求めるなら、商用の疲労専用ソフトの投資価値は高い。

低サイクル疲労解析で典型的な問題だ。まず、収束しないタイムステップ付近で、部品のどこかに「過度の塑性変形」が発生していないか確認せよ。局部で材料が破断に近い状態になると、剛性マトリクスが特異になり収束しない。対策は、(1) 荷重増分を小さくする（Abaqusでは*STATICオプションのINITIALやMINIMUMを調整）、(2) 問題箇所のメッシュを細かくする（塑性ひずみの勾配を捉えるため）、(3) 材料モデルの塑性硬化パラメータを見直す（特に降伏後の軟化挙動を定義していないか）。また、接触条件が急激に変化する場合も原因となる。

絶対に違う。これは「計算上の応力/ひずみ振幅が、材料の疲労限界（エンドランス）以下だった」ことを意味する。最も多い原因は、弾性解析の応力をそのまま使った場合だ。先述したように、応力集中部では弾性計算値が実際より高く出るが、逆に局部塑性による応力再配分を考慮しないため、平均応力が過大評価され、結果として応力振幅が小さく見積もられてしまう。対策は、弾塑性解析を実施して真のひずみ振幅を求めること。あるいは、弾性解析結果を使う場合は、Neuber法や線形化法といった「弾塑性修正」を必ず適用することだ。nCodeでは「Strain-Life (EN)」メソッドでこの修正オプションを選択できる。

これが低サイクル疲労解析の難しい点だ。寿命がひずみ振幅