Matlab Code FEM Beam Example_9_3_2

Matlab Code FEM Beam Example_9_3_2 - Handouts - MME 412:...

Info iconThis preview shows pages 1–3. Sign up to view the full content.

View Full Document Right Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: Handouts - MME 412: Mechanics of Materials (Fall-2007) - 12 - MATLAB Simulation: Finite Element Formulation: Beam (Example 1) Clear the workspace clear all; close all; clc; Define the number of elements ne=2; Given Physical Parameters (Only valid for ne=2, if more elements the code needs to be modified) Length of element 1 L1=12 L1 = 12 Length of element 2 L2=24 L2 = 24 Total length of the beam L=L1+L2 L = 36 Moment of inertia of element 1 I1=0.2 I1 = 0.2000 Moment of inertia of element 2 I2=0.2 I2 = 0.2000 Modulus of elasticity E=30e6 E = 30000000 hh=[L1 L2]'; In=[I1 I2]'; Intensity of the distributed load on each element q=[-20 -10]' q = -20 -10 Intensity of the Point load acting on node 2 P=200 P = 200 Built Elemental stiffness matrix and load vectors for i=1:ne h=hh(i,1); Handouts - MME 412: Mechanics of Materials (Fall-2007) - 13 - q0=q(i,1); Ke(:,:,i)=E*In(i,1)/h^3*[12 6*h -12 6*h ; 6*h 4*h^2 -6*h 2*h^2; -12 -6*h 12 -6*h; 6*h 2*h^2 -6*h 4*h^2]; % Element Load vector Need to modify the load vector based...
View Full Document

This document was uploaded on 11/04/2011 for the course MME 512 at Miami University.

Page1 / 3

Matlab Code FEM Beam Example_9_3_2 - Handouts - MME 412:...

This preview shows document pages 1 - 3. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online