Assignment 1: Finite Element Analysis Essay example

2560 Words Jan 24th, 2016 11 Pages
[pic]

ENGD3016 Solid Mechanics
Assignment 1: Finite Element Analysis

Name: Wei Zhang ID: P14021978 Date: Dec17th 2015
Abstract
1.0 Introduction
2.0 Objectives
3.0 Matlab
4.0 Solidworks 4.1 Model of truss 1 4.2 Model of truss 2
5.0 Comparison of the two trusses
6.0 Comparison between MATLAB and SOLIDWORKS 6.1 Comparison of results 6.2 Advantages and disadvantages
7.0 Conclusion
Appendix
Abstract

The purpose of this report is to based on two different 2D pictures of trusses to finite element analysis from two ways. First carried
…show more content…
MATLAB coding
The information has been accessed
E = 2.05E11
A=0.000144 m²
The length of each bar
L1=0.6m
L2=0.8m
L3=1m
Deal with bar① angle1=input('Enter the value of angle1(degree):' ) s1=sind(angle1); c1=cosd(angle1);

Here need to input the angle of the bar ① and X
K1= E*A/L1 *[c1^2 c1*s1 -c1^2 -c1*s1 0 0;c1*s1 s1^2 -c1*s1 -s1^2 0 0;-c1^2 -c1*s1 c1^2 c1*s1 0 0;-c1*s1 -s1^2 c1*s1 s1^2 0 0;0 0 0 0 0 0;0 0 0 0 0 0]

K1 is a matrix of 6 * 6
This change is based on the formula
[pic]
The same code the bar ②③ angle2=input('Enter the value of angle2(degree):' ) s2=sind(angle2); c2=cosd(angle2);
K2=E*A/L2 *[c2^2 c2*s2 0 0 -c2^2 -c2*s2;c2*s2 s2^2 0 0 -c2*s2 -s2^2;0 0 0 0 0 0;0 0 0 0 0 0;-c2^2 -c2*s2 0 0 c2^2 c2*s2;-c2*s2 -s2^2 0 0 c2*s2 s2^2]

angle3=input('Enter the value of angle3(degree):' ) s3=sind(angle3); c3=cosd(angle3);
K3=E*A/L3 *[0 0 0 0 0 0;0 0 0 0 0 0;0 0 c3^2 c3*s3 -c3^2 -c3*s3;0 0 c3*s3 s3^2 -c3*s3 -s3^2;0 0 -c3^2 -c3*s3 c3^2 c3*s3;0 0 -c3*s3 -s3^2 c3*s3 s3^2]
Give an overall stiffness matrix
K=K1+K2+K3
K =

1.0e+07 *

3.6900 0 0 0 -3.6900 0 0 4.9200 0 -4.9200 0 0 0 0 1.8828 -1.4188 -1.8828 1.4188 0 -4.9200 -1.4188 5.9892 1.4188 -1.0692 -3.6900 0 -1.8828 1.4188 5.5728 -1.4188 0 0 1.4188 -1.0692 -1.4188 1.0692

Related Documents