delta-theory 6.9.0__py3-none-any.whl

validation/fatigue_redis_api.py

@@ -0,0 +1,334 @@
+#!/usr/bin/env python3
+"""
+FatigueData-AM2022 Upstash Redis API
+=====================================
+
+δ理論検証用 疲労データベースAPI
+
+Usage:
+    from fatigue_redis_api import FatigueDB
+    
+    db = FatigueDB()
+    
+    # 材料一覧
+    materials = db.list_materials()
+    
+    # S-Nデータ取得
+    ti64_sn = db.get_sn('Ti-6Al-4V')
+    
+    # R=-1 のデータのみ
+    ti64_r1 = db.get_sn('Ti-6Al-4V', R=-1.0)
+    
+    # σ_y付きデータのみ
+    ti64_sy = db.get_sn('Ti-6Al-4V', with_sigma_y=True)
+    
+    # da/dN データ
+    dadn = db.get_dadn('Ti-6Al-4V')
+
+Author: 環 & ご主人さま
+Date: 2026-02-02
+"""
+
+import json
+from typing import Dict, List, Optional, Any
+from dataclasses import dataclass
+
+try:
+    from upstash_redis import Redis
+except ImportError:
+    raise ImportError("pip install upstash-redis --break-system-packages")
+
+
+# =============================================================================
+# Configuration
+# =============================================================================
+UPSTASH_URL = "https://casual-crayfish-21486.upstash.io"
+UPSTASH_TOKEN = "AVPuAAIncDFiNWEyN2ExYmEzZGU0N2I5ODlkNmUxMjRkM2UzNWMwY3AxMjE0ODY"
+
+
+# =============================================================================
+# Data Classes
+# =============================================================================
+@dataclass
+class SNPoint:
+    """S-N データポイント"""
+    N: float           # 破断サイクル数
+    S: float           # 応力振幅 [MPa]
+    runout: bool       # ランアウト（未破断）フラグ
+    R: Optional[float] # 応力比
+    sigma_y: Optional[float]   # 降伏応力 [MPa]
+    sigma_uts: Optional[float] # 引張強さ [MPa]
+    doi: Optional[str] # 出典DOI
+
+
+@dataclass  
+class ENPoint:
+    """ε-N データポイント"""
+    N: float           # 破断サイクル数
+    e: float           # ひずみ振幅
+    runout: bool       # ランアウト
+    R: Optional[float]
+    sigma_y: Optional[float]
+    E: Optional[float] # ヤング率 [GPa]
+    doi: Optional[str]
+
+
+@dataclass
+class DaDNPoint:
+    """da/dN-ΔK データポイント"""
+    dK: float          # 応力拡大係数範囲 [MPa√m]
+    dadn: float        # き裂進展速度 [m/cycle]
+    R: Optional[float]
+    doi: Optional[str]
+
+
+# =============================================================================
+# Main API Class
+# =============================================================================
+class FatigueDB:
+    """FatigueData-AM2022 Upstash Redis API"""
+    
+    def __init__(self, url: str = UPSTASH_URL, token: str = UPSTASH_TOKEN):
+        """Initialize connection to Upstash Redis"""
+        self.redis = Redis(url=url, token=token)
+        self._meta = None
+        self._materials = None
+    
+    # -------------------------------------------------------------------------
+    # Meta & Materials
+    # -------------------------------------------------------------------------
+    
+    def get_meta(self) -> Dict[str, Any]:
+        """データベースメタ情報"""
+        if self._meta is None:
+            self._meta = json.loads(self.redis.get('fatigue:meta'))
+        return self._meta
+    
+    def list_materials(self, sort_by: str = 'sn_count') -> List[Dict]:
+        """材料一覧を取得
+        
+        Args:
+            sort_by: ソートキー ('sn_count', 'en_count', 'dadn_count', 'name')
+        
+        Returns:
+            材料リスト [{name, sn_count, en_count, dadn_count, sigma_y_min, sigma_y_max}, ...]
+        """
+        if self._materials is None:
+            self._materials = json.loads(self.redis.get('fatigue:materials'))
+        
+        result = []
+        for name, info in self._materials.items():
+            result.append({'name': name, **info})
+        
+        if sort_by == 'name':
+            result.sort(key=lambda x: x['name'])
+        elif sort_by in ['sn_count', 'en_count', 'dadn_count']:
+            result.sort(key=lambda x: x.get(sort_by, 0), reverse=True)
+        
+        return result
+    
+    def get_material_info(self, material: str) -> Optional[Dict]:
+        """特定材料の情報"""
+        if self._materials is None:
+            self.list_materials()
+        return self._materials.get(material)
+    
+    # -------------------------------------------------------------------------
+    # S-N Data
+    # -------------------------------------------------------------------------
+    
+    def get_sn(
+        self, 
+        material: str,
+        R: Optional[float] = None,
+        with_sigma_y: bool = False,
+        as_dataclass: bool = False
+    ) -> List[Dict | SNPoint]:
+        """S-Nデータを取得
+        
+        Args:
+            material: 材料名 (e.g., 'Ti-6Al-4V', '316L')
+            R: 応力比でフィルタ (e.g., -1.0, 0.1)
+            with_sigma_y: σ_yが存在するデータのみ
+            as_dataclass: SNPointクラスで返す
+        
+        Returns:
+            S-Nデータリスト
+        """
+        key = f'fatigue:sn:{material}'
+        raw = self.redis.get(key)
+        if raw is None:
+            return []
+        
+        data = json.loads(raw)
+        
+        # フィルタ
+        if R is not None:
+            data = [d for d in data if d.get('R') == R]
+        if with_sigma_y:
+            data = [d for d in data if d.get('sigma_y') is not None]
+        
+        if as_dataclass:
+            return [SNPoint(
+                N=d['N'], S=d['S'], runout=bool(d.get('runout', 0)),
+                R=d.get('R'), sigma_y=d.get('sigma_y'), 
+                sigma_uts=d.get('sigma_uts'), doi=d.get('doi')
+            ) for d in data]
+        
+        return data
+    
+    def get_sn_for_delta(
+        self,
+        material: str,
+        R: float = -1.0
+    ) -> List[Dict]:
+        """δ理論検証用にr = σ_a/σ_y を計算済みのデータを取得
+        
+        Args:
+            material: 材料名
+            R: 応力比 (default: -1.0 for fully reversed)
+        
+        Returns:
+            [{N, S, sigma_y, r, runout}, ...]
+        """
+        data = self.get_sn(material, R=R, with_sigma_y=True)
+        
+        result = []
+        for d in data:
+            sigma_y = d['sigma_y']
+            # R=-1 の場合、S = σ_max = 2 * σ_a なので σ_a = S/2... ではない
+            # 実際は S が応力振幅として記録されていることが多い
+            # データセットの定義を確認: S は stress amplitude or max stress?
+            # AM2022では Sは応力振幅（amplitude）として扱う
+            sigma_a = d['S']  # stress amplitude [MPa]
+            r = sigma_a / sigma_y
+            
+            result.append({
+                'N': d['N'],
+                'S': sigma_a,
+                'sigma_y': sigma_y,
+                'r': r,
+                'runout': d.get('runout', 0),
+                'doi': d.get('doi')
+            })
+        
+        return result
+    
+    # -------------------------------------------------------------------------
+    # ε-N Data
+    # -------------------------------------------------------------------------
+    
+    def get_en(
+        self,
+        material: str,
+        R: Optional[float] = None,
+        as_dataclass: bool = False
+    ) -> List[Dict | ENPoint]:
+        """ε-Nデータを取得"""
+        key = f'fatigue:en:{material}'
+        raw = self.redis.get(key)
+        if raw is None:
+            return []
+        
+        data = json.loads(raw)
+        
+        if R is not None:
+            data = [d for d in data if d.get('R') == R]
+        
+        if as_dataclass:
+            return [ENPoint(
+                N=d['N'], e=d['e'], runout=bool(d.get('runout', 0)),
+                R=d.get('R'), sigma_y=d.get('sigma_y'),
+                E=d.get('E'), doi=d.get('doi')
+            ) for d in data]
+        
+        return data
+    
+    # -------------------------------------------------------------------------
+    # da/dN-ΔK Data
+    # -------------------------------------------------------------------------
+    
+    def get_dadn(
+        self,
+        material: str,
+        R: Optional[float] = None,
+        as_dataclass: bool = False
+    ) -> List[Dict | DaDNPoint]:
+        """da/dN-ΔKデータを取得（チャンク対応）"""
+        # チャンクがあるか確認
+        chunks_key = f'fatigue:dadn:{material}:chunks'
+        n_chunks = self.redis.get(chunks_key)
+        
+        if n_chunks:
+            # チャンクから読み込み
+            data = []
+            for i in range(int(n_chunks)):
+                chunk = json.loads(self.redis.get(f'fatigue:dadn:{material}:chunk{i}'))
+                data.extend(chunk)
+        else:
+            # 単一キー
+            key = f'fatigue:dadn:{material}'
+            raw = self.redis.get(key)
+            if raw is None:
+                return []
+            data = json.loads(raw)
+        
+        if R is not None:
+            data = [d for d in data if d.get('R') == R]
+        
+        if as_dataclass:
+            return [DaDNPoint(
+                dK=d['dK'], dadn=d['dadn'],
+                R=d.get('R'), doi=d.get('doi')
+            ) for d in data]
+        
+        return data
+    
+    # -------------------------------------------------------------------------
+    # Utilities
+    # -------------------------------------------------------------------------
+    
+    def search_materials(self, query: str) -> List[str]:
+        """材料名を検索"""
+        if self._materials is None:
+            self.list_materials()
+        
+        query_lower = query.lower()
+        return [name for name in self._materials.keys() 
+                if query_lower in name.lower()]
+    
+    def summary(self) -> str:
+        """データベースサマリー"""
+        meta = self.get_meta()
+        return f"""
+FatigueData-AM2022 @ Upstash Redis
+===================================
+Source: {meta['source']}
+DOI: {meta['doi']}
+License: {meta['license']}
+
+Materials: {meta['total_materials']}
+S-N points: {meta['total_sn_points']:,}
+ε-N points: {meta['total_en_points']:,}
+da/dN points: {meta['total_dadn_points']:,}
+TOTAL: {meta['total_sn_points'] + meta['total_en_points'] + meta['total_dadn_points']:,}
+
+Updated: {meta['updated']}
+"""
+
+
+# =============================================================================
+# CLI
+# =============================================================================
+if __name__ == '__main__':
+    db = FatigueDB()
+    print(db.summary())
+    
+    print("\n📦 Top 5 materials by S-N count:")
+    for mat in db.list_materials()[:5]:
+        print(f"   {mat['name']}: {mat['sn_count']} S-N, σ_y={mat.get('sigma_y_min')}-{mat.get('sigma_y_max')} MPa")
+    
+    print("\n🔍 Ti-6Al-4V S-N (R=-1, with σ_y):")
+    ti64 = db.get_sn_for_delta('Ti-6Al-4V', R=-1.0)
+    print(f"   {len(ti64)} points")
+    print(f"   r range: {min(d['r'] for d in ti64):.4f} - {max(d['r'] for d in ti64):.4f}")