fixed indeces of stringers

aseredesign/aseredesign.py

@@ -147,19 +147,28 @@ def RF_panel(panelNo, Loadcase):
     return(result)
 
 
-def GetLambda(stringerNo):
+def GetLambda(i): # checked and correct
     I_comb = GetIcomb(i)
     r_gyr = math.sqrt(I_comb/(areaStringer + 0.5 * pThickness[i] * pWidth + 0.5 * pThickness[i+1] * pWidth))
     c = 1
-    return(1*pLength/r_gyr)
+    sLambda = c*pLength/r_gyr
+    return(sLambda)
 
 
 def GetSigCRIT(stringerNo, Loadcase):
-    sLambda_crit = math.sqrt((2*math.pi*math.pi*Ebasis)/GetSigCrip(i,0))
+    sLambda_crit = math.sqrt((2*math.pi*math.pi*Ebasis)/GetSigCrip(stringerNo,Loadcase))
-    if(GetLambda(stringerNo)<sLambda_crit):
+    # Lambda_crit correct
-        return(GetSigCrip(stringerNo, Loadcase) - ((1/Ebasis)*((GetSigCrip(stringerNo, Loadcase)/(2*math.pi))**2)*(GetLambda(stringerNo)**2)))
+    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)))
     else:
-        return((math.pi**2)*Ebasis / (GetLambda(stringerNo) ** 2))
+        return((math.pi**2)*Ebasis / (hLambda ** 2))
 
 
 def GetRFstringer(stringerNo, Loadcase):