Concrete Mix Design Calculator

def concrete_mix_design(grade, wcr, slump=40, agg_size=20, fm=2.8):
    """
    Calculate concrete mix proportions using IS 10262 method
    
    grade: Concrete grade (MPa) - 25, 30, 35, 40
    wcr: Water-cement ratio (0.4 to 0.8)
    slump: Slump value (mm) - 20, 40, 60
    agg_size: Nominal aggregate size (mm) - 20, 40
    fm: Fineness modulus of sand (2.4 to 3.0)
    
    Returns:
        Dictionary with cement, water, fine sand, coarse aggregate (kg/m³)
    """
    
    water_req = {
        30: {'20mm': 195, '40mm': 160, '60mm': 140},
        35: {'20mm': 205, '40mm': 170, '60mm': 150},
        40: {'20mm': 215, '40mm': 180, '60mm': 160},
    }
    
    closest_slump = 60 if slump > 40 else (40 if slump > 20 else 20)
    slump_key = f'{closest_slump}mm'
    
    water_content = water_req[grade][slump_key]
    cement_content = water_content / wcr
    
    # Simplified aggregate calculation
    fine_agg_pct = 30 + (2.8 - fm) * 2
    coarse_agg_pct = 100 - fine_agg_pct - 2
    
    fine_sand = (fine_agg_pct / 100) * 2650 * 0.5
    coarse_agg = (coarse_agg_pct / 100) * 2750 * 0.5
    
    return {
        'Cement': round(cement_content),
        'Water': round(water_content),
        'Fine Sand': round(fine_sand),
        'Coarse Aggregate': round(coarse_agg),
        'W/C Ratio': wcr
    }

# Example
grade = 30
wcr = 0.55
slump = 40

result = concrete_mix_design(grade, wcr, slump)
print(f"Concrete Mix Design (M{grade})")
for key, value in result.items():
    print(f"{key}: {value}")

function mix_design = concrete_mix_design(grade, wcr, varargin)
    % Concrete Mix Design Calculator (IS 10262 Method)
    % mix_design = concrete_mix_design(grade, wcr)
    % mix_design = concrete_mix_design(grade, wcr, slump, agg_size, fm)
    
    % Default parameters
    slump = 40;
    agg_size = 20;
    fm = 2.8;
    
    if nargin > 2
        slump = varargin{1};
    end
    if nargin > 3
        agg_size = varargin{2};
    end
    if nargin > 4
        fm = varargin{3};
    end
    
    % Water requirement table
    water_data = containers.Map(...
        [25, 30, 35, 40], ...
        {[186, 150, 130], [195, 160, 140], [205, 170, 150], [215, 180, 160]});
    
    % Find water content
    slump_idx = (slump > 40) * 3 + (slump <= 40 && slump > 20) * 2 + 1;
    water_content = water_data(grade);
    water_content = water_content(slump_idx);
    
    % Calculate cement content
    cement_content = water_content / wcr;
    
    % Aggregate calculations
    fine_agg_pct = 30 + (2.8 - fm) * 2;
    coarse_agg_pct = 100 - fine_agg_pct - 2;
    
    % Material quantities (kg/m³)
    fine_sand = (fine_agg_pct / 100) * 2650 * 0.5;
    coarse_agg = (coarse_agg_pct / 100) * 2750 * 0.5;
    
    % Return structure
    mix_design.cement = round(cement_content);
    mix_design.water = round(water_content);
    mix_design.fine_sand = round(fine_sand);
    mix_design.coarse_agg = round(coarse_agg);
    mix_design.wcr = wcr;
end

% Example usage
mix = concrete_mix_design(30, 0.55, 40);
fprintf('Concrete Mix Design (M%d)\n', 30);
fprintf('Cement: %d kg/m³\n', mix.cement);
fprintf('Water: %d kg/m³\n', mix.water);
fprintf('Fine Sand: %d kg/m³\n', mix.fine_sand);
fprintf('Coarse Aggregate: %d kg/m³\n', mix.coarse_agg);