Suraj Kumar
ed24s029
Question 4
disp('Suraj')
disp('ed24s029')
% Material properties: [Density (kg/m^3), Yield Strength (MPa), Elastic
Modulus (GPa)]
materials = {
'MB - Steel AISI 4130', 7850, 951, 205;
'MC - Steel A36', 7800, 250, 200;
'MD - Carbon fiber', 1750, 1600, 40;
};
% Convert table to array for calculations
numMaterials = size(materials, 1);
results = cell(numMaterials, 4); % Store [Material, Failure Index, Crack
Width, Remarks]
% Assume same sigma and crack size a
% We can define a normalized factor for comparison
for i = 1:numMaterials
name = materials{i, 1};
density = materials{i, 2};
sigma_y = materials{i, 3}; % in MPa
E = materials{i, 4} * 1e3; % Convert GPa to MPa for consistency
% Failure Index = sigma_applied / sigma_yield
% Use a dummy applied stress, say sigma = 500 MPa (same for all)
sigma_applied = 500; % MPa
failure_index = sigma_applied / sigma_y;
% Crack width estimate: delta ~ K^2 / (E * sigma_y)
% Assume K^2 is same for all, so delta ~ 1 / (E * sigma_y)
crack_width = 1 / (E * sigma_y);
% Remarks based on failure index
if failure_index > 1
remark = 'Likely to fail';
else
remark = 'Safe under load';
end
% Store results
results{i, 1} = name;
results{i, 2} = failure_index;
results{i, 3} = crack_width;
results{i, 4} = remark;
end
% Display results
disp('Material Analysis:');
fprintf('%-25s | Failure Index | Crack Width (normalized) | Remark\n',
'Material');
fprintf('%s\n', repmat('-', 1, 70));
for i = 1:numMaterials
fprintf('%-25s | %13.3f | %23.3e | %s\n', ...
results{i,1}, results{i,2}, results{i,3}, results{i,4});
end
Suraj
ed24s029
Material Analysis:
Material | Failure Index | Crack Width (normalized) |
Remark
----------------------------------------------------------------------
MB - Steel AISI 4130 | 0.526 | 5.129e-09 | Safe
under load
MC - Steel A36 | 2.000 | 2.000e-08 |
Likely to fail
MD - Carbon fiber | 0.312 | 1.562e-08 | Safe
under load