Scribd
Upload
35 views1 page

06/10/16 14:02 C:/Users/g1020522/Documents/MATLAB/LESM 10-06/elem.m 1 of 1

This document contains MATLAB code functions for calculating member axial and transversal fixed end forces (FEF) given different load cases. It includes functions to calculate the axial FEF for a linearly distributed load, axial FEF for a temperature variation load, and transversal FEF for a linearly distributed load assuming simple support on both ends. Each function takes in relevant parameters, performs calculations based on structural mechanics equations, and returns the 2x1 FEF vector.

Uploaded by

Catarina Campelo
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
35 views1 page

06/10/16 14:02 C:/Users/g1020522/Documents/MATLAB/LESM 10-06/elem.m 1 of 1

This document contains MATLAB code functions for calculating member axial and transversal fixed end forces (FEF) given different load cases. It includes functions to calculate the axial FEF for a linearly distributed load, axial FEF for a temperature variation load, and transversal FEF for a linearly distributed load assuming simple support on both ends. Each function takes in relevant parameters, performs calculations based on structural mechanics equations, and returns the 2x1 FEF vector.

Uploaded by

Catarina Campelo
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 1

06/10/16 14:02 C:\Users\g1020522\Documents\MATLAB\LESM 10-06\Elem.

m 1 of 1

%-----------------------------------------------------------------------
% Generate member axial fixed end force vector for a linearly distributed
% load applied to given element.
% Output:
% fea: a 2x1 vector with axial fixed end force components
% Input arguments:
% qxi: axial load value at init node
% qxf: axial load value at end (final) node
function fea = axialLinearLoadFEF(obj,qxi,qxf)
L = obj.length; % member length

% *** COMPLETE HERE - Elem-07 ***


fea = [ qxi*L/3+qxf*L/8;
qxi*L/6+qxf*L/3];
% *** COMPLETE HERE - Elem-07 ***
end

%-----------------------------------------------------------------------
% Generate member axial fixed end force vector for a temperature variation
% of a given element.
% Output:
% fea: a 2x1 vector with axial fixed end force components
function fea = axialThermalLoadFEF(obj)
E = obj.material.elasticity;
alpha = obj.material.thermexp;
A = obj.section.area_X;
deltaT = obj.tempVar_GC;
EA = E*A;

% *** COMPLETE HERE - Elem-08 ***


fea = [ EA*alpha*deltaT ;
-EA*alpha*deltaT ];
% *** COMPLETE HERE - Elem-08 ***
end

%-----------------------------------------------------------------------
% Generate member transversal fixed end force vector for a linearly distributed
% load applied to given element assuming that it is simply support on both ends.
% Output:
% fes: a 2x1 vector with transversal simply supported fixed end forces:
% fes(1,1) -> transversal force at init node
% fes(2,1) -> transversal force at end node
% Input arguments:
% qyi: transversal load value at init node
% qyf: transversal load value at end (final) node
function fes = simplySuppLinearLoadFEF(obj,qyi,qyf)
L = obj.length; % member length

q0 = qyi; % uniform contribution of load


q1 = qyf - qyi; % linear contribution of load

% *** COMPLETE HERE - Elem-09 ***


fes = [ -q0*L/2-q1*L/6 ;
-q0*L/2-q1*L/3 ];
% *** COMPLETE HERE - Elem-09 ***
end

You might also like