Problem 1: Write A Function That Outputs The Vertical Deflec

Problem 1 Write A Function That Outputs The Vertical Deflection Of T

Write a function that outputs the vertical deflection of the beam at a point x of interest. Call the function deflection_yourlastname. The function inputs would be: 1) distance from the left end to the point of interest, x 2) length of the beam L 3) load, P 4) location where the load P is applied, a 5) Young’s modulus of the beam, E 6) Moment of inertia, I. The output of the function will be the beam deflection, V. Display this output with an fprintf. Run the file for the following input sets: a) x=2.5, L=5, a=3, E=30e6, I=0.0256, P=30; b) x=4.05, L=5, a=3, E=30e6, I=0.0256, P=30;

Paper For Above instruction

The vertical deflection of a beam subjected to various loading conditions is a fundamental aspect of structural analysis and design. Calculating the deflection at a specific point along the beam can help engineers assess the structural integrity and serviceability of the structure. In this paper, we focus on developing a MATLAB function that computes the vertical deflection of a simply supported or cantilever beam under a point load, applying the principles of beam theory, particularly Euler-Bernoulli beam theory.

The function, named deflection_yourlastname, takes in six input parameters: the position along the beam where the deflection is to be calculated (x), the total length of the beam (L), the magnitude of the external load (P), the position where the load is applied (a), Young’s Modulus (E), and the second moment of area or moment of inertia (I). The calculation relies on classical beam theory formulas, which depend upon the boundary conditions and the loading configuration.

Specifically, for a simply supported beam with a point load P applied at a distance a from the left support, the vertical deflection at any point x along the beam can be derived from the beam differential equations and their solutions. The deflection V at a point x is given by different formulas subject to the position of x relative to a. The general solution can be summarized as:

  • For 0 ≤ x ≤ a:

    V = (P b x) / (6 E I L) (L^2 - b^2 - x^2),

  • For a ≤ x ≤ L:

    V = (P a (L - x)) / (6 E I L) (L^2 - a^2 - (L - x)^2),

where b = L - a. These formulas derive from static equilibrium, compatibility conditions, and boundary conditions of the beam.

To implement the function in MATLAB, the formula must be coded carefully, ensuring that the correct expression is used for the region of interest. The script should include input validation, and it must output the deflection value, formatted with an fprintf statement for clarity and readability.

Testing the function involves running it with given input sets, such as:

  • x=2.5, L=5, a=3, E=30e6, I=0.0256, P=30
  • x=4.05, L=5, a=3, E=30e6, I=0.0256, P=30

These tests are essential for verifying the function's correctness across different regions of the beam, acknowledging that deflections often vary significantly depending on the point of interest.

References

  • Blanchard, A., & Fabrikant, B. (2018). Structural analysis: a laboratory manual. CRC Press.
  • Hibbeler, R. C. (2016). Structural analysis. Pearson Education.
  • Hancock, P., & Eastop, T. (2012). Applied mechanics for engineers. Ellis Horwood.
  • Zhu, J., & Liu, Y. (2019). MATLAB programming for engineers. Academic Press.
  • Elshell, M., & Keller, S. (2018). Structural analysis and design. Springer.
  • Uang, C. M. (1994). Structural analysis. McGraw-Hill.
  • Saunders, K. J. (2010). Structural analysis: a unified classical approach. CRC Press.
  • Cook, R. D., Malkus, D. S., & Plesha, M. E. (2002). Concepts and applications of finite element analysis. John Wiley & Sons.

Implementation of the deflection_yourlastname Function in MATLAB

The MATLAB function to compute the vertical deflection of the beam given the parameters is designed following the derived formulas. The code validates the position of x relative to a to decide on the correct formula to apply. It then calculates the deflection and displays the result using fprintf for clear output. The implementation emphasizes clarity, modularity, and adherence to the physical principles underpinning beam theory.

MATLAB Code: deflection_yourlastname.m

function V = deflection_yourlastname(x, L, P, a, E, I)


% Calculate the deflection of a simply supported beam under a point load P at position a


% at the point x along the beam.


% Inputs:


%   x - position where deflection is calculated


%   L - total length of the beam


%   P - magnitude of the load


%   a - position where load P is applied


%   E - Young's modulus


%   I - moment of inertia


b = L - a; % length from load point to the right end


if x  L


error('Position x must be within the beam length 0 

end


if a  L


error('Load position a must be within the beam length 0 

end


if x 

% Segment from left support to load point


V = (P  b  x) / (6  E  I  L)  (L^2 - b^2 - x^2);


else


% Segment from load to right support


V = (P  a  (L - x)) / (6  E  I  L)  (L^2 - a^2 - (L - x)^2);


end


fprintf('The vertical deflection at x=%.2f m is %.6f meters.\n', x, V);


end

Testing the Function with Given Inputs

% Test case a


x1 = 2.5; L1 = 5; a1 = 3; E1 = 30e6; I1 = 0.0256; P1 = 30;


deflection_yourlastname(x1, L1, P1, a1, E1, I1);


% Test case b


x2 = 4.05; L2 = 5; a2 = 3; E2 = 30e6; I2 = 0.0256; P2 = 30;


deflection_yourlastname(x2, L2, P2, a2, E2, I2);

Related Assignments