.

Part F/F_sol.asv

@@ -0,0 +1,6 @@
+ n=20:5:40;
+ P=0.001; %MPa
+ T = zeros(1,length(n));
+ for i = 1:length(n)
+     [T(i) ea(i)] = tension_sol(P,n(i));
+ end

Part F/F_sol.m

@@ -0,0 +1,7 @@
+ n=20:5:40;
+ P=0.001; %MPa
+ T = zeros(1,length(n));
+ ea = zeros(1,length(n));
+ for i = 1:length(n)
+     [T(i), ea(i)] = tension_sol(P,n(i));
+ end

Part F/SE_diff.m

@@ -0,0 +1,26 @@
+function [pw_se,w]=SE_diff(T,P,n)
+E = 1000000; %TPa Units may need to be changed
+v = .31; %Poissons ratio
+t = .0003; %um
+h = 10/(n+1); %nm
+w = membrane_solution(T,P,n);
+z = zeros(n+2);
+z(2:end-1,2:end-1) = reshape(w,[n n]);
+num = n + 1;
+wbar = zeros(num);
+for i = 1:num
+    for j = 1:num
+        wbar(i,j) = mean([z(i,j),z(i+1,j),z(i,j+1),z(i+1,j+1)]);
+    end
+end
+pw = sum(sum(wbar.*h^2.*P));
+dwdx = zeros(num);
+dwdy = zeros(num);
+for i = 1:num
+    for j = 1:num
+        dwdx(i,j) = mean([z(i+1,j)-z(i,j),z(i+1,j+1)-z(i,j+1)]);
+        dwdy(i,j) = mean([z(i,j+1)-z(i,j),z(i+1,j+1)-z(i+1,j)]);
+    end
+end
+se = E*t*h^2/(2*(1-v^2))*sum(sum(0.25.*dwdx.^4+.25.*dwdy.^4+0.5.*(dwdx.*dwdy).^2));
+pw_se = abs(pw-se);

Part F/tension_sol.asv

@@ -1,7 +1,6 @@
 function T = tension_sol(P,n)
-num = length(n)
-for i = 1:nu
-od = ones(n^2-1,1);
-od(n:n:end) = 0;
-k = -4*diag(ones(n^2,1))+diag(ones((n^2)-n,1),n)+diag(ones((n^2)-n,1),-n)+diag(od,1)+diag(od,-1);
-y = -(10/(n+1))^2*(P/T)*ones(n^2,1);
+num = length(n);
+for i = 1:num
+    y =@(T) SE_diff(
+    [root,fx,ea,iter]=bisect(y,0,.006)
+end

Part F/tension_sol.m

@@ -1,6 +1,3 @@
-function T = tension_sol(P,n)
-num = length(n);
-for i = 1:num
-    y =@(T) (P./T)+3.6534e+04;
-    [root,fx,ea,iter]=bisect(y,0,.006)
-end
+function [T,ea] = tension_sol(P,n)
+y =@(T) SE_diff(T,P,n);
+[T,fx,ea,iter]=bisect(y,.01,1);

Part G/part_g.asv

@@ -0,0 +1,5 @@
+P = linspace(.001,.01,10);
+n = 20;
+for i = 1:length(P)
+    T(i) = tension_sol(P(i),n);
+

Part G/part_g.m

@@ -0,0 +1,10 @@
+P = linspace(.001,.01,10);
+n = 20;
+T = zeros(1,length(P));
+wmax = zeros(1,length(P));
+for i = 1:length(P)
+    T(i) = tension_sol(P(i),n);
+    w = membrane_solution(T(i),P(i),n);
+    wmax(i) = max(w);
+end
+plot(wmax,P,'o')