From 3217bec2baefdb9797b60fd7e29f12b838227cd0 Mon Sep 17 00:00:00 2001 From: Matthew Maliniak Date: Thu, 14 Dec 2017 20:13:47 -0500 Subject: [PATCH] Update README.md --- README.md | 16 +++++++++++++++- 1 file changed, 15 insertions(+), 1 deletion(-) diff --git a/README.md b/README.md index 822ba54..0c9f8c6 100644 --- a/README.md +++ b/README.md @@ -125,6 +125,12 @@ 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 = pw-se; ``` +### Approach +- Using the membrane_solution function, a vector of the displacements is formed +- Next, the average displacement for each element is calculated. For each elements, the dispalcement at all four courners is taken and then averaged +- Using these values, the work done by pressure can be calculated +- For the change in dispalcement over the x and y coordinate system, the values of the change on the left and right (y-axis) or top and bottom (x-axis) are taken and averaged +- Using these values, the strain enegery can be calculated # Part F ```matlab @@ -214,7 +220,9 @@ end |35|0.0602|0.09%| |40|0.0603|0.06%| - +### Approach +- This problem uses the bisect method for locating roots and the SE_diff function for calculating the difference in work and strain energy of the membrane +- The script runs through all iterations of different amounts of nodes, zeroing the SE_diff function output and saving the values for tension in the T variable as a vector # Part G ```matlab P = linspace(.001,.01,10); @@ -240,3 +248,9 @@ ylabel('Pressure (MPa)') print('Part g','-dpng') ``` ![](https://github.uconn.edu/ltd13002/ME3255_Final_Project/blob/master/Part%20G/Part%20g.png) + +### Approach +- This script uses the tension_sol function as described in the part above, running though all iterations of pressure +- Additionally, the max deflection for each pressure and tension is calculated at each iteration +- From there, a general linear regression is calculated with the formula P(x) = A*w^3 +- The results are plotted