Slab Load Calculator

def slab_load(thickness_mm, finishes_kN_m2, live_load_kN_m2, density_kN_m3=25.0, trib_width_m=0.0, standard="ASCE_ACI"):
    """
    Calculate area and tributary beam line loads for RC slabs.
    """
    if thickness_mm <= 0 or density_kN_m3 <= 0 or finishes_kN_m2 < 0 or live_load_kN_m2 < 0:
        raise ValueError("Invalid negative or zero parameters.")

    # Self-weight
    self_weight = (thickness_mm / 1000.0) * density_kN_m3
    dead_load = self_weight + finishes_kN_m2
    
    # Factored Ultimate Load (ULS)
    if standard == "EUROCODE":
        factored_load = 1.35 * dead_load + 1.50 * live_load_kN_m2
    else: # ASCE 7 / ACI 318
        factored_load = 1.20 * dead_load + 1.60 * live_load_kN_m2

    linear_factored_load = factored_load * trib_width_m if trib_width_m > 0 else 0.0

    return {
        "self_weight_kN_m2": round(self_weight, 2),
        "total_dead_load_kN_m2": round(dead_load, 2),
        "factored_area_load_kN_m2": round(factored_load, 2),
        "factored_beam_line_load_kN_m": round(linear_factored_load, 2)
    }

# Example
res = slab_load(
    thickness_mm=180,
    finishes_kN_m2=1.5,
    live_load_kN_m2=2.0,
    trib_width_m=3.5,
    standard="ASCE_ACI"
)

print(f"Factored Area Load: {res['factored_area_load_kN_m2']} kN/m²")
print(f"Beam Line Load (ULS): {res['factored_beam_line_load_kN_m']} kN/m")

function [q_factored, w_beam] = slab_load(t_mm, finishes, live_load, density, trib_w, is_eurocode)
    % t_mm: slab thickness (mm), finishes & live_load: kN/m^2
    % density: kN/m^3 (default 25), trib_w: tributary width (m)
    if nargin < 4, density = 25; end
    if nargin < 5, trib_w = 0; end
    if nargin < 6, is_eurocode = false; end

    sw = (t_mm / 1000) * density;
    dl = sw + finishes;

    if is_eurocode
        q_factored = 1.35 * dl + 1.50 * live_load;
    else
        q_factored = 1.20 * dl + 1.60 * live_load;
    end

    w_beam = q_factored * trib_w;
end

% Example
[q_ult, w_line] = slab_load(180, 1.5, 2.0, 25, 3.5, false);
fprintf('Factored Slab Area Load: %.2f kN/m^2\n', q_ult);
fprintf('Factored Line Load on Beam: %.2f kN/m\n', w_line);