Finding Stress Function and Evaluating Stresses Applying Prandtl Membrane Anology, Using Fem

Download as pdf or txt
Download as pdf or txt
You are on page 1of 5

FINDING STRESS FUNCTION AND

EVALUATING STRESSES APPLYING PRANDTL


MEMBRANE ANOLOGY, USING FEM
AS A PART OF FINITE ELEMENT METHOD AND APPLICATIONS
(ME6612D)
PROJECT WORK
SUBMIT TED BY

KISHAN NATH SIDH ULLAS U


M190064ME M190197ME
INTRODUCTION
To find Torsion of a complex cross section, it is necessary to find the stress function of the section.
 There are three equations of equilibrium, which are to be satisfied. Assuming pure torsion, the first
two conditions are automatically satisfied. The third condition is satisfied assuming a function called
as stress function.

 Membrane Analogy given by Prandtl is a simplified approach to find the stress function.
The deformation of the membrane covering a pipe of the required cross section, under unit normal
pressure is obtained, and is used to obtain shear stresses and twisting moment.
A membrane element will be used for the finite element model.
A membrane element is a 2D element, which has 3 or 4 nodes. The element can still have a
deformation in third dimension.

2
STEPS OF WORK AND TIMELINE
GANTT CHART SHOWING PROJECT TIMELINE Problem Identification

PARTICULARS Mar-20 Apr-20


W1 W2 W3 W4 W1 W2 Mathematical
formulation
Literature study Colour code
Mathematical formulation and modelling Finite element
Start or In progress
meshing
Theoritical calculations
Completed MATLAB Coding Analytical calculation

Finite element meshing


Not started
Results
Coding on Matlab

PARTICULARS Mar-20 Apr-20 Check error and


Obtaining results
W1 W2 W3 W4 W1 W2 Iterate
Literature study Colour code
Verifying and Iterating
Reports and
Mathematical formulation and modelling
Start or In progress presentation
Report and Presentation
Theoritical calculations
Completed

Finite element meshing


Not started

Coding on Matlab
3
MATHEMATICAL MODELLING
A thin homogenous membrane is stretched over a cross section having the same shape as the
shaft under test.

Internal induced forces + pressure x area = 0


After balancing vertical forces on element, we get

The equation can be compared to the differential equation,


Where,
∅ is the torsional stress function.

4
MATHEMATICAL MODELLING (CONTD…)
𝑝
According to Prandtl membrane analogy, if = 2𝐺𝜃, then height ‘z’ can be replaced by
𝐹
torsion stress function. 2 2
𝜕 ∅ 𝜕 ∅ 𝑝
+ + =0
𝜕𝑥 2 𝜕𝑦 2 𝐹
If the above relation of p/F is satisfied, z becomes numerically equal top torsional stress
function, the slope equal to the shear stresses, and the twice the volume under the
membrane is equal to the twisting moment.
𝜕2 ∅ 𝜕2 ∅ 𝑝
Residue over the domain = 𝑅Ω = 𝜕𝑥 2
+ 𝜕𝑦 2+𝐹=0

ഥ = 0 and ∅
Dirichlet Boundary condition, ∅ - ∅ ഥ = 0 on Γ𝜙
Primary variable ∅ (Torsional stress function)
Secondary variable p (Applied pressure)
5

You might also like