Compression and Bending of a Symmetric Wedge by Concentrated Forces Applied to Its Vertex (Michell’s Problem)

test422

Objective: Determination of the stress state of a symmetric wedge of unit thickness in polar coordinates subjected to compression and bending by concentrated forces applied to its vertex.

Initial data file:  4_22.spr

Problem formulation: The compressive force Px1 acting along the symmetry axis of the wedge OX1 and the bending force Px2, which is a skew-symmetric load with respect to the symmetry axis of the wedge OX1, are applied to the vertex of the wedge of unit thickness. Determine the stress tensor components in polar coordinates σrr, σθθ, σ at a radial distance r = 5.0 m from the vertex of the wedge.

References: S.P. Demidov, Theory of Elasticity. — Moscow: High school, 1979.

Initial data:

E = 3.0·107 kPa - elastic modulus;
μ = 0.2 - Poisson’s ratio;
h = 1.0 m - thickness of the wedge;
2·α = 30º - apex angle of the wedge;
R = 15.0 m - radius of the fixed end of the wedge;
Px1 = -5.0 kN - concentrated force compressing the wedge (horizontal);
Px2 = 5.0 kN - concentrated force bending the wedge (vertical).
 

 

Finite element model: Design model – general type system, wedge elements– 280 eight-node elements of type 50. The spacing of the finite element mesh in the radial direction is 0.5 m, and in the tangential direction is 3º. The direction of the output of internal forces is radial tangential. Since in the case of a cylindrical surface of a small radius a the force Px1 at the vertex of the wedge cannot be represented as a resultant of stresses distributed according to the law of the analytical solution given below, the edge of the wedge is modeled by a rigid body with a master node at the vertex of the wedge and the slave nodes at the radial distance of a = 1.0 m from the vertex of the wedge (member type – 100). Since there are no forces distributed according to the law of the analytical solution at the fixed end of the wedge, in order to obtain an exact solution at the radial distance r = 5.0 m from the action of the force Px2 the radial distance to the fixed end is taken as R = 15.0 m. Number of nodes in the design model – 918.

 

Results in SCAD

ScreenShot059
Design model

ScreenShot060ScreenShot061ScreenShot074ScreenShot075
Values of stresses σrr (kN/m2) under the compressive force Px1 and the bending force Px2

ScreenShot062ScreenShot063ScreenShot076ScreenShot077
Values of stresses σθθ (kN/m2) under the compressive force Px1 and the bending force Px2

ScreenShot064ScreenShot065ScreenShot078ScreenShot079
Values of stresses σrθ (kN/m2) under the compressive force Px1 and the bending force Px2

Comparison of solutions:

Stress tensor components at a radial distance r = 5.0 m from the vertex of the wedge under the compressive force Px1

Angle θ

Stresses σrr (kN/m2)

Theory

SCAD

Deviations, %

-15º

-1.8873

-1.8845

0.15

-1.9539

-1.9508

0.16

+15º

-1.8873

-1.8845

0.15

 

Angle θ

Stresses σθθ (kN/m2)

Theory

SCAD

Deviations, %

-15º

0.0000

-0.0005

0.0000

-0.0007

+15º

0.0000

-0.0005

 

Angle θ

Stresses σ (kN/m2)

Theory

SCAD

Deviations, %

-15º

0.0000

-0.0498

0.0000

0.0000

+15º

0.0000

0.0498

 

Stress tensor components at a radial distance r = 5.0 m from the vertex of the wedge under the bending force Px2

Angle θ

Stresses σrr (kN/m2)

Theory

SCAD

Deviations, %

-15º

-21.9350

-21.9098

0.11

0.0000

0.0000

+15º

21.9350

21.9098

0.11

 

Angle θ

Stresses σθθ (kN/m2)

Theory

SCAD

Deviations, %

-15º

0.0000

-0.0086

0.0000

0.0000

+15º

0.0000

0.0086

 

Angle θ

Stresses σ (kN/m2)

Theory

SCAD

Deviations, %

-15º

0.0000

0.5314

0.0000

0.0494

+15º

0.0000

-05314

 

Notes: In the analytical solution the stresses σrr, σθθ, σ in the body of the wedge subjected to the compressive force Px1 are determined according to the following formulas (S.P. Demidov, Theory of Elasticity. — Moscow: High school, 1979, p. 273):

\[ \sigma_{rr} =\frac{2\cdot P\cdot \cos \theta }{r\cdot \left( {2\cdot \alpha +\sin \left( {2\cdot \alpha } \right)} \right)}; \quad \sigma_{\theta \theta } =0; \quad \sigma_{r\theta } =0. \]

In the analytical solution the stresses σrr, σθθ, σ in the body of the wedge subjected to the bending force Px2 are determined according to the following formulas (S.P. Demidov, Theory of Elasticity. — Moscow: High school, 1979, p. 275):

\[ \sigma_{rr} =\frac{2\cdot P\cdot \sin \theta }{r\cdot \left( {2\cdot \alpha -\sin \left( {2\cdot \alpha } \right)} \right)}; \quad \sigma_{\theta \theta } =0; \quad \sigma_{r\theta } =0. \]