第1章：積層造形の基礎

flowchart LR
    A[1986
SLA発明
Chuck Hull] --> B[1988
SLS登場
Carl Deckard]
    B --> C[1992
FDM特許
Stratasys社]
    C --> D[2005
RepRap
オープンソース化]
    D --> E[2012
金属AM普及
EBM/SLM]
    E --> F[2023
産業化加速
大型・高速化]

    style A fill:#e3f2fd
    style B fill:#fff3e0
    style C fill:#e8f5e9
    style D fill:#f3e5f5
    style E fill:#fce4ec
    style F fill:#fff9c4
        

flowchart TD
    AM[積層造形
7つのプロセス] --> MEX[Material Extrusion
材料押出]
    AM --> VPP[Vat Photopolymerization
液槽光重合]
    AM --> PBF[Powder Bed Fusion
粉末床溶融結合]
    AM --> MJ[Material Jetting
材料噴射]
    AM --> BJ[Binder Jetting
結合剤噴射]
    AM --> SL[Sheet Lamination
シート積層]
    AM --> DED[Directed Energy Deposition
指向性エネルギー堆積]

    MEX --> MEX_EX[FDM/FFF
低コスト・普及型]
    VPP --> VPP_EX[SLA/DLP
高精度・高表面品質]
    PBF --> PBF_EX[SLS/SLM/EBM
高強度・金属対応]

    style AM fill:#f093fb
    style MEX fill:#e3f2fd
    style VPP fill:#fff3e0
    style PBF fill:#e8f5e9
    style MJ fill:#f3e5f5
    style BJ fill:#fce4ec
    style SL fill:#fff9c4
    style DED fill:#fce4ec
        

1.3 STLファイル形式とデータ処理

1.3.1 STLファイルの構造

STL（STereoLithography）は、AMで最も広く使用される3Dモデルファイル形式で、1987年に3D Systems社が開発しました。STLファイルは物体表面を三角形メッシュ（Triangle Mesh）の集合として表現します。

STLファイルの基本構造

ASCII STL形式の例：

solid cube
  facet normal 0 0 1
    outer loop
      vertex 0 0 10
      vertex 10 0 10
      vertex 10 10 10
    endloop
  endfacet
  facet normal 0 0 1
    outer loop
      vertex 0 0 10
      vertex 10 10 10
      vertex 0 10 10
    endloop
  endfacet
  ...
endsolid cube

STLフォーマットの2つの種類：

1. ASCII STL: 人間が読めるテキスト形式。ファイルサイズ大（同じモデルでBinaryの10-20倍）。デバッグ・検証に有用。
2. Binary STL: バイナリ形式、ファイルサイズ小、処理高速。産業用途で標準。構造：80バイトヘッダー + 4バイト（三角形数） + 各三角形50バイト（法線12B + 頂点36B + 属性2B）。

1.3.2 STLファイルの重要概念

1. 法線ベクトル（Normal Vector）

各三角形面には法線ベクトル（外向き方向）が定義され、物体の「内側」と「外側」を区別します。法線方向は右手の法則で決定されます：

頂点順序ルール： 頂点v1, v2, v3は反時計回り（CCW: Counter-ClockWise）に配置され、外から見て反時計回りの順序で法線が外向きになります。

2. 多様体（Manifold）条件

STLメッシュが3Dプリント可能であるためには、多様体（Manifold）でなければなりません：

• エッジ共有: すべてのエッジ（辺）は正確に2つの三角形に共有される
• 頂点共有: すべての頂点は連続した三角形扇（fan）に属する
• 閉じた表面: 穴や開口部がなく、完全に閉じた表面を形成
• 自己交差なし: 三角形が互いに交差・貫通していない

1.3.3 STLファイルの品質指標

STLメッシュの品質は以下の指標で評価されます：

