Spur Gear Mesh Calculator

import math

def calculate_gear_mesh(m, z1, z2, alpha_deg=20):
    """
    Calculate properties of a meshing spur gear pair.
    
    m: Module (mm)
    z1: Pinion teeth
    z2: Gear teeth
    alpha_deg: Pressure angle (degrees)
    """
    if m <= 0 or z1 < 1 or z2 < 1:
        raise ValueError("Invalid gear inputs")

    d1 = m * z1
    d2 = m * z2
    a = (d1 + d2) / 2
    i = z2 / z1
    
    alpha_rad = math.radians(alpha_deg)
    db1 = d1 * math.cos(alpha_rad)
    db2 = d2 * math.cos(alpha_rad)

    return {
        "center_distance": a,
        "ratio": i,
        "pinion_pitch": d1,
        "gear_pitch": d2,
        "pinion_base": db1,
        "gear_base": db2
    }

# Example
mesh = calculate_gear_mesh(3.0, 18, 54, 20)
print(f"Center Distance: {mesh['center_distance']:.2f} mm")
print(f"Pinion Base Dia: {mesh['pinion_base']:.2f} mm")

function mesh = spur_gear_mesh(m, z1, z2, alpha_deg)
    % Calculate properties of meshing spur gear pair
    
    if m <= 0 || z1 < 1 || z2 < 1
        error('Invalid gear inputs');
    end

    mesh.d1 = m * z1;
    mesh.d2 = m * z2;
    mesh.a = (mesh.d1 + mesh.d2) / 2;
    mesh.ratio = z2 / z1;
    
    alpha_rad = alpha_deg * (pi / 180);
    mesh.db1 = mesh.d1 * cos(alpha_rad);
    mesh.db2 = mesh.d2 * cos(alpha_rad);
end

% Example
mesh = spur_gear_mesh(3.0, 18, 54, 20);
fprintf('Center Distance: %.2f mm\n', mesh.a);
fprintf('Pinion Base Dia: %.2f mm\n', mesh.db1);