Elastic flexion

Storyboard

When a torque is applied to a body, it can be bent or flexed. The way this happens depends on both the geometry of the body and how the torque is applied. Additionally, it is possible to estimate the elastic energy absorbed by the body based on the deformation it undergoes.

>Model

ID:(2062, 0)



Mechanisms

Iframe

>Top



Code
Concept

Mechanisms

ID:(15570, 0)



Bending of an element

Concept

>Top


ID:(15884, 0)



Forces on a volume element

Concept

>Top


In general, the torque on beam ($M_y$) must be applied to both ends to induce bending:



The difference between both torques is equal to the torque generated by the vertical force to the axis of the beam ($Q_z$), with an interval corresponding to the variation of the position along beam ($x$):

$\displaystyle\frac{d M_y }{d x } = Q_z $



In turn, the axial force on beam ($N_x$) can vary along the axis:



This variation corresponds to the axial force along beam length ($n$) along the axis:

$\displaystyle\frac{d N_x }{d x } = - n $



Finally, the vertical force to the axis of the beam ($Q_z$) varies as a function of the load per length on the beam ($q_z$), which may correspond to the beam's own weight:



Thus, the variation of the vertical force to the axis of the beam ($Q_z$) along the position along beam ($x$) corresponds to the load per length on the beam ($q_z$) as follows:

$\displaystyle\frac{d Q_z }{dx} = - q_z $



By combining the first and last equations, we obtain the deformation equation:

$\displaystyle\frac{d^2 M_y }{d x ^2} = - q $

ID:(15885, 0)



Condition along beam

Concept

>Top


ID:(15886, 0)



Boundary condition

Concept

>Top


ID:(15887, 0)



Model

Top

>Top



Parameters

Symbol
Text
Variable
Value
Units
Calculate
MKS Value
MKS Units
$E$
E
Modulus of Elasticity
Pa
$I_b$
I_b
Moment of inertia of the section
m^4

Variables

Symbol
Text
Variable
Value
Units
Calculate
MKS Value
MKS Units
$n$
n
Axial force along beam length
N/m
$N_x$
N_x
Axial force on beam
N
$u_z$
u_z
Displacement in z
m
$q_z$
q_z
Load per length on the beam
N/m
$x$
x
Position along beam
m
$z$
z
Position at height in beam
m
$y$
y
Position within width of beam
m
$M_y$
M_y
Torque on beam
N m
$Q_z$
Q_z
Vertical force to the axis of the beam
N

Calculations


First, select the equation: to , then, select the variable: to

Calculations

Symbol
Equation
Solved
Translated

Calculations

Symbol
Equation
Solved
Translated

Variable Given Calculate Target : Equation To be used




Equations

#
Equation

$ E I_b \displaystyle\frac{d^4 u_z }{d x ^4} = q_z $

E * I_b * DIFF( u_z , x , 4 ) = q_z


$ I_b =\displaystyle\int dz dy z^2 $

I_b = @INT2( z ^2 , y , z )


$ M_y = - E I_b \displaystyle\frac{d^2 u_z }{d x ^2}$

M_y = - E * I_b * @DIF( u_z , x , 2)


$\displaystyle\frac{d M_y }{d x } = Q_z $

@DIFF( M_y , x ) = Q_z


$\displaystyle\frac{d^2 M_y }{d x ^2} = - q $

@DIFF( M_y , x , 2 ) = q


$\displaystyle\frac{d N_x }{d x } = - n $

@DIFF( N_x , x ) = - n


$\displaystyle\frac{d Q_z }{dx} = - q_z $

@DIFF( Q_z , x ) = - q_z

ID:(15571, 0)



Vertical force balance

Equation

>Top, >Model


The load per length on the beam ($q_z$) can be calculated based on how the vertical force to the axis of the beam ($Q_z$) varies along the position along beam ($x$), therefore:

$\displaystyle\frac{d Q_z }{dx} = - q_z $

$q$
Load per length on the beam
$N/m$
10432
$x$
Position along beam
$m$
10433
$Q$
Vertical force to the axis of the beam
$N$
10434

If we observe how the vertical force to the axis of the beam ($Q_z$) varies as a function of the load per length on the beam ($q_z$) over an element of the position along beam ($x$), we have:

$Q - qdx - (Q + dQ) = 0$



therefore:

$-qdx - dQ = 0$



which leads to:

$\displaystyle\frac{d Q_z }{dx} = - q_z $

ID:(15888, 0)



Axial force balance

Equation

>Top, >Model


The axial force along beam length ($n$) can be calculated based on how the axial force on beam ($N_x$) varies along the position along beam ($x$), therefore:

$\displaystyle\frac{d N_x }{d x } = - n $

