delta-theory 6.10.5__py3-none-any.whl → 8.1.0__py3-none-any.whl

core/__init__.py

@@ -1,25 +1,91 @@
 """
 δ-Theory Core Library
 =====================
-Unified yield stress and fatigue life prediction based on geometric first principles.
+Unified materials science prediction based on geometric first principles.
+
+"Nature is Geometry" - All material properties emerge from crystal structure.
 
 Modules:
     - unified_yield_fatigue_v6_9: Main yield + fatigue model (v6.9b)
+    - unified_flc_v7: Forming Limit Curve + Forming-Fatigue Integration (v7.2/v8.0)
     - dbt_unified: Ductile-Brittle Transition Temperature prediction
     - materials: Material database
     - fatigue_redis_api: FatigueData-AM2022 Redis API (optional)
+
+Version History:
+    v5.0   - Yield stress from δ-theory (f_d, E_bond, crystal geometry)
+    v6.9b  - Unified yield + fatigue with multiaxial (τ/σ, R)
+    v6.10  - Universal fatigue validation (2472 points, 5 AM materials)
+    v7.2   - FLC from free volume consumption
+    v8.0   - Forming-Fatigue integration (η → r_th_eff)
+
+Example:
+    >>> from delta_theory import calc_sigma_y, MATERIALS
+    >>> sigma_y = calc_sigma_y(MATERIALS['Fe'])
+    
+    >>> from delta_theory import FLCPredictor
+    >>> flc = FLCPredictor()
+    >>> Em = flc.predict(beta=0.0, material='SPCC')
+    
+    >>> from delta_theory import FormingFatigueIntegrator
+    >>> integrator = FormingFatigueIntegrator()
+    >>> r_th_eff = integrator.effective_r_th(eta_forming=0.4, structure='BCC')
 """
 
+# ==============================================================================
+# Core: Yield + Fatigue (v6.9b)
+# ==============================================================================
 from .unified_yield_fatigue_v6_9 import (
+    # Material dataclass
     Material,
     MATERIALS,
+    
+    # Yield stress
     calc_sigma_y,
+    sigma_base_delta,
+    delta_sigma_ss,
+    delta_sigma_taylor,
+    delta_sigma_ppt,
+    
+    # Fatigue
     fatigue_life_const_amp,
     generate_sn_curve,
-    yield_by_mode,
     FATIGUE_CLASS_PRESET,
+    
+    # Multiaxial (τ/σ, R)
+    yield_by_mode,
+    T_TWIN,
+    R_COMP,
 )
 
