Propped Cantilever with an end Moment first created 09/11/09 - last modified 09/11/09 Page Author: Ty Harness
A silicon beam 1[mm] long and a square cross section 10um x 10um fixed one end and propped the other is shown in figure 1. The free end is subjected to moment,M of 1[uNm]. Find the maximum deflection and the slope at the propped end using an analytic method and compare with any available finite element analysis program. Work together in groups of 2 or 3 and present your solutions at the next tutorial meet.

text here
Figure 1 - Propped Cantilever



Analytical Soloution

$E*I*(d^2y)/dx^2 = M - F*(L-x)$

$E*I*(dy)/dx = M*x - F*(L*x-x^2/2) + A$

$E*I*y = M*x^2/2 - F*(L*x^2/2-x^3/6) + A*x + B$

$y = 0$ when $x=0$ therefore $B=0$

$dy/dx=0$ when $x = 0$ therefore $A = 0$

$y = 0$ when $x=L$ therefore:

$0 = M*L^2/2 - F*(L*L^2/2-L^3/6)$

$F = (3*M)/(2*L)$

The slope equations is then:

$ E*I*dy/dx = -(M/2) * x + (3*M)/(4*L)*x^2 $

The slope at the propped end x=L:

$ dy/dx = (M*L)/(4*E*I) $

The deflection equation is then:

$ y(x) = (1/(E*I))* (-M/4*x^2 + M/(4*L)*x^3 )$

The max/min deflection will be given when dy/dx = 0

$0 = -M/2 * x + (3*M)/(4*L)*x^2 $

Using the quadratic formula you can obtain the values of x:

$x = 0$ and $x=2/3*L$

Substituting 2/3L into the deflection equation gives us the max deflection:

$y_max = y(2/3*L) = -(M*L^2)/27 $
Applying the analytical solution to the above problem

Use the following material properties: A typical value of 1.08e11 can be used for the Young's Modulus for silicon[1].

The second moment of area:

$I = (b*d^3)/12 = (10*10^-6 (10*10^-6)^3) / 12 = 8.333*10^-22 [m]^4 $

End Slope at the propped end:

$ dy/dx = (M*L)/(4*E*I) = (1*10^-6*1*10^-3)/(4*1.08*10^11 * 8.333*10^-22) = 2.778$

And the maximum y deflection is given by:

$y_max = -(M*L^2)/27 = (-1*10^-6*(1*10^-3)^2)/27 =4.115*10^-4 [m]$
Finte Element Solution

I've chosen to use Prof. Rieg's Z88 Finite Elements Program because it's free and unlimited and ideal for student use. The theory behind FEM is beyond the scope of this tutorial but you can learn a great deal from the Z88 manual. The input files Z88I1.txt and Z88I2.txt are below if you wish to repeat the analysis using the Z88 software.
Z88i1.txt

2 10 9 30 1 0 1 0 Z88I1.TXT,typed in by Ty
1 3 +0.00000E+000 +0.00000E+000 +0.00000E+000 node #1
2 3 +1.11111E-004 +0.00000E+000 +0.00000E+000 node #2
3 3 +2.22222E-004 +0.00000E+000 +0.00000E+000 node #3
4 3 +3.33333E-004 +0.00000E+000 +0.00000E+000 node #4
5 3 +4.44444E-004 +0.00000E+000 +0.00000E+000 node #5
6 3 +5.55556E-004 +0.00000E+000 +0.00000E+000 node #6
7 3 +6.66667E-004 +0.00000E+000 +0.00000E+000 node #7
8 3 +7.77778E-004 +0.00000E+000 +0.00000E+000 node #8
9 3 +8.88889E-004 +0.00000E+000 +0.00000E+000 node #9
10 3 +1.00000E-003 +0.00000E+000 +0.00000E+000 node #10
1 13 element #1
1 2
2 13 element #2
2 3
3 13 element #3
3 4
4 13 element #4
4 5
5 13 element #5
5 6
6 13 element #6
6 7
7 13 element #7
7 8
8 13 element #8
8 9
9 13 element #9
9 10
1 9 +1.08000E+011 +3.00000E-001 1 +1.00000E-010 0 0 +8.33333E-022 +5.00000E-006 0 0

5 , typed in from ty
1 1 2 0.00000E+000
1 2 2 0.00000E+000
1 3 2 0.00000E+000
10 2 2 0.00000E+000
10 3 1 1.00000E-006


text here
Figure XX - Sample picture description here



output file Z88O2.TXT : displacements, computed by Z88F V10
*************

Knoten U(1) U(2) U(3) U(4) U(5) U(6)

1 +0.0000000E+000 +0.0000000E+000 +0.0000000E+000
2 +0.0000000E+000 -3.0483113E-005 -5.1440309E-001
3 +0.0000000E+000 -1.0669091E-004 -8.2304519E-001
4 +0.0000000E+000 -2.0576109E-004 -9.2592630E-001
5 +0.0000000E+000 -3.0483134E-004 -8.2304642E-001
6 +0.0000000E+000 -3.8103986E-004 -5.1440185E-001
7 +0.0000000E+000 -4.1152280E-004 +1.8518535E-006
8 +0.0000000E+000 -3.7341868E-004 +7.2016654E-001
9 +0.0000000E+000 -2.4386518E-004 +1.6460922E+000
10 +0.0000000E+000 +0.0000000E+000 +2.7777789E+000
References


[1] Zhang L M, Uttamcham i D., "Measurement of the Mechanical Properties of Silicon Microresonators", Sensors and Actuators, A 29 79-84.

Z88 Finite Element Program

ASCII to MathML used in this page: ASCIItoMathML homepage