Introduction to Finite Element Method
|
|
- Brice Bennett
- 6 years ago
- Views:
Transcription
1 Introduction to Finite Element Method Dr. Rakesh K Kapania Aerospace and Ocean Engineering Department Virginia Polytechnic Institute and State University, Blacksburg, VA AOE 524, Vehicle Structures Summer, 216 c 216 Rakesh K. Kapania, Mitchell Professor, Aerospace and Ocean Engineering, Virginia Polytechnic Institute and State University, Blacksburg, VA,
2 Introduction to Finite Element Method Complex structures under different types of loads Solution is often governed by ordinary or partial differential equations Simple domains or simple forcing functions analytical techniques closed form solution Complicated problems approximate solutions c Rakesh K. Kapania AOE 524, Vehicle Structures 2
3 Approximate methods Approximate methods 1. Semi-Analytical Methods 2. Numerical Methods Semi-Analytical Methods Series Expansion using Fourier, Bessel and other orthogonal functions Rayleigh-Ritz Method based on a variational principle Methods of Weighted Residual (MWR) (Galerkin Method and the Collocation Method) The Green s Functions approach (Integral Equations Approach) c Rakesh K. Kapania AOE 524, Vehicle Structures 3
4 Approximate methods (contd...) Numerical Methods Finite difference method Finite element method Spectral method Boundary element method Numerical methods (except the finite difference method) are essentially derived from one or the other semi-analytical methods. For example Finite element method discrete implementation of the Rayleigh-Ritz Method or the Method of Weighted Residuals Spectral method Series expansion method with use of FFT (coefficients) Boundary element method Integral equation approach. c Rakesh K. Kapania AOE 524, Vehicle Structures 4
5 Approximate methods (contd...) A key feature of both of these categories is to convert an ordinary or partial differential equations into a set of linear (or nonlinear depending on the problem at hand) algebraic equations, i.e. Ka = Q. K is called the stiffness matrix. a is the vector of generalized co-ordinates Q is a vector of generalized forces. This set of equations can be solved numerically. c Rakesh K. Kapania AOE 524, Vehicle Structures 5
6 Semi- Analytical Methods In the semi-analytic methods, the solution u(x, y) is generally expressed as u(x, y) = φ o (x, y) + N n=1 a n φ n (x, y) where a n are called the generalized co-ordinates and φ n (x, y) are called the trial functions and are defined over the entire domain Ω. The function φ o (x, y) is chosen such that it satisfies the non-homogeneous part of the boundary conditions. Note that this function is completely known and there is no undetermined coefficient in front of this function. c Rakesh K. Kapania AOE 524, Vehicle Structures 6
7 Trial Functions The trial functions i.e. functions φ n (x, y) have global support. φ n (x, y) satisfy the homogeneous part of the boundary conditions and must satisfy certain requirements of completeness and linear independence (the trial functions should be linearly independent) Selection of the trial functions is non- trivial Accuracy of the solution generally increases as the number of terms N in the expansion are increased Semi-analytical methods differ from each other depending upon the manner in which the coefficients a n are determined c Rakesh K. Kapania AOE 524, Vehicle Structures 7
8 Methods of Weighted Residuals In the Methods of Weighted Residuals the above approximation to the solution is substituted in the governing equations and the resulting error or the residual R(x, y) is minimized in various ways. In general, the following condition is used to obtain the values a n. R(x, y)ψ j (x, y)dω =, j = 1, 2, 3,...N Ω where the ψ j (x, y) are called the weight (hence the name Methods of Weighted Residuals) or the test functions, Ω, represents the domain. N is the number of terms used in the original expansion. c Rakesh K. Kapania AOE 524, Vehicle Structures 8
9 Methods of Weighted Residuals (contd...) The above condition is thus applied N different times and leads to a set of N algebraic equations in the coefficients, a n. The weight functions may or may not be same as the trial functions. For example, in the collocation scheme, the weight functions are the Dirac s Delta functions and in the Galerkin Method the weight functions and the trial functions are the same. c Rakesh K. Kapania AOE 524, Vehicle Structures 9
10 Ritz Method (Principle of Stationary Potential Energy) The unknown coefficients a n s are determined by using the underlying variational principle. In structures, for example, the coefficients are determined from the condition that at equilibrium the total potential energy Π of a structure under given loading must be minimum, i.e. the first variation of the total potential energy be equal to zero. In this method the N algebraic equations are obtained by satisfying the following N conditions Π a n = ; 1 n N Note that the total Potential energy is defined as Π = U + V, where U is the strain energy and V is the potential of the applied loads. c Rakesh K. Kapania AOE 524, Vehicle Structures 1
11 Principle of Virtual Work The principle of minimum potential energy on which the Ritz method is based, is restricted to conservative external forces i.e. to forces that can be derived from a potential. For non conservative forces, one has to resort to the more general Principle of Virtual Work which can be stated as: If a structure is in equilibrium and remains in equilibrium while it is subjected to a virtual distortion, the external virtual work δw vir done by the external forces acting on the structure is equal to the internal virtual work δw vir done by the internal stresses (δw vir = δw int ). c Rakesh K. Kapania AOE 524, Vehicle Structures 11
12 Principle of Virtual Work (contd..) The term virtual displacement implies an imaginary infinitesimal displacement that is superimposed on the deformed configuration while keeping the external loads unchanged and which is consistent with the boundary conditions (i.e. it is zero on that part of the boundary where the essential boundary conditions are specified). The virtual displacements, which should be piecewise continuous, to satisfy the displacement compatibility inside the domain, can be considered as the delta (δ) operator (a variation of the displacement field) or they may be considered as the test functions that are used in the method of weighted residuals (e.g. in Galerkin Method). c Rakesh K. Kapania AOE 524, Vehicle Structures 12
13 Principle of Virtual Work (contd...) Approximate methods using the principle of virtual work can be as easily derived as those using the principle of minimum potential energy. Also note that the principle of virtual work is both a sufficient and necessary condition for satisfying the equilibrium of a structure. Indeed, you can derive the equations of equilibrium and the associated boundary conditions from the principle of virtual work. c Rakesh K. Kapania AOE 524, Vehicle Structures 13
14 Principle of Virtual Work (contd...) For a general three dimensional structure, the expressions for external and internal virtual work can be written as : δw vir = B i δu i dvol + T (v) i δu i ds Vol δw int = Vol S 1 τ ij δɛ ij dvol where δu i is the virtual displacement field, B i is the body force field, T (v) i is the external surface traction acting on S 1 part of the boundary, v is the unit normal to the body. τ i j are the nine internal stress components (i, j = 1,2,3) and δɛ i j are the virtual strain components. c Rakesh K. Kapania AOE 524, Vehicle Structures 14
15 Principle of Virtual Work (contd...) Also note that a repeated index indicates a summation over that index (i.e. u i v i = u 1 v 1 + u 2 v 2 + u 3 v 3 ). Following this convention, the right hand side in the second equation is a sum of nine components. c Rakesh K. Kapania AOE 524, Vehicle Structures 15
16 Ritz Method For conservative forces, the principle of virtual work and the principle of stationary potential energy lead to the same equation, namely δ(u + V ) =. Here U is the strain energy and V is the potential of the applied loads. For a general 3-D structure the Strain Energy is given as U = dudvol Where du is the strain energy density and is given as du = Vol ɛij τ ij dɛ ij For a linear elastic material, du = 1 2 τ ijɛ ij Also recall that following the summation convention, the right hand side is a sum of nine terms. c Rakesh K. Kapania AOE 524, Vehicle Structures 16
17 Ritz Method (contd...) The potential V of the applied loads for a general 3-D structure is given as V = B i u i dvol + u i ds Vol Note that δv = δw vir S 1 T (v) i For a discrete load, Q, with a corresponding displacement, q, V can be written as: V = Qq c Rakesh K. Kapania AOE 524, Vehicle Structures 17
18 Example A Cantilever Beam Restrained By a Spring. The expression for the strain energy for a beam supported by a spring at the free end (x = L) is given as U = 1 2 EI ( d 2 ) 2 w dx 2 dx + 1 ( ) 2 2 k w (x=l) c Rakesh K. Kapania AOE 524, Vehicle Structures 18
19 Example (contd..) The expression for the potential of the applied load V is given as V = p(x)w(x)dx where w(x) is the transverse deflection of the beam. The approximate solution is obtained by using the principle of minimum total potential energy (Π): where Π = U + V. δπ = c Rakesh K. Kapania AOE 524, Vehicle Structures 19
20 Example (contd...) To obtain the approximate solution, we assume the transverse deflection as: N w(x) = C i φ i (x) i=1 where φ i (x) are the trial functions. Since there are no prescribed non homogeneous displacements and slope boundary conditions, we do not need φ (x). c Rakesh K. Kapania AOE 524, Vehicle Structures 2
21 Example (contd...) The trial functions must satisfy some requirements: (i) functions should be linearly independent. This means that the following equation will hold if and only if all C i s are zero. N C i φ i (x) = i=1 (ii) Functions must satisfy the essential boundary conditions. For the problem at hand, it means that both φ i (x) = and φ i (x) = at x =. Note that in the Ritz method the force (natural) boundary conditions are not explicitly satisfied. As the number of terms in the expansion for the dependent variable is increased, the natural boundary conditions tend to be satisfied naturally (Hence the name Natural Boundary Conditions). c Rakesh K. Kapania AOE 524, Vehicle Structures 21
22 Example (contd...) (iii) The trial functions must form a complete set. This means that the solution of the problem can be approximated arbitrarily close by a linear combination of the given trial functions as the number of terms in the expansion are increased. If we are using sine (or cosine) series, we must use sin(ix)(cosine(ix)) before we use sin(i + 1)x(cosine(i + 1)x). Similarly as we use (x) i+1, we must use (x) i. c Rakesh K. Kapania AOE 524, Vehicle Structures 22
23 Example (contd...) Substituting the expansion for w(x) in the expression for total potential energy, and using the principle of minimum potential energy, we get: N δπ δπ = δc i = i=1 δc i Since the variation in each of the generalized displacements C i s arbitrary and independent from the variation in the rest of the generalized displacements C i s, it means that Π C i = ; 1 i N c Rakesh K. Kapania AOE 524, Vehicle Structures 23
24 Example (contd...) This will yield a set of N equations. For the present beam problem, the expression for Π can be written as: C i Π L ( = 2 ) w 2 w EI C i x 2 C i x 2 dx + kw w L L L p(x) w dx C i C i [ N ] [ N ] = C j φ j (x) φ i (x)dx + k C j φ j (L) φ i (L) j=1 j=1 p(x)φ i (x)dx ] = N [ EI (x)φ i (x)φ j (x)dx + kφ i (L)φ j (L) j=1 p(x)φ i (x)dx = N j=1 C j K ij C j Q i c Rakesh K. Kapania AOE 524, Vehicle Structures 24
25 Example (contd...) Where the stiffness coefficient K ij and the generalized force Q i are given as: K ij = EI (x)φ i (x)φ j (x)dx + kφ i (L)φ j (L) Q i = p(x)φ i (x)dx For the present problem, we can choose where i varies from 1 to N. ( x ) i+1 φ i (x) = L c Rakesh K. Kapania AOE 524, Vehicle Structures 25
26 Example (contd...) Note that these trial functions satisfy all the essential boundary conditions, namely, φ() = and φ () =. Also the φ s are linearly independent and note that to satisfy the essential boundary conditions, we are starting our expansion from ( x ) 2 ( x ), (for i = l) and not from i.e. a constant term as L L required by the completeness requirement. The solution to our problem does not contain any constant or linear terms in the expansion. The present expansion will thus satisfy the completeness requirement. c Rakesh K. Kapania AOE 524, Vehicle Structures 26
27 Example (contd...) Substituting the chosen trial functions in the expressions for the stiffness matrix and the generalized forces (note (φ i (L) = 1), we obtain: ( ) x i+j 2 K ij = EI (x)(i)(j)(i + 1)(j + 1) dx + k = Q i = = ( (i)(j)(i + 1)(j + 1) (i + j 1) ( x ) ( x ) i+1 p dx L L p L (i + 3) ) ( EI L 3 L i+j+2 ) + α ( EI L 3 Using the above expressions we can obtain the approximate solutions using various values of N. c Rakesh K. Kapania AOE 524, Vehicle Structures 27 )
28 Example (contd...) Case (i) (n = 1) The set of linear equations for this case can be written as: EI L 3 (4 + α)c 1 = p L 4 p L C 1 = 4(4 + α)ei c Rakesh K. Kapania AOE 524, Vehicle Structures 28
29 Example (contd...) Case (ii) (n = 2) The set of linear equations for this case can be written as: [ EI 4 + α 6 + α L α 12 + α C 2 ] { C1 C 2 } = p L { Solution of this set of two equations will give us the values of the undetermined constants C 1 and C 2. { } C1 p L = 4 { } 36 + α 2EI (12 + 4α) (14 + α) } c Rakesh K. Kapania AOE 524, Vehicle Structures 29
30 Example (contd...) Case (iii) N = 3 The set of three linear equations is: EI C 1 L 3 C 2 = p L C α 6 + α 8 + α 6 + α 12 + α 18 + α 8 + α 18 + α α Using Gauss elimination, the solution is C 1 = (8α + 123)p L 4 24EI (α + 3) C 2 = (13α + 72)p L 4 24EI (α + 3) C 3 = p L 4 48EI c Rakesh K. Kapania AOE 524, Vehicle Structures 3
31 Example (contd...) Case (iv) N = 4 The set of four equations is: EI L α 6 + α 8 + α 1 + α 6 + α 12 + α 18 + α 24 + α 8 + α 18 + α α 4 + α α 24 + α 4 + α 7 + α C 1 C 2 C 3 C 4 = p L Solution of these four equations using Gauss elimination yields: C 1 = (7α + 2)p L 4 24EI (α + 3) C 3 = C 2 = (3α + 2)p L 4 8EI (α + 3) C 4 = p L 4 12EI It can be shown that this is also the exact solution to this problem. As a result, we do not have to proceed any further. c Rakesh K. Kapania AOE 524, Vehicle Structures
32 Example (contd...) For Complex problem exact solution is not known We have to decide how many terms in the expansion will yield an accurate solution We stop when we see that addition of more terms does not make any significant change in the desired response quantity (generally stresses) For a given number of terms, the accuracy of the stresses will be less than that for the displacements Differentiation of an approximation is almost always less accurate than the function c Rakesh K. Kapania AOE 524, Vehicle Structures 32
33 Penalty Approach Solution for various values of the spring constant can be easily determined by using appropriate value of α. As a special case, consider the case of a simply-supported end at x = L. This case is achieved by letting the value of the spring constant to be infinite (α infinity). The limiting value of the solution is: w(x) = p L 4 [ 7 ( x ) 2 3 ( x ) 3 1 ( x ) ] 5 + EI 24 L 8 L 12 L c Rakesh K. Kapania AOE 524, Vehicle Structures 33
34 Penalty Approach (contd...) Note that in obtaining this solution we did not satisfy the essential boundary condition at x = L explicitly. This suggests a way around satisfying the essential boundary conditions explicitly, by placing a spring with very large stiffness at the point where the displacement is zero. This approach is called the penalty approach. Instead of looking for trial functions that are zero at x =, and, x = L and have slope = at x =, we only satisfied the essential boundary conditions at x = exactly and will satisfy the essential boundary condition at x = L only approximately. c Rakesh K. Kapania AOE 524, Vehicle Structures 34
35 Limitations of the Ritz Method In the above approach, difficulty is that the selection of trial functions is very time consuming and for complicated domains with complex set of restraints, it may well be impossible to come up with a set of trial functions that will satisfy all the essential boundary conditions. c Rakesh K. Kapania AOE 524, Vehicle Structures 35
36 Finite Element Method In finite element method, in which the domain is divided ( discretized ) into a large number of smaller domains ( called finite elements) and simple polynomials that are defined only locally (i.e. only over a given finite element) are used as trial functions. Using either a weighted residual method or Ritz method, a set of algebraic equations are derived for every finite element and are then assembled (in a procedure called assemblage ) together to form a global set of equations. Another advantage the finite element method has over the global methods is that the matrix K is banded (sparse). This is due to the fact that a given trial function is defined only locally. It can easily solve linear and nonlinear problems and static and dynamic problems. c Rakesh K. Kapania AOE 524, Vehicle Structures 36
37 Finite Element Method (contd...) c Rakesh K. Kapania AOE 524, Vehicle Structures 37
38 Finite Element Method (contd...) Consider the solution of the ordinary differential equation governing the axial response of a bar with axial rigidity EA(x), and subjected to an axial force p(x): d dx ( EA(x) du dx ) = p(x) Let the boundary conditions be : u() = and u(l) = 1. An approximate solution of this equation can be found by satisfying the stationary condition for the total potential potential energy Π, which for an axial bar is given as: Π = 1 ( ) du 2 EA(x) dx p(x)u(x)dx 2 dx In the finite element method, we assume u(x) to be : N u(x) = a i φ i (x) i=1 c Rakesh K. Kapania AOE 524, Vehicle Structures 38
39 Finite Element Method (contd...) In the finite element method, the generalized coefficients are the values of the dependent variables (the axial displacement u in the present case) at the nodes (i.e a i = u i ). This is an important difference between the finite element method and the Ritz method. In the Ritz method, the generalized co-ordinates a i s, in general, do not have any physical significance. c Rakesh K. Kapania AOE 524, Vehicle Structures 39
40 Finite Element Method (contd...) By taking the first variation of the potential energy with respect to the nodal displacement values, we obtain: Π δπ = u i = u i N i=1 Since the variation of each of the generalized coordinates (nodal displacements u i s) is arbitrary and independent from the variation of all the other generalized co-ordinates, we get a set of N equations: Π =, i = 1, 2...N u i c Rakesh K. Kapania AOE 524, Vehicle Structures 4
41 Finite Element Method (contd...) The left hand side in the above equation can also be written as: Π L ( ) ( ) u u L = EA(x) dx p(x) (u(x))dx u i x u i x) u i = = N j=1 [ EA(x) [ ( N ) u j φ j (x) j=1 φ i (x) p(x)φ i (x) ] dx ] EA(x)φ i (x)φ j (x)dx u j p(x)φ i (x)dx where a prime indicates the derivative with respect to x. c Rakesh K. Kapania AOE 524, Vehicle Structures 41
42 Finite Element Method (contd...) For a linear problem, the resulting set of equations will be linear and can be written as: Ka = Q. Here, K the stiffness matrix a the vector of nodal displacements Q the nodal load vector. The elements of the stiffness matrix K, and the load vector Q are given as: K ij = Q i = EA(x)φ i (x)φ j (x)dx p(x)φ i (x)dx c Rakesh K. Kapania AOE 524, Vehicle Structures 42
43 Finite Element Method (contd...) The stiffness coefficient K ij have a physical meaning: it is the force required at node i due to a unit displacement at node j provided all other displacements are zero. Note that the stiffness matrix is symmetric but sparse. For example, the stiffness coefficients K 13, K 14,... =. To illustrate this, let us examine the stiffness co-efficient K 13. K 13 = φ 1(x)φ 3(x)dx Product is zero because the trial function φ 1 (x) is zero on that part of the domain on which the trial function φ 3 (x) is non-zero and vice-versa. Similarly we can see that for the second-order ordinary differential equations, the stiffness matrix will have a bandwidth of 3. c Rakesh K. Kapania AOE 524, Vehicle Structures 43
44 Finite Element Method (contd...) For example, in the third row, the three nonzero coefficients will be : K 32, K 33, and K 34, and in the i th row, the three non-zero coefficients will be K i,i 1, K i,i, and K i,i+1. Let P i be an axial force acting at node i. In that case, the ith generalized force will be modified as Q i = P i + Q i c Rakesh K. Kapania AOE 524, Vehicle Structures 44
45 Finite Element Method (contd...) This can be easily seen from the fact that if a point force P k is acting at the k th node, the expression for the total potential energy will be modified to contain an additional term, P k u k. When the variation of this term is taken, it leads to the additional force term P i in the expression for the generalized force Q i. If for the given problem, we divide the domain (, L) using four elements (or five nodes), the resulting set of equations will be as follows: K 11 K 12 K 21 K 22 K 23 K 32 K 33 K 34 K 43 K 44 K 45 K 54 K 55 u 1 = u 2 u 3 u 4 u 5 = 1 = Q 1 Q 2 Q 3 Q 4 Q 5 c Rakesh K. Kapania AOE 524, Vehicle Structures 45
46 Finite Element Method (contd...) Note that in the above equation, there are five unknowns including: the two generalized forces Q 1 and Q 5 and the three generalized (nodal) displacements u 2, u 3, and u 4. Also note that the 5 by 5 global stiffness matrix K in the above set of linear equations is singular. This is due to the fact that the axial bar is not restrained from moving freely in the space. If we determine the eigenvalues of this global stiffness matrix, one of the eigenvalue will be found to be zero and all the elements in the corresponding eigenvector will be unity (i.e. a rigid body motion in the x direction). c Rakesh K. Kapania AOE 524, Vehicle Structures 46
47 Finite Element Method (contd...) For numerical results to converge, it is very important that the global stiffness matrix be singular and contain appropriate rigid body modes. Lack of this singularity will imply that a rigid body motion (in our case, all u i = 1) will result in a nonzero force vector, i.e. a free body motion will result in a non-zero axial strain which is of course not true. To obtain a unique solution for a given problem we remove this singularity by using the essential boundary conditions. In the present problem, this will be accomplished by substituting u 1 = and u 5 = 1. c Rakesh K. Kapania AOE 524, Vehicle Structures 47
48 Finite Element Method (contd...) The three displacements can be found by solving the following reduced set of equations: K 22 K 23 u 2 Q 2 K 32 K 33 K 34 u 3 = Q 3 u 5 K 43 K 44 u 4 Q 4 K 45 The first vector on right hand side can be obtained by performing the required integrals in the expressions for generalized forces and the second vector contains products of some stiffness coefficients with the generalized displacement u 5 = 1. c Rakesh K. Kapania AOE 524, Vehicle Structures 48
49 Finite Element Method (contd...) Once the three unknown displacements are determined from the above equations, we can obtain the unknown generalized forces Q 1 and Q 5 by substituting the values of all the nodal displacements in the full stiffness matrix. Note that Q 1 represents the reaction at the node 1 and the Q 5 represents the force required to produce u 5 = 1. We can easily change the boundary conditions in the finite element method. For example, suppose that we have applied a force P at the free edge (x = L). In that case, the displacement at the free edge (u 5 if we divide the domain in four elements) will be unknown but the force Q 5 = Q 5 + P 5 will be known. c Rakesh K. Kapania AOE 524, Vehicle Structures 49
50 Finite Element Method (contd...) As another example, suppose we have a spring of stiffness k at the free end (x = L). This is equivalent to applying a force P 5 = ku 5 at the free end. Since this term contains the unknown u 5 we can take it to the left hand side in the fifth (last) equation. This will change the stiffness coefficient K 55 to K 55 + k. c Rakesh K. Kapania AOE 524, Vehicle Structures 5
51 Example EI 4 w x 4 = p (p is uniform) I (w) = w() = w () = EIw x=l = M EI (w ) = [ ] 1 2 EI (w ) 2 pw dx M w x x=l c Rakesh K. Kapania AOE 524, Vehicle Structures 51
52 Example (contd...) Assume w(x) = N C i φ i (x) (N-parameter Ritz method) j=1 φ i must satisfy essential boundary conditions. i.e. φ(x) = at x = ; φ () = Furthermore, the φ i s should be linearly i.e. independent a 1 φ 1 + a 2 φ a n φ n = if a 1 = a 2 = = a n = c Rakesh K. Kapania AOE 524, Vehicle Structures 52
53 Example (contd...) i.e. if we chose φ 1 = x 2, we cannot choose another φ j = ax 2, where a is a constant. Let φ j (x) = x j+1, j = 1, 2, 3,, N since x j+1 = at x = at x j+1 x = (j + 1)x j = x = c Rakesh K. Kapania AOE 524, Vehicle Structures 53
54 Example (contd...) Furthermore, x j+1 is linearly independent of x j+2, x j+3, w = C 1 x 2 + C 2 x 3 + C 3 x 4 + C 4 x C N x N+1 if N = 2, we stop at x 3 if N = 3, we stop at x 4.. w = w = N C j x j+1 j=1 N C j (j + 1)(j)x j 1 j=1 c Rakesh K. Kapania AOE 524, Vehicle Structures 54
55 Example (contd...) I = 1 2 EI ( N C j φ j j=1 M [ N j=1 ) 2 p C j φ j x=l ( N j=1 C j φ j ) ] dx = I [C 1, C 2, C 3,, C N ] I C 1 = I I = = C 2 C N which are a set of linear equations. c Rakesh K. Kapania AOE 524, Vehicle Structures 55
56 Example (contd...) I C i = 2 EI 2 [ N C j φ j j=1 ] φ i dx l pφ i dx M φ i x=l = for i = 1, 2,, N N [ ] EI φ i φ j dx C j = j=1 for i=1,2,,n l pφ i dx + M φ i x=l c Rakesh K. Kapania AOE 524, Vehicle Structures 56
57 Example (contd...) Let K ij = f i = EI φ i φ j dx = B(φ i, φ j ) pφ i dx + M φ i x=l = l(φ i ) N K ij C j j=1 for i=1,2,3,,n K 11 K 12 K 1N K 21 K 22 K 2N K N1 K N2 K 1N = f i c Rakesh K. Kapania AOE 524, Vehicle Structures 57 C 1 C 2. C N = f 1 f 2. f N
58 Example (contd...) Note K ij = K ji [K ]{C } = {f } where [K] is the stiffness matrix, { C } is the generalized displacement vector, and { f } is the generalized load vector. Also it can be shown that [K] is positive definite. So, {C } = [K ] 1 {f } c Rakesh K. Kapania AOE 524, Vehicle Structures 58
59 Example (contd...) Let N=2: w = C 1 x 2 + C 2 x 3 i.e. φ 1 = x 2, φ 2 = x 3, φ 1 = 2, φ 2 = 6x K 11 = = = 4EIL EI φ 1 φ 1 dx EI 2 2 dx K 12 = K 21 = = = 6EIL 2 EI 2 6x dx EI φ 1 φ 2 dx K 22 = = 12EIL 3 EI φ 2 φ 2 dx = EI 6 6 dx c Rakesh K. Kapania AOE 524, Vehicle Structures 59
60 Example (contd...) f 1 = = pφ 1 dx + M (φ 1) x=l px 2 dx + M (2x) x=l = pl3 + 2M L 3 [ ] { } 4L 6L 2 C1 EI 6L 2 12L 3 C 2 f 2 = = = pφ 2 dx + M (φ 2) x=l px 3 dx + M (3x 2 ) x=l = pl4 + 3M L 2 4 { } pl M L pl M L 2 c Rakesh K. Kapania AOE 524, Vehicle Structures 6
61 Example (contd...) { C1 C 2 } = 1 EI [ 4L 6L 2 6L 2 12L 3 ] { 1 pl M L pl M L 2 } = [ 1 EI 1 12L 3 6L 2 12L 4 6L 2 4L ] { pl M L pl M L 2 } = { 1 4pL M L EIL 4 4 pl6 18M L 4 2pL 5 12M L 3 + pl 4 L + 12M L 3 } = 1 { 5 2 pl6 + 6M L 4 12EIL 4 pl 5 c Rakesh K. Kapania AOE 524, Vehicle Structures 61 }
62 Example (contd...) w(x) = C 1 x 2 + C 2 x 3 ( 5pL = 2 ) + 12M x EI ( ) pl x 3 12EI Now, for N = 3 w(x) = C 1 x 2 + C 2 x 3 + C 3 x 4 A better approximation Objective : To determine C 1, C 2, C 3 φ 1 = x 2 ; φ 1 = 2x; φ 1 = 2 φ 2 = x 3 ; φ 2 = 3x 2 ; φ 2 = 6x φ 3 = x 4 ; φ 3 = 4x 3 ; φ 3 = 12x 2 c Rakesh K. Kapania AOE 524, Vehicle Structures 62
63 Example (contd...) K 13 = K 31 K 23 = K 32 = EI φ 1 φ 3 dx = EI φ 2 φ 3 dx = EI 2 12x 2 dx = EI 6x 12x 2 dx K 13 = 8EIL 3 K 23 = 18EIL 4 c Rakesh K. Kapania AOE 524, Vehicle Structures 63
64 Example (contd...) K 22 = = EI φ 3 φ 3 dx EI (12x 2 ) 2 dx f 3 = = pφ 3 dx + M (φ 1) x=l p x 4 dx + M 4x 3 x=l K 22 = 144EIL5 5 f 3 = pl M L 3 c Rakesh K. Kapania AOE 524, Vehicle Structures 64
65 Example (contd...) EIL 4 6L 8l 2 6L 12L 2 18L 3 8L 2 18L L4 SOLVING GIVES: C 1 C 2 C 3 = PL M L PL M L 2 PL M L 3 C 1 = pl2 + 2M 4EI C 2 = pl 6EI p C 3 = 24EI c Rakesh K. Kapania AOE 524, Vehicle Structures 65
Chapter 3 Variational Formulation & the Galerkin Method
Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation
More informationCOPYRIGHTED MATERIAL. Index
Index A Admissible function, 163 Amplification factor, 36 Amplitude, 1, 22 Amplitude-modulated carrier, 630 Amplitude ratio, 36 Antinodes, 612 Approximate analytical methods, 647 Assumed modes method,
More informationQuintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationAdvanced Vibrations. Distributed-Parameter Systems: Approximate Methods Lecture 20. By: H. Ahmadian
Advanced Vibrations Distributed-Parameter Systems: Approximate Methods Lecture 20 By: H. Ahmadian ahmadian@iust.ac.ir Distributed-Parameter Systems: Approximate Methods Rayleigh's Principle The Rayleigh-Ritz
More informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
More informationVirtual Work and Variational Principles
Virtual Work and Principles Mathematically, the structural analysis problem is a boundary value problem (BVP). Forces, displacements, stresses, and strains are connected and computed within the framework
More informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More informationRitz Method. Introductory Course on Multiphysics Modelling
Ritz Method Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy of Sciences Warsaw Poland
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 06
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 06 In the last lecture, we have seen a boundary value problem, using the formal
More informationStress analysis of a stepped bar
Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has cross-sectional areas of A ) and A ) over the lengths l ) and l ), respectively.
More informationChapter 12 Elastic Stability of Columns
Chapter 12 Elastic Stability of Columns Axial compressive loads can cause a sudden lateral deflection (Buckling) For columns made of elastic-perfectly plastic materials, P cr Depends primarily on E and
More informationLecture 15: Revisiting bars and beams.
3.10 Potential Energy Approach to Derive Bar Element Equations. We place assumptions on the stresses developed inside the bar. The spatial stress and strain matrices become very sparse. We add (ad-hoc)
More informationContents. Prologue Introduction. Classical Approximation... 19
Contents Prologue........................................................................ 15 1 Introduction. Classical Approximation.................................. 19 1.1 Introduction................................................................
More informationPart D: Frames and Plates
Part D: Frames and Plates Plane Frames and Thin Plates A Beam with General Boundary Conditions The Stiffness Method Thin Plates Initial Imperfections The Ritz and Finite Element Approaches A Beam with
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture 1-20 September, 2017
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Lecture 1-20 September, 2017 Institute of Structural Engineering Method of Finite Elements II 1 Course
More informationComputational Stiffness Method
Computational Stiffness Method Hand calculations are central in the classical stiffness method. In that approach, the stiffness matrix is established column-by-column by setting the degrees of freedom
More informationLevel 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method
9210-203 Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method You should have the following for this examination one answer book No additional data is attached
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationBasic Energy Principles in Stiffness Analysis
Basic Energy Principles in Stiffness Analysis Stress-Strain Relations The application of any theory requires knowledge of the physical properties of the material(s) comprising the structure. We are limiting
More informationMethods of Analysis. Force or Flexibility Method
INTRODUCTION: The structural analysis is a mathematical process by which the response of a structure to specified loads is determined. This response is measured by determining the internal forces or stresses
More informationA The Finite Element Method
A The Finite Element Method Most of this book has been concerned with exact analytical methods for solving problems in mechanics of materials. Real engineering applications seldom involve geometries or
More informationAdvanced Vibrations. Distributed-Parameter Systems: Exact Solutions (Lecture 10) By: H. Ahmadian
Advanced Vibrations Distributed-Parameter Systems: Exact Solutions (Lecture 10) By: H. Ahmadian ahmadian@iust.ac.ir Distributed-Parameter Systems: Exact Solutions Relation between Discrete and Distributed
More informationDHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS PART A (2 MARKS)
DHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS UNIT I : FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PART A (2 MARKS) 1. Write the types
More informationStructural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian
Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationPrepared by M. GUNASHANKAR AP/MECH DEPARTMENT OF MECHANICAL ENGINEERING
CHETTINAD COLLEGE OF ENGINEERING AND TECHNOLOGY-KARUR FINITE ELEMENT ANALYSIS 2 MARKS QUESTIONS WITH ANSWER Prepared by M. GUNASHANKAR AP/MECH DEPARTMENT OF MECHANICAL ENGINEERING FINITE ELEMENT ANALYSIS
More informationBHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I
BHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I 635 8 54. Third Year M E C H A NICAL VI S E M ES TER QUE S T I ON B ANK Subject: ME 6 603 FIN I T E E LE ME N T A N A L YSIS UNI T - I INTRODUCTION
More informationa x Questions on Classical Solutions 1. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress
Questions on Classical Solutions. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress σ xy = T is applied at infinity. Determine the value of the stress σ θθ on the edge
More informationStructural Dynamics A Graduate Course in Aerospace Engineering
Structural Dynamics A Graduate Course in Aerospace Engineering By: H. Ahmadian ahmadian@iust.ac.ir The Science and Art of Structural Dynamics What do all the followings have in common? > A sport-utility
More informationCHAPTER 5. Beam Theory
CHPTER 5. Beam Theory SangJoon Shin School of Mechanical and erospace Engineering Seoul National University ctive eroelasticity and Rotorcraft Lab. 5. The Euler-Bernoulli assumptions One of its dimensions
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 11
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Module - 01 Lecture - 11 Last class, what we did is, we looked at a method called superposition
More informationCIVL 8/7117 Chapter 12 - Structural Dynamics 1/75. To discuss the dynamics of a single-degree-of freedom springmass
CIV 8/77 Chapter - /75 Introduction To discuss the dynamics of a single-degree-of freedom springmass system. To derive the finite element equations for the time-dependent stress analysis of the one-dimensional
More informationFinal Exam Solution Dynamics :45 12:15. Problem 1 Bateau
Final Exam Solution Dynamics 2 191157140 31-01-2013 8:45 12:15 Problem 1 Bateau Bateau is a trapeze act by Cirque du Soleil in which artists perform aerial maneuvers on a boat shaped structure. The boat
More informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationInstitute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I
Institute of Structural Engineering Page 1 Chapter 2 The Direct Stiffness Method Institute of Structural Engineering Page 2 Direct Stiffness Method (DSM) Computational method for structural analysis Matrix
More informationA HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS
A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras (Refer Slide Time: 00:25) Module - 01 Lecture - 13 In the last class, we have seen how
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationFINITE ELEMENT METHOD: APPROXIMATE SOLUTIONS
FINITE ELEMENT METHOD: APPROXIMATE SOLUTIONS I Introduction: Most engineering problems which are expressed in a differential form can only be solved in an approximate manner due to their complexity. The
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 17
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil engineering Indian Institute of Technology, Madras Module - 01 Lecture - 17 In the last class, we were looking at this general one-dimensional
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Variational Problems of the Dirichlet BVP of the Poisson Equation 1 For the homogeneous
More information(L, t) = 0, t > 0. (iii)
. Sturm Liouville Boundary Value Problems 69 where E is Young s modulus and ρ is the mass per unit volume. If the end x = isfixed, then the boundary condition there is u(, t) =, t >. (ii) Suppose that
More informationLecture 27: Structural Dynamics - Beams.
Chapter #16: Structural Dynamics and Time Dependent Heat Transfer. Lectures #1-6 have discussed only steady systems. There has been no time dependence in any problems. We will investigate beam dynamics
More informationEigenvalues of Trusses and Beams Using the Accurate Element Method
Eigenvalues of russes and Beams Using the Accurate Element Method Maty Blumenfeld Department of Strength of Materials Universitatea Politehnica Bucharest, Romania Paul Cizmas Department of Aerospace Engineering
More informationLecture Slides. Chapter 4. Deflection and Stiffness. The McGraw-Hill Companies 2012
Lecture Slides Chapter 4 Deflection and Stiffness The McGraw-Hill Companies 2012 Chapter Outline Force vs Deflection Elasticity property of a material that enables it to regain its original configuration
More informationME 1401 FINITE ELEMENT ANALYSIS UNIT I PART -A. 2. Why polynomial type of interpolation functions is mostly used in FEM?
SHRI ANGALAMMAN COLLEGE OF ENGINEERING AND TECHNOLOGY (An ISO 9001:2008 Certified Institution) SIRUGANOOR, TIRUCHIRAPPALLI 621 105 Department of Mechanical Engineering ME 1401 FINITE ELEMENT ANALYSIS 1.
More informationBAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS
Journal of Computational and Applied Mechanics, Vol.., No. 1., (2005), pp. 83 94 BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Vladimír Kutiš and Justín Murín Department
More informationCIVL 7/8117 Chapter 4 - Development of Beam Equations - Part 2 1/34. Chapter 4b Development of Beam Equations. Learning Objectives
CIV 7/87 Chapter 4 - Development of Beam Equations - Part /4 Chapter 4b Development of Beam Equations earning Objectives To introduce the work-equivalence method for replacing distributed loading by a
More informationIntroduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.
Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.
More information2 marks Questions and Answers
1. Define the term strain energy. A: Strain Energy of the elastic body is defined as the internal work done by the external load in deforming or straining the body. 2. Define the terms: Resilience and
More informationUnit 18 Other Issues In Buckling/Structural Instability
Unit 18 Other Issues In Buckling/Structural Instability Readings: Rivello Timoshenko Jones 14.3, 14.5, 14.6, 14.7 (read these at least, others at your leisure ) Ch. 15, Ch. 16 Theory of Elastic Stability
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationAssumptions: beam is initially straight, is elastically deformed by the loads, such that the slope and deflection of the elastic curve are
*12.4 SLOPE & DISPLACEMENT BY THE MOMENT-AREA METHOD Assumptions: beam is initially straight, is elastically deformed by the loads, such that the slope and deflection of the elastic curve are very small,
More informationWeighted Residual Methods
Weighted Residual Methods Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy of Sciences
More information1 Nonlinear deformation
NONLINEAR TRUSS 1 Nonlinear deformation When deformation and/or rotation of the truss are large, various strains and stresses can be defined and related by material laws. The material behavior can be expected
More informationGeneralized Single Degree of Freedom Systems
Single Degree of Freedom http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March, 8 Outline Until now our were described
More informationVIBRATION PROBLEMS IN ENGINEERING
VIBRATION PROBLEMS IN ENGINEERING FIFTH EDITION W. WEAVER, JR. Professor Emeritus of Structural Engineering The Late S. P. TIMOSHENKO Professor Emeritus of Engineering Mechanics The Late D. H. YOUNG Professor
More informationStiffness Matrices, Spring and Bar Elements
CHAPTER Stiffness Matrices, Spring and Bar Elements. INTRODUCTION The primary characteristics of a finite element are embodied in the element stiffness matrix. For a structural finite element, the stiffness
More informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More informationDr. D. Dinev, Department of Structural Mechanics, UACEG
Lecture 6 Energy principles Energy methods and variational principles Print version Lecture on Theory of Elasticity and Plasticity of Dr. D. Dinev, Department of Structural Mechanics, UACEG 6.1 Contents
More informationSteps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen
Steps in the Finite Element Method Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen General Idea Engineers are interested in evaluating effects such as deformations, stresses,
More informationSection 12.6: Non-homogeneous Problems
Section 12.6: Non-homogeneous Problems 1 Introduction Up to this point all the problems we have considered are we what we call homogeneous problems. This means that for an interval < x < l the problems
More information2C9 Design for seismic and climate changes. Jiří Máca
2C9 Design for seismic and climate changes Jiří Máca List of lectures 1. Elements of seismology and seismicity I 2. Elements of seismology and seismicity II 3. Dynamic analysis of single-degree-of-freedom
More informationEML4507 Finite Element Analysis and Design EXAM 1
2-17-15 Name (underline last name): EML4507 Finite Element Analysis and Design EXAM 1 In this exam you may not use any materials except a pencil or a pen, an 8.5x11 formula sheet, and a calculator. Whenever
More informationWeak form of Boundary Value Problems. Simulation Methods in Acoustics
Weak form of Boundary Value Problems Simulation Methods in Acoustics Note on finite dimensional description of functions Approximation: N f (x) ˆf (x) = q j φ j (x) j=1 Residual function: r(x) = f (x)
More informationAdvanced Vibrations. Elements of Analytical Dynamics. By: H. Ahmadian Lecture One
Advanced Vibrations Lecture One Elements of Analytical Dynamics By: H. Ahmadian ahmadian@iust.ac.ir Elements of Analytical Dynamics Newton's laws were formulated for a single particle Can be extended to
More informationPost Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method
9210-220 Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method You should have the following for this examination one answer book scientific calculator No
More informationChapter 4 Deflection and Stiffness
Chapter 4 Deflection and Stiffness Asst. Prof. Dr. Supakit Rooppakhun Chapter Outline Deflection and Stiffness 4-1 Spring Rates 4-2 Tension, Compression, and Torsion 4-3 Deflection Due to Bending 4-4 Beam
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationMECh300H Introduction to Finite Element Methods. Finite Element Analysis (F.E.A.) of 1-D Problems
MECh300H Introduction to Finite Element Methods Finite Element Analysis (F.E.A.) of -D Problems Historical Background Hrenikoff, 94 frame work method Courant, 943 piecewise polynomial interpolation Turner,
More informationDr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum
STRUCTURAL DYNAMICS Dr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum Overview of Structural Dynamics Structure Members, joints, strength, stiffness, ductility Structure
More information16.21 Techniques of Structural Analysis and Design Spring 2005 Unit #8 Principle of Virtual Displacements
6. Techniques of Structural Analysis and Design Spring 005 Unit #8 rinciple of irtual Displacements Raúl Radovitzky March 3, 005 rinciple of irtual Displacements Consider a body in equilibrium. We know
More informationMechanical Design in Optical Engineering
OPTI Buckling Buckling and Stability: As we learned in the previous lectures, structures may fail in a variety of ways, depending on the materials, load and support conditions. We had two primary concerns:
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 03 Quiz # April 8, 03 Name: SOUTION ID#: PS.: A the
More information14. *14.8 CASTIGLIANO S THEOREM
*14.8 CASTIGLIANO S THEOREM Consider a body of arbitrary shape subjected to a series of n forces P 1, P 2, P n. Since external work done by forces is equal to internal strain energy stored in body, by
More informationMore Examples Of Generalized Coordinates
Slides of ecture 8 Today s Class: Review Of Homework From ecture 7 Hamilton s Principle More Examples Of Generalized Coordinates Calculating Generalized Forces Via Virtual Work /3/98 /home/djsegal/unm/vibcourse/slides/ecture8.frm
More informationConcepts. 3.1 Numerical Analysis. Chapter Numerical Analysis Scheme
Chapter 3 Concepts The objective of this work is to create a framework to implement multi-disciplinary finite element applications. Before starting, it is necessary to explain some basic concepts of the
More informationA study of the critical condition of a battened column and a frame by classical methods
University of South Florida Scholar Commons Graduate Theses and Dissertations Graduate School 003 A study of the critical condition of a battened column and a frame by classical methods Jamal A.H Bekdache
More informationMethod of Finite Elements I
Method of Finite Elements I PhD Candidate - Charilaos Mylonas HIL H33.1 and Boundary Conditions, 26 March, 2018 Institute of Structural Engineering Method of Finite Elements I 1 Outline 1 2 Penalty method
More informationChapter 2: Deflections of Structures
Chapter 2: Deflections of Structures Fig. 4.1. (Fig. 2.1.) ASTU, Dept. of C Eng., Prepared by: Melkamu E. Page 1 (2.1) (4.1) (2.2) Fig.4.2 Fig.2.2 ASTU, Dept. of C Eng., Prepared by: Melkamu E. Page 2
More informationInterpolation Functions for General Element Formulation
CHPTER 6 Interpolation Functions 6.1 INTRODUCTION The structural elements introduced in the previous chapters were formulated on the basis of known principles from elementary strength of materials theory.
More informationDevelopment of Truss Equations
CIVL 7/87 Chapter 3 - Truss Equations - Part /53 Chapter 3a Development of Truss Equations Learning Objectives To derive the stiffness matri for a bar element. To illustrate how to solve a bar assemblage
More informationBeam Models. Wenbin Yu Utah State University, Logan, Utah April 13, 2012
Beam Models Wenbin Yu Utah State University, Logan, Utah 843-4130 April 13, 01 1 Introduction If a structure has one of its dimensions much larger than the other two, such as slender wings, rotor blades,
More informationMathematical Properties of Stiffness Matrices
Mathematical Properties of Stiffness Matrices CEE 4L. Matrix Structural Analysis Department of Civil and Environmental Engineering Duke University Henri P. Gavin Fall, 0 These notes describe some of the
More informationFormulation of the displacement-based Finite Element Method and General Convergence Results
Formulation of the displacement-based Finite Element Method and General Convergence Results z Basics of Elasticity Theory strain e: measure of relative distortions u r r' y for small displacements : x
More informationLecture 7: The Beam Element Equations.
4.1 Beam Stiffness. A Beam: A long slender structural component generally subjected to transverse loading that produces significant bending effects as opposed to twisting or axial effects. MECH 40: Finite
More informationNON LINEAR BUCKLING OF COLUMNS Dr. Mereen Hassan Fahmi Technical College of Erbil
Abstract: NON LINEAR BUCKLING OF COLUMNS Dr. Mereen Hassan Fahmi Technical College of Erbil The geometric non-linear total potential energy equation is developed and extended to study the behavior of buckling
More informationAlgorithms for Scientific Computing
Algorithms for Scientific Computing Finite Element Methods Michael Bader Technical University of Munich Summer 2016 Part I Looking Back: Discrete Models for Heat Transfer and the Poisson Equation Modelling
More informationUsing MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,
Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business
More informationLecture 11: The Stiffness Method. Introduction
Introduction Although the mathematical formulation of the flexibility and stiffness methods are similar, the physical concepts involved are different. We found that in the flexibility method, the unknowns
More informationDiscrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method
131 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 21 (2008) Published online (http://hdl.handle.net/10114/1532) Discrete Analysis for Plate Bending Problems by Using
More informationMechanical Design in Optical Engineering. For a prismatic bar of length L in tension by axial forces P we have determined:
Deformation of Axial Members For a prismatic bar of length L in tension by axial forces P we have determined: σ = P A δ ε = L It is important to recall that the load P must act on the centroid of the cross
More informationENGI Second Order Linear ODEs Page Second Order Linear Ordinary Differential Equations
ENGI 344 - Second Order Linear ODEs age -01. Second Order Linear Ordinary Differential Equations The general second order linear ordinary differential equation is of the form d y dy x Q x y Rx dx dx Of
More informationWeighted Residual Methods
Weighted Residual Methods Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Table of Contents Problem definition. oundary-value Problem..................
More informationFig. 1. Circular fiber and interphase between the fiber and the matrix.
Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In
More informationM5 Simple Beam Theory (continued)
M5 Simple Beam Theory (continued) Reading: Crandall, Dahl and Lardner 7.-7.6 In the previous lecture we had reached the point of obtaining 5 equations, 5 unknowns by application of equations of elasticity
More informationAA242B: MECHANICAL VIBRATIONS
AA242B: MECHANICAL VIBRATIONS 1 / 57 AA242B: MECHANICAL VIBRATIONS Dynamics of Continuous Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations: Theory
More information