From 2aa8abc1fb9fa094a81fd545ec272ccba66d9c40 Mon Sep 17 00:00:00 2001
From: Sasha <sasha@noreply.localhost>
Date: Thu, 3 Jul 2025 19:05:38 +0100
Subject: [PATCH] fixed indeces of stringers

---
 aseredesign/aseredesign.py | 21 +++++++++++++++------
 1 file changed, 15 insertions(+), 6 deletions(-)

diff --git a/aseredesign/aseredesign.py b/aseredesign/aseredesign.py
index 8816d25..98734c3 100644
--- a/aseredesign/aseredesign.py
+++ b/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))
-    if(GetLambda(stringerNo)<sLambda_crit):
-        return(GetSigCrip(stringerNo, Loadcase) - ((1/Ebasis)*((GetSigCrip(stringerNo, Loadcase)/(2*math.pi))**2)*(GetLambda(stringerNo)**2)))
+    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)))
     else:
-        return((math.pi**2)*Ebasis / (GetLambda(stringerNo) ** 2))
+        return((math.pi**2)*Ebasis / (hLambda ** 2))
 
 
 def GetRFstringer(stringerNo, Loadcase):