aseredesign/aseredesign.py

@@ -147,28 +147,19 @@ def RF_panel(panelNo, Loadcase):
     return(result)
 
 
-def GetLambda(i): # checked and correct
+def GetLambda(stringerNo):
     I_comb = GetIcomb(i)
     r_gyr = math.sqrt(I_comb/(areaStringer + 0.5 * pThickness[i] * pWidth + 0.5 * pThickness[i+1] * pWidth))
     c = 1
-    sLambda = c*pLength/r_gyr
-    return(sLambda)
+    return(1*pLength/r_gyr)
 
 
 def GetSigCRIT(stringerNo, Loadcase):
-    sLambda_crit = math.sqrt((2*math.pi*math.pi*Ebasis)/GetSigCrip(stringerNo,Loadcase))
-    # Lambda_crit correct
-    hSigCrip = GetSigCrip(stringerNo, Loadcase)
-    # SigCrip  correct!!
-    hLambda = GetLambda(stringerNo-1)
-    #print("Lmabda: " + str(hLambda))
-    # ISSUE HERE - Lambda off from results
-    # no double checking now
-        
-    if(hLambda<sLambda_crit):
-        return(hSigCrip - ((1/Ebasis)*((hSigCrip/(2*math.pi))**2)*(hLambda**2)))
+    sLambda_crit = math.sqrt((2*math.pi*math.pi*Ebasis)/GetSigCrip(i,0))
+    if(GetLambda(stringerNo)<sLambda_crit):
+        return(GetSigCrip(stringerNo, Loadcase) - ((1/Ebasis)*((GetSigCrip(stringerNo, Loadcase)/(2*math.pi))**2)*(GetLambda(stringerNo)**2)))
     else:
-        return((math.pi**2)*Ebasis / (hLambda ** 2))
+        return((math.pi**2)*Ebasis / (GetLambda(stringerNo) ** 2))
 
 
 def GetRFstringer(stringerNo, Loadcase):