۳۱) ادامه‌ی تشریح بخش های یک کد بهینه سازی بخش ۳ و ۴: تابع هدف و قیود… — EQ SOFT — TG.ME

class FrameAnalyzer:
    def evaluate_design(self, section_indices: List[int]) -> Dict:

total_weight = 0.0
for idx in beam_indices:
    section = SECTION_LIST[idx]
    weight = section.A * length * 7850  # چگالی فولاد
    total_weight += weight

w_u = 1.2 * w_DL + 1.6 * w_LL

seismic_load = self.calculate_seismic_load(total_weight)  # V = Cs × W

# ظرفیت خمشی: Mn = Fy × Zx
phi_Mn = inputs.phi_b * inputs.Fy * section.Zx
bending_ratio = Mu / phi_Mn  # باید ≤ 1.0

# ظرفیت برشی: Vn = 0.6 × Fy × Aw
phi_Vn = inputs.phi_v * 0.6 * inputs.Fy * section.d * section.tw
shear_ratio = Vu / phi_Vn  # باید ≤ 1.0

# تنش کمانش بحرانی (Euler buckling)
lambda_c = (KL / section.ry) * sqrt(Fy / (π² × E))
if lambda_c <= 4.71 * sqrt(E/Fy):
    Fcr = (0.658^(Fy/Fe)) × Fy  # کمانش غیرالاستیک
else:
    Fcr = 0.877 × Fe            # کمانش الاستیک
Pn = Fcr × A
compression_ratio = axial_load / (phi_c × Pn)  # باید ≤ 1.0

if not constraints_ok:
    excess = sum(max(0, r - 1.0) for r in all_ratios)
    penalty = 1000 * excess
objective_value = total_weight + penalty