ffseval 0.3.4__tar.gz

ffseval-0.3.4/FFSeval/FFS.py

@@ -0,0 +1,1765 @@
+import numpy as np
+from scipy.optimize import fsolve
+import matplotlib.pyplot as plt
+class Base:
+    def __init__(self):
+        self.title=None
+        self.data=None
+        self.res=None
+        self.refer=None
+    def SetTitle(self,title:str)->None:
+        self.title=title
+    def Title(self)->str:
+        return self.title
+    def SetData(self,data:dict):
+        self.data=data
+    def GetData(self):
+        return self.data
+    def SetRes(self,res:dict):
+        self.res=res
+    def GetRes(self):
+        return self.res
+    def SetRefer(self,refer:str):
+        self.refer=refer
+    def GetRefer(self):
+        return self.refer
+    def CalcKc(self):
+        nu=self.data['Nu']
+        E=self.data['E']
+        j1c=self.data['J1c']
+        Ed=E/(1-nu*nu)
+        return np.sqrt(Ed*j1c)
+    def Option1(self,Kr,Lr):
+        kk=(1-0.14*Lr*Lr)*(0.3+0.7*np.exp(-0.65*Lr**6))
+        flag=False
+        if Kr<kk:
+            flag=True
+        return flag,kk
+    def DrawOption1(self,Lr0,Kr0,Su,Sy):
+        """
+        R6法-Rev.3のOption1破壊評価曲線の描画
+        Lr0,Kr0描画点
+        Su:引張強さ
+        Sy:降伏強さ
+        """
+        compute_Kr = lambda Lr: (1 - 0.14 * Lr**2) * (0.3 + 0.7 * np.exp(-0.65 * Lr**6))
+        Sf=(Sy+Su)/2
+        Lrmax=Sf/Sy
+        Krmax=compute_Kr(Lrmax)
+        # Lrの範囲を生成
+        Lr_values = np.linspace(0, Lrmax, 500)
+        Kr_values = compute_Kr(Lr_values)
+
+        # グラフ描画
+        plt.figure(figsize=(8, 5))
+        plt.plot(Lr_values, Kr_values,color='blue')
+        plt.plot(Lr0, Kr0, 'ro')  # 赤い点をプロット
+        plt.plot([Lrmax,Lrmax],[0,Krmax],color='blue')
+        plt.ylim(0,1.2)
+        plt.xlabel('Lr',fontsize=16)
+        plt.ylabel('Kr',fontsize=16)
+        plt.grid(True)
+        plt.legend()
+        plt.tight_layout()
+        plt.show()        
+    def Margin(self,K:float,L:float)->dict:
+        '''
+        評価点(L,K)について，Option1曲線に対する安全裕度を評価する
+        計算する値
+        (L0,K0): 原点と評価点の延長線とOption1曲線との交点
+        margin: 安全裕度　1以下であれば安全
+        '''
+
+        # Krの定義
+        Kr1 = lambda Lr: (1 - 0.14 * Lr**2) * (0.3 + 0.7 * np.exp(-0.65 * Lr**6))
+        Kr2 = lambda Lr:(K / L) * Lr
+        equation=lambda Lr:Kr1(Lr) - Kr2(Lr)
+        # 初期値の推測（範囲によって複数回試すと良い）
+        initial_guess = 0.5
+        solution = fsolve(equation, initial_guess)
+        res={}
+        res['L0']=solution[0]
+        res['K0']=Kr1(solution[0])
+        res['margin']=K/res['K0']
+        return res
+    def RolfeBarsom(self,Cv,Sy)->float:
+        #Cv=self.data['Cv']
+        #Sy=self.data['Sy']
+        cc=np.array(Cv)
+        c=0.6478*(cc/Sy-0.0098)
+        K1c=np.sqrt(c)*Sy
+        return K1c
+    def JR(self,C:float,m:float,da:float)->float:
+        return C*da**m
+class Fatigue:
+    '''
+    JSME維持規格における炭素鋼および低合金鋼の大気中における疲労亀裂進展特性
+    '''
+    def __init__(self,cls,data,pfm=False,cov=0.1):
+        self.cls=cls
+        self.cls.SetData(data)
+        self.data=data
+        self.pfm=pfm #PFM計算のときTrue
+        self.cov=cov #PFM計算のとき，係数Cのcov
+    def dadN(self,a,c,Pmin,Pmax):
+        self.data['a']=a
+        self.data['c']=c
+        self.data['P']=Pmin
+        self.cls.SetData(self.data)
+        self.cls.Calc()
+        resMin=self.cls.GetRes()
+        self.data['P']=Pmax
+        self.cls.SetData(self.data)
+        self.cls.Calc()
+        resMax=self.cls.GetRes()
+        dKA=resMax['KA']-resMin['KA']
+        dKB=resMax['KB']-resMin['KB']
+        da=self.FatigueSteel(dKA)
+        dc=self.FatigueSteel(dKB)
+        return da,dc,resMax
+    def FatigueSteel(self,dK):
+        n=3.07
+        da=self.C*dK**n
+        return da
+    def EvalAC(self,a0,c0,Pmin,Pmax,R,n):
+        S=25.72*(2.88-R)**(-3.07)
+        C=3.88e-9*S
+        if self.pfm:#PFM計算のときには，正規乱数を発生してCに割り当てる
+            mean=C
+            std_dev=C*self.cov
+            C=np.random.normal(mean,std_dev)
+        self.C=C
+        self.data['a']=a0
+        self.data['c']=c0
+        self.data['P']=Pmax
+        self.cls.SetData(self.data)
+        self.cls.Calc()
+        res0=self.cls.GetRes()
+        a=a0
+        c=c0
+        for i in range(n):
+            da,dc,resMax=self.dadN(a,c,Pmin,Pmax)
+            a += da/1000
+            c += dc/1000
+        crack={'a':a,
+               'c':c}
+        res1=resMax
+        return res0,res1,crack
+        
+        
+class Treat:
+    def Set(self,spec:str):
+        '''
+        対象とする解析記号名を文字列でセット
+        '''
+        spec2=spec.replace("-","_")
+        df=self.Registered()
+        dd=df[spec2[0]]
+        if spec2 not in dd:
+            print(spec+' is not registered yet!')
+            return
+        cls=globals()[spec2]
+        instance=cls()
+        return instance
+    def Registered(self):
+        df={'J':[
+            'J_2_k_a',
+            'J_2_k_b',
+            'J_7_a'
+            ],
+            'K':[
+            'K_1_a_1',
+            'K_1_a_2',
+            'K_1_a_3',
+            'K_1_b_1',
+            'K_1_c_1',
+            'K_1_d_1',
+            'K_1_d_2',
+            'K_1_e_1',
+            'K_1_e_2',
+            'K_1_e_3',
+            'K_2_a_1',
+            'K_2_a_2',
+            'K_2_a_3',
+            'K_2_b_1',
+            'K_2_b_2',
+            'K_2_c_1',
+            'K_2_c_2',
+            'K_2_d',
+            'K_2_e_1',
+            'K_2_e_2',
+            'K_2_e_3',
+            'K_2_f_1',
+            'K_2_f_2',
+            'K_2_f_3',
+            'K_2_g',
+            'K_2_k_2'
+        ],'L':[
+            'L_1_a',
+            'L_2_b'
+        ]
+        }
+        return df
+import Kriging as kr
+import csv
+import matplotlib.pyplot as plt
+from mpl_toolkits.mplot3d import Axes3D
+from typing import Tuple
+import importlib.resources as pkg_resources
+from importlib import resources
+from FFSeval import data
+class dmanage:
+    '''
+    係数が表形式で与えられるときに，表データをcsvファイルから読み取り，
+    Kriging法で内挿した上で，評価値を返す
+    [使用法]二次元の場合*****************************
+    dm=dmanage()
+    data=dm.Finput('J-2-k-b.csv')
+    X,W=dm.CSV2XW(data,1)
+    r2_score=dm.KrigCalc(X,W)
+    dm.DrawRes(X,W) #鳥観図による確認
+    x_val=0.5; y_val=2.0 #評価したい点
+    target_point=np.array([[x_val,y_val]])
+    w_pred,sigma=dm.Eval(target_point)# w_predが予測値
+    [使用法]三次元の場合*****************************
+    dm=dmanage()
+    data=dm.Finput('J-2-k-b.csv')
+    z=[0.0625,0.125,0.25,0.37,0.50]#三次元目のデータリスト(表の数だけ存在)
+    X,W=dm.Fconv3D(data,z)
+    r2_score=dm.KrigCalc(X,W)
+    target_point=np.array([[0.5,2,0.18]])
+    w_pred,sigma=dm.Eval(target_point)# w_predが予測値
+    '''
+    def __init__(self):
+        self.df=None
+    def Fconv3D(self,data:list,z:list)->Tuple[np.ndarray, np.ndarray]:
+        '''
+        3次元の入力テーブルに対する対応
+        dataの中から，全データの処理を行いX,Wを構成して戻す
+        三次元目の情報はz:listで与える
+        '''
+        # STARTを含む要素の数をカウント
+        count = 0
+        for sublist in data:
+            for item in sublist:
+                if "START" in item:
+                    count += 1
+        nz=len(z)
+        if count!=nz:
+            print(f'STARTの数{count}とzのサイズ{nz}が一致していない')
+            return
+        W=np.array([])
+        X=np.empty((0,3))
+        for i in range(count):
+            ith=i+1
+            df=self.dRead(data,ith)
+            xval=df['xval']
+            yval=df['yval']
+            xn=len(xval)
+            yn=len(yval)
+            for iy in range(yn):
+                for ix in range(xn):
+                    d=[yval[iy],xval[ix],z[i]]
+                    X=np.append(X,[d], axis=0)
+                    W=np.append(W,df['arr'][iy,ix])
+        return X,W
+    def Fconv4D(self,data:list,z1:list,z2:list)->Tuple[np.ndarray, np.ndarray]:
+        '''
+        4次元の入力テーブルに対する対応
+        dataの中から，全データの処理を行いX,Wを構成して戻す
+        三次元目，四次元目の情報はz1:list,z2:listで与える
+        '''
+        # STARTを含む要素の数をカウント
+        count = 0
+        for sublist in data:
+            for item in sublist:
+                if "START" in item:
+                    count += 1
+        nz1=len(z1)
+        nz2=len(z2)
+        if count!=nz1 or count!=nz2:
+            print(f'STARTの数{count}とzのサイズ{nz1,nz2}が一致していない')
+            return
+        W=np.array([])
+        X=np.empty((0,4))
+        for i in range(count):
+            ith=i+1
+            df=self.dRead(data,ith)
+            xval=df['xval']
+            yval=df['yval']
+            xn=len(xval)
+            yn=len(yval)
+            for iy in range(yn):
+                for ix in range(xn):
+                    d=[yval[iy],xval[ix],z1[i],z2[i]]
+                    X=np.append(X,[d], axis=0)
+                    W=np.append(W,df['arr'][iy,ix])
+        return X,W            
+    def Finput(self,fname:str)->list:
+        '''
+        csvファイルを読み取りリストに格納する
+        '''
+#        with resources.files("FFSeval.data").joinpath("J-2-k-a-2.csv").open("r", encoding="utf-8", newline='') as csvfile:
+        with resources.files("FFSeval.data").joinpath(fname).open("r", encoding="utf-8", newline='') as csvfile:
+
+            reader = csv.reader(csvfile)
+
+            data=[]
+            for row in reader:
+                data.append(row)
+        return data        
+    def dRead(self,data:list,ith:int)->dict:
+        '''
+        dataのith番目のテーブルを辞書として返す
+        '''
+        n=len(data)
+        ii=0
+        flag=False
+        res=[]
+        l=0
+        for ll in data:
+            if ll[0]=='START':
+                ii+=1
+                if ii==ith:
+                    flag=True
+                    l+=1
+                    continue
+            if flag and l!=0:
+                if ll[0]=='END':
+                    break
+                if l==1:
+                    nx=int(ll[0])
+                    ny=int(ll[1])
+                    l+=1
+                    continue
+                if l==2:
+                    numlist=[float(x) for x in ll[:nx]]
+                    xval=numlist
+                    l+=1
+                    continue
+                if l==3:
+                    numlist=[float(x) for x in ll[:ny]]
+                    yval=numlist
+                    l+=1
+                    continue
+                numlist=[float(x) for x in ll[:nx]]
+                res.append(numlist)
+                l+=1
+        arr=np.array(res)
+        df={}
+        df['xval']=xval
+        df['yval']=yval
+        df['arr']=arr
+        return df
+    def MakeInp(self,df:dict)->Tuple[np.ndarray, np.ndarray]:
+        '''
+        辞書型データを，Kriging入力用のnp.arrayに変換する
+        '''
+        xval=df['xval']
+        yval=df['yval']
+        xn=len(xval)
+        yn=len(yval)
+        W=np.array([])
+        X=np.empty((0,2))
+        for iy in range(yn):
+            for ix in range(xn):
+                d=[yval[iy],xval[ix]]
+                X=np.append(X,[d], axis=0)
+                W=np.append(W,df['arr'][iy,ix])
+        return X,W
+    def CSV2XW(self,data:list,ith:int)->Tuple[np.ndarray, np.ndarray]:
+        df=self.dRead(data,ith)
+        X,W=self.MakeInp(df)
+        self.df=df
+        return X,W
+    def GetDf(self):
+        return self.df
+    def KrigCalc(self,X,W,alpha=5e-4):
+        self.krig=kr.Kriging()
+        self.krig.setData(X,W)
+        r2_score=self.krig.Fit(alpha=alpha)
+        return r2_score
+    def Eval(self,target_point:np.array)->float:
+        #target_point = np.array([[x, y]])  
+        w_pred,sigma=self.krig.Predict(target_point)
+        return w_pred[0],sigma        
+    def DrawRes(self,X,W)->None:
+        # 1. 描画用のメッシュグリッドを作成
+        x = np.linspace(np.min(X[:,0]), np.max(X[:,0]), 50)
+        y = np.linspace(np.min(X[:,1]), np.max(X[:,1]), 50)
+        X_grid, Y_grid = np.meshgrid(x, y)
+
+        # 2. メッシュ座標を一次元にして予測
+        XY = np.vstack([X_grid.ravel(), Y_grid.ravel()]).T
+
+        # Kriging モデルで予測（タプルで返る）
+        Z_pred, _ = self.krig.Predict(XY)
+
+        # 予測値をグリッド形状に整形
+        Z_grid = Z_pred.reshape(X_grid.shape)
+
+        # 4. 描画
+        fig = plt.figure(figsize=(10, 7))
+        ax = fig.add_subplot(111, projection='3d')
+
+        # 予測面
+        ax.plot_surface(X_grid, Y_grid, Z_grid, cmap='viridis', alpha=0.7)
+
+        # 元データ点も重ねる
+        #ax.scatter(X[:,0], X[:,1], W, color='r', s=30, label='Data points')
+        # 予測に使ったデータ点を赤い球で表示
+        ax.scatter(X[:, 0], X[:, 1], W, color='black', s=30, marker='^', label='Training data')
+        # 軸ラベル
+        ax.set_xlabel('X')
+        ax.set_ylabel('Y')
+        ax.set_zlabel('W')
+
+        # タイトル・凡例
+        ax.set_title('Kriging Model Surface')
+        ax.legend()
+        plt.show()               
+class J_2_k_a(Base):
+    #by S.Sakai
+    def __init__(self):
+        super().SetTitle('周方向貫通亀裂 Zahoorの解')
+        super().SetRefer('Zahoor, A.:Ductile Fracture Handbook Volume 1,EPRI NP-6301-D,1989')
+        # Kriging法による応答曲面を計算
+        self.dm_P=dmanage()
+        data=self.dm_P.Finput('J-2-k-a-1.csv')
+        z=[5.0,10.0,20.0]
+        X,W=self.dm_P.Fconv3D(data,z)
+        r2_score=self.dm_P.KrigCalc(X,W)
+        self.dm_M=dmanage()
+        data=self.dm_M.Finput('J-2-k-a-2.csv')
+        z=[5.0,10.0,20.0]
+        X,W=self.dm_M.Fconv3D(data,z)
+        r2_score=self.dm_M.KrigCalc(X,W)                  
+    def Calc(self):
+        df=super().GetData()
+        th=df['th']
+        plane=df['plane']
+        if plane=='stress': beta=2
+        if plane=='strain': beta=6
+        M=df['M']
+        R=df['R']
+        t=df['t']
+        P=df['P']
+        S0=df['S0']
+        alpha=df['alpha']
+        e0=df['e0']
+        n=df['n']
+        E=df['E']
+        A=0.0
+        P0=2.0*S0*R*t*(np.pi-th-2.0*np.arcsin(0.5*np.sin(th))) 
+        M0=4.0*S0*R*R*t*(np.cos(th/2.0)-0.5*np.sin(th))
+        if df['Case']=='Collapse': #塑性崩壊値の計算
+            res={
+                'P0':P0,
+                'M0':M0
+            }
+            super().SetRes(res)
+            return
+        if R/t >= 5.0 and R/t<10.0:
+            A=(0.125*(R/t)-0.25)**0.25
+        if R/t>=10.0 and R/t<=20.0:
+            A=(0.4*R/t-3.0)**0.25
+        if plane=='stress': beta=2
+        if plane=='strain': beta=6
+        if df['Case']=='PR': #塑性崩壊強度の計算
+            pass #将来開発すること
+        if df['Case']=='PJ':
+            target_point=np.array([[th/np.pi,n,R/t]])
+            H1,sigma=self.dm_P.Eval(target_point)           
+            Ft=1.0+A*(5.3303*(th/np.pi)**1.5+18.773*(th/np.pi)**4.24)
+            St=P/(2.0*np.pi*R*t)
+            the=th*(1.0+(Ft*Ft/beta)*(n-1)/(n+1)*(St/S0)**2/(1+(P/P0)**2))
+            ft=(the/np.pi)*(1.0+A*(5.3303*(the/np.pi)**1.5+18.773*(the/np.pi)**4.24))**2
+            J=ft*P*P/(4.0*np.pi*R*t*t*E)+alpha*S0*e0*(np.pi-th)*H1*(P/P0)**(n+1)
+            res={'J':J}
+            super().SetRes(res)
+            return
+        if df['Case']=='MR': #塑性崩壊強度の計算
+            target_point=np.array([[th/np.pi,n,R/t]])
+            H1,sigma=self.dm_M.Eval(target_point)
+            JR=df['JR']
+            MR=M0*(JR/(alpha*S0*e0*np.pi*R*(1.0-th/np.pi)**2*H1))**(1./(n+1.))
+            res={'MR':MR,
+                    'M0':M0,
+                    'H1':H1}
+            super().SetRes(res)
+            return
+        if df['Case']=='MJ':
+            target_point=np.array([[th/np.pi,n,R/t]])
+            H1,sigma=self.dm_M.Eval(target_point)
+            Sb=M/(np.pi*R*R*t)
+
+            Fb=1.0+A*(4.5967*(th/np.pi)**1.5+2.6422*(th/np.pi)**4.24)
+            the=th*(1.0+Fb*Fb/beta*(n-1)/(n+1)*(Sb/S0)**2/(1+(M/M0)**2))
+            fb=(the/np.pi)*(1.0+A*(4.5967*(the/np.pi)**1.5+2.6422*(the/np.pi)**4.24))**2
+            J=fb*M*M/(R*R*R*t*t*E)+alpha*S0*e0*np.pi*R*(1-th/np.pi)**2*H1*(M/M0)**(n+1) 
+            res={'J':J,
+                'M0':M0,
+                'H1':H1}
+            super().SetRes(res)       
+class J_2_k_b(Base):
+    #by S.Sakai
+    def __init__(self):
+        super().SetTitle('周方向貫通亀裂　Zahoorの解')
+        super().SetRefer('Zahoor, A.:Ductile Fracture Handbook Volume 1,EPRI NP-6301-D,1989')
+
+             
+             
+    def Calc(self):
+        df=super().GetData()
+        th=df['Th']/180.0*np.pi
+        plane=df['plane']
+        if plane=='stress': beta=2
+        if plane=='strain': beta=6
+        M=df['M']
+        R=df['R']
+        t=df['t']
+        P=df['P']
+        St=P/(2.*np.pi*R*t)
+        Sb=M/(np.pi*R*R*t)
+        S0=df['S0']
+        P0=2*S0*R*t*(np.pi-th-2*np.arcsin(0.5*np.sin(th)))
+        M0=4*S0*R*R*t*(np.cos(th/2)-0.5*np.sin(th))
+        lam=M/P/R
+        P0d=0.5*(-lam*R*P0*P0/M0+np.sqrt((lam*R*P0*P0/M0)**2+4*P0*P0))
+        n=df['n']
+        dm=dmanage()
+        x=lam/(1.+lam)
+        data=dm.Finput('J-2-k-b.csv')
+        z=[0.0625,0.125,0.25,0.37,0.50]#三次元目のデータリスト(表の数だけ存在)
+        X,W=dm.Fconv3D(data,z)
+        r2_score=dm.KrigCalc(X,W)
+        target_point=np.array([[x,n,th/np.pi]])
+        h1,sigma=dm.Eval(target_point)
+        E=df['E']
+        alpha=df['alpha']
+        e0=df['e0']
+        the=th*(1+(1/beta)*((n-1)/(n+1))*((St*Ft+Sb*Fb)**2/S0/S0)/(1+(P/P0d)**2))
+        ft=(the/np.pi)*(1+A*(5.3303*(the/np.pi)**1.5+18.773*(the/np.pi)**4.24))**2
+        fb=(the/np.pi)*(1+A*(4.5967*(the/np.pi)**1.5+2.6422*(the/np.pi)**4.24))**2
+        J=ft*P*P/(4*R*t*t*E)+fb*M*M/(R*R*R*t*t*E)+alpha*S0*e0*R*(np.pi-th)*(th/np.pi)*h1*(P/P0d)**(n+1)
+        res={'J':J}
+        super().SetRes(res)
+class J_7_a(Base):
+    #by S.Sakai
+    def __init__(self):
+        super().SetTitle('円孔縁のき裂 片側貫通亀裂 Zahoorの解')
+        super().SetRefer('Zahoor, A.:Ductile Fracture Handbook Volume 3,EPRI NP-6301-D,1991')  
+    def Calc(self):
+        dm=dmanage()
+        df=super().GetData()
+
+        a=df['a']
+        R=df['R']
+        alpha=df['alpha']
+        n=df['n']
+        S=df['sigma'] # σ 
+        S0=df['sigma0'] # σ0
+        E = df['E']
+        e0=S0/E # ε0
+
+        data=dm.Finput('J-7-a.csv')
+        ith=1
+        X,W=dm.CSV2XW(data,ith)
+        r2_score=dm.KrigCalc(X,W)
+        target_point=np.array([[a/R,n]]) 
+        H1,sigma=dm.Eval(target_point) 
+
+        a_over_R=a/R #代入してみましたが、代入するのとしないのとどちらが良いですか？  
+        B0=(4.0/3.0)/(1-0.08696*(1+0.5*a_over_R))**2
+        F=(2.8041-4.9327*a_over_R+7.986*a_over_R**2-6.9783*a_over_R**3+2.4132*a_over_R**4)
+        ae=a*(1+0.5*F**2*((n-1)/(n+1))*((S/S0)**2/(1+B0*(S/S0)**2)))
+        ae_over_R=ae/R #代入してみましたが、代入するのとしないのとどちらが良いですか？
+        f=np.pi*ae_over_R*(2.8041-4.9327*ae_over_R+7.986*ae_over_R**2-6.9783*ae_over_R**3+2.4132*ae_over_R**4)**2
+        #J=f*R*S**2/E+alpha*S0*e0*H1*(S/S0)**(n+1)
+        J=f*R*S**2/E+alpha*S0*e0*R*H1*(S/S0)**(n+1) #R*が抜けていなした
+
+        #res={'J':J,
+            #'H1':H1}
+        res={'J':J}#H1は戻す必要はありません
+        super().SetRes(res) 
+
+
+class K_1_a_1(Base):
+    #by S.Sakai
+    def __init__(self):
+        super().SetTitle('平板の半楕円表面亀裂，Raju-Newmanの解')
+        super().SetRefer('Newman,J.C>Jr., and Raju,I.S.:Stress-Intensity Factor Equations for Cracks in Three-Dimentional Finite Bodies Subjected to Tension and Bending Loads, NASA Technical Memorandum, 85793, NASA,1984')
+    def Calc(self):
+        df=super().GetData()
+        a=df['a']
+        c=df['c']
+        b=df['b']
+        t=df['t']
+        P=df['P']
+        M=df['M']
+        Sm=P/(2*b*t)
+        Sb=3*M/(b*t*t)
+        if a/c <=1.0:
+            Q=1+1.464*(a/c)**1.65
+            g=1
+            fphai=1
+            fw=np.sqrt(1.0/np.cos(np.pi*c/(2*b)*np.sqrt(a/t)))
+            H=1+(-1.22-0.12*a/c)*a/t+(0.55-1.05*(a/c)**0.75+0.47*(a/c)**1.5)*(a/t)**2
+            FA=(1.13-0.09*a/c+(-0.54+0.89/(0.2+a/c))*(a/t)**2+(0.5-1/(0.65+a/c)+14*(1-a/c)**24)*(a/t)**4)*g*fphai*fw
+        else:
+            Q=1+1.464*(c/a)**1.65
+            g=1
+            fphai=np.sqrt(c/a)
+            fw=np.sqrt(1.0/np.cos(np.pi*c/(2*b)*np.sqrt(a/t)))
+            FA=(np.sqrt(c/a)*(1+0.04*c/a)+0.2*(c/a)**4*(a/t)**2-0.11*(c/a)**4*(a/t)**4)*g*fphai*fw
+            H=1+(-2.11+0.77*c/a)*a/t+(0.55-0.72*(c/a)**0.75+0.14*(c/a)**1.5)*(a/t)**2
+        KA=FA*(Sm+H*Sb)*np.sqrt(np.pi*a/Q)
+        if a/c <=1.0:
+            Q=1+1.464*(a/c)**1.65
+            g=1.1+0.35*(a/t)**2
+            fphai=np.sqrt(a/c)
+            fw=np.sqrt(1.0/np.cos(np.pi*c/(2*b)*np.sqrt(a/t)))
+            FB=(1.13-0.09*a/c+(-0.54+0.89/(0.2+a/c))*(a/t)**2+(0.5-1/(0.65+a/c)+14*(1-a/c)**24)*(a/t)**4)*g*fphai*fw
+            H=1-0.34*a/t-0.11*a/c*a/t
+        else:
+            Q=1+1.464*(c/a)**1.65
+            g=1.1+0.35*c/a*(a/t)**2
+            fphai=1
+            fw=np.sqrt(1.0/np.cos(np.pi*c/(2*b)*np.sqrt(a/t)))
+            FB=(np.sqrt(c/a)*(1+0.04*c/a)+0.2*(c/a)**4*(a/t)**2-0.11*(c/a)**4*(a/t)**4)*g*fphai*fw
+            H=1+(-0.04-0.41*c/a)*a/t+(0.55-1.93*(c/a)**0.75+1.38*(c/a)**1.5)*(a/t)**2
+        KB=FB*(Sm+H*Sb)*np.sqrt(np.pi*a/Q)
+        res={
+            'KA':KA,
+            'KB':KB
+        }
+        super().SetRes(res)
+class K_1_a_2(Base):
+    #by S.Sakai
+    def __init__(self):
+        super().SetTitle('半楕円表面き裂 ASME Section XI, Appendix A の解')
+
+        super().SetRefer(
+            "ASME Boiler and Pressure Vessel Code, Section XI, Rules for Inservice Inspection of Nuclear Power Plant Components, 2004\n"
+            "Cipolla, R. C.: Technical Basis for the Residual Stress Intensity Factor Equation for Surface Flaws in ASME Section XI Appendix A, ASME PVP, Vol. 313-1, p. 105, 1995"
+        )
+
+        # Applicable range: 0 < a/t ≤ 0.8, 0 < a/c ≤ 1 
+
+    def Calc(self):
+        dm = dmanage()
+        df = super().GetData()
+
+        a = df['a']          # crack depth
+        t = df['t']          # thickness
+        c = df['c']          # half surface length
+        sigma0 = df['sigma0']
+        sigma1 = df['sigma1']
+        sigma2 = df['sigma2']
+        sigma3 = df['sigma3']
+        sigmaY = df['Sy']
+
+        data = dm.Finput('K-1-a-2.csv')
+        target_point = np.array([[a / t, a / c]])
+
+        FA = np.zeros(4, dtype=float)           # F0A..F3A
+        for ith in range(1, 5):                 # 1,2,3,4  (→ F0A..F3A)
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FA[ith - 1] = float(Fi)
+
+        FB = np.zeros(4, dtype=float)           # F0B..F3B
+        for ith in range(5, 9):                 # 5,6,7,8  (→ F0B..F3B)
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FB[ith - 5] = float(Fi)
+
+        SA = FA[0]*sigma0 + FA[1]*sigma1 + FA[2]*sigma2 + FA[3]*sigma3
+        qyA = (SA**2) / (6.0 * sigmaY**2)
+        QA = 1.0 + 1.464 * (a / c)**1.65 - qyA
+        KA = SA * np.sqrt(np.pi * a / QA)
+
+        SB = FB[0]*sigma0 + FB[1]*sigma1 + FB[2]*sigma2 + FB[3]*sigma3
+        qyB = (SB**2) / (6.0 * sigmaY**2)
+        QB = 1.0 + 1.464 * (a / c)**1.65 - qyB
+        KB = SB * np.sqrt(np.pi * a / QB)
+
+        res = {
+            'KA': KA, 'KB': KB
+        }
+        super().SetRes(res)
+class K_1_a_3(Base):
+    #by S.Sakai
+    def __init__(self):
+        super().SetTitle('半楕円表面き裂 白鳥の解')
+
+        super().SetRefer(
+            "白鳥: 影響関数法による応力拡大係数の解析, 日本機械学会講習会教材, 表面き裂—その解析と評価—, No.900-2, p. 1, 1990"
+        )
+
+        # Applicable range: 0.1 ≤ a/t ≤ 0.8, 0.2 ≤ a/c ≤ 2
+
+    def Calc(self):
+        dm = dmanage()
+        df = super().GetData()
+
+        a = df['a']          # crack depth
+        t = df['t']          # thickness
+        c = df['c']          # half surface length
+        sigma0 = df['sigma0']
+        sigma1 = df['sigma1']
+        sigma2 = df['sigma2']
+        sigma3 = df['sigma3']
+        # sigma0, sigma1, sigma2, sigma3 は、き裂深さ方向の応力分布
+        #   σ = σ0 + σ1·ξ + σ2·ξ² + σ3·ξ³
+        # における多項式展開の係数を表す
+        # ここで ξ = 1 - u/a （u：表面からの距離，a：き裂深さ）
+        data = dm.Finput('K-1-a-3.csv')
+        target_point = np.array([[a / t, a / c]])
+
+        FA = np.zeros(4, dtype=float)           # F0A..F3A
+        for ith in range(1, 5):                 # 1,2,3,4  (→ F0A..F3A)
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FA[ith - 1] = float(Fi)
+
+        FB = np.zeros(4, dtype=float)           # F0B..F3B
+        for ith in range(5, 9):                 # 5,6,7,8  (→ F0B..F3B)
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FB[ith - 5] = float(Fi)
+
+        Q = 1.0 + 1.464 * (a / c)**1.65
+
+        # 最深点の応力拡大係数
+        SA = FA[0]*sigma0 + FA[1]*sigma1 + FA[2]*sigma2 + FA[3]*sigma3
+        KA = SA * np.sqrt(np.pi * a / Q)
+
+        # 表面点の応力拡大係数
+        SB = FB[0]*sigma0 + FB[1]*sigma1 + FB[2]*sigma2 + FB[3]*sigma3
+        KB = SB * np.sqrt(np.pi * a / Q)
+
+        res = {
+            'KA': KA, 'KB': KB
+        }
+        super().SetRes(res)
+class K_1_b_1(Base):
+    def __init__(self):
+        super().SetTitle('長い表面き裂(片側) Tada らの解')
+        super().SetRefer('Tada, H., Paris, P. C.,and Irwin, G. R.: The Stress Analysis of Cracks Handbook, Third edition, ASME, 2000')
+
+        # Applicable range: 0 < a/t < 1 
+
+    def Calc(self):
+        df = super().GetData()
+
+        a = df['a']          # crack depth
+        t = df['t']          # thickness
+        b = df['b']          # half width
+        P = df['P']          # axial load
+        M = df['M']          # bending moment
+
+        sigma_m = P / (2.0 * b * t) 
+        sigma_b = 3.0 * M / (b * t**2)
+        theta = np.pi * a / (2.0 * t)
+        root_term = np.sqrt((2.0 * t) / (np.pi * a) * np.tan(theta))
+        cos_term = np.cos(theta)
+        Fm = root_term * (0.752 + 2.02 * (a / t) + 0.37 * (1.0 - np.sin(theta))**3) / cos_term
+        Fb = root_term * (0.923 + 0.199 * (1.0 - np.sin(theta))**4) / cos_term
+
+        K = (Fm * sigma_m + Fb * sigma_b) * np.sqrt(np.pi * a)
+
+        res = {'K': K}
+        super().SetRes(res)
+class K_1_c_1(Base):
+    def __init__(self):
+        super().SetTitle('長い表面き裂(両側) Tada らの解')
+        super().SetRefer('Tada, H., Paris, P. C., and Irwin, G. R.: The Stress Analysis of Cracks Handbook, Third edition, ASME, 2000')
+
+        # Applicable range: 0 < a/t < 0.5
+
+    def Calc(self):
+        df = super().GetData()
+
+        a = df['a']          # crack depth
+        t = df['t']          # thickness
+        b = df['b']          # half width
+        P = df['P']          # axial load
+        M = df['M']          # bending moment
+
+        sigma_m = P / (2.0 * b * t)
+        theta = np.pi * a / t
+        Fm = (1.0 + 0.122 * np.cos(theta)**4) * np.sqrt((t / (np.pi * a)) * np.tan(theta))
+        
+        K = Fm * sigma_m * np.sqrt(np.pi * a)
+
+        res = {'K': K}
+        super().SetRes(res)
+
+class K_1_d_1(Base):
+    def __init__(self):
+        super().SetTitle('中央貫通き裂 Shih らの解（き裂が短い場合）')
+        super().SetRefer('Shih, G. C., Paris, P. C., and Erdogan, F.: Stress Intensity Factors for Plane Extension and Plate Bending Problems, Trans. ASME, J. of Applied Mechanics, 29, p. 306, 1962')
+
+        # Applicable range: 0 < c/b ≤ 0.15
+
+    def Calc(self):
+        df = super().GetData()
+
+        c = df['c']          # crack half-length
+        t = df['t']          # thickness
+        b = df['b']          # half width
+        P = df['P']          # axial load
+        M = df['M']          # bending moment
+
+        sigma_m = P/(2.0*b*t)
+        sigma_b = 3.0*M/(b*t**2)
+
+        #開口側の応力拡大係数
+
+        FmA, FbA = 1.0, 1.0
+
+        KA = (FmA*sigma_m + FbA*sigma_b)*np.sqrt(np.pi*c)
+
+        #閉口側の応力拡大係数
+
+        FmB, FbB = 1.0, -1.0
+        
+        KB = (FmB*sigma_m + FbB*sigma_b)*np.sqrt(np.pi*c)
+
+        res = {'KA': KA, 'KB': KB}
+        super().SetRes(res)
+class K_1_d_2(Base):
+    def __init__(self):
+        super().SetTitle('中央貫通き裂 Tada らの解（き裂が長い場合）')
+        super().SetRefer('Tada, H., Paris, P. C., and Irwin, G. R.: The Stress Analysis of Cracks Handbook, Third edition, ASME, 2000')
+
+        # Applicable range: 0 < c/b < 1
+
+    def Calc(self):
+        df = super().GetData()
+
+        c = df['c']          # crack half-length
+        t = df['t']          # thickness
+        b = df['b']          # half width
+        P = df['P']          # axial load
+        M = df['M']          # bending moment
+
+        sigma_m = P/(2.0*b*t)
+        ratio = c/b
+        theta = np.pi*c/(2.0*b)
+        Fm = (1.0 - 0.025*ratio**2 + 0.06*ratio**4) / np.sqrt(np.cos(theta))
+        
+        K = Fm * sigma_m * np.sqrt(np.pi*c)
+
+        res = {'K': K}
+        super().SetRes(res)
+class K_1_e_1(Base):
+    def __init__(self):
+        super().SetTitle('楕円内部き裂 Ovchinnikov–Vasiltchenko の解')
+        super().SetRefer("Ovchinnikov, A. V., and Vasiltchenko, G. S.: The Defect Schematization and SIF Determination for Assessment of the Vessel and Piping Integrity, Final Report, CHIITMASH Project 125-01-90, p. 1, 1990")
+
+        # Applicable range: 0 < a/t ≤ 0.45(1 - 2e/t), 0 < a/c ≤ 1, 0 ≤ e/t ≤ 0.5, 0 < c/b ≤ 0.15
+
+    def Calc(self):
+        df = super().GetData()
+
+        a = df['a']          # crack depth (semi-minor)
+        c = df['c']          # crack half-length (semi-major)
+        t = df['t']          # thickness
+        b = df['b']          # half width
+        e = df['e']          # offset from mid-plane
+        P = df['P']          # axial load
+        M = df['M']          # bending moment
+
+        sigma_m = P/(2.0*b*t)
+        sigma_b = 3.0*M/(b*t**2)
+
+        ac = a/c
+        at = a/t
+        et = e/t
+        L = (2.0*at)/(1.0 - 2.0*et)
+
+        #評価点Aにおける応力拡大係数
+
+        denom_A = (1.0 - (L**1.8)*(1.0 - 0.4*ac - et**2))**0.54
+        FmA = (1.01 - 0.37*ac)/denom_A
+        FbA = (1.01 - 0.37*ac)*(2.0*et + at + 0.34*ac*at)/denom_A
+
+        KA = (FmA*sigma_m + FbA*sigma_b)*np.sqrt(np.pi*a)
+
+        #評価点Bにおける応力拡大係数
+
+        denom_B = (1.0 - (L**1.8)*(1.0 - 0.4*ac - 0.8*et**2))**0.54
+        FmB = (1.01 - 0.37*ac)/denom_B
+        FbB = (1.01 - 0.37*ac)*(2.0*et + at + 0.34*ac*at)/denom_B
+
+        KB = (FmB*sigma_m + FbB*sigma_b)*np.sqrt(np.pi*a)
+
+        res = {'KA': KA, 'KB': KB}
+        super().SetRes(res)
+class K_1_e_2(Base):
+    def __init__(self):
+        super().SetTitle('楕円内部き裂 ASME Section XI, Appendix A の解')
+
+        super().SetRefer(
+            "ASME Boiler and Pressure Vessel Code, Section XI, Rules for Inservice Inspection of Nuclear Power Plant Components, 2004"
+        )
+
+        # Applicable range: 0.075 < a/t ≤ 0.325, 0 ≤ e/t ≤ 0.35
+
+    def Calc(self):
+        dm = dmanage()
+        df = super().GetData()
+
+        a = df['a']          # crack half-depth (semi-minor axis)
+        b = df['b']          # plate half-width
+        c = df['c']          # crack half-length (semi-major axis)
+        t = df['t']          # thickness
+        e = df['e']          # offset from mid-thickness
+        P = df['P']          # axial force 
+        M = df['M']          # bending moment
+        Sy = df['Sy']        # yield strength
+
+        sigma_m = P / (2.0 * b * t)
+        sigma_b = 3.0 * M / (b * t**2)
+
+        data = dm.Finput('K-1-e-2.csv')
+        target_point = np.array([[a / t, e / t]])
+
+        # 1 -> FmA, 2 -> FmB, 3 -> FbA, 4 -> FbB
+        
+        X1, W1 = dm.CSV2XW(data, 1)
+        r2_score = dm.KrigCalc(X1, W1)
+        FmA, _sigma = dm.Eval(target_point)
+        FmA = float(FmA)
+
+        X2, W2 = dm.CSV2XW(data, 2)
+        r2_score = dm.KrigCalc(X2, W2)
+        FmB, _sigma = dm.Eval(target_point)
+        FmB = float(FmB)
+        
+        X3, W3 = dm.CSV2XW(data, 3)
+        r2_score = dm.KrigCalc(X3, W3)
+        FbA, _sigma = dm.Eval(target_point)
+        FbA = float(FbA)
+        
+        X4, W4 = dm.CSV2XW(data, 4)
+        r2_score = dm.KrigCalc(X4, W4)
+        FbB, _sigma = dm.Eval(target_point)
+        FbB = float(FbB)
+
+        # 評価点Aにおける応力拡大係数
+        SA = FmA * sigma_m + FbA * sigma_b
+        qyA = (SA**2) / (6.0 * Sy**2)
+        QA = 1.0 + 1.464 * (a / c)**1.65 - qyA
+        KA = SA * np.sqrt(np.pi * a / QA)
+
+        # 評価点Bにおける応力拡大係数
+        SB = FmB * sigma_m + FbB * sigma_b
+        qyB = (SB**2) / (6.0 * Sy**2)
+        QB = 1.0 + 1.464 * (a / c)**1.65 - qyB
+        KB = SB * np.sqrt(np.pi * a / QB)
+
+        super().SetRes({'KA': KA, 'KB': KB})
+class K_1_e_3(Base):
+    def __init__(self):
+        super().SetTitle('楕円内部き裂 Raju–Newman の解')
+        super().SetRefer('Newman, J. C. Jr., and Raju, I. S.: Stress-Intensity Factor Equations for Cracks in Three-Dimensional Finite Bodies, NASA Technical Memorandum 83200, NASA, 1981')
+
+        # Applicable range: c/b < 0.5
+
+    def Calc(self):
+        df = super().GetData()
+
+        a = df['a']          # crack semi-minor (depth)
+        c = df['c']          # crack semi-major (half-length)
+        t = df['t']          # thickness
+        b = df['b']          # half width
+        P = df['P']          # axial load
+
+        sigma_m = P/(2.0*b*t)
+
+        ac = a/c
+        at = a/t
+        theta_w = (np.pi*c/(2.0*b))*np.sqrt(at)
+        f_w = (1.0/np.cos(theta_w))**0.5
+
+        if ac <= 1.0:
+            Q = 1.0 + 1.464*(ac**1.65)
+            base = 1.0 + (0.05/(0.11 + ac**1.5))*(at**2) + (0.29/(0.23 + ac**1.5))*(at**4)
+            # Point A
+            gA = 1.0
+            fphiA = 1.0
+            FA = base*gA*fphiA*f_w
+            # Point C
+            gC = 1.0 - (at**4)/(1.0 + 4.0*(ac))
+            fphiC = np.sqrt(ac)
+            FC = base*gC*fphiC*f_w
+        else:
+            Q = 1.0 + 1.464*((c/a)**1.65)
+            base = np.sqrt(c/a) + (0.05/(0.11 + ac**1.5))*(at**2) + (0.29/(0.23 + ac**1.5))*(at**4)
+            # Point A
+            gA = 1.0
+            fphiA = np.sqrt(c/a)
+            FA = base*gA*fphiA*f_w
+            # Point C
+            gC = 1.0 - (at**4)/(1.0 + 4.0*(ac))
+            fphiC = 1.0
+            FC = base*gC*fphiC*f_w
+
+        #評価点Aにおける応力拡大係数
+
+        KA = FA*sigma_m*np.sqrt(np.pi*a/Q)
+
+        #評価点Cにおける応力拡大係数
+
+        KC = FC*sigma_m*np.sqrt(np.pi*a/Q)
+
+        res = {'KA': KA, 'KC': KC}
+        super().SetRes(res)
+
+class K_2_a_1(Base):
+    def __init__(self):
+        super().SetTitle('軸方向内表面半楕円表面き裂　Fett らの解')
+
+        super().SetRefer(
+            "Fett, T., Munz, D., and Neuman, J.: Local Stress Intensity Factors for Surface Cracks in Plates under Power-Shaped Stress Distributions, Engineering Fracture Mechanics, 36, 4, p. 647, 1990\n"
+            "Raju, I. S. and Newman, J. C.: Stress Intensity Factor Influence Coefficients for Internal and External Surface Cracks in Cylindrical Vessels, ASME PVP, 58, p. 37, 1978"
+        )
+
+        # Applicable range: 0 < a/t ≤ 0.8, 0.2 ≤ a/c ≤ 1, 4 ≤ Ri/t ≤ 10
+
+    def Calc(self):
+        dm = dmanage()
+        df = super().GetData()
+
+        a = df['a']          # crack depth
+        t = df['t']          # wall thickness
+        c = df['c']          # half surface length
+        Ri = df['Ri']        # inner radius
+        sigma0 = df['sigma0']
+        sigma1 = df['sigma1']
+        sigma2 = df['sigma2']
+        sigma3 = df['sigma3']
+
+        data = dm.Finput('K-2-a-1.csv')
+        target_point = np.array([[a / t, a / c, Ri / t]])
+
+        FA = np.zeros(4, dtype=float)           # F0A..F3A
+        for ith in range(1, 5):                 # 1,2,3,4  (→ F0A..F3A)
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FA[ith - 1] = float(Fi)
+
+        FB = np.zeros(4, dtype=float)           # F0B..F3B
+        for ith in range(5, 9):                 # 5,6,7,8  (→ F0B..F3B)
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FB[ith - 5] = float(Fi)
+
+        SA = FA[0]*sigma0 + FA[1]*sigma1 + FA[2]*sigma2 + FA[3]*sigma3
+        SB = FB[0]*sigma0 + FB[1]*sigma1 + FB[2]*sigma2 + FB[3]*sigma3
+
+        # 最深点の応力拡大係数
+        KA = SA * np.sqrt(np.pi * a)
+
+        # 表面点の応力拡大係数
+        KB = SB * np.sqrt(np.pi * a)
+
+        res = {
+            'KA': KA, 'KB': KB
+        }
+        super().SetRes(res)
+class K_2_a_2(Base):
+    def __init__(self):
+        super().SetTitle('軸方向内表面半楕円表面き裂　白鳥の解')
+
+        super().SetRefer(
+            "白鳥: 影響関数法による応力拡大係数の解析, 日本機械学会講習会教材, 表面き裂—その解析と評価—, No.900-2, p. 1, 1990"
+        )
+
+        # Applicable range: 0.1 ≤ a/t ≤ 0.8, 0.2 ≤ a/c ≤ 1, 1/9 ≤ Ri/t ≤ 10
+
+    def Calc(self):
+        dm = dmanage()
+        df = super().GetData()
+
+        a = df['a']          # crack depth
+        t = df['t']          # wall thickness
+        c = df['c']          # half surface crack length
+        Ri = df['Ri']        # inner radius
+        sigma0 = df['sigma0']
+        sigma1 = df['sigma1']
+        sigma2 = df['sigma2']
+        sigma3 = df['sigma3']
+
+        data = dm.Finput('K-2-a-2.csv')
+        target_point = np.array([[a / t, a / c, Ri / t]])
+
+        FA = np.zeros(4, dtype=float)           # F0A..F3A
+        for ith in range(1, 5):                 # 1,2,3,4  (→ F0A..F3A)
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FA[ith - 1] = float(Fi)
+
+        FB = np.zeros(4, dtype=float)           # F0B..F3B
+        for ith in range(5, 9):                 # 5,6,7,8  (→ F0B..F3B)
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FB[ith - 5] = float(Fi)
+
+        SA = FA[0]*sigma0 + FA[1]*sigma1 + FA[2]*sigma2 + FA[3]*sigma3
+        SB = FB[0]*sigma0 + FB[1]*sigma1 + FB[2]*sigma2 + FB[3]*sigma3
+
+        Q = 1.0 + 1.464 * (a / c)**1.65
+
+        # 最深点の応力拡大係数
+        KA = SA * np.sqrt(np.pi * a / Q)
+
+        # 表面点の応力拡大係数
+        KB = SB * np.sqrt(np.pi * a / Q)
+
+        res = {
+            'KA': KA, 'KB': KB
+        }
+        super().SetRes(res)
+
+class K_2_a_3(Base):
+    # by S.Sakai
+    def __init__(self):
+        super().SetTitle('軸方向無い表面半楕円表面亀裂，Zahoorの解')
+        super().SetRefer('Zahoor,A.:Ductile Fracture Handbook Volume 3,EPRI NP-6301-D,1991')
+    def Calc(self):
+        df=super().GetData()
+        Ro=df['Ro']
+        Ri=df['Ri']
+        p=df['p']
+        Sm=(Ro*Ro+Ri*Ri)/(Ro*Ro-Ri*Ri)*p
+        a=df['a']
+        c=df['c']
+        t=df['t']
+        ar=(a/t)/(a/c)**0.58
+        FA=0.25+(0.4759*ar+0.1262*ar*ar)/(0.102*(Ri/t)-0.02)**0.1
+        KA=FA*Sm*np.sqrt(np.pi*t)
+        FB=FA*(1.06+0.28*(a/t)**2)*(a/c)**0.41
+        KB=FB*Sm*np.sqrt(np.pi*t)
+        res={'KA':KA,
+             'KB':KB}
+        super().SetRes(res)
+class K_2_a_3(Base):
+    def __init__(self):
+        super().SetTitle('軸方向内表面半楕円表面き裂 Zahoor の解')
+        super().SetRefer('Zahoor, A.: Ductile Fracture Handbook Volume 3, EPRI NP-6301-D, 1991')
+
+        # Applicable range: 0.05 ≤ a/t ≤ 0.85, 0.1 ≤ a/c ≤ 1, 0.2 ≤ α, 1 ≤ Ri/t ≤ 10
+
+    def Calc(self):
+        df = super().GetData()
+
+        a  = df['a']          # crack depth
+        c  = df['c']          # half surface crack length
+        t  = df['t']          # wall thickness
+        Ri = df['Ri']         # inner radius
+        Ro = df['Ro']         # outer radius
+        p  = df['p']          # internal pressure
+
+        sigma_m = ((Ro**2 + Ri**2) / (Ro**2 - Ri**2)) * p
+        alpha   = (a/t) / ((a/c)**0.58)
+
+        # 最深点の応力拡大係数
+        denomA = (0.102*(Ri/t) - 0.02)**0.1
+        FA = 0.25 + (0.4759*alpha + 0.1262*alpha**2) / denomA
+        KA = FA * sigma_m * np.sqrt(np.pi * t)
+
+        # 表面点の応力拡大係数
+        FB = FA * (1.06 + 0.28*(a/t)**2) * (a/c)**0.41
+        KB = FB * sigma_m * np.sqrt(np.pi * t)
+
+        res = {'KA': KA, 'KB': KB}
+        super().SetRes(res)
+class K_2_b_1(Base):
+    def __init__(self):
+        super().SetTitle('軸方向内表面長い表面き裂 Fuhley-Osage の解')
+
+        super().SetRefer(
+            "American Petroleum Institute: Recommended Practice for Fitness-for-Service, "
+            "API RP 579, 2000"
+        )
+
+        # Applicable range: 0 ≤ a/t ≤ 0.8, 2 ≤ Ri/t ≤ 1000
+
+    def Calc(self):
+        dm = dmanage()
+        df = super().GetData()
+
+        a = df['a']          # crack depth
+        t = df['t']          # thickness
+        Ri = df['Ri']        # inner radius
+        sigma0 = df['sigma0']
+        sigma1 = df['sigma1']
+        sigma2 = df['sigma2']
+        sigma3 = df['sigma3']
+        sigma4 = df['sigma4']
+
+        data = dm.Finput('K-2-b-1.csv')
+        target_point = np.array([[a / t, Ri / t]])
+
+        F = np.zeros(5, dtype=float)            # F0..F4
+        for ith in range(1, 6):                 # 1,2,3,4,5 (→ F0..F4)
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            F[ith - 1] = float(Fi)
+
+        # 応力拡大係数
+        S = F[0]*sigma0 + F[1]*sigma1*(a/t) + F[2]*sigma2*(a/t)**2 + F[3]*sigma3*(a/t)**3 + F[4]*sigma4*(a/t)**4
+        K = S * np.sqrt(np.pi * a)
+
+        res = {
+            'K': K
+        }
+        super().SetRes(res)
+
+class K_2_b_2(Base):
+    #by S.Sakai
+    def __init__(self):
+        super().SetTitle('軸方向内表面長い表面亀裂，Zahoorの解')
+        super().SetRefer('Zahoor,A.:Closed Form Expressions for Fracture Mechanics Analysis of Cracked Pipes, Trans.ASME, J. of Pressure Vessel Technology,107,p.203,1987')
+    def Calc(self):
+        df=super().GetData()
+        Ro=df['Ro']
+        Ri=df['Ri']
+        p=df['p']
+        Sm=2*Ro*Ro/(Ro*Ro-Ri*Ri)*p
+        a=df['a']
+        t=Ro-Ri
+        if Ri/t >=5.0 and Ri/t<10.0:
+            A=(0.125*(Ri/t)-0.25)**0.25
+        if Ri/t >=10.0 and Ri/t<=20.0:
+            A=(0.2*(Ri/t)-1)**0.25
+        F=1.1+A*(4.951*(a/t)**2+1.092*(a/t)**4)
+        K=F*Sm*np.sqrt(np.pi*a)
+        res={'K':K}
+        super().SetRes(res)
+class K_2_c_1(Base):
+    def __init__(self):
+        super().SetTitle('軸方向外表面半楕円表面き裂　Fett らの解')
+
+        super().SetRefer(
+            "Fett, T., Munz, D., and Neuman, J.: Local Stress Intensity Factors for Surface Cracks in Plates under Power-Shaped Stress Distributions, Engineering Fracture Mechanics, 36, 4, p. 647, 1990\n"
+            "Raju, I. S. and Newman, J. C.: Stress Intensity Factor Influence Coefficients for Internal and External Surface Cracks in Cylindrical Vessels, ASME PVP, 58, p. 37, 1978"
+        )
+
+        # Applicable range: 0 < a/t ≤ 0.8, 0.2 ≤ a/c ≤ 1, 4 ≤ Ri/t ≤ 10
+
+    def Calc(self):
+        dm = dmanage()
+        df = super().GetData()
+
+        a = df['a']          # crack depth
+        t = df['t']          # wall thickness
+        c = df['c']          # half surface crack length
+        Ri = df['Ri']        # inner radius
+        sigma0 = df['sigma0']
+        sigma1 = df['sigma1']
+        sigma2 = df['sigma2']
+        sigma3 = df['sigma3']
+
+        data = dm.Finput('K-2-c-1.csv')
+        target_point = np.array([[a / t, a / c, Ri / t]])
+
+        FA = np.zeros(4, dtype=float)           # F0A..F3A
+        for ith in range(1, 5):                 # 1,2,3,4  (→ F0A..F3A)
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FA[ith - 1] = float(Fi)
+
+        FB = np.zeros(4, dtype=float)           # F0B..F3B
+        for ith in range(5, 9):                 # 5,6,7,8  (→ F0B..F3B)
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FB[ith - 5] = float(Fi)
+
+        SA = FA[0]*sigma0 + FA[1]*sigma1 + FA[2]*sigma2 + FA[3]*sigma3
+        SB = FB[0]*sigma0 + FB[1]*sigma1 + FB[2]*sigma2 + FB[3]*sigma3
+
+        # 最深点の応力拡大係数
+        KA = SA * np.sqrt(np.pi * a)
+
+        # 表面点の応力拡大係数
+        KB = SB * np.sqrt(np.pi * a)
+
+        res = {
+            'KA': KA, 'KB': KB
+        }
+        super().SetRes(res)
+class K_2_c_2(Base):
+    def __init__(self):
+        super().SetTitle('軸方向外表面半楕円表面き裂 Zahoor の解')
+        super().SetRefer('Zahoor, A.: Ductile Fracture Handbook Volume 3, EPRI NP-6301-D, 1991')
+
+        # Applicable range: 0.05 ≤ a/t ≤ 0.85, 0.1 ≤ a/c ≤ 1, 0.2 ≤ α, 2 ≤ Ri/t ≤ 10
+
+    def Calc(self):
+        df = super().GetData()
+
+        a  = df['a']    # crack depth
+        c  = df['c']    #half surface crack length
+        t  = df['t']    # wall thickness
+        Ri = df['Ri']   # inner radius
+        Ro = df['Ro']   # outer radius
+        p  = df['p']    # internal pressure
+
+        # 最深点の応力拡大係数
+        sigma_m = (2.0 * Ri**2 / (Ro**2 - Ri**2)) * p
+        alpha = (a/t) / ((a/c)**0.58)
+        denomA = (0.11*(Ri/t) - 0.1)**0.1
+        FA = 0.25 + (0.42*alpha + 0.21*alpha**2) / denomA
+        KA = FA * sigma_m * np.sqrt(np.pi * t)
+
+        # 表面点の応力拡大係数
+        FB = FA * (1.0 + 0.33*(a/t)**2) * (a/c)**0.47
+        KB = FB * sigma_m * np.sqrt(np.pi * t)
+
+        res = {'KA': KA, 'KB': KB}
+        super().SetRes(res)
+class K_2_d(Base):
+    def __init__(self):
+        super().SetTitle('軸方向外表面長い表面き裂 Fuhley-Osage の解')
+
+        super().SetRefer(
+            "American Petroleum Institute: Recommended Practice for Fitness-for-Service, "
+            "API RP 579, 2000"
+        )
+
+        # Applicable range: 0 ≤ a/t ≤ 0.8, 2 ≤ Ri/t ≤ 1000
+
+    def Calc(self):
+        dm = dmanage()
+        df = super().GetData()
+
+        a = df['a']          # crack depth
+        t = df['t']          # wall thickness
+        Ri = df['Ri']        # inner radius
+        sigma0 = df['sigma0']
+        sigma1 = df['sigma1']
+        sigma2 = df['sigma2']
+        sigma3 = df['sigma3']
+        sigma4 = df['sigma4']
+
+        data = dm.Finput('K-2-d.csv')
+        target_point = np.array([[a / t, Ri / t]])
+
+        F = np.zeros(5, dtype=float)            # F0..F4
+        for ith in range(1, 6):                 # 1,2,3,4,5 (→ F0..F4)
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            F[ith - 1] = float(Fi)
+
+        S = F[0]*sigma0 + F[1]*sigma1*(a/t) + F[2]*sigma2*(a/t)**2 + F[3]*sigma3*(a/t)**3 + F[4]*sigma4*(a/t)**4
+        K = S * np.sqrt(np.pi * a)
+
+        res = {
+            'K': K
+        }
+        super().SetRes(res)
+class K_2_e_1(Base):
+    def __init__(self):
+        super().SetTitle('軸方向貫通き裂　Erdogan-Kibler の解')
+
+        super().SetRefer(
+            "Erdogan, F., and Kibler, J. J.: Cylindrical and Spherical Shells with Cracks, International Journal of Fracture Mechanics, 5, p. 229, 1969"
+        )
+
+        # Applicable range: 0 < c/t ≤ 12.5, 10 ≤ Ri/t ≤ 20
+
+    def Calc(self):
+        dm = dmanage()
+        df = super().GetData()
+
+        c = df['c']          # half crack length
+        t = df['t']          # wall thickness
+        Ri = df['Ri']        # inner radius
+        sigma_m = df['sigma_m']
+        sigma_b = df['sigma_b']
+
+        data = dm.Finput('K-2-e-1.csv')
+        target_point = np.array([[c / t, Ri / t]])
+
+        FA = np.zeros(2, dtype=float)           # FmA, FbA
+        for ith in range(1, 3):                 # 1,2  (→ FmA, FbA)
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FA[ith - 1] = float(Fi)
+
+        FB = np.zeros(2, dtype=float)           # FmB, FbB
+        for ith in range(3, 5):                 # 3,4  (→ FmB, FbB)
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FB[ith - 3] = float(Fi)
+
+        # 内表面における応力拡大係数
+        KA = (FA[0]*sigma_m + FA[1]*sigma_b) * np.sqrt(np.pi * c)
+
+        # 外表面における応力拡大係数
+        KB = (FB[0]*sigma_m + FB[1]*sigma_b) * np.sqrt(np.pi * c)
+
+        res = {
+            'KA': KA, 'KB': KB
+        }
+        super().SetRes(res)
+
+class K_2_e_2(Base):
+    #by S.Sakai
+    def __init__(self):
+        super().SetTitle('軸方向貫通亀裂，ASME Code Case N-513の解')
+        super().SetRefer('ASME Boiler and Pressure Vessel Code, Code Case N-513, Evaluation Criteria for Temporary Acceptance of Flaws in Calss 3 Piping, 1997')
+    def Calc(self):
+        df=super().GetData()
+        p=df['p']
+        R=df['R']
+        t=df['t']
+        Sm=p*R/t
+        c=df['c']
+        l=c/np.sqrt(R*t)
+        F=1+0.072449*l+0.64856*l*l-0.2327*l*l*l+0.038154*l**4-0.0023487*l**5
+        K=F*Sm*np.sqrt(np.pi*c)
+        res={'K':K}
+        super().SetRes(res)
+class K_2_e_3(Base):
+    def __init__(self):
+        super().SetTitle('軸方向貫通き裂　Zang の解')
+
+        super().SetRefer(
+            "Zang, W.: Stress Intensity Factor Solutions for Axial and Circumferential Through-Wall Cracks in Cylinders, "
+            "SINTAP/SAQ/02, SAQ Control AB, 1997"
+        )
+
+        # Applicable range: 0.5 ≤ c/t ≤ 25, 5 ≤ Ri/t ≤ 100
+
+    def Calc(self):
+        dm = dmanage()
+        df = super().GetData()
+
+        c = df['c']          # half crack length
+        t = df['t']          # wall thickness
+        Ri = df['Ri']        # inner radius
+        sigma0 = df['sigma0']
+        sigma1 = df['sigma1']
+        sigma2 = df['sigma2']
+        sigma3 = df['sigma3']
+        sigma4 = df['sigma4']
+
+        data = dm.Finput('K-2-e-3.csv')
+        target_point = np.array([[c / t, Ri / t]])
+
+        FA = np.zeros(5, dtype=float)           # F0A..F4A
+        for ith in range(1, 6):                 # 1..5 -> F0A..F4A
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FA[ith - 1] = float(Fi)
+
+        FB = np.zeros(5, dtype=float)           # F0B..F4B
+        for ith in range(6, 11):                # 6..10 -> F0B..F4B
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FB[ith - 6] = float(Fi)
+
+        SA = FA[0]*sigma0 + FA[1]*sigma1 + FA[2]*sigma2 + FA[3]*sigma3 + FA[4]*sigma4
+        SB = FB[0]*sigma0 + FB[1]*sigma1 + FB[2]*sigma2 + FB[3]*sigma3 + FB[4]*sigma4
+
+        # 内表面における応力拡大係数
+        KA = SA * np.sqrt(np.pi * c)
+
+        # 外表面における応力拡大係数
+        KB = SB * np.sqrt(np.pi * c)
+
+        res = {
+            'KA': KA, 'KB': KB
+        }
+        super().SetRes(res)
+class K_2_f_1(Base):
+    def __init__(self):
+        super().SetTitle('周方向内表面半楕円表面き裂　Chapuliot らの解')
+        super().SetRefer(
+            "Chapuliot, S.: Formulaire de KI Pour les Tubes Comportant un Defaut de Surface Semi-elliptique "
+            "Longitudinal ou Circonférentiel, interne ou externe, Rapport CEA-R-5900, 2000"
+        )
+        # Applicable range: 0 < a/t ≤ 0.8, 0 < a/c ≤ 1, 1 ≤ Ri/t < ∞
+
+    def Calc(self):
+        dm = dmanage()
+        df = super().GetData()
+
+        a = df['a']          # crack depth
+        t = df['t']          # wall thickness
+        c = df['c']          # half surface crack length
+        Ri = df['Ri']        # inner radius
+        sigma0 = df['sigma0']
+        sigma1 = df['sigma1']
+        sigma2 = df['sigma2']
+        sigma3 = df['sigma3']
+        sigma_bg = df['sigma_bg']
+
+        data = dm.Finput('K-2-f-1.csv')
+        target_point = np.array([[a / c, a / t, Ri / t]])
+
+        FA = np.zeros(5, dtype=float)
+        for ith in range(1, 5 + 1):
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FA[ith - 1] = float(Fi)
+
+        FB = np.zeros(5, dtype=float)
+        for ith in range(6, 10 + 1):
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FB[ith - 6] = float(Fi)
+
+        SA = FA[0]*sigma0 + FA[1]*sigma1*(a/t) + FA[2]*sigma2*(a/t)**2 + FA[3]*sigma3*(a/t)**3 + FA[4]*sigma_bg
+        SB = FB[0]*sigma0 + FB[1]*sigma1*(a/t) + FB[2]*sigma2*(a/t)**2 + FB[3]*sigma3*(a/t)**3 + FB[4]*sigma_bg
+
+        # 最深点の応力拡大係数
+        KA = SA * np.sqrt(np.pi * a)
+
+        # 表面点の応力拡大係数
+        KB = SB * np.sqrt(np.pi * a)
+
+        res = {'KA': KA, 'KB': KB}
+        super().SetRes(res)
+class K_2_f_2(Base):
+    def __init__(self):
+        super().SetTitle('周方向内表面半楕円表面き裂 白鳥の解')
+
+        super().SetRefer(
+            "白鳥: 影響関数法による応力拡大係数の解析, 日本機械学会講習会教材, "
+            "表面き裂―その解析と評価―, No. 900-2, p. 1, 1990"
+        )
+
+        # Applicable range: 0.1 ≤ a/t ≤ 0.8, 0.2 ≤ a/c ≤ 1, 1.25 ≤ Ri/t ≤ 10
+
+    def Calc(self):
+        dm = dmanage()
+        df = super().GetData()
+
+        a = df['a']          # crack depth
+        t = df['t']          # wall thickness
+        c = df['c']          # half surface crack length
+        Ri = df['Ri']        # inner radius
+        sigma_m  = df['sigma_m']
+        sigma_bg = df['sigma_bg']
+
+        data = dm.Finput('K-2-f-2.csv')
+        target_point = np.array([[a / t, a / c, Ri / t]])
+
+        FA = np.zeros(2, dtype=float)           # FmA, FbgA
+        for ith in range(1, 3):                 # 1,2 (→ FmA, FbgA)
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FA[ith - 1] = float(Fi)
+
+        FB = np.zeros(2, dtype=float)           # FmB, FbgB
+        for ith in range(3, 5):                 # 3,4 (→ FmB, FbgB)
+            X_ith, W_ith = dm.CSV2XW(data, ith)
+            r2_score = dm.KrigCalc(X_ith, W_ith)
+            Fi, _sigma = dm.Eval(target_point)
+            FB[ith - 3] = float(Fi)
+
+        Q = 1.0 + 1.464 * (a / c)**1.65
+
+        # 最深点の応力拡大係数
+        SA = FA[0]*sigma_m + FA[1]*sigma_bg
+        KA = SA * np.sqrt(np.pi * a / Q)
+
+        # 表面点の応力拡大係数
+        SB = FB[0]*sigma_m + FB[1]*sigma_bg
+        KB = SB * np.sqrt(np.pi * a / Q)
+
+        res = {
+            'KA': KA, 'KB': KB
+        }
+        super().SetRes(res)
+class K_2_f_3(Base):
+    def __init__(self):
+        super().SetTitle('周方向内表面半楕円表面き裂 Zahoor の解')
+
+        super().SetRefer(
+            "Zahoor, A.: Closed Form Expressions for Fracture Mechanics Analysis of Cracked Pipes, "
+            "Trans. ASME, J. of Pressure Vessel Technology, 107, p. 203, 1985\n"
+            "Zahoor, A.: Ductile Fracture Handbook Volume 2, EPRI NP-6301-D, 1990"
+        )
+
+        # Applicable range (軸方向荷重に対して):
+        # 0.05 ≤ a/t ≤ 0.8, 1/6 ≤ a/c ≤ 2/3, 5 ≤ R/t ≤ 20
+        # (曲げモーメントに対して):
+        # 0.05 ≤ a/t ≤ 0.8, a/c ≤ 2/3, 5 ≤ R/t ≤ 160
+
+    def Calc(self):
+        dm = dmanage()
+        df = super().GetData()
+
+        a = df['a']          # crack depth
+        t = df['t']          # wall thickness
+        c = df['c']          # half surface crack length
+        Ri = df['Ri']        # inner radius
+        P  = df['P']         # axial load
+        M  = df['M']         # bending moment
+
+        data = dm.Finput('K-2-f-3.csv')
+        target_point = np.array([[a / t, a / c, Ri / t]])
+
+        X_1, W_1 = dm.CSV2XW(data, 1)          # 1列目 → F_bg^A
+        r2_score = dm.KrigCalc(X_1, W_1)
+        FbgA, _sigma = dm.Eval(target_point)
+        FbgA = float(FbgA)
+
+        sigma_m = P/(2.0*np.pi*Ri*t)
+        sigma_bg = M/(np.pi*Ri**2*t)
+        Q = 1.0 + 1.464 * (a / c)**1.65
+        alpha = 2.0 * c / t
+        FmA = (1.0 + (0.02 + 0.0103*alpha + 0.00617*alpha**2 + 0.0035*(1.0 + 0.7*alpha)*(Ri/t - 5.0)**0.7) * Q**2) / Q**0.5
+     
+        # 最深点の応力拡大係数
+        K = FmA*sigma_m * np.sqrt(np.pi * a) + FbgA*sigma_bg * np.sqrt(np.pi * t)
+
+        res = {
+            'K': K
+        }
+        super().SetRes(res)
+class K_2_g(Base):
+    def __init__(self):
+        super().SetTitle('周方向内表面扇形表面き裂 ASME Section XI, Appendix C の解')
+        super().SetRefer('ASME Boiler and Pressure Vessel Code, Section XI, Rules for Inservice Inspection of Nuclear Power Plant Components, 2004')
+
+        # Applicable range: 0.08 ≤ a/t ≤ 0.8, 0 < a/c ≤ 1, 0.05 ≤ c/(πRi) ≤ 1
+
+    def Calc(self):
+        df = super().GetData()
+
+        a  = df['a']   # crack depth
+        c  = df['c']   # half surface crack length
+        t  = df['t']   # wall thickness
+        Ri = df['Ri']  # inner radius
+        R  = df['R']   # mean radius
+        P  = df['P']   # axial load
+        M  = df['M']   # bending moment
+
+        sigma_m = P/(2.0*np.pi*R*t)
+        sigma_bg = M/(np.pi*R**2*t)
+
+        cp = c/(np.pi*Ri)
+        if 0.5 <= cp <= 1.0:
+            cp = 0.5
+
+        x = (a/t)*cp
+
+        # 最深点の応力拡大係数
+        FmA  = 1.1 + (a/t)*(0.15241 + 16.722*x**0.855 - 14.944*x)
+        FbgA = 1.1 + (a/t)*(-0.09967 + 5.0057*x**0.565 - 2.8329*x)
+        KA = (FmA*sigma_m + FbgA*sigma_bg)*np.sqrt(np.pi*a)
+
+        res = {'KA': KA}
+        super().SetRes(res)
+
+class K_2_k_2(Base):
+    def __init__(self):
+        super().SetTitle('周方向貫通亀裂，ASME Code Case N-513の解')
+        super().SetRefer('ASME Boiler and Pressure Vessel Code, Code Cae N-513,Evaluation Criteria for Temporary Acceptance of Flaws in Class 3 Piping,1997')
+    def Calc(self):
+        df=super().GetData()
+        R=df['R']
+        c=df['c']
+        P=df['P']
+        M=df['M']
+        t=df['t']
+        Sm=P/(2*np.pi*R*t)
+        Sbg=M/(np.pi*R*R*t)
+        evaluate_cubic = lambda x, c: c[0] + c[1]*x + c[2]*x**2 + c[3]*x**3
+        x=R/t
+        coeffs=[-2.02917,1.67763,-0.07987,0.00176]
+        Am=evaluate_cubic(x,coeffs)
+        coeffs=[7.09987,-4.42394,0.21036,-0.00463]
+        Bm=evaluate_cubic(x,coeffs)
+        coeffs=[7.79661,5.16676,-0.24577,0.00541]
+        Cm=evaluate_cubic(x,coeffs)
+        coeffs=[-3.26543,1.52784,-0.072698,0.0016011]
+        Abg=evaluate_cubic(x,coeffs)
+        coeffs=[11.36322,-3.91412,0.18619,-0.004099]
+        Bbg=evaluate_cubic(x,coeffs)
+        coeffs=[-3.18609,3.84763,-0.18304,0.00403]
+        Cbg=evaluate_cubic(x,coeffs)
+        evaluate_F = lambda x, c: c[0] + c[1]*x**1.5 + c[2]*x**2.5 + c[3]*x**3.5
+        x=c/(np.pi*R)
+        coeffs=[1,Am,Bm,Cm]
+        Fm=evaluate_F(x,coeffs)
+        coeffs=[1,Abg,Bbg,Cbg]
+        Fbg=evaluate_F(x,coeffs)
+        K=(Fm*Sm+Fbg*Sbg)*np.sqrt(np.pi*c)
+        res={'K':K
+            }
+        super().SetRes(res)
+
+
+class L_1_a(Base):
+    def __init__(self):
+        super().SetTitle('平板の亀裂，半楕円表面亀裂')
+        super().SetRefer('Dillstrom,P.andSattari-Far,I.:Limit Load Solutions for Surface Cracks in Plates and Cylinders, RSE R & D Report,No.2002/01,Det Norske Veritas AB 2002')
+    def Calc(self):
+        df=super().GetData()
+        w=df['b']
+        a=df['a']
+        t=df['t']
+        P=df['P']
+        M=df['M']
+        l=df['c']*2
+        Sy=df['Sy']
+        Sm=P/(2*w*t)
+        Sb=3*M/(w*t*t)
+        z=a*l/(t*(l+2*t))
+        Lr=((1-z)**1.58*Sb/3+np.sqrt((1-z)**3.16*Sb*Sb/9+(1-z)**3.14*Sm*Sm))/((1-z)**2*Sy)
+        res={'Lr':Lr}
+        super().SetRes(res)
+
+        
+        
+class L_2_b(Base):
+    def __init__(self):
+        super().SetTitle('軸方向内表面長い表面亀裂')
+        super().SetRefer('Kumar V.,German M.D. and Shih C.F.:EPRI NP-1931,Electric Power Research Institute, Palo Alto,CA,July 1981')
+    def Calc(self):
+        df=super().GetData()
+        a=df['a']
+        t=df['t']
+        Sy=df['Sy']
+        Ri=df['Ri'] 
+        p=df['p']
+        z=a/t
+        p0=(2/np.sqrt(3))*Sy*(t/Ri)
+        pc=p0*((1-z)/(1+z/(Ri/t)))
+        Lr=p/pc
+        res={'Lr':Lr,
+             'pc':pc}
+        super().SetRes(res)
+
+
+
+