×

Close

- Finite Element Methods - FEM
- Note
**6908 Views**- 312 Offline Downloads
- Uploaded 11 months ago

Misrimal Navajee Munoth Jain Engineering College Department of Mechanical Engineering ME 6603 – FINITE ELEMENT ANALYSIS Assignment-1 B.E. Mechanical 3rd Year Due Date:02.01.2018 1. Solve the differential equation for the physical problem expressed as with the boundary d2y 500x 2 0, 0 x 1 conditions as dx2 y(0)=0 and y(1)=0 using a two-term trial function by (i) point collocation method, (ii) subdomain collocation method, (iii) Least square method and (iv) Galerkin’s method. d2y 300x 2 0, 0 x 1 with 2. Consider the differential equation for a problem as dx2 the boundary conditions y(0)=0 and y(1)=0. Find the solution by using trial function as y a1x(1 x3 ) . Use (i) point collocation method, (ii) subdomain collocation method, (iii) Least square method and (iv) Galerkin’s method. 3. Solve the differential equation for the physical problem expressed as with the boundary d2y 100 0, 0 x 10 conditions as dx2 y(0)=0 and y(10)=0 using (i) point collocation method, (ii) subdomain collocation method, (iii) Least square method and (iv) Galerkin’s method 4. A simply supported beam (span L and flexural rigidity EI) subjected to a point load at the centre of the span. Calculate the deflection at the mid-span using Galerkin’s weighted residual method and compare with the exact solution. Assume trial solution x as y aSin . L 5. A simply supported beam (span L and flexural rigidity EI) subjected to a point load at the centre of the span. Calculate the deflection at the mid-span using Rayleigh-Ritz method and compare with the exact solution. 6. A cantilever beam is subjected to a uniformly distributed load throughout its length. Find the deflection at free end using Rayleigh-Ritz method. 7. A cantilever beam of length L is loaded with a point load at the free end. Find the maximum deflection and maximum bending moment using Rayleigh-Ritz method x using the function y a 1 Cos . Given EI is constant. 2 L d2y 400x 2 0, 0 x 1 subject to the 2 dx boundary conditions y(0)=0 and y(1)=0. The functional corresponding to the problem 2 1 dy 2 I 0 . 5 400 x y to be extremised is given by dx . Find the solution using 0 dx Rayleigh-Ritz method by considering a two-term trial solution as y( x) a1x(1 x) a2 x(1 x3 ) 8. Consider the differential equation

Misrimal Navajee Munoth Jain Engineering College Department of Mechanical Engineering Formulae for ME 6603 - Finite Element Analysis (6th Semester) Unit -1 Introduction Weighted Residual methods: 1. Point collocation method The residual(R) is forced to zero at discrete number of points. The number of points is equal to the number of undetermined parameters in the assumed trial solution. 2. Sub-domain method The average residual in the sub-domain is forced to zero. The number of sub-domains is equal to the number of undetermined parameters in the assumed trial solution. 3. Least-square method The integral of the weighted square of the residual over the domain is minimized. max R I 2 dx min Suppose if there are two undetermined coefficients C1 and C2, then I R R dx 0 C1 min C1 max I R R dx 0 C2 min C2 max 4. Galerkin Weighted residual method Galerkin introduced the idea of letting w(x) to be the same as trial functions max w ( x) R( x)dx 0 i min Rayleigh-Ritz method: 1. Assume a trial solution (generally polynomial) u ( x) c1 c2 x c3 x 2 ... 2. Evaluate the total potential U V Where U-Strain energy and V-Work done by external forces 3. Setup and solve the system of equations 0 , where i=1,2,3,…No. of undetermined parameters in trial solution ci Prepared by Dr. B. Janarthanan and Mr. A. Prakash Faculty members in Mechanical Engineering Page 1

Unit-2 One Dimensional problems 1. The general finite element equation is {F e } [ K e ]{u e } Where{Fe}-nodal force vector [Ke]-stiffness matrix {ue}-nodal field variable vector 2. For linear bar element AE 1 1 [K e ] L 1 1 3. The field variable in terms of shape functions (interpolation functions) u( x) N1 ( x)u1 N 2 ( x)u2 N3 ( x)u3 ... 4. Strain in the element, u1 dN 3 dN 2 du dN1 u 2 dx dx dx dx u3 Bu Where B-Strain displacement matrix 5. The element stiffness matrix L K B T DBdx , for 1-D case D=E, Young’s modulus e 0 6. Element stresses D 7. Body force vector { f e } [ N T qAdx , where q-body force per unit length 8. For Truss element l2 lm l 2 lm m 2 lm m 2 AE lm e [K ] Le l 2 lm l 2 lm 2 lm m 2 lm m x x y y1 Sin where l 2 1 Cos and m 2 Le Le 9. One-Dimensional heat transfer k 1 1 PhL 2 1 T1 Ph L 2 Q0 q L 1 1 6 A 1 2 T 0 A T L Q 2 c c 2 l Where Q0 and Ql represent the heat flux at the ends of the element (nodes) 10. Plane beam element 6 L 12 6 L 12 6 L 4 L2 6 L 2 L2 EI [K e ] 3 L 12 6 L 12 6 L 2 2 6L 2L 6L 4L Prepared by Dr. B. Janarthanan and Mr. A. Prakash Faculty members in Mechanical Engineering Page 2

11. Dynamic analysis K 2 m u F Where [K]-stiffness matrix [m]-mass matrix ω – natural frequency 12. Longitudinal vibration 2u 2u Equation of motion is AE 2 A 2 x t AL 2 1 Consistent mass matrix is m 6 1 2 AL 1 0 Lumped mass matrix is m 2 0 1 13. Transverse vibration 4v 2v Equation of motion is EI 4 A 2 0 dx t 156 4 L2 AL 22 L Consistent mass matrix is m 13L 420 54 2 13L 3L 1 AL 0 Lumped mass matrix is m 2 0 0 0 0 0 0 0 0 1 0 symmetric 156 22 L 4 L2 0 0 0 0 Unit-3 Two-Dimensional Scalar Variable Problems 3 ( x3,y3) 1. For a three noded triangular element Field variable, u(x,y) is given as u1 v 1 u N1 0 N 2 0 N 3 0 u 2 1 2 v 0 N1 0 N 2 0 N 3 v 2 (x1,y1) (x2,y2) u3 v3 3 x 3 1 x 1 2 x 2 Where N1 1 , N2 2 and N 3 3 2A 2A 2A 1 x2 y3 x3 y2 , 2 x3 y1 x1 y3 and 3 x1 y2 x2 y1 1 y2 y3 , 2 y3 y1 and 3 y1 y2 1 x3 x2 , 2 x1 x3 and 3 x2 x1 Prepared by Dr. B. Janarthanan and Mr. A. Prakash Faculty members in Mechanical Engineering Page 3

## Leave your Comments