MATLAB Script for residual stress modeling validation

clc;
clear all;

Average Normal Stress Before Release

1️⃣10x10x10

x10b=0:2.54:25.4;
y10b=[-116.675 37.3138 53.9142 -104.139 -198.5 -223.202 -185.361... -109.179 4.73441 189.466 302.992];
area10b=trapz(x10b,y10b);
avg10b=area10b/25.4;

2️⃣20x20x20

x20b=0:1.27:25.4;
y20b=[-123.046 153.332 286.696 142.477 59.0121 -64.4223 -200.745... -305.516 -359.947 -367.224 -340.813 -294.851 -239.42 -179.41 -116.162... -48.4717 26.4926 114.329 229.632 382.93 468.404];
area20b=trapz(x20b,y20b);
avg20b=area20b/25.4;

3️⃣40x40x40

x40b=0:.6350:25.4;
y40b=[-109.872 327.533 675.345 472.736 359.383 328.183 314.33 267.602... 165.266 20.0981 -145.266 -320.089 -491.852 -561.195 -558.552 -587.847... -609.893 -607.23 -585.05 -547.869 -500.893 -447.508 -389.736 -330.545... -272.239 -216.58 -164.434 -115.151 -67.076 -20.9095 23.8104 67.1678... 111.216 158.702 214.452 236.996 295.826 376.009 478.65 604.114 669.595];
area40b=trapz(x40b,y40b);
avg40b=area40b/25.4;

Average Normal Stress After Release

y??=[]（第二行）和before release不一样

1️⃣10x10x10

x10=0:2.54:25.4;
y10=[76.7958 185.584 137.523 -63.2302 -183.398 -221.85 -1