$n$
Axial force along beam length
$N/m$
10437
$N$
Axial force on beam
$N$
10436
$x$
Position along beam
$m$
10433

If we observe how the axial force on beam ($N_x$) varies as a function of the axial force along beam length ($n$) over an element of the position along beam ($x$), we get:

$N - ndx - (N + dN) = 0$



therefore:

$-ndx - dN = 0$



which leads us to:

$\displaystyle\frac{d N_x }{d x } = - n $

ID:(15891, 0)



Torque balance

Equation

>Top, >Model


The vertical force to the axis of the beam ($Q_z$) can be calculated based on how the torque on beam ($M_y$) varies along the position along beam ($x$), therefore:

$\displaystyle\frac{d M_y }{d x } = Q_z $

$x$
Position along beam
$m$
10433
$M$
Torque on beam
$N m$
10435
$Q$
Vertical force to the axis of the beam
$N$
10434

The variation of the torque on beam ($M_y$) along the position along beam ($x$) is of the order of the length of the arm of element the position along beam ($x$) multiplied by the vertical force to the axis of the beam ($Q_z$), therefore:

$-M - Qdx + (M + dM) = 0$



which implies:

$-Qdx + dM = 0$



in other words:

$\displaystyle\frac{d M_y }{d x } = Q_z $

ID:(15889, 0)



Torque and load equation

Equation

>Top, >Model


The curvature of the torque on beam ($M_y$) in the position along beam ($x$) is equal to negative the load per length on the beam ($q_z$):

$\displaystyle\frac{d^2 M_y }{d x ^2} = - q $

$q$
Load per length on the beam
$N/m$
10432
$x$
Position along beam
$m$
10433
$M$
Torque on beam
$N m$
10435

Since the torque on beam ($M_y$) derived with respect to the position along beam ($x$) gives the vertical force to the axis of the beam ($Q_z$):

$\displaystyle\frac{d M_y }{d x } = Q_z $



and the derivative of the vertical force to the axis of the beam ($Q_z$) is equal to negative the load per length on the beam ($q_z$):

$\displaystyle\frac{d Q_z }{dx} = - q_z $



therefore, the derivative of the first leads us to:

$\displaystyle\frac{d^2 M_y }{d x ^2} = - q $

ID:(15890, 0)



Moment of inertia

Equation

>Top, >Model


The moment of inertia of the section ($I_b$) for a section of a bar is calculated by integrating over the section in the plane the position within width of beam ($y$) and the position at height in beam ($z$):

$ I_b =\displaystyle\int dz dy z^2 $

$I_b$
Moment of inertia of the section
0
$m^4$
10026
$z$
Position at height in beam
$m$
10438
$y$
Position within width of beam
$m$
10439

ID:(15881, 0)



Bending moment

Equation

>Top, >Model


El momento de flexión se calcula sumando todas las fuerzas que genera el sólido en la sección multiplicado por la distancia de esta al eje en que rota la sección. La suma de estos genera el momento de tensión que tiene la forma

$ M_y = - E I_b \displaystyle\frac{d^2 u_z }{d x ^2}$

$u_z$
Displacement in z
$m$
10440
$E$
Modulus of Elasticity
$Pa$
5357
$I_b$
Moment of inertia of the section
0
$m^4$
10026
$x$
Position along beam
$m$
10433
$M_y$
Torque on beam
$N m$
10435

ID:(14195, 0)



Deformation equation

Equation

>Top, >Model


With the equation for the torque on beam ($M_y$) as a function of the load per length on the beam ($q_z$) and the position along beam ($x$), and the equation for flexion involving the modulus of Elasticity ($E$), the moment of inertia of the section ($I_b$), and the displacement in z ($u_z$), we get:

$ E I_b \displaystyle\frac{d^4 u_z }{d x ^4} = q_z $

$u_z$
Displacement in z
$m$
10440
$q_z$
Load per length on the beam
$N/m$
10432
$E$
Modulus of Elasticity
$Pa$
5357
$I_b$
Moment of inertia of the section
0
$m^4$
10026
$x$
Position along beam
$m$
10433

With the torque on beam ($M_y$), the load per length on the beam ($q_z$), and the position along beam ($x$), the torque equation is established:

$\displaystyle\frac{d^2 M_y }{d x ^2} = - q $



which, combined with the displacement equation involving the modulus of Elasticity ($E$), the moment of inertia of the section ($I_b$), and the displacement in z ($u_z$):

$ M_y = - E I_b \displaystyle\frac{d^2 u_z }{d x ^2}$



results in:

$ E I_b \displaystyle\frac{d^4 u_z }{d x ^4} = q_z $

ID:(15892, 0)