1. 三角形数（Triangle Count）: 通常10,000-500,000個。過少（粗いモデル）または過多（ファイルサイズ大・処理遅延）は避ける。
2. エッジ長の一様性: 極端に大小の三角形が混在すると造形品質低下。理想的には0.1-1.0 mm範囲。
3. アスペクト比（Aspect Ratio）: 細長い三角形（高アスペクト比）は数値誤差の原因。理想的にはアスペクト比 < 10。
4. 法線の一貫性: すべての法線が外向き統一。反転法線が混在すると内外判定エラー。

1.3.4 Pythonライブラリによる STL処理

PythonでSTLファイルを扱うための主要ライブラリ：

1. numpy-stl: 高速STL読込・書込、体積・表面積計算、法線ベクトル操作。シンプルで軽量。
2. trimesh: 包括的な3Dメッシュ処理ライブラリ。メッシュ修復、ブーリアン演算、レイキャスト、衝突検出。多機能だが依存関係多い。
3. PyMesh: 高度なメッシュ処理（リメッシュ、サブディビジョン、フィーチャー抽出）。インストールやや複雑。

numpy-stlの基本的な使用法：

from stl import mesh
import numpy as np

# STLファイルを読み込み
your_mesh = mesh.Mesh.from_file('model.stl')

# 基本的な幾何情報
volume, cog, inertia = your_mesh.get_mass_properties()
print(f"Volume: {volume:.2f} mm³")
print(f"Center of Gravity: {cog}")
print(f"Surface Area: {your_mesh.areas.sum():.2f} mm²")

# 三角形数
print(f"Number of Triangles: {len(your_mesh.vectors)}")

1.4 スライシングとツールパス生成

STLファイルを3Dプリンタが理解できる指令（G-code）に変換するプロセスをスライシング（Slicing）といいます。このセクションでは、スライシングの基本原理、ツールパス戦略、そしてG-codeの基礎を学びます。

1.4.1 スライシングの基本原理

スライシングは、3Dモデルを一定の高さ（レイヤー高さ）で水平に切断し、各層の輪郭を抽出するプロセスです：

flowchart TD
    A[3Dモデル
STLファイル] --> B[Z軸方向に
層状にスライス]
    B --> C[各層の輪郭抽出
Contour Detection]
    C --> D[シェル生成
Perimeter Path]
    D --> E[インフィル生成
Infill Path]
    E --> F[サポート追加
Support Structure]
    F --> G[ツールパス最適化
Retraction/Travel]
    G --> H[G-code出力]

    style A fill:#e3f2fd
    style H fill:#e8f5e9
        

レイヤー高さ（Layer Height）の選択

レイヤー高さは造形品質と造形時間のトレードオフを決定する最重要パラメータです：

1.4.2 シェルとインフィル戦略

シェル（外殻）の生成

シェル（Shell/Perimeter）は、各層の外周部を形成する経路です：

• シェル数（Perimeter Count）: 通常2-4本。外部品質と強度に影響。
  • 1本: 非常に弱い、透明性高い、装飾用のみ
  • 2本: 標準（バランス良好）
  • 3-4本: 高強度、表面品質向上、気密性向上
• シェル順序: 内側→外側（Inside-Out）が一般的。外側→内側は表面品質重視時に使用。

インフィル（内部充填）パターン

インフィル（Infill）は内部構造を形成し、強度と材料使用量を制御します：

1.4.3 サポート構造の生成

オーバーハング角度が45度を超える部分は、サポート構造（Support Structure）が必要です：

サポートのタイプ

• Linear Support（直線サポート）: 垂直な柱状サポート。シンプルで除去しやすいが、材料使用量多い。
• Tree Support（ツリーサポート）: 樹木状に分岐するサポート。材料使用量30-50%削減、除去しやすい。CuraやPrusaSlicerで標準サポート。
• Interface Layers（接合層）: サポート上面に薄い接合層を設ける。除去しやすく、表面品質向上。通常2-4層。

サポート設定の重要パラメータ

1.4.4 G-codeの基礎

G-codeは、3DプリンタやCNCマシンを制御する標準的な数値制御言語です。各行が1つのコマンドを表します：

