diff --git a/PartG.m b/PartG.m new file mode 100644 index 0000000..252162d --- /dev/null +++ b/PartG.m @@ -0,0 +1,22 @@ +No = linspace(.001,.01,10); +n = 20; + +T = zeros(1,length(No)); +wmax = zeros(1,length(No)); + +for i = 1:length(No) + T(i) = tension_sol(No(i),n); + w = membrane_solution(T(i),No(i),n); + wmax(i) = max(w); +end + +clf +x = wmax'; +y = No'; +Z=x.^3; +a=Z\y; +x_fcn=linspace(min(x),max(x)); +plot(x,y,'o',x_fcn,a*x_fcn.^3) +title('Pressure vs Maximum Deflection') +xlabel('Maximum Deflection (\muM)') +ylabel('Pressure (MPa)') \ No newline at end of file diff --git a/README.md b/README.md index ce17b4b..22f28dd 100644 --- a/README.md +++ b/README.md @@ -154,3 +154,30 @@ ea = zeros(1,length(n)); |30|0.0602|0.16%| |35|0.0602|0.09%| |40|0.0603|0.05%| + +## Part G +``` +No = linspace(.001,.01,10); +n = 20; + +T = zeros(1,length(No)); +wmax = zeros(1,length(No)); + +for i = 1:length(No) + T(i) = tension_sol(No(i),n); + w = membrane_solution(T(i),No(i),n); + wmax(i) = max(w); +end + +clf +x = wmax'; +y = No'; +Z=x.^3; +a=Z\y; +x_fcn=linspace(min(x),max(x)); +plot(x,y,'o',x_fcn,a*x_fcn.^3) +title('Pressure vs Maximum Deflection') +xlabel('Maximum Deflection (\muM)') +ylabel('Pressure (MPa)') +``` +![](assets/README-bbfc19b7.png) diff --git a/assets/PartG.png b/assets/PartG.png new file mode 100644 index 0000000..4b744b6 Binary files /dev/null and b/assets/PartG.png differ diff --git a/assets/README-bbfc19b7.png b/assets/README-bbfc19b7.png new file mode 100644 index 0000000..c48c03d Binary files /dev/null and b/assets/README-bbfc19b7.png differ