Калкулатор за монтаж на амортисьор за врата
==================================================
Минимална дължина на амортисьора (L_min): 180 mm
Максимална дължина на амортисьора (L_max): 260 mm
Маса на вратата: 20 kg ( не знам колко е точно ... просто кг )
Височина на вратата: 500 mm

РЕЗУЛТАТИ ЗА МОНТАЖ:
------------------------------
Разстояние от пантите до точка на вратата (PA): 70.3 mm
Разстояние от пантите до точка на рамката (PB): 250.3 mm

ПРОВЕРКА НА ГЕОМЕТРИЯТА:
------------------------------
Затворена врата: 180.0 mm (очаквано: 180 mm)
Отворена врата: 260.0 mm (очаквано: 260 mm)

ИЗЧИСЛЯВАНЕ НА СИЛАТА
----------------------------------------
Момент от тежестта:        49.05 Nm
Сила при 90° (без резерв): 697.6 N
Препоръчителна сила:       907.0 N
Изберете амортисьор:       900–950 N

Момент от тежестта:        19.62 Nm
Сила при 90° (без резерв): 279.0 N
Препоръчителна сила:       363.0 N
Изберете амортисьор:       350–400 N

import math

# Входни параметри
L_min = 180  # mm
L_max = 260  # mm
door_mass = 8  # kg
door_height = 500  # mm

def calculate_required_force(
    door_mass_kg: float,
    R_mm: float, # разстояние от пантите до точката на вратата (PA)
    center_of_gravity_mm: float = None, # ако не се подаде – приема се door_height/2
    door_height_mm: float = None,
    safety_factor: float = 1.3, # 30% резерв – стандартна практика
    gravity: float = 9.81
) -> dict:
    """
    Изчислява необходимата сила на газовия амортисьор (в Нютони)
    за врата, която се отваря хоризонтално на 90°.

    Parameters:
        door_mass_kg         – маса на вратата в kg
        R_mm                 – разстояние от оста на пантите до точката на вратата (mm)
        center_of_gravity_mm – разстояние от пантите до центъра на тежестта (mm)
        door_height_mm       – използва се ако не е подаден center_of_gravity
        safety_factor        – коефициент на безопасност (обикн. 1.2–1.5)
        gravity              – 9.81 m/s²

    Returns:
        dict с резултати
    """
    if center_of_gravity_mm is None:
        if door_height_mm is None:
            raise ValueError("Трябва да се подаде или center_of_gravity_mm, или door_height_mm!")
        center_of_gravity_mm = door_height_mm / 2

    # Връщане в метри за правилни единици
    cg_m = center_of_gravity_mm / 1000
    R_m  = R_mm / 1000

    # Момент на силата от тежестта на вратата (Nm)
    torque_Nm = door_mass_kg * gravity * cg_m

    # При 90° отворена врата – ръката на силата на амортисьора е точно R
    # Силата е перпендикулярна на вратата → F * R = Torque
    force_at_90_deg_N = torque_Nm / R_m

    # С реален резерв (препоръчително 20–50%)
    recommended_force_N = force_at_90_deg_N * safety_factor

    # Най-близки стандартни стойности (обикн. на стъпки по 50 N)
    step = 50
    lower  = (recommended_force_N // step) * step
    upper  = lower + step
    if recommended_force_N > lower + step / 2:
        lower = upper

    return {
        "torque_Nm":                round(torque_Nm, 2),
        "force_at_90_deg_N":        round(force_at_90_deg_N, 1),
        "recommended_force_N":      round(recommended_force_N, 0),
        "suggested_spring_range":   f"{int(lower)}–{int(lower + 50)} N",
        "R_mm":                     R_mm,
        "center_of_gravity_mm":     center_of_gravity_mm,
        "safety_factor":            safety_factor
    }

def calculate_gas_spring_mounting_points():
    """
    Изчислява точките на монтаж за амортисьор на врата.
    """
    print("Калкулатор за монтаж на амортисьор за врата")
    print("=" * 50)

    print(f"Минимална дължина на амортисьора (L_min): {L_min} mm")
    print(f"Максимална дължина на амортисьора (L_max): {L_max} mm")
    print(f"Маса на вратата: {door_mass} kg")
    print(f"Височина на вратата: {door_height} mm")
    print()

    # Проверка за валидност на входните данни
    if L_min >= L_max:
        print("ГРЕШКА: Минималната дължина трябва да е по-малка от максималната!")
        return

    # Решаване на системата уравнения
    a = 1
    b = L_min
    c = (L_min**2 - L_max**2) / 2

    discriminant = b**2 - 4*a*c

    if discriminant < 0:
        print("ГРЕШКА: Няма реални решения за дадените параметри!")
        return

    R1 = (-b + math.sqrt(discriminant)) / (2*a)
    R2 = (-b - math.sqrt(discriminant)) / (2*a)
    R = max(R1, R2)
    D = R + L_min

    print("РЕЗУЛТАТИ ЗА МОНТАЖ:")
    print("-" * 30)
    print(f"Разстояние от пантите до точка на вратата (PA): {R:.1f} mm")
    print(f"Разстояние от пантите до точка на рамката (PB): {D:.1f} mm")
    print()

    # Проверка на резултатите
    print("ПРОВЕРКА НА ГЕОМЕТРИЯТА:")
    print("-" * 30)
    calculated_L_min = D - R
    calculated_L_max = math.sqrt(R**2 + D**2)
    print(f"Затворена врата: {calculated_L_min:.1f} mm (очаквано: {L_min} mm)")
    print(f"Отворена врата: {calculated_L_max:.1f} mm (очаквано: {L_max} mm)")

    return R

if __name__ == "__main__":
    R = calculate_gas_spring_mounting_points()

    result = calculate_required_force(
        door_mass_kg=door_mass,
        R_mm=R,
        door_height_mm=door_height
    )

    print("\nИЗЧИСЛЯВАНЕ НА СИЛАТА")
    print("-" * 40)
    print(f"Момент от тежестта:        {result['torque_Nm']} Nm")
    print(f"Сила при 90° (без резерв): {result['force_at_90_deg_N']} N")
    print(f"Препоръчителна сила:       {result['recommended_force_N']} N")
    print(f"Изберете амортисьор:       {result['suggested_spring_range']}")