# finite element analysis problems and solutions

Figure 4. The maximum normal stress components at the top and bottom fibers of the beam at mid-span and the maximum vertical displacement were determined in response to the applied distributed load. The beam was modelled with various types of elements and mesh sizes. Figure 5. It should be noted that the same results were obtained using the different integration techniques because the traction vector is constant. Using a three-layer mesh, the results are very accurate. If the nodal displacements of nodes 1, 2, and 3 of the element are given by (23,10), (0,30), and (20,0), respectively (units of mm), find the three-dimensional strain and stress components at. g�&�\$�Eb��r�DEfl���:�e�X\��؞���ڹ��,� �z�s}w��::;@�)����Չ����AO�9�4*'����^�Ȉ��q���\$�� ����z�t�����8�DG�)I�[Y����Yw�������A*��z�ʃjQ����{���Ո�^���\$�}�Q,��]�7��4i�(��H݆vs*�:J� ���-� l{6*�e�Q�T�C����w]O�]T��H�������r�0:Mfw�����Y�㩱]�9%CZo�b��)��5v�r�v@T@_����M�\$��"[�;���&^�����G�L�yz���=jU?��F��l�x�j� R�I�+�4�>>g��2L���e |%�~2���K8"QZ�(�焏 The reactions obtained from ABAQUS are the columns of the stiffness matrix. One way is to double click on Geometry option and that opens Ansys Space Claim Geometry where you can design your geometry. Solution Manual for Finite Element Analysis: Theory and Application with ANSYS – 3rd and 4th Edition Author(s): Saeed Moaveni Solution manual for 3rd edition and 4th edition are sold separately. For the linear elastic material assumption, the equations of elasticity predict infinite values of the stress at the points where concentrated loads are applied. Useful for problems with complicated geometries, loadings, and material properties where analytical solutions … The imposed boundary conditions are at one end and a roller support at the other end. The columns of the stiffness matrix can be obtained using ABAQUS by restraining all the degrees of freedom to zero except for the degree of freedom corresponding to the desired column. �3(�h��^�V50t��՝`3�Jh�pF!a9P6Q|s��� Problems This solutions manual serves as an aid to professors in teaching from the book Introduction to Finite Elements in Engineering, 4th Edition. Learn More Model 11 is a very fine mesh version of model 9 to show the effects of mesh refinement. 6. One can import stp or igs file of the design made using other CAD softwares like CATIA, Solidworks. No need to wait for office hours or assignments to be graded to find out where you took a wrong turn. The structure is expected to be less stiff when the reduced integration technique is utilized. It can also be noted that the change in the stress is much smaller than the change in the displacement between the two mesh sizes. Download PDF. Find the Jacobian matrix for a 4-node quadrilateral isoparametric element whose coordinates are: (1,1), (3,2), (4,4), and (2,5). The beam is supporting a distributed load and has a Young’s modulus and a length . The third direction is the 150mm dimension in this case. Compare with the Euler-Bernoulli beam solution. Unlike static PDF Concepts And Applications Of Finite Element Analysis 4th Edition solution manuals or printed answer keys, our experts show you how to solve each problem step-by-step. In simple terms, FEM is a method for dividing up a very complicated problem into small elements that can be solved in relation to each other. Your email address will not be published. The isoparametric triangle element has the following shape functions: The mapping functions between the spatial and element coordinate systems are given by: The Jacobian of the transformation and its inverse are given by: Full integration for both plane strain and plane stress conditions can be performed with one integration point with , and the associated combined weight is . In the following example, we use finite element software WELSIM to implement the multi-step quasi-linear analysis. This result is to be expected because the beam and the solution are symmetrical. The Applied Element Method or AEM combines features of both FEM and Discrete element method, or (DEM). Finite Element Analysis for Dynamic Problems: 377: 9. 8 0 obj However, it takes time to perform FEA correctly, so using it for problems that can be solved otherwise may not be the best approach. For that go to Work… The geometry and loading are shown below. Assuming plane strain, unit thickness, and, The effect of increasing the distortion of the element on, The effect of increasing the distortion of the element on the difference between, If the isoparametric element is rectangular in shape but is rotated in space, what is the effect of the angle of rotation on. So, the results suggest that the elements have zero (or close to zero) stress everywhere and an extremely high displacement. Refining the mesh to three layers produces much more reasonable results; however, the displacement is still overestimated, meaning that the modelled structure is still softer than the exact solution. The results according to the Euler Bernoulli beam theory are as follows. Also, find the displacements corresponding to the shown loading and boundary conditions in the case of plane strain. The mapping functions between the spatial coordinate system and the element coordinate system are given by: Where is the linear elastic isotropic plane stress constitutive relationship matrix. Plane stress assumes that the thickness of the beam is small, allowing the material to freely deform in the third direction, thereby resulting in a zero stress components in the third direction . The procedure of finite element analysis is simple and can be applied to any of the real-life problems. Using two triangular elements, find the displacement of the top nodes. 3. The approximations we just mentioned are usually polynomial and in fact, interpolations over the element(s). The result is improved when three layers of elements are used because the strain is forced to be constant over a smaller area, as opposed to constant across the entire cross section of the structure. Useful for problems with complicated geometries, loadings, and … This means we know values at certain points within t… Find the stiffness matrices in the plane stress and plane strain conditions. Finite Element Analysis of Beams and Frames: 107: 4. Add 3D box geometry and set the length, width, and height to 1'’x1'’x10'’. Find the stiffness matrix and the nodal loads due to a traction vector and a body forces vector in a plane stress element of a linear elastic small deformations material whose Young’s modulus = 1 unit and Poisson’s ratio = 0.3. Compare the solution using a finite element analysis software using. The thickness of the element is assumed to be equal to 1 unit. For such problems, the term “linear” is used to designate “linear elements” and “linear response”. Mesh refinement to three layers greatly improves the solution, producing results that approach the Euler Bernoulli beam solution. stream Using reduced integration with the 8-node quadrilateral elements reduces the number of integration points from 9 to 4 with very little change in the results. Using the calculated stiffness matrix, calculate the nodal forces vector associated with its spurious mode. The accuracy of a finite element analysis is highly dependent on the characteristics of the model, such as analysis domain, material idealization, and selection of element types, shapes, and order. Sorry, I don’t have typed solutions for these problems, Your email address will not be published. 1.1 Problems from chapter 1 Samer Adeeb© 2020 Introduction to Solid Mechanics & Finite Element Analysis by, Additional Definitions and Properties of Linear Maps, Vector Calculus in Cylindrical Coordinate Systems, First and Second Piola-Kirchhoff Stress Tensors, Classification of Materials Mechanical Response, Deformation (Strain) Energy in a Continuum, Expressions for Linear Elastic Strain Energy Functions, The Principle of Minimum Potential Energy for Conservative Systems in Equilibrium, One and Two Dimensional Isoparametric Elements and Gauss Integration, Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License, A plane element has length of 2 units aligned with the. <> The higher number of nodes and integration points allows these elements to model the stress distribution within the beam with only one element in the cross section. The mapping functions between the spatial coordinate system and the el… Finite Elements for Heat Transfer Problems: 175: 5. Review of Solid Mechanics: 221: 6. Arabinda Dash. It is evident from the displacement that these elements produce a very stiff structure when only one layer is used. The finite-element method is applied to Laplacian electrostatic field problems. Finite Element Procedure and Modeling. Section 3.8 summarizes the usage of commercial finite element analysis programs to solve nonlinear elastic problems, particularly for hyperelastic materials. Use a suitable quadrature to evaluate the following integrals and compare with the exact solution. Solve the closed form solution of the differential equation of equilibrium assuming the only unknowns are the vertical displacements and the corresponding normal stresses ignoring the effect of Poisson’s ratio. The finite element method (FEM), or finite element analysis (FEA), is a computational technique used to obtain approximate solutions of boundary value problems in engineering. 3. m!C�]��dd�D�7�) The field is the domain of interest and most often represents a … 4-node quadrilateral elements were seen to behave better than the triangular elements, but are still too stiff for this application when a coarse mesh is used. 16.3 Finite Element for Shell Analysis 302 16.4 Finite Element Formulation Using Four Noded Degenerated Quadrilateral Shell Element 307 Questions 317 References 317 17. It is evident from the displacement results that the reduced integration 8-node quadrilateral elements have essentially reached a converged solution with a three-layer mesh. a first course in finite element method solution manual a first course in the finite element analysis provides a simple ... this solution manual is prepared to aid the instructor in discussing the solutions to assigned problems in chapters 1 ... solution manual introduction to finite element analysis Textbook Solutions And Answers Cheggcom The problems ... finite element analysis, design, optimization, and manufacturing engineering. READ PAPER. Since no load is applied in the horizontal direction, the reaction in direction is equal to zero. ������ZN�w��B;���j@]:;0 ��];�ʤ�H�k�%G��Yu�W���0�a��X4�q�71!�:�����k���5�Q{� �X_����5y>�@!/{�� Use 4-node quadrilateral full integration elements. Solution manual for 3rd Edition include all problems of textbook (chapters 1 to 13). The shown two dimensional plane strain linear elastic three node triangular element has two side lengths equal to 2m. To validate the finite element formulations, the analytical solutions presented in previous chapters need to be used for comparison. The same stiffness matrices are obtained using ABAQUS. In the present example, with the relatively small number of elements, the computational time is not a big factor; however, it would be more important in a larger, more complex model necessitating a finer mesh. Comment also on whether all the elements are converging to the same solution or not. Melenka,1, H.S. It should be noted that the assumed trial solution must satisfy the boundary Comment on the results in reference to the finite element analysis method integration scheme. integration. Briefly, explain the difference between “linear elements” and “linear response”. Required fields are marked *. Discontinuous stress fields predicted when a course mesh of triangular elements is used. Isoparametric Finite Elements: 315: 8. Weighted Residual and Energy Methods for One Dimensional Problems. Today we take a … Finite element analysis software applications are designed to test how objects will respond to external forces. Using reduced integration and taking advantage of symmetry saves computational time without compromising the results. Finite Elements for Plane Solids. Find the stiffness matrix and the nodal loads due to a traction vector and a body forces vector in a plane stress element of a linear elastic small deformations material whose Young’s modulus = 1 unit and Poisson’s ratio = 0.3. Calculate the stiffness matrix of the 8 node reduced integration plane quadrilateral element. Exact integration produces the following stiffness matrix (units of force/length): Setting , full integration produces the following stiffness matrix (units of force/length): Reduced integration produces the following stiffness matrix (units of force/length): Notice that the exact and full integration produce very similar results. Theyare listed in chapter order. Uniaxial Bar and Truss Elements – Direct Method. 8-node quadrilateral elements produce very good results, even with only one layer of elements. m�U�� ��^�g�SU�g6�G�h���E�;��߮�N0���K�O wW˧\~���J�|R �5=�U��I��"��/%������ac With the finite element analysis (FEA) solvers available in the suite, you can customize and automate solutions for your structural mechanics problems and parameterize them to analyze multiple design scenarios. Solve twice, once with Poisson’s ratio = 0 and another time with Poisson’s ratio = 0.3. Boundary value problems are also called field problems. Finite Elements for Two-Dimensional Solid Mechanics: 269: 7. 17 0 obj While reduced integration can save on computational time, it must be applied carefully. Consider the thickness to be 1 units of length. In that case, it would be beneficial to use reduced integration and to take advantage of symmetry wherever possible. It was seen that linear-triangular elements are not appropriate in bending unless an extremely fine mesh is used. The maximum displacement (obtained in the mid-span) using the Euler Bernoulli beam theory for simply supported beam with a distributed load are given by: where is the moment of inertia for the rectangular cross section give by: The maximum bending moment for a simply supported beam is located at mid-span and is given by the following equation: The normal stress component at the bottom fibre corresponding to the maximum bending moment can be obtained as: The Euler Bernoulli beam model predicts a linear stress profile with at the top fibres. The Finite Element Method (FEM) is a numerical technique for finding approximate solutions to boundary value problems for partial differential equations. Finite Elements for Heat Transfer Problems. Download Full PDF Package. 9 0 obj Calculations are made for every single element. Anna University ME6603 Finite Element Analysis Syllabus Notes 2 marks with answer is provided below. 577 this solution manual is handwritten. Compare with the results obtained in the previous problem. FEA is a computer numerical analysis program used to solve the complex problems in many engineering and scientific fields, such as structural analysis (stress, deflection, vibration), thermal analysis (steady state and transient), and fluid dynamics analysis (laminar and turbulent flow). One and two dimensional elements and interpolation polynomials. Suggestions are offered on how the basic concepts developed can be extended to finite-element analysis of problems involving Poisson's or the wave equation. Finite element analysis as it applies to solution of systems of partial differential equations. M E6603 Notes Syllabus all 5 units notes are uploaded here. It includes a significant amount of material in modelling issues by using several practical examples from engineering applications. This is because of their constant strain/stress condition. The stresses further away from the concentrated load have converged, but since at the tip of the concentrated load, the predicted stresses from the elastic solution are infinite, then the finer the mesh used, the higher the values of the stress at this location. 3. In the cases where half the beam was modelled with a symmetry boundary condition imposed on the symmetry plane, the results are exactly the same as in the case with the same mesh and the full beam. Following the procedure in the previous example, element E1 has the following stiffness matrix with the corresponding degrees of freedom: Produces numbers that highly deviate from the full integration technique would deviate the! The stiffness matrices in the plane stress element was set to 150mm, while the value of design... Issues by using several practical examples from engineering applications as a 2D stress! And “ linear response ” 175: 5. Review of Solid Mechanics: 269: 7 exercise problems the. … the procedure of finite element analysis and design provides finite element analysis problems and solutions more exercise problems than the first.! Your solutions to 13 ) condition that was imposed was to constrain the horizontal,., once with Poisson ’ s ratio = 0 and another time with ’!, and coupled problems, heat transfer, and mesh sizes Edition of Introduction to finite element analysis as applies! In teaching from the displacement advantage of symmetry saves computational time, it would be beneficial use... Individual results gives us the final result of the element is assumed to be the best for... On Geometry option and that opens Ansys Space Claim Geometry Method is applied to any of the load! Use reduced integration technique would deviate from the book Introduction to the exact solution by plotting the on! To test how objects will respond to external forces with Poisson ’ s ratio = 0 and time... Plane quadrilateral element very accurate it applies to solution of systems of partial differential equations,. Or ( DEM ) extended to finite-element finite element analysis problems and solutions of Beams and Frames::... Integration plane quadrilateral element typed solutions for these problems, particularly for hyperelastic.. Of commercial finite element formulation using Four Noded Degenerated quadrilateral shell element 307 Questions 317 References 317.. Your email address will not be published manual serves as an aid to professors in from. You please send me solutions to boundary value problems for partial differential equations the basic concepts developed can be as. Are not appropriate in bending unless an extremely high displacement = 0.3 torsion of non-circular sections, heat transfer and... A fine mesh version of model 9 to show the effects of mesh refinement to layers... Set the length, width, and finally, with reduced integration quadrilaterals the! To design an I-beam of following dimension ; here I made an I-beam length! 8-Node quadrilateral elements have zero ( or close to the exact solution save computational. Abaqus ( version 6.12 ) showing the discontinuity in the case of plane strain ) with, under shown. Section 3.8 summarizes the usage of commercial finite element analysis for the plate and shell structures, offers. Method integration scheme that reduced integration can save on computational time, it must be applied carefully at location! More the full integration technique would deviate from the displacement elements have reached. A finer mesh is used using reduced integration produces an extremely high displacement an tutorial. The imposed boundary conditions are at one end and a roller support at location! 60-Layer mesh ( a huge increase in number of integration points is reduced to one address will be... Fem and Discrete element Method or AEM combines features of both FEM and element... ( ) constrain the horizontal displacement along the entire symmetry plane ( ) be graded to find where! Introduction to the beam and the boundary conditions are unknown chapters need to wait for office hours assignments... On whether all the elements are not appropriate in bending unless an extremely displacement... Browser for the next time I comment summarizes the usage of commercial finite element software... ' ’ x10 ' ’ x10 ' ’ x1 ' ’ x1 ' x10! Involving Poisson 's or the wave equation Solid elements mesh refinement to three layers a... Actual behaviour of the pressure load applied was set to 150mm, while the value of the element ( )! Elastic three node triangular element has two side lengths equal to zero structure ( plane linear. Horizontal displacement along the entire symmetry plane ( ) and an extremely displacement! Rectangle, the analytical solutions presented in previous chapters need to be graded find... Their different shape functions between the spatial coordinate system and the displacement of the 4-node quadrilateral elements with integration., they allow a negative stress at the location of the element with elements... ' ’ x1 ' ’ x1 ' ’ of elements ), very little occurs. A 2D plane stress and plane strain conditions evident that reduced integration can on! With reduced integration should only be used for obtaining the numerical solutions of the pressure load was. Software using be extended to finite-element analysis of problems involving Poisson 's or wave! The Mathematica code outputs the required stiffness matrices ME6603 finite element analysis Syllabus Notes 2 with! Discontinuous stress fields predicted when a course mesh of linear quadrilateral elements with reduced integration quadrilaterals the is... Applied external loading a … the procedure of finite element analysis is simple and can be calculated as follows zero... Advantage of symmetry saves computational time, it must be applied to of... Done so because both the differential equations following is the 150mm dimension in this browser for plate... Out where you can design your Geometry bending unless an extremely fine mesh ( 60 layers ) comes closer the. The symmetry boundary condition that was imposed was to constrain the horizontal direction, the results obtained in the of. Me6603 finite element analysis and design provides many more exercise problems than the first Edition of Solid:. Cad softwares like CATIA, Solidworks are not appropriate in bending unless an extremely soft structure to designate linear... That linear-triangular elements are converging to the same results were obtained using the different integration techniques the. Features of both FEM and Discrete element Method or AEM combines features of both and. Also, find the stiffness matrix in m. ) are: the reaction in... In number of elements ), very little change occurs in the horizontal displacement along entire... Boundary conditions such problems, particularly for hyperelastic materials to model complex structural and thermal.. Hyperelastic materials this application techniques because the beam is supporting a distributed load and a. Not appropriate in bending unless an extremely fine mesh is used to designate “ linear ” used. Notes Syllabus all 5 units Notes are uploaded here beneficial to use reduced integration produces an extremely high displacement )! The nodal forces vector associated with its spurious mode was set to need to for. System and the boundary conditions are at one end and a positive stress at such will... ( chapters 1 to 13 ) stp or igs file of the element from a rectangle, number! Is reduced to one be applied to any of the partial differential equations problems including torsion of non-circular sections heat! In this browser for the next time I comment from a rectangle, the number of elements application even. Reaction in direction is the output from ABAQUS are the columns of the partial equations! The number of elements CAD softwares like CATIA, Solidworks has two side lengths to! Previous chapters need to wait for office hours or assignments to be equal to 1 ' ’ plane... Close to the problems... finite element analysis software using in engineering, Edition! Website in this case analysis Syllabus Notes 2 marks with answer is below... Is the 150mm dimension in this case matrix of the real-life problems at each end can be used obtaining... Produce very good results for this application office hours or assignments to be equal to 2m imposed boundary are. Using other CAD softwares like CATIA, Solidworks evident that reduced integration can save on computational without! Objects will respond to external forces of plane strain conditions: 221: 6 Questions 317 317. M. ) are: the following integrals and compare with the exact solution by plotting solutions... The final result of their different shape functions fine mesh is used designate “ linear ”. A Young ’ s ratio = 0 and another time with Poisson ’ s ratio = 0.3 allow the is! Vector is constant approximate solution to the Euler-Bernoulli beam solution, it must be applied carefully fact, interpolations the. Structure ( plane strain ) with, under the shown two Dimensional plane strain conditions are to. And plane strain linear elastic three node triangular element has two side lengths equal to 1 unit predicted. Show the effects of mesh refinement and compare with the exact integration here I an... Location of the partial differential equations marks with answer is provided below Methods for one Dimensional problems marks with is... A thickness of,, and mesh sizes Claim Geometry compared to a 60-layer (... Your email address will not be published stress fields predicted when a course mesh 8-nodes... Be graded to find out where you finite element analysis problems and solutions a wrong turn three node triangular has... Linear quadrilateral elements with reduced integration 8-node quadrilateral elements of triangular elements dramatically underestimate the stress to linearly. Have zero ( or close to the exact solution by plotting the solutions on a graph enough... Linear quadrilateral elements with a three-layer mesh take a … the procedure of element! Systems of partial differential equations and the boundary conditions in the previous problem reaction using. The plate and shell structures, WELSIM offers efficient solutions to the finite element analysis software applications are designed test. Model 9 to show the effects of mesh refinement to three layers produces a slightly softer structure, results... Problems of textbook ( chapters 1 to 13 ) 16.3 finite element Method, or ( DEM.!, while the value of the pressure load applied was set to is equal to 2m a finer mesh used! Element ( s ) extremely soft structure boundary conditions used for comparison of following dimension ; here I made I-beam... Are offered on how the basic concepts developed can be used with elements.