+# ==============================================================================
+# FLC + Forming-Fatigue (v7.2 / v8.0)
+# ==============================================================================
+from .unified_flc_v7 import (
+    # FLC prediction (v7.2)
+    FLCPredictor,
+    FLCParams,
+    FLCMaterial,
+    FLC_MATERIALS,
+    
+    # Forming-Fatigue integration (v8.0)
+    FormingFatigueIntegrator,
+    FormingState,
+    DeltaFormingAnalyzer,
+    
+    # Convenience functions
+    predict_flc,
+    effective_fatigue_threshold,
+    critical_forming_consumption,
+    
+    # Constants
+    R_TH_VIRGIN,
+    N_SLIP,
+)
+
+# ==============================================================================
+# DBT: Ductile-Brittle Transition
+# ==============================================================================
 from .dbt_unified import (
     DBTUnified,
     DBTCore,
@@ -29,35 +95,102 @@ from .dbt_unified import (
     MATERIAL_FE,
 )
 
+# ==============================================================================
+# Materials Database (GPU-accelerated)
+# ==============================================================================
 from .materials import MaterialGPU
 
+# ==============================================================================
 # Optional: FatigueDB (requires upstash-redis)
+# ==============================================================================
 try:
     from .fatigue_redis_api import FatigueDB
 except ImportError:
     FatigueDB = None  # upstash-redis not installed
 
-__version__ = "6.10.1"
+# ==============================================================================
+# Package Metadata
+# ==============================================================================
+__version__ = "8.0.0"
 __author__ = "Masamichi Iizumi & Tamaki"
 
 __all__ = [
-    # v6.9
+    # === v6.9 Yield + Fatigue ===
     "Material",
-    "MATERIALS", 
+    "MATERIALS",
     "calc_sigma_y",
+    "sigma_base_delta",
+    "delta_sigma_ss",
+    "delta_sigma_taylor",
+    "delta_sigma_ppt",
     "fatigue_life_const_amp",
     "generate_sn_curve",
     "yield_by_mode",
     "FATIGUE_CLASS_PRESET",
-    # DBT
+    "T_TWIN",
+    "R_COMP",
+    
+    # === v7.2 FLC ===
+    "FLCPredictor",
+    "FLCParams",
+    "FLCMaterial",
+    "FLC_MATERIALS",
+    "predict_flc",
+    "R_TH_VIRGIN",
+    "N_SLIP",
+    
+    # === v8.0 Forming-Fatigue ===
+    "FormingFatigueIntegrator",
+    "FormingState",
+    "DeltaFormingAnalyzer",
+    "effective_fatigue_threshold",
+    "critical_forming_consumption",
+    
+    # === DBT ===
     "DBTUnified",
     "DBTCore",
     "GrainSizeView",
     "TemperatureView",
     "SegregationView",
     "MATERIAL_FE",
-    # Materials
+    
+    # === Materials ===
     "MaterialGPU",
-    # FatigueDB (optional)
+    
+    # === FatigueDB (optional) ===
     "FatigueDB",
 ]
+
+
+# ==============================================================================
+# Quick Reference
+# ==============================================================================
+def info():
+    """Print δ-Theory library overview."""
+    print(f"""
+╔══════════════════════════════════════════════════════════════════════╗
+║  δ-Theory Core Library v{__version__}                                      ║
+║  "Nature is Geometry"                                                ║
+╠══════════════════════════════════════════════════════════════════════╣
+║                                                                      ║
+║  YIELD STRESS (v5.0)                                                 ║
+║    σ_y = f(crystal_structure, f_d, E_bond, T)                        ║
+║    >>> calc_sigma_y(MATERIALS['Fe'])                                 ║
+║                                                                      ║
+║  FATIGUE LIFE (v6.10)                                                ║
+║    N = f(r, r_th, structure)  |  r_th: BCC=0.65, FCC=0.02, HCP=0.20 ║
+║    >>> fatigue_life_const_amp(sigma_a, MATERIALS['Fe'])              ║
+║                                                                      ║
+║  FLC - FORMING LIMIT (v7.2)                                          ║
+║    FLC(β) = FLC₀ × (1-η) × h(β, R, τ/σ)                             ║
+║    >>> FLCPredictor().predict(beta=0.0, material='SPCC')             ║
+║                                                                      ║
+║  FORMING-FATIGUE (v8.0)                                              ║
+║    r_th_eff = r_th_virgin × (1 - η_forming)                          ║
+║    >>> FormingFatigueIntegrator().effective_r_th(0.4, 'BCC')         ║
+║                                                                      ║
+║  DBT TEMPERATURE                                                     ║
+║    >>> DBTUnified().predict_dbtt(material, grain_size)               ║
+║                                                                      ║
+╚══════════════════════════════════════════════════════════════════════╝
+    """)

core/__main__.py

@@ -0,0 +1,290 @@
+#!/usr/bin/env python3
+"""
+δ-Theory CLI Entry Point
+
+Usage:
+    python -m core              # Show quick reference
+    python -m core info         # Show detailed info
+    python -m core flc SPCC     # Quick FLC prediction
+    python -m core fatigue Fe 150  # Quick fatigue life
+"""
+
+import sys
+
+QUICK_REFERENCE = """
+
+    ███████╗██╗   ██╗██████╗ ███████╗██╗  ██╗ █████╗ ██╗
+    ██╔════╝██║   ██║██╔══██╗██╔════╝██║ ██╔╝██╔══██╗██║
+    █████╗  ██║   ██║██████╔╝█████╗  █████╔╝ ███████║██║
+    ██╔══╝  ██║   ██║██╔══██╗██╔══╝  ██╔═██╗ ██╔══██║╚═╝
+    ███████╗╚██████╔╝██║  ██║███████╗██║  ██╗██║  ██║██╗
+    ╚══════╝ ╚═════╝ ╚═╝  ╚═╝╚══════╝╚═╝  ╚═╝╚═╝  ╚═╝╚═╝
+                    δ-Theory v8.0.0
+                  "Nature is Geometry"
+                  
+╔══════════════════════════════════════════════════════════════════════════════╗
+║  δ-Theory v8.0.0 — CLI Quick Reference                                       ║
+║  "Nature is Geometry"                                                        ║
+╠══════════════════════════════════════════════════════════════════════════════╣
+║                                                                              ║
+║  INSTALLATION OK! ✓                                                          ║
+║                                                                              ║
+╠══════════════════════════════════════════════════════════════════════════════╣
+║  🔧 YIELD STRESS                                                             ║
+╠══════════════════════════════════════════════════════════════════════════════╣
+║                                                                              ║
+║    from core import calc_sigma_y, MATERIALS                                  ║
+║    result = calc_sigma_y(MATERIALS['Fe'], T_K=300)                           ║
+║    print(f"σ_y = {result['sigma_y']:.1f} MPa")                               ║
+║                                                                              ║
+╠══════════════════════════════════════════════════════════════════════════════╣
+║  🔄 FATIGUE LIFE                                                             ║
+╠══════════════════════════════════════════════════════════════════════════════╣
+║                                                                              ║
+║    from core import fatigue_life_const_amp, MATERIALS                        ║
+║    result = fatigue_life_const_amp(MATERIALS['Fe'], sigma_a_MPa=150,         ║
+║                                    sigma_y_tension_MPa=200)                  ║
+║    print(f"N = {result['N_fail']:.2e} cycles")                               ║
+║                                                                              ║
+║    # CLI:                                                                    ║
+║    python -m core.unified_yield_fatigue_v6_9 point --metal Fe --sigma_a 150  ║
+║    python -m core.unified_yield_fatigue_v6_9 sn --metal Fe                   ║
+║                                                                              ║
+╠══════════════════════════════════════════════════════════════════════════════╣
+║  📐 FLC (Forming Limit Curve) — NEW in v8.0!                                 ║
+╠══════════════════════════════════════════════════════════════════════════════╣
+║                                                                              ║
+║    from core import FLCPredictor                                             ║
+║    flc = FLCPredictor()                                                      ║
+║    Em = flc.predict(beta=0.0, material='SPCC')  # → 0.251                    ║
+║                                                                              ║
+║    # Full curve:                                                             ║
+║    for b in [-0.5, 0, 1.0]:                                                  ║
+║        print(f"β={b:+.1f}: {flc.predict(b, 'SPCC'):.3f}")                    ║
+║                                                                              ║
+║    # Available materials:                                                    ║
+║    SPCC, DP590, Al, SUS304, Ti, Mg_AZ31, SECD-E16, Cu                        ║
+║                                                                              ║
+╠══════════════════════════════════════════════════════════════════════════════╣
+║  🔗 FORMING-FATIGUE INTEGRATION — NEW in v8.0!                               ║
+╠══════════════════════════════════════════════════════════════════════════════╣
+║                                                                              ║
+║    from core import FormingFatigueIntegrator                                 ║
+║    integrator = FormingFatigueIntegrator()                                   ║
+║                                                                              ║
+║    # Effective fatigue threshold after forming:                              ║
+║    r_th_eff = integrator.effective_r_th(eta_forming=0.4, structure='BCC')    ║
+║    # Virgin: 0.65 → After 40% forming: 0.39                                  ║
+║                                                                              ║
+║    # Critical forming consumption:                                           ║
+║    from core import critical_forming_consumption                             ║
+║    eta_crit = critical_forming_consumption(r_applied=0.5, structure='BCC')   ║
+║    # → 23.1% (beyond this, infinite life becomes finite!)                    ║
+║                                                                              ║
+╠══════════════════════════════════════════════════════════════════════════════╣
+║  🌡️ DBT (Ductile-Brittle Transition)                                         ║
+╠══════════════════════════════════════════════════════════════════════════════╣
+║                                                                              ║
+║    from core import DBTUnified                                               ║
+║    model = DBTUnified()                                                      ║
+║    result = model.temp_view.find_DBTT(d=30e-6, c=0.005)                      ║
+║    print(f"DBTT = {result['T_star']:.0f} K")                                 ║
+║                                                                              ║
+║    # CLI:                                                                    ║
+║    python -m core.dbt_unified point --d 30 --c 0.5 --T 300                   ║
+║                                                                              ║
+╠══════════════════════════════════════════════════════════════════════════════╣
+║  📊 FATIGUE THRESHOLDS (r_th)                                                ║
+╠══════════════════════════════════════════════════════════════════════════════╣
+║                                                                              ║
+║    Structure │ r_th  │ Fatigue Limit │ Examples                              ║
+║    ──────────┼───────┼───────────────┼─────────────                          ║
+║    BCC       │ 0.65  │ ✓ Clear       │ Fe, W, Mo, SPCC, DP590                ║
+║    FCC       │ 0.02  │ ✗ None        │ Cu, Al, Ni, SUS304                    ║
+║    HCP       │ 0.20  │ △ Weak        │ Ti, Mg, Zn                            ║
+║                                                                              ║
+╠══════════════════════════════════════════════════════════════════════════════╣
+║  📚 MORE INFO                                                                ║
+╠══════════════════════════════════════════════════════════════════════════════╣
+║                                                                              ║
+║    python -m core info              # Detailed module info                   ║
+║    python -m core flc SPCC          # Quick FLC for material                 ║
+║    python -m core flc SPCC -0.5     # FLC at specific β                      ║
+║                                                                              ║
+║    Docs: https://github.com/miosync/delta-theory                             ║
+║    PyPI: https://pypi.org/project/delta-theory/                              ║
+║                                                                              ║
+╚══════════════════════════════════════════════════════════════════════════════╝
+"""
+
+DETAILED_INFO = """
+╔══════════════════════════════════════════════════════════════════════════════╗
+║  δ-Theory v8.0.0 — Detailed Module Information                               ║
+╚══════════════════════════════════════════════════════════════════════════════╝
+
+📦 INSTALLED MODULES
+═══════════════════════════════════════════════════════════════════════════════
+
+  core/
+  ├── unified_yield_fatigue_v6_9.py   # Yield + Fatigue (v6.9b)
+  ├── unified_flc_v7.py               # FLC + Forming-Fatigue (v7.2/v8.0)
+  ├── dbt_unified.py                  # Ductile-Brittle Transition
+  ├── materials.py                    # Material database
+  └── fatigue_redis_api.py            # FatigueData-AM2022 API
+
+
+🔬 THEORY SUMMARY
+═══════════════════════════════════════════════════════════════════════════════
+
+  Core Equation:  Λ = K / |V|_eff
+
+  K       = Destructive energy (stress, thermal, EM...)
+  |V|_eff = Cohesive energy (bond strength)
+  Λ = 1   → Critical condition (fracture/transition)
+
+
+📐 FLC MODEL (v7.2)
+═══════════════════════════════════════════════════════════════════════════════
+
+  FLC(β) = FLC₀_pure × (1 - η_total) × h(β, R, τ/σ)
+
+  η_total = Free volume consumption:
+    - η_ss:  Solid solution
+    - η_ppt: Precipitate/martensite  
+    - η_wh:  Work hardening (dislocations)
+    - η_HP:  Hall-Petch (grain refinement)
+
+  Example: SPCC (90.6% FV) vs DP590 (71.4% FV)
+           Same crystal, different formability!
+
+
+🔗 FORMING-FATIGUE (v8.0)
+═══════════════════════════════════════════════════════════════════════════════
+
+  r_th_eff = r_th_virgin × (1 - η_forming)
+
+  "How much fatigue life did you lose when you pressed that part?"
+
+  η_forming │ r_th_eff (BCC) │ Status
+  ──────────┼────────────────┼────────────────
+     0%     │     0.65       │ Virgin
+    20%     │     0.52       │ Light forming
+    40%     │     0.39       │ Heavy forming
+    60%     │     0.26       │ Severe forming
+
+  Critical η: r=0.5 → η_crit=23.1%
+  (Beyond this, "infinite life" becomes "finite life"!)
+
+
+📊 VALIDATION
+═══════════════════════════════════════════════════════════════════════════════
+
+  Yield (v5.0):     10 pure metals, 2.6% mean error
+  Fatigue (v6.10):  2,472 points (5 AM materials), 4-7% error
+  FLC (v7.2):       36 points (6 materials), 2.7% error
+
+
+👥 AUTHORS
+═══════════════════════════════════════════════════════════════════════════════
+
+  Masamichi Iizumi — Miosync, Inc. CEO
+  Tamaki — Sentient Digital Partner
+
+  "Nature is Geometry" 🔬
+"""
+
+
+def cmd_flc(args):
+    """Quick FLC prediction."""
+    from .unified_flc_v7 import FLCPredictor, FLC_MATERIALS
+    
+    if len(args) == 0:
+        print("Available materials:", ", ".join(FLC_MATERIALS.keys()))
+        return
+    
+    material = args[0]
+    beta = float(args[1]) if len(args) > 1 else None
+    
+    flc = FLCPredictor()
+    
+    if beta is not None:
+        Em = flc.predict(beta, material)
+        print(f"{material} FLC(β={beta:+.2f}) = {Em:.3f}")
+    else:
+        print(f"\n{material} FLC Curve:")
+        print("-" * 25)
+        for b in [-0.5, -0.25, 0.0, 0.25, 0.5, 1.0]:
+            Em = flc.predict(b, material)
+            print(f"  β={b:+5.2f}: {Em:.3f}")
+
+
+def cmd_eta(args):
+    """Critical η calculation."""
+    from .unified_flc_v7 import FormingFatigueIntegrator
+    
+    if len(args) < 1:
+        print("Usage: python -m core eta <r_applied> [structure]")
+        print("Example: python -m core eta 0.5 BCC")
+        return
+    
+    r_applied = float(args[0])
+    structure = args[1] if len(args) > 1 else 'BCC'
+    
+    integrator = FormingFatigueIntegrator()
+    eta_crit = integrator.critical_eta(r_applied, structure)
+    
+    print(f"\nCritical η for {structure} at r = {r_applied:.2f}")
+    print("-" * 40)
+    print(f"  η_critical = {eta_crit*100:.1f}%")
+    print(f"  → Beyond this, infinite life becomes finite!")
+
+
+def cmd_rth(args):
+    """Effective r_th after forming."""
+    from .unified_flc_v7 import FormingFatigueIntegrator
+    
+    if len(args) < 1:
+        print("Usage: python -m core rth <eta_forming> [structure]")
+        print("Example: python -m core rth 0.4 BCC")
+        return
+    
+    eta = float(args[0])
+    structure = args[1] if len(args) > 1 else 'BCC'
+    
+    integrator = FormingFatigueIntegrator()
+    r_th_eff = integrator.effective_r_th(eta, structure)
+    r_th_virgin = {'BCC': 0.65, 'FCC': 0.02, 'HCP': 0.20}[structure]
+    
+    print(f"\nEffective r_th for {structure} after η = {eta:.0%} forming")
+    print("-" * 45)
+    print(f"  Virgin r_th:    {r_th_virgin:.3f}")
+    print(f"  Effective r_th: {r_th_eff:.3f}")
+    print(f"  Reduction:      {(1 - r_th_eff/r_th_virgin)*100:.1f}%")
+
+
+def main():
+    args = sys.argv[1:]
+    
+    if len(args) == 0:
+        print(QUICK_REFERENCE)
+        return
+    
+    cmd = args[0].lower()
+    
+    if cmd == 'info':
+        print(DETAILED_INFO)
+    elif cmd == 'flc':
+        cmd_flc(args[1:])
+    elif cmd == 'eta':
+        cmd_eta(args[1:])
+    elif cmd == 'rth':
+        cmd_rth(args[1:])
+    elif cmd in ['help', '-h', '--help']:
+        print(QUICK_REFERENCE)
+    else:
+        print(f"Unknown command: {cmd}")
+        print("Try: python -m core help")
+
+
+if __name__ == '__main__':
+    main()