Skip to content

Commit

Permalink
Browse files Browse the repository at this point in the history
Update README.md
  • Loading branch information
ltd13002 committed Dec 15, 2017
1 parent fa1b75c commit 9fec7f1
Showing 1 changed file with 8 additions and 0 deletions.
8 changes: 8 additions & 0 deletions README.md
Expand Up @@ -113,6 +113,8 @@ Use membranesolution function to solve for w and plot the result with surf(X,Y,W


# Part E
### Problem
Create a function that calculates the difference in strain energy and work done by pressure for n-by-n elements '[pw_se,w]=SE_diff(T,P,n)'.
```matlab
function [pw_se,w]=SE_diff(T,P,n)
% SE_diff: calculates difference between strain energy and work done by pressure in
Expand Down Expand Up @@ -161,6 +163,8 @@ pw_se = pw-se;
- Using these values, the strain energy can be calculated

# Part F
### Problem
Use a root-finding method to calculate the tension in the membrane given a pressure, P=0.001 MPa, and n=[20:5:40] interior nodes. Show the decrease in tension error in a table.
```matlab
% Part F script
% From Problem Statement:
Expand Down Expand Up @@ -257,6 +261,8 @@ end
- 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
### Problem
Plot the pressure vs maximum deflection for pressure range 0.001-0.01 MPa with 10 steps using a root-finding method to determine tension, T, at each pressure.
```matlab
% Part G Script
% From Problem Statement
Expand Down Expand Up @@ -301,6 +307,8 @@ print('Part g','-dpng')
- The results are plotted

# Part H
### Problem
Show that the constant A is converging as the number of nodes is increased.
```matlab
% Part G Script
% From Problem Statement
Expand Down

0 comments on commit 9fec7f1

Please sign in to comment.