diff --git a/README.md b/README.md index 7363dee..87f3abc 100644 --- a/README.md +++ b/README.md @@ -186,6 +186,35 @@ For a given pressure range of 0.001 to 0.01 MPa with 10 steps, determine the mem ### Approach Using the given pressure range, the tension was found at each different pressure value using the bisect method. The best fit line was found +```matlab +clear +P = linspace(0.001,0.01,10); + +for i = 1:length(P) + func = @(T) SE_diff(T,P(i),10); + [root(i),fx,ea,iter] = bisect_final_project(func,0.001,1,.1); +end + +for i = 1:length(root) + w = membrane_solution(root(i),P(i),10); + w1(:,i) = w; +end + +w_max = max(w1); +coefficients = polyfit(w_max,P,3); +Y = polyval(coefficients,w_max); + +plot(w_max,P,w_max,Y,'or') +xlabel('Max Deflection (micron)') +ylabel('Pressure (MPa)') +title('Pressure vs. Maximum Deflection') +``` + +Output - Based on 10 interior nodes: +- A = 0.4443 + +![Pressure vs. Max Deflection](./figures/figure3.png) + ## Part H ### Problem Statement diff --git a/figures/figure3.png b/figures/figure3.png new file mode 100644 index 0000000..4b8b2df Binary files /dev/null and b/figures/figure3.png differ