はじめに：なぜパラメータ最適化は堂々巡りになるのか
金属レーザ粉末床溶融結合（LPBF）および高分子焼結プロセスの現場調整では、次のような場面がよく見られる。エンジニアがレーザ出力を250 Wから280 Wへ変更すると、密度が0.3%上昇したため、さらに出力を上げる。300 Wまで上げるとキーホールが発生し、再び出力を下げて調整する。このように十数個の試験片を作製しても、得られるのは断片的な経験曲線に過ぎず、装置や粉末ロットが変われば最初からやり直しになる。問題はエンジニアの努力不足ではなく、単一因子を順番に変える手法そのものが交互作用を扱えないことにある。出力と走査速度は体積エネルギー密度を共に決定し、粉末層厚さはさらに溶融池の放熱条件を変化させる。3つの因子は相互に結合しており、一度に1変数だけを動かしても全体像は決して見えない。
実験計画法（Design of Experiments、DoE）は、こうした多因子結合の課題を解決するためのエンジニアリング手法である。最小限の試験片数で、因子の主効果、交互作用および応答曲面を一括して推定し、信頼区間を伴う承認可能な工程ウィンドウを提示できる。本稿では3Dプリンティング工程開発を対象に、直交スクリーニング、応答曲面モデリング、ウィンドウ検証までの実行可能な一連の方法を示す。
1. 3種類の実験を区別する：スクリーニング、モデリング、検証
工程開発は1回の実験だけで完了するものではない。通常は3層で進められ、各層で目的、因子数、試験片数が異なる。
第1層はスクリーニング実験であり、6～10個の候補因子から真に有意な2～4因子を選び出すことが目的である。代表的な候補因子には、レーザ出力、走査速度、ハッチ間隔、層厚、予熱温度、走査戦略（ストライプ幅）、層間回転角、保護ガス流量がある。スクリーニング段階では部分因子実験またはPlackett-Burman設計を用い、8～12回の実験で完了できる。応答は密度と表面粗さという2つの概略指標のみを確認する。
第2層はモデリング実験であり、スクリーニングで選定された2～3個の重要因子について応答曲面モデルを構築することが目的である。一般に中心複合計画法（CCD）またはBox-Behnken計画を用いる。3因子CCDでは、6個の軸点、8個の立方点、6回の中心点繰返しを含み、合計20回の実験となる。応答は密度、引張強さ、寸法偏差、表面粗さRaの4項目まで拡張する。
第3層は検証実験である。モデルが予測する最適領域内で3～5個の確認点を選び、繰り返し造形を行う。モデル予測値と実測値の差が±5%以内かを確認する。検証に合格した後、量産用工程票のパラメータ範囲として、ウィンドウ境界を10%狭める。
2. 直交実験：最少の試験片で最大の空間をカバーする
因子数が少なく、各因子の水準数が等しい場合、直交表は最も経済的な選択肢である。316Lステンレス鋼のLPBFを例に、レーザ出力P、走査速度v、ハッチ間隔hの3因子・3水準を評価する場合、完全因子実験には27回が必要となるが、L9(3^4)直交表を使用すれば9回で済む。
具体的には、Pを200/240/280 W、vを800/1000/1200 mm/s、hを0.09/0.11/0.13 mmに設定する。L9に従って実験を割り付けた後、各水準組合せで10×10×10 mmの立方体試験片を1組造形し、アルキメデス法で密度を測定すると同時に、表面粗さRaを記録する。
データ処理では、範囲分析により迅速に順位付けできる。各因子について各水準における平均応答を算出し、範囲R = max（平均応答）− min（平均応答）とする。Rが大きいほど、その因子の影響は大きい。R(P)=2.1%、R(v)=1.4%、R(h)=0.6%であれば、影響度の順序は出力 > 速度 > 間隔となる。さらに分散分析（ANOVA）を実施し、F検定によって有意性を判定する。F値が臨界値を上回り、かつp値が0.05未満の因子を次段階のモデリングに進める。
直交実験の限界は、因子間に交互作用がないことを前提としている点である。出力と速度の間に強い結合が疑われる場合、すなわち体積エネルギー密度E = P/(v·h·t)自体が積の関係である場合は、直交表を無理に継続使用するのではなく、モデリング段階で応答曲面計画へ移行すべきである。
3. 応答曲面法：離散点から連続曲面へ
応答曲面法（RSM）の本質は、二次多項式によって応答と因子の関係を近似することである。
Y = β0 + Σβi·xi + Σβii·xi² + Σβij·xi·xj + ε
ここで、交差項βijが交互作用を定量的に表す。CCDを例にすると、3因子計画の重要なパラメータは軸距離αである。回転可能計画では、α = 2^(k/4) = 2^(3/4) ≈ 1.682（kは因子数）となる。軸点は符号化空間で±1.682、立方点は±1に位置し、中心点は純粋誤差を推定するために6回繰り返す。
フィッティング完了後には、必ず3つの検証を行う。第1に、適合不足検定（lack-of-fit）であり、p値が0.05より大きければ、モデル形式とデータの間に有意な乖離がないことを示す。第2に、決定係数R²と調整済みR²であり、両者の差は0.2未満であるべきで、そうでなければ過学習の可能性がある。第3に、残差の正規性および等分散性を確認する。実測では、R²が0.95に達している一方で調整済みR²が0.86にとどまる場合がある。この場合は、有意でない三次項を除外して再フィッティングするべきである。
多応答最適化は3Dプリンティングにおける難所である。密度は最大化し、粗さは最小化し、寸法偏差はゼロに近づける必要があるが、この3つの目的はしばしば相反する。実務では望ましさ関数（desirability function）を用いて各応答を0～1の範囲に写像し、その幾何平均から総合望ましさDを算出する。Dが最大となる点でパラメータ組合せを決定する。重みは顧客の受入優先順位に応じて配分し、例えば航空宇宙部品では密度を0.5、粗さを0.3、寸法偏差を0.2と設定する。
4. 工程ウィンドウの設定：最適点を承認可能な範囲に変える
最適点を見つけることは第一歩に過ぎず、量産に必要なのは範囲である。実行可能な手法は、重ね合わせ等高線図（overlay contour plot）を用いることである。各応答の合格閾値を等高線として描き（密度 ≥ 99.5%、Ra ≤ 12 μm、寸法偏差 ≤ ±0.1 mm）、すべての等高線が重なる領域を工程ウィンドウとする。
ウィンドウ決定後には、必ず安全マージンを確保しなければならない。推奨される方法は、ウィンドウ境界を内側へ10%縮小することである。例えば、モデルが出力の実行可能範囲を235～265 Wと示した場合、工程票には238～262 Wと記載する。その理由は、粉末ロットのばらつき、保護ガス純度の変動、光学系の出力低下によって、実際の工程点が徐々に移動するためである。マージンを残すことで、工程能力指数Cpkを1.0から1.33以上へ高めることができる。
ウィンドウ検証では3ロットによる確認を行う。各ロット5個とし、ウィンドウ中心、上限境界の内側、下限境界の内側でそれぞれ1組を造形する。連続する3ロットがすべて合格し、ロット間の平均値ドリフトが応答標準偏差の1/3未満であれば、ウィンドウは安定していると判定できる。検証データは同時に工程台帳へ記録し、後続のトレーサビリティにおけるベースラインとする。
5. よくある誤解とチェックリスト
誤解その1：試験片数だけをコストと考えること。実際には、9回の直交実験と20回のCCD実験を合わせた試験片コストは、目的なく30回の単一因子調整を行うよりはるかに低く、しかも再利用可能なモデルが得られる。
誤解その2：ブロック効果を無視すること。20回のCCD実験を3日間に分けて実施する場合、装置状態、粉末ロット、環境温湿度の変化が残差に混入する。正しい方法は、日付をブロック因子として計画に組み込むか、ブロック内に中心点の繰返しを配置してドリフトを監視することである。
誤解その3：密度だけを見ること。密度が99.5%に達する工程点でも、機械特性が合格するとは限らない。未溶融欠陥は密度測定では過小評価される可能性がある。少なくとも引張強さと伸びを同時に評価し、必要に応じて顕微鏡組織観察および産業用CTの抜取り検査を追加することが望ましい。
実行チェックリスト：因子水準は装置が安定して運転できる範囲をカバーしているか。純粋誤差を推定するための中心点繰返しを設定したか。適合不足検定を実施したか。多応答の重みは顧客の受入基準と一致しているか。ウィンドウには10%のマージンを確保したか。検証ロットはウィンドウ境界をカバーしているか。
まとめ
実験計画法は、3Dプリンティング工程開発を経験的な試行錯誤から、定量化・再現可能なエンジニアリング活動へ変える。直交実験では9～12回で因子スクリーニングを完了し、応答曲面法では約20回で交互作用を含む二次モデルを構築できる。望ましさ関数は多目的間の競合を解決し、重ね合わせ等高線は安全マージンを備えた承認可能な工程ウィンドウを示す。設計データから実製品の納入までを扱うBlueprint 3Dのようなライフサイクル管理プラットフォームにとって、DoEの価値はパラメータ調整にとどまらない。工程知識を検索可能、再利用可能、検証可能なパラメータモデルとして蓄積し、あらゆる工程開発を次回の出発点に変え、毎回ゼロから始めることを防ぐ。