1

cpl11001 · Dec 14, 2017 · d2ac2d2 · d2ac2d2
1 parent f61bb53
commit d2ac2d2
Show file tree

Hide file tree

Showing 2 changed files with 29 additions and 0 deletions.
diff --git 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