主要なG-codeコマンド

G-codeの例（造形開始部分）

; === Start G-code ===
M140 S60       ; ベッドを60°Cに加熱開始（非待機）
M104 S210      ; ノズルを210°Cに加熱開始（非待機）
G28            ; 全軸ホーミング
G29            ; オートレベリング（ベッドメッシュ計測）
M190 S60       ; ベッド温度到達を待機
M109 S210      ; ノズル温度到達を待機
G92 E0         ; 押出量をゼロリセット
G1 Z2.0 F3000  ; Z軸を2mm上昇（安全確保）
G1 X10 Y10 F5000  ; プライム位置へ移動
G1 Z0.3 F3000  ; Z軸を0.3mmへ降下（初層高さ）
G1 X100 E10 F1500 ; プライムライン描画（ノズル詰まり除去）
G92 E0         ; 押出量を再度ゼロリセット
; === 造形開始 ===

1.4.5 主要スライシングソフトウェア

1.4.6 ツールパス最適化戦略

効率的なツールパスは、造形時間・品質・材料使用量を改善します：

• リトラクション（Retraction）: 移動時にフィラメントを引き戻してストリング（糸引き）を防止。
  • 距離: 1-6mm（ボーデンチューブ式は4-6mm、ダイレクト式は1-2mm）
  • 速度: 25-45 mm/s
  • 過度なリトラクションはノズル詰まりの原因
• Z-hop（Z軸跳躍）: 移動時にノズルを上昇させて造形物との衝突を回避。0.2-0.5mm上昇。造形時間微増だが表面品質向上。
• コーミング（Combing）: 移動経路をインフィル上に制限し、表面への移動痕を低減。外観重視時に有効。
• シーム位置（Seam Position）: 各層の開始/終了点を揃える戦略。
  • Random: ランダム配置（目立たない）
  • Aligned: 一直線に配置（後加工でシームを除去しやすい）
  • Sharpest Corner: 最も鋭角なコーナーに配置（目立ちにくい）

Example 1: STLファイルの読み込みと基本情報取得

# ===================================
# Example 1: STLファイルの読み込みと基本情報取得
# ===================================

import numpy as np
from stl import mesh

# STLファイルを読み込む
your_mesh = mesh.Mesh.from_file('model.stl')

# 基本的な幾何情報を取得
volume, cog, inertia = your_mesh.get_mass_properties()

print("=== STLファイル基本情報 ===")
print(f"Volume: {volume:.2f} mm³")
print(f"Surface Area: {your_mesh.areas.sum():.2f} mm²")
print(f"Center of Gravity: [{cog[0]:.2f}, {cog[1]:.2f}, {cog[2]:.2f}] mm")
print(f"Number of Triangles: {len(your_mesh.vectors)}")

# バウンディングボックス（最小包含直方体）を計算
min_coords = your_mesh.vectors.min(axis=(0, 1))
max_coords = your_mesh.vectors.max(axis=(0, 1))
dimensions = max_coords - min_coords

print(f"\n=== バウンディングボックス ===")
print(f"X: {min_coords[0]:.2f} to {max_coords[0]:.2f} mm (幅: {dimensions[0]:.2f} mm)")
print(f"Y: {min_coords[1]:.2f} to {max_coords[1]:.2f} mm (奥行: {dimensions[1]:.2f} mm)")
print(f"Z: {min_coords[2]:.2f} to {max_coords[2]:.2f} mm (高さ: {dimensions[2]:.2f} mm)")

# 造形時間の簡易推定（レイヤー高さ0.2mm、速度50mm/sと仮定）
layer_height = 0.2  # mm
print_speed = 50    # mm/s
num_layers = int(dimensions[2] / layer_height)
# 簡易計算: 表面積に基づく推定
estimated_path_length = your_mesh.areas.sum() / layer_height  # mm
estimated_time_seconds = estimated_path_length / print_speed
estimated_time_minutes = estimated_time_seconds / 60

print(f"\n=== 造形推定 ===")
print(f"レイヤー数（0.2mm/層）: {num_layers} 層")
print(f"推定造形時間: {estimated_time_minutes:.1f} 分 ({estimated_time_minutes/60:.2f} 時間)")

# 出力例:
# === STLファイル基本情報 ===
# Volume: 12450.75 mm³
# Surface Area: 5832.42 mm²
# Center of Gravity: [25.34, 18.92, 15.67] mm
# Number of Triangles: 2456
#
# === バウンディングボックス ===
# X: 0.00 to 50.00 mm (幅: 50.00 mm)
# Y: 0.00 to 40.00 mm (奥行: 40.00 mm)
# Z: 0.00 to 30.00 mm (高さ: 30.00 mm)
#
# === 造形推定 ===
# レイヤー数（0.2mm/層）: 150 層
# 推定造形時間: 97.2 分 (1.62 時間)

Example 2: メッシュの法線ベクトル検証

# ===================================
# Example 2: メッシュの法線ベクトル検証
# ===================================

import numpy as np
from stl import mesh

def check_normals(mesh_data):
    """STLメッシュの法線ベクトルの整合性をチェック

    Args:
        mesh_data: numpy-stlのMeshオブジェクト

    Returns:
        tuple: (flipped_count, total_count, percentage)
    """
    # 右手系ルールで法線方向を確認
    flipped_count = 0
    total_count = len(mesh_data.vectors)

    for i, facet in enumerate(mesh_data.vectors):
        v0, v1, v2 = facet

        # エッジベクトルを計算
        edge1 = v1 - v0
        edge2 = v2 - v0

        # 外積で法線を計算（右手系）
        calculated_normal = np.cross(edge1, edge2)

        # 正規化
        norm = np.linalg.norm(calculated_normal)
        if norm > 1e-10:  # ゼロベクトルでないことを確認
            calculated_normal = calculated_normal / norm
        else:
            continue  # 縮退三角形をスキップ

        # ファイルに保存されている法線と比較
        stored_normal = mesh_data.normals[i]
        stored_norm = np.linalg.norm(stored_normal)

        if stored_norm > 1e-10:
            stored_normal = stored_normal / stored_norm

        # 内積で方向の一致をチェック
        dot_product = np.dot(calculated_normal, stored_normal)

        # 内積が負なら逆向き
        if dot_product < 0:
            flipped_count += 1

    percentage = (flipped_count / total_count) * 100 if total_count > 0 else 0

    return flipped_count, total_count, percentage

# STLファイルを読み込み
your_mesh = mesh.Mesh.from_file('model.stl')

# 法線チェックを実行
flipped, total, percent = check_normals(your_mesh)

print("=== 法線ベクトル検証結果 ===")
print(f"総三角形数: {total}")
print(f"反転法線数: {flipped}")
print(f"反転率: {percent:.2f}%")

if flipped == 0:
    print("\n✅ すべての法線が正しい方向を向いています")
    print("   このメッシュは3Dプリント可能です")
elif percent < 5:
    print("\n⚠️ 一部の法線が反転しています（軽微）")
    print("   スライサーが自動修正する可能性が高い")
else:
    print("\n❌ 多数の法線が反転しています（重大）")
    print("   メッシュ修復ツール（Meshmixer, netfabb）での修正を推奨")

# 出力例:
# === 法線ベクトル検証結果 ===
# 総三角形数: 2456
# 反転法線数: 0
# 反転率: 0.00%
#
# ✅ すべての法線が正しい方向を向いています
#    このメッシュは3Dプリント可能です

Example 3: マニフォールド性のチェック

# ===================================
# Example 3: マニフォールド性（Watertight）のチェック
# ===================================

import trimesh

# STLファイルを読み込み（trimeshは自動で修復を試みる）
mesh = trimesh.load('model.stl')

print("=== メッシュ品質診断 ===")

# 基本情報
print(f"Vertex count: {len(mesh.vertices)}")
print(f"Face count: {len(mesh.faces)}")
print(f"Volume: {mesh.volume:.2f} mm³")

# マニフォールド性をチェック
print(f"\n=== 3Dプリント適性チェック ===")
print(f"Is watertight (密閉性): {mesh.is_watertight}")
print(f"Is winding consistent (法線一致性): {mesh.is_winding_consistent}")
print(f"Is valid (幾何的妥当性): {mesh.is_valid}")

# 問題の詳細を診断
if not mesh.is_watertight:
    # 穴（hole）の数を検出
    try:
        edges = mesh.edges_unique
        edges_sorted = mesh.edges_sorted
        duplicate_edges = len(edges_sorted) - len(edges)
        print(f"\n⚠️ 問題検出:")
        print(f"   - メッシュに穴があります")
        print(f"   - 重複エッジ数: {duplicate_edges}")
    except:
        print(f"\n⚠️ メッシュ構造に問題があります")

# 修復を試みる
if not mesh.is_watertight or not mesh.is_winding_consistent:
    print(f"\n🔧 自動修復を実行中...")

    # 法線を修正
    trimesh.repair.fix_normals(mesh)
    print("   ✓ 法線ベクトルを修正")

    # 穴を埋める
    trimesh.repair.fill_holes(mesh)
    print("   ✓ 穴を充填")

    # 縮退三角形を削除
    mesh.remove_degenerate_faces()
    print("   ✓ 縮退面を削除")

    # 重複頂点を結合
    mesh.merge_vertices()
    print("   ✓ 重複頂点を結合")

    # 修復後の状態を確認
    print(f"\n=== 修復後の状態 ===")
    print(f"Is watertight: {mesh.is_watertight}")
    print(f"Is winding consistent: {mesh.is_winding_consistent}")

    # 修復したメッシュを保存
    if mesh.is_watertight:
        mesh.export('model_repaired.stl')
        print(f"\n✅ 修復完了！ model_repaired.stl として保存しました")
    else:
        print(f"\n❌ 自動修復失敗。Meshmixer等の専用ツールを推奨")
else:
    print(f"\n✅ このメッシュは3Dプリント可能です")

# 出力例:
# === メッシュ品質診断 ===
# Vertex count: 1534
# Face count: 2456
# Volume: 12450.75 mm³
#
# === 3Dプリント適性チェック ===
# Is watertight (密閉性): True
# Is winding consistent (法線一致性): True
# Is valid (幾何的妥当性): True
#
# ✅ このメッシュは3Dプリント可能です

学習目標の確認

この章を通じて、以下を説明できるようになったか確認してください。

基本理解

• ✅ 積層造形（AM）の定義とISO/ASTM 52900規格の基本概念を説明できる
• ✅ 7つのAMプロセスカテゴリ（MEX, VPP, PBF, MJ, BJ, SL, DED）の原理と特徴を説明できる
• ✅ STLファイル形式の構造（三角形メッシュ、法線ベクトル、頂点順序）を説明できる
• ✅ スライシングの基本原理（レイヤー高さ、シェル、インフィル、サポート）を説明できる

実践スキル

• ✅ numpy-stlでSTLファイルを読み込み、体積・表面積・バウンディングボックスを計算できる
• ✅ 法線ベクトルの整合性を検証し、反転法線を検出できる
• ✅ trimeshでメッシュの多様体性（watertight）を診断し、自動修復できる
• ✅ 主要なG-codeコマンド（G0/G1/G28/M104など）を読み解ける

応用力

• ✅ 用途要求（精度・強度・コスト・材料）に応じて最適なAMプロセスを選択できる
• ✅ レイヤー高さ・インフィル・サポートを目的に応じて最適化できる
• ✅ STLメッシュの品質を評価し、プリント適性を判断できる

演習問題

Easy（基礎確認）

Medium（応用）

import trimesh

mesh = trimesh.load('model.stl')

# ここにコードを追加: マニフォールド性をチェックし、
# 問題があれば自動修復を行い、修復後のメッシュを
# 'model_fixed.stl'として保存してください

import trimesh

mesh = trimesh.load('model.stl')

# マニフォールド性をチェック
print(f"Is watertight: {mesh.is_watertight}")
print(f"Is winding consistent: {mesh.is_winding_consistent}")

# 問題がある場合は修復
if not mesh.is_watertight or not mesh.is_winding_consistent:
    print("メッシュ修復を実行中...")

    # 法線を修正
    trimesh.repair.fix_normals(mesh)

    # 穴を埋める
    trimesh.repair.fill_holes(mesh)

    # 縮退三角形を削除
    mesh.remove_degenerate_faces()

    # 重複頂点を結合
    mesh.merge_vertices()

    # 修復結果を確認
    print(f"修復後 watertight: {mesh.is_watertight}")

    # 修復したメッシュを保存
    if mesh.is_watertight:
        mesh.export('model_fixed.stl')
        print("修復完了: model_fixed.stl として保存")
    else:
        print("⚠️ 自動修復失敗。Meshmixer等を使用してください")
else:
    print("✓ メッシュは3Dプリント可能です")

次のステップ

第1章では積層造形（AM）の基礎として、ISO/ASTM 52900による7つのプロセス分類、STLファイル形式の構造、スライシングとG-codeの基本を学びました。次の第2章では、材料押出（FDM/FFF）の詳細な造形プロセス、材料特性、プロセスパラメータ最適化について学びます。

参考文献

1. Gibson, I., Rosen, D., & Stucker, B. (2015). Additive Manufacturing Technologies: 3D Printing, Rapid Prototyping, and Direct Digital Manufacturing (2nd ed.). Springer. pp. 1-35, 89-145, 287-334. - AM技術の包括的教科書、7つのプロセスカテゴリとSTLデータ処理の詳細解説
2. ISO/ASTM 52900:2021. Additive manufacturing — General principles — Fundamentals and vocabulary. International Organization for Standardization. - AM用語とプロセス分類の国際標準規格、産業界で広く参照される
3. Kruth, J.P., Leu, M.C., & Nakagawa, T. (1998). "Progress in Additive Manufacturing and Rapid Prototyping." CIRP Annals - Manufacturing Technology, 47(2), 525-540. - 選択的レーザー焼結とバインディング機構の理論的基礎
4. Hull, C.W. (1986). Apparatus for production of three-dimensional objects by stereolithography. US Patent 4,575,330. - 世界初のAM技術（SLA）の特許、AM産業の起源となる重要文献
5. Wohlers, T. (2023). Wohlers Report 2023: 3D Printing and Additive Manufacturing Global State of the Industry. Wohlers Associates, Inc. pp. 15-89, 156-234. - AM市場動向と産業応用の最新統計レポート、年次更新される業界標準資料
6. 3D Systems, Inc. (1988). StereoLithography Interface Specification. - STLファイル形式の公式仕様書、ASCII/Binary STL構造の定義
7. numpy-stl Documentation. (2024). Python library for working with STL files. https://numpy-stl.readthedocs.io/ - STLファイル読込・体積計算のためのPythonライブラリ
8. trimesh Documentation. (2024). Python library for loading and using triangular meshes. https://trimsh.org/ - メッシュ修復・ブーリアン演算・品質評価の包括的ライブラリ

使用ツールとライブラリ

• NumPy (v1.24+): 数値計算ライブラリ - https://numpy.org/
• numpy-stl (v3.0+): STLファイル処理ライブラリ - https://numpy-stl.readthedocs.io/
• trimesh (v4.0+): 3Dメッシュ処理ライブラリ（修復、検証、ブーリアン演算） - https://trimsh.org/
• Matplotlib (v3.7+): データ可視化ライブラリ - https://matplotlib.org/
• SciPy (v1.10+): 科学技術計算ライブラリ（最適化、補間） - https://scipy.org/