Reference no: EM133350771
Computational Method for Aerostructures
Problem Statement
You are required to create a programming algorithm to calculate the necessary displacement of nodes for one-dimensional elements. Your program should consist of the following functions that are written as a separate function file:
Table 1: Functions for one-dimensional elements
Function file name
|
Description
|
Input
|
Output
|
Gauss_Elim.m
|
To perform gauss elimination on
a square matrix and output the solution
|
Assuming [A][x] = [b], inputs are [A] and [b]
|
[x]
|
Bar_Elem.m
|
To determine the stiffness matrix and forcing matrix by considering one-dimensional bar elements and calculate the nodal displacement of each
node
|
A, E, nodal matrix, element connectivity matrix
|
Nodal displacements
|
Truss_Elem.m
|
To determine the stiffness matrix and forcing matrix by considering one-dimensional truss elements and calculate the nodal displacement of each
node
|
A, E, nodal matrix, element connectivity matrix
|
Nodal displacements
|
Beam_Elem.m
|
To determine the stiffness matrix and forcing matrix by considering one-dimensional beam elements and calculate the nodal displacement of each
node
|
E, I, nodal matrix, element connectivity matrix
|
Nodal displacements / angle
|
Your Task
You will need to create the functions listed in Table 1. You will need to create the Gauss_Elim.m function first, which solves matrices using the Gauss Elimination method. You will need to call this function every time you are required to solve a matrix. DO NOT use the in-built functions to solve matrices.
For the other three functions, based on the given inputs, your functions should:
1. Be able to accept the inputs listed in Table 1 for any size matrix
2. Determine the stiffness and forcing matrix of each element
3. Form the global stiffness and forcing matrix
4. Based on the boundary conditions, reduce the matrix by removing the rows and columns corresponding to a solution of zero displacements or zero angles.
5. Calculate the nodal displacements or angles and the corresponding deformed coordinates for each node
6. Output the nodal displacements and plot the deformed structure based on the determined deformed coordinates
Note that your program should also have an error detection algorithm to detect incorrect/invalid inputs. For instance, it should be able to detect if a node is missing.
Try to complete the program for one type of element first. For the other types, the code will be somewhat similar. You just have to make a small tweak to it only.
Attachment:- Computational Method for Aerostructures.rar