1 d0=12;
2 L0=50;
3 A0=pi*(d0/2)^2;
4
5 Lf = 72.8;
6 df = 6.3;
7 Af = pi*(df/2)^2;
8 BreakingLoad = 27.2;
9
10 data = [
11 2 0.0041;
12 4 0.0082;
13 6 0.0132;
14 8 0.0183;
15 10 0.0226;
16 12 0.0267;
17 14 0.0310;
18 16 0.0351;
19 18 0.0391;
20 20 0.0445;
21 22 0.0485;
22 24 0.0518;
23 25.2 0.51;
24 28 1.52;
25 30 2.03;
26 34 3.05;
27 38.4 4.57;
28 40 6.60;
29 40.4 7.62;
30 40.8 12.7;
31 40.2 14.7;
32 38.6 15.7;
33 36.4 17.8;
34 32.4 19.3
35 ];
36
37
38 P_kN=data(:,1);
39 elong_mm=data(:,2);
40
41 P_N=P_kN*1000;
42 % Calculate the stress and strain
43 stress = P_N / A0; % Stress in Pascals
44 strain = elong_mm / L0; % Strain (dimensionless)
45
46 % Plot stress vs strain
47 figure;
48 plot(strain, stress, '-o');
49 xlabel('Strain (dimensionless)');
50 ylabel('Stress (Pa)');
51 title('Stress-Strain Curve');
52 grid on;
53
54
55 % young modulus
56 youngModulus = stress ./ strain; % Calculate Young's modulus
57
58 % modulus of elasticity
59 modulusOfElasticity = mean(youngModulus); % Calculate the average modulus of elasticity
60
61 % Proportional limit
62 proportionalLimit = max(stress(strain <= 0.01)); % Calculate the proportional limit for strain
<= 0.01
63
64 % offset yield stress
65 offsetYieldStress = proportionalLimit + (BreakingLoad / A0); % Calculate the offset yield stress
66
67 % tensile strength
68 tensileStrength = max(stress); % Calculate the tensile strength as the maximum stress
69
70 % fracture stress
71 fractureStress = max(stress(strain > 0.1)); % Calculate the fracture stress for strain > 0.1
72
73 % percentage elongation
74 percentageElongation = (elong_mm(end) - elong_mm(1)) / elong_mm(1) * 100; % Calculate percentage
75
76 % reduction of area
77 reductionOfArea = (A0 - Af) / A0 * 100; % Calculate the reduction of area as a percentage
78
79 % toughness
80 toughness = trapz(strain, stress); % Calculate toughness as the area under the stress-strain
curve
81
82 % resiliance
83 resilience = trapz(strain(stress < offsetYieldStress), stress(stress < offsetYieldStress)); %
Calculate resilience as the area under the elastic portion of the stress-strain curve
84
85 % true stress strain curve
86 % Calculate true stress and true strain
87 trueStrain = log(1 + strain); % True strain
88 trueStress = stress .* (1 + strain); % True stress
89 % Plot true stress vs true strain
90 figure;
91 plot(trueStrain, trueStress, '-o');
92 xlabel('True Strain (dimensionless)');\\\\\\\\\\\\\\\\\\\\\
93 ylabel('True Stress (Pa)');
94 title('True Stress-Strain Curve');
95 grid on;
96