Nonlinear analysis in ADINA Structures
|
|
- Marcus Robertson
- 6 years ago
- Views:
Transcription
1 Nonlinear analysis in ADINA Structures Theodore Sussman, Ph.D. ADINA R&D, Inc,
2 Topics presented Types of nonlinearities Materially nonlinear only Geometrically nonlinear analysis Deformation-dependent loads Theory of nonlinear iterations full Newton iterations line searches limiting displacements ATS low-speed dynamics Stiffness stabilization TLA and TLA-S methods Birth-death Prescribed displacements Geometric imperfections ADINA R&D, Inc,
3 Nonlinear analysis in ADINA Structures, theory Much of the theory used by ADINA Structures can be found in the following reference: K.J. Bathe, Finite Element Procedures, 2 nd ed, Cambridge, MA, Klaus-Jürgen Bathe, 2014 and also in the papers by Bathe and co-workers. Much of the material presented in this lecture is taken from Chapter 6 of the textbook. ADINA R&D, Inc,
4 Materially nonlinear only Materially nonlinear only (MNO) Elastic-plastic material, viscoelastic material, etc. u nodal point displacements e strains e E elastic strains e P plastic strains P e change in plastic strains τ Cauchy stresses F nodal point forces Symbolically τ f E ( e ) ADINA R&D, Inc,
5 Geometric nonlinearities Geometric nonlinearities include Euler buckling P effect stress stiffening large rotations large strains 3 HL 3EI All of these effects are included by the use of a Lagrangian formulation (for example, total Lagrangian), with appropriate stress and strain measures, and with an appropriate material model. The details depend upon the element type and material model. We will discuss the nonlinearities due to deformation-dependent loading later. ADINA R&D, Inc,
6 Geometric nonlinearities - large disp, small strains Large displacements / small strains with an elastic orthotropic material: ε Green-Lagrange strains S 2nd Piola-Kirchhoff stresses S11 C11 C12 C S C C S 33 C S12 C S 23 symmetric C S C Symbolically τ f () e is replaced by S f () ε same stress-strain matrix as in linear analysis ADINA R&D, Inc,
7 Geometric nonlinearities - large disp, small strains This procedure includes all of the geometric nonlinearities mentioned above, with the exception of nonlinearities due to large strains. This procedure gives correct results provided that the strains are small, because, for small strains ε e and S τ. In fact, the definition of small strains is based on these conditions. Numerically, small strains are less than about 2%. If the strains are larger, the results will start to deviate from the results that you might expect. For example, the response in compression will not be the same as the response in tension. There is no real material for which S is linearly related to ε, when the strains are large. ADINA R&D, Inc,
8 Geometric nonlinearities - large disp, small strains Large displacements / small strains with an elastic-plastic material: ε E elastic Green-Lagrange strains ε P plastic Green-Lagrange strains Symbolically E τ f ( e ) is replaced by E S f ( ε ) As in the elastic-orthotropic case, this procedure gives correct results provided that the strains are small. We will give more details in the lecture on plasticity. ADINA R&D, Inc,
9 Geometric nonlinearities - large disp, large strains Large displacements / large strains with an elastic-plastic material, ULH formulation: X deformation gradient X P plastic deformation gradient X E elastic deformation gradient R E elastic rotation tensor U E elastic stretch tensor E ER elastic Hencky strains E Symbolically τ f ( e ) ER is replaced by τ f ( E ) τ rotated Cauchy stresses ADINA R&D, Inc,
10 Geometric nonlinearities - large disp, large strains The Hencky strains are a three-dimensional extension of the logarithmic (true) strains, and the Cauchy stresses are the true stresses (force per unit deformed area). So the material relationship is a good approximation for real materials. Therefore the ULH procedure can be used in both small and large strain situations. In addition, large rotations are accounted for. For some materials, there is also the possibility of using the ULJ formulation (updated Lagrangian Jaumann), but we will not discuss this further. ADINA R&D, Inc,
11 Geometric nonlinearities - large disp, large strains Large displacements / large strains with a rubberlike material. W S ε W strain energy density Because the material model is directly applicable to large strains, the procedure also is directly applicable to large strains. ADINA R&D, Inc,
12 Deformation-dependent loads When the model undergoes large displacements, the loads can also be geometrically nonlinear. If the model deforms significantly, there are two options: ADINA R&D, Inc,
13 Deformation-dependent loads By default, a pressure load is deformation-dependent if it is applied to a geometrically nonlinear element (large displacement kinematics). Deformation-dependent pressure loads can slow down the rate of convergence of the equilibrium iterations. Concentrated loads can also be deformation-dependent (follower loads). Centrifugal loads can also be deformation-dependent (load softening). ADINA R&D, Inc,
14 Equations of motion, statics The finite element equations of motion for ADINA Structures are written in terms of nodal point forces. The left superscript t+t means evaluated at time t+t. tt tt 0 R F Vector of external nodal point forces (from applied loads) Vector of internal nodal point forces (from element stresses) ADINA R&D, Inc,
15 Equations of motion The unknowns in ADINA Structures are the nodal point displacements and rotations. t t U F F U tt tt Lots of theory contained in this equation! R R U tt tt In many cases, the external forces are independent of the nodal point displacements (deformation-independent loads). ADINA R&D, Inc,
16 Equations of motion, linear special case If R is deformation-independent and F is a linear function of U, then the analysis is linear, and we can then write F K U tt tt and tt tt K U R But frequently F is a nonlinear function of U, or R is deformationdependent, so the equations of motion are nonlinear and we must solve them iteratively. tt tt 0 R F ADINA R&D, Inc,
17 Equations of motion, nonlinear case We suppose that the solution at time t is known, and we seek the solution at time t+t. The first approximation to the the solution at time t+t is U tt (0) t where (i) is an iteration counter (starting from 0). U t t (0) U Using we compute t t (0) F, t t (0) R and also t t (0) K where tt K (0) tt F U (0) ideally. ADINA R&D, Inc,
18 Equations of motion, nonlinear case Next, we compute the out-of-balance loads vector and a correction to the displacement vector K ΔU R F t t (0) t t (0) t t (0) to obtain a new trial displacement vector U = U +ΔU tt (1) tt (0) R F tt (0) tt (0) U tt ( i1) In general, given, ADINA Structures iterates as follows: Full Newton iterations without line searches Compute F tt ( i1), R tt ( i1) K U R F tt ( i1) tt ( i1) tt ( i1) U U U tt () i tt ( i1), K tt ( i1) ADINA R&D, Inc,
19 One DOF visualization of iterations The solution process can be visualized for one DOF in graphical form: Force t t R Slope t K Slope t t (1) K t t (1) F Iteration 0 Iteration 1 t R t U t t (1) U t t U Displacement ADINA R&D, Inc,
20 Convergence norms Norm of the out-of-balance force vector: R F tt ( i1) tt ( i1) 2 Norm of the incremental displacement vector: t t () i U 2 (for each of these quantities, there are two norms, one for the translational degrees of freedom and one for the rotational degrees of freedom) Out-of-balance energy norm U R F () it tt ( i1) tt ( i1) This norm has units of energy, but is not related to the total strain energy of the model or any other physical energy. ADINA R&D, Inc,
21 Convergence tolerances Force/moment Displacement/rotation R tt ( i1) tt ( i1) RNORM F 2 RTOL Energy t t () i U DNORM 2 DTOL U R F U R F () it tt ( i1) tt ( i1) (1) T tt (0) tt (0) ETOL ADINA R&D, Inc,
22 Line searches Instead of we use U U U tt () i tt ( i1) U U U tt () i tt ( i1) where is iteratively chosen to satisfy U R F U R F T tt () i tt () i T tt ( i1) tt ( i1) STOL The default is no line searches. ADINA R&D, Inc,
23 Limiting displacements Instead of U U U tt () i tt ( i1) FAC U we use in which FAC is chosen to make the largest incremental displacement component smaller than a certain value (MAXDISP in the TOLERANCES ITERATION command). The default value of MAXDISP (selected by the choice MAXDISP=0.0) depends on the model: - If the model is static and includes contact, MAXDISP is 1% of the largest model dimension. - Otherwise MAXDISP is set very large so the limiting displacement feature is not used. ADINA R&D, Inc,
24 Example convergence history Here is sample printout from the.out file (.log file for FSI). OUT-OF- NORM OF BALANCE OUT-OF-BALANCE NORM OF INCREMENTAL ENERGY FORCE MOMENT DISP. ROTN. CFORCE NODE-DOF NODE-DOF NODE-DOF NODE-DOF CFNORM MAX VALUE MAX VALUE MAX VALUE MAX VALUE ITE= E E E E E E Y 25-X 250-Z 121-Y 0.00E E E E E-05 ITE= E E E E E E X 45-Z 1-X 24-Y 0.00E E E E E-05 ITE= E E E E E E X 45-Y 56-Z 111-X 0.00E E E E E-05 ITE= E E E E E E Y 47-Z 1-X 114-Z 0.00E E E E E-06 ITE= E E E E E E Y 47-Z 1-X 114-Z 0.00E E E E E-06 ITE= E E E E E E Y 47-Z 1-X 114-Z 0.00E E E E E-06 CONVERGENCE RATIOS CONVERGENCE RATIOS OUT-OF-BALANCE LOAD FOR OUT-OF-BALANCE FOR INCREMENTAL VECTOR CALCULATION ENERGY FORCE DISP. CFORCE BETA RATIO MOMENT ROTN. (ITERNS) COMPARE WITH COMPARE WITH ETOL RTOL DTOL RCTOL (NOT USED)(NOT USED) 1.00E E E E E E E E E E E E E E+00 ( 1) 1.70E E E E E E E E+00 ( 1) 7.60E E E E E E E E+00 ( 3) 7.20E E E E E E E E+00 ( 3) 2.20E E E E E E E E+00 ( 1) ADINA R&D, Inc,
25 Convergence history view (not the same data as shown in the previous slide) ADINA R&D, Inc,
26 Iteration patterns, quadratic convergence Example: ITE= E E E E E E X Y E E E E E+00 ITE= E E E E E E Y 0 21-Z E E E E E+00 ITE= E E E E E E Z 0 21-Z E E E E E+00 ITE= E E E E E E X Y E E E E E+00 ITE= E E E E E E X Z E E E E E+00 This is the fastest convergence. Near convergence, the current force norm is roughly the square of the previous force norm. Quadratic convergence is only observed if the stiffness matrix is truly tangent. Observing quadratic convergence gives some assurance that the solution has truly converged. ADINA R&D, Inc,
27 Iteration patterns, linear convergence Example: ITE= E E E E E E X 22-Y 34-X 16-Z 4.60E E E E E-05 Near convergence, the current force norm is a factor times the previous force norm. ITE= E E E E E E Z 55-Y 66-X 23-X 4.60E E E E E-06 ITE= E E E E E E X 56-Y 17-X 45-Z 4.60E E E E E-07 ITE= E E E E E E X 56-Y 17-X 45-Z 4.60E+01) 2.50E E E E-08 Linear convergence is observed if the stiffness matrix is not truly tangent. Deformationdependent loads Certain material models ADINA R&D, Inc,
28 Iteration patterns, divergence Example: ITE= E E E E E E X 2314-Y 1232-X 3421-X 0.00E E E E E-03 ITE= E E E E E E X 3238-Y 1523-X 1962-Y 1.20E E E E E-02 ITE= E E E E E E X 456-X 132-Y 429-Z 3.60E E E E E-02 ITE= E E E E E E X 962-Y 434-Z 347-X 6.60E E E E E-02 ITE= E E E E E E X 562-Y 983-X 4713-Z 1.20E E E E E-03 ADINA R&D, Inc,
29 What can we do when the solution doesn t converge? Use more iterations (only if the solution was converging already) Use line searches Use smaller load increments Manually reduce the time step Use the ATS method (with or without low-speed dynamics) Use physical intuition to determine and resolve the problem. ADINA R&D, Inc,
30 ATS (automatic time-stepping method) When ADINA Structures uses the ATS method, ADINA Structures cuts the step size increment and retries the step, when there is no convergence. For a 1 DOF system Force t t R 1 2 t t R t R t t U 1 t 2 U tt U Displacement ADINA R&D, Inc,
31 ATS example, graphical representation t t Original time step, t NO CONVERGENCE ADINA R&D, Inc,
32 ATS example, graphical representation t t/2 t First subdivided time step NO CONVERGENCE Original time step, t NO CONVERGENCE ADINA R&D, Inc,
33 ATS example, graphical representation t t /4 t/2 t First subdivided time step NO CONVERGENCE Second subdivided time step CONVERGED Original time step, t NO CONVERGENCE ADINA R&D, Inc,
34 ATS example, graphical representation t t/2 t t /4 CONVERGED t Same as the Time Step Prior to Subdivision ADINA R&D, Inc,
35 ATS example, graphical representation t t/2 t t /4 CONVERGED t /4 t Same as the Time Step that Gave Convergence Same as the Time Step Prior to Subdivision ADINA R&D, Inc,
36 ATS example, graphical representation t t/2 t t /4 CONVERGED t /4 t Same as the Time Step that Gave Convergence (t t /4) Same as the Time Step Prior to Subdivision Return to Original Time Step Specified ADINA R&D, Inc,
37 ATS notes The ATS method only saves the solution for the user-specified time steps, and for the last converged solution if this solution does not correspond to a user-specified time step. After the solution is obtained for a user-specified time step, the program chooses the next time step according to one of the following options: Same as time step prior to subdivision (RESTORE=YES) Same as time step that gave convergence (RESTORE=NO) Solution time matches original next solution time (always used for iterative FSI) (RESTORE=ORIGINAL) Automatic, based on problem characteristics (RESTORE=AUTOMATIC) If the ATS method does not converge, the program can optionally turn on low-speed dynamics. An additional option specifies the time duration during which lowspeed dynamics is activated; after this period of time, the program deactivates low-speed dynamics. ADINA R&D, Inc,
38 Low-speed dynamics The ATS method includes a special low-speed dynamics option in static analysis. When low-speed dynamics is active, ADINA Structures includes mass and damping effects in the (otherwise static) analysis. M U C U K U R F tt () i tt () i tt ( i1) () i tt tt ( i1) where C M K,, user-specified Reasons to use low-speed dynamics: Allow for solution of models with rigid-body modes. Allow for solutions of models with local snap-through or buckling instabilities Dynamics smooths the response, hence convergence is easier. ADINA R&D, Inc,
39 Low-speed dynamics, continued The mass matrix is evaluated from the density and the low-speed inertia factor (default = 1). The damping matrix is evaluated using Rayleigh damping: C M K where and are user-specified parameters (defaults 0 and 10-4 ). Two basic options for damping specification: Damping factors applied to entire model (in which case = 0). Damping factors applied to individual element groups, using the Rayleigh damping options. ADINA R&D, Inc,
40 Low-speed dynamics, continued When low-speed dynamics is used, the time step size will influence the results. can be interpreted as a decay time (time constant of exponential decay). Increasing increases the decay time. It is recommended that either The time step size be at least 10, or After the solution is obtained, the load be held constant for a period of time at least 10 to allow the dynamic solution to decay away. ADINA R&D, Inc,
41 Low-speed dynamics, continued The solution indicators can be used to assess the magnitudes of the inertia and damping forces: External force indicator Damping force indicator Inertia force indicator I I R U EF I R U DF I R I D U ==================================================================================== S O L U T I O N A C C U R A C Y I N D I C A T O R S FOR LOW SPEED DYNAMICS, CONTACT DAMPING, SHELL DRILLING AND STIFFNESS STABILIZATION EXTERNAL DRILLING DAMPING INERTIA CONTACT STIFFNESS FORCES FORCES FORCES FORCES DAMP.FORCES STABIL E E E E E % of ext.forces ====================================================================================== ADINA R&D, Inc,
42 Stiffness stabilization Sometimes the structure is unstable (has rigid body modes) at the start of the analysis. The equation solver cannot give a solution when there are rigid-body modes. Rigid-body modes cause zero pivots in the factorized stiffness matrix. ADINA R&D, Inc,
43 Zero pivots can be removed using stiffness stabilization Scale all diagonal stiffness terms (excluding contact diagonals) without modifying the right-hand-side load vector K 1 K ii stab ii Physically, attaches weak springs to all degrees of freedom, but in such a way to not affect the solution in nonlinear analysis. ADINA R&D, Inc,
44 Stiffness stabilization does not affect the solution in nonlinear analysis Stiffness stabilization only modifies K, element nodal forces F are not modified i1 i i1 i1 K U R F tt tt tt converges when out-of-balance load is sufficiently small Hence, the converged solution is same as without stabilization. However, as K is modified, rate of convergence might be worsened, so that more iterations are required. ADINA R&D, Inc,
45 TLA, TLA-S methods Goal: apply a load in nonlinear static analysis without excessive user input. User applies total load level (TLA = Total Load Application). ADINA Structures automatically applies the load using a ramp time function. ADINA Structures chooses the number of time steps. ADINA Structures increases/decreases the time step size depending upon the number of equilibrium iterations in the previous time step. By default: 1 st time step has size 1/50 th of the total time. Max number of equilibrium iterations = 30. Max number of time step subdivisions used in ATS method = 64 Max limiting incremental displacement = 5% model dimension ADINA R&D, Inc,
46 TLA-S = TLA with stabilization low-speed dynamics contact damping stiffness stabilization TLA, TLA-S methods At the end of the equilibrium iterations for each step, ADINA Structures prints the solution indicators: ================================================================================ S O L U T I O N A C C U R A C Y I N D I C A T O R S FOR LOW SPEED DYNAMICS (LSD), CONTACT DAMPING, AND SHELL DRILLING STIFFNESS EXTERNAL DRILLING LSD DAMPING LSD INERTIA CONTACT FORCES FORCES FORCES FORCES DAMP.FORCES E E E E % of ext.forces ================================================================================ If indicators are < 1%, solution is accurate ADINA R&D, Inc,
47 Element birth/death Elements can be born, or can die, during the solution. You can specify the birth and death times as part of the model definition: Birth Death Birth/death Elements can die as a result of material rupture. ADINA R&D, Inc,
48 Element birth/death You should not specify the birth and death times to be exactly equal to solution times. Solution times 1, 2, etc. If you want the element to be born at time 2, birth time can be between and If you want the element to die at time 3, death time can be between and The element is assumed to be stress-free when it is born. Birth time = to 1.999, element is stress free at time 1. ADINA R&D, Inc,
49 The dying element mass and stiffness can be removed from the model over a period of time. The element starts to die at the userspecified time, and is completely removed from the model after an additional period of time. Element death decay time Element death decay time ADINA R&D, Inc,
50 Prescribed displacements Displacements can be prescribed as a function of time (similar to other load types). When the arrival time option is used, the prescribed displacement becomes active only after the arrival time. Arrival time is 10 Relative prescribed displacement option. Total prescribed displacement = preexisting displacement + relative prescribed displacement. Arrival time is 10 ADINA R&D, Inc,
51 Prescribed displacements - time unloading Time unloading with zero user-specified unloading force: Unloading time is time 2.5 Time unloading with nonzero user-specified unloading force R u : Unloading time is time 2.5 ADINA R&D, Inc,
52 Prescribed displacements - force unloading Unloading starts when the specified force is exceeded: User-specified unloading force = 1000 ADINA R&D, Inc,
53 Geometric imperfections Before the analysis begins, the nodal coordinates can be updated by imperfection displacements. The "imperfect" nodal coordinates are used during the analysis, and the computed displacements are based on the imperfect nodal coordinates. The imperfect mesh is stress-free at the start of the analysis. ADINA R&D, Inc,
54 Geometric imperfections in collapse analysis Buckling modes come from a linearized buckling analysis. Each buckling mode can be independently scaled. ADINA R&D, Inc,
Contact analysis - theory and concepts
Contact analysis - theory and concepts Theodore Sussman, Ph.D. ADINA R&D, Inc, 2016 1 Overview Review of contact concepts segments, surfaces, groups, pairs Interaction of contactor nodes and target segments
More informationA Demonstrative Computer Session Using ADINA- Nonlinear Analysis
Topic 22 A Demonstrative Computer Session Using ADINA- Nonlinear Analysis Contents: Use of ADINA for elastic-plastic analysis of a plate with a hole Computer laboratory demonstration-part II Selection
More informationCourse in. Geometric nonlinearity. Nonlinear FEM. Computational Mechanics, AAU, Esbjerg
Course in Nonlinear FEM Geometric nonlinearity Nonlinear FEM Outline Lecture 1 Introduction Lecture 2 Geometric nonlinearity Lecture 3 Material nonlinearity Lecture 4 Material nonlinearity it continued
More informationGlossary. Glossary of Symbols. Glossary of Roman Symbols Glossary of Greek Symbols. Contents:
Glossary Glossary of Symbols Contents: Glossary of Roman Symbols Glossary of Greek Symbols Glossary G-l Glossary of Roman Symbols The Euclidean norm or "two-norm." For a vector a The Mooney-Rivlin material
More informationMITOCW MITRES2_002S10nonlinear_lec15_300k-mp4
MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
More informationCAEFEM v9.5 Information
CAEFEM v9.5 Information Concurrent Analysis Corporation, 50 Via Ricardo, Thousand Oaks, CA 91320 USA Tel. (805) 375 1060, Fax (805) 375 1061 email: info@caefem.com or support@caefem.com Web: http://www.caefem.com
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 informationNon-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 17, 2017, Lesson 5
Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 17, 2017, Lesson 5 1 Politecnico di Milano, February 17, 2017, Lesson 5 2 Outline
More informationNonlinear Theory of Elasticity. Dr.-Ing. Martin Ruess
Nonlinear Theory of Elasticity Dr.-Ing. Martin Ruess geometry description Cartesian global coordinate system with base vectors of the Euclidian space orthonormal basis origin O point P domain of a deformable
More informationThe Finite Element Method II
[ 1 The Finite Element Method II Non-Linear finite element Use of Constitutive Relations Xinghong LIU Phd student 02.11.2007 [ 2 Finite element equilibrium equations: kinematic variables Displacement Strain-displacement
More informationNon-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises
Non-linear and time-dependent material models in Mentat & MARC Tutorial with Background and Exercises Eindhoven University of Technology Department of Mechanical Engineering Piet Schreurs July 7, 2009
More informationDemonstrative Exam~le Solutions
Topic 12 Demonstrative Exam~le Solutions in Static Analysis Contents: Analysis of various problems to demonstrate, study, and evaluate solution methods in statics Example analysis: Snap-through of an arch
More informationSolution of Nonlinear Dynamic Response-Part I
Topic 13 Solution of Nonlinear Dynamic Response-Part I Contents: Basic procedure of direct integration The explicit central difference method, basic equations, details of computations performed, stability
More informationAlvaro F. M. Azevedo A. Adão da Fonseca
SECOND-ORDER SHAPE OPTIMIZATION OF A STEEL BRIDGE Alvaro F. M. Azevedo A. Adão da Fonseca Faculty of Engineering University of Porto Portugal 16-18 March 1999 OPTI 99 Orlando - Florida - USA 1 PROBLEM
More informationLaboratory 4 Topic: Buckling
Laboratory 4 Topic: Buckling Objectives: To record the load-deflection response of a clamped-clamped column. To identify, from the recorded response, the collapse load of the column. Introduction: Buckling
More informationFinite Element Modelling with Plastic Hinges
01/02/2016 Marco Donà Finite Element Modelling with Plastic Hinges 1 Plastic hinge approach A plastic hinge represents a concentrated post-yield behaviour in one or more degrees of freedom. Hinges only
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationUsing Energy History Data to Obtain Load vs. Deflection Curves from Quasi-Static Abaqus/Explicit Analyses
Using Energy History Data to Obtain Load vs. Deflection Curves from Quasi-Static Abaqus/Explicit Analyses Brian Baillargeon, Ramesh Marrey, Randy Grishaber 1, and David B. Woyak 2 1 Cordis Corporation,
More informationUse of Elastic Constitutive Relations in Total Lagrangian Formulation
Topic 15 Use of Elastic Constitutive Relations in Total Lagrangian Formulation Contents: Basic considerations in modeling material response Linear and nonlinear elasticity Isotropic and orthotropic materials
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationTheoretical Manual Theoretical background to the Strand7 finite element analysis system
Theoretical Manual Theoretical background to the Strand7 finite element analysis system Edition 1 January 2005 Strand7 Release 2.3 2004-2005 Strand7 Pty Limited All rights reserved Contents Preface Chapter
More informationLeaf Spring (Material, Contact, geometric nonlinearity)
00 Summary Summary Nonlinear Static Analysis - Unit: N, mm - Geometric model: Leaf Spring.x_t Leaf Spring (Material, Contact, geometric nonlinearity) Nonlinear Material configuration - Stress - Strain
More informationOn Nonlinear Buckling and Collapse Analysis using Riks Method
Visit the SIMULIA Resource Center for more customer examples. On Nonlinear Buckling and Collapse Analysis using Riks Method Mingxin Zhao, Ph.D. UOP, A Honeywell Company, 50 East Algonquin Road, Des Plaines,
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015 Institute of Structural Engineering Method of Finite Elements II 1 Constitutive
More informationGEO E1050 Finite Element Method Autumn Lecture. 9. Nonlinear Finite Element Method & Summary
GEO E1050 Finite Element Method Autumn 2016 Lecture. 9. Nonlinear Finite Element Method & Summary To learn today The lecture should give you overview of how non-linear problems in Finite Element Method
More informationON EFFECTIVE IMPLICIT TIME INTEGRATION IN ANALYSIS OF FLUID-STRUCTURE PROBLEMS
SHORT COMMUNICATIONS 943 ON EFFECTIVE IMPLICIT TIME INTEGRATION IN ANALYSIS OF FLUID-STRUCTURE PROBLEMS KLAUS-JURGEN BATHE? AND VUAY SONNADS Dcpaflment of Mechanical Engineering, Massachusetts Institute
More informationBiomechanics. Soft Tissue Biomechanics
Biomechanics cross-bridges 3-D myocardium ventricles circulation Image Research Machines plc R* off k n k b Ca 2+ 0 R off Ca 2+ * k on R* on g f Ca 2+ R0 on Ca 2+ g Ca 2+ A* 1 A0 1 Ca 2+ Myofilament kinetic
More informationMITOCW MITRES2_002S10nonlinear_lec05_300k-mp4
MITOCW MITRES2_002S10nonlinear_lec05_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
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 informationNonlinear Buckling Prediction in ANSYS. August 2009
Nonlinear Buckling Prediction in ANSYS August 2009 Buckling Overview Prediction of buckling of engineering structures is a challenging problem for several reasons: A real structure contains imperfections
More informationMSC Nastran N is for NonLinear as in SOL400. Shekhar Kanetkar, PhD
MSC Nastran N is for NonLinear as in SOL400 Shekhar Kanetkar, PhD AGENDA What is SOL400? Types of Nonlinearities Contact Defining Contact Moving Rigid Bodies Friction in Contact S2S Contact CASI Solver
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 informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationFLEXIBILITY METHOD FOR INDETERMINATE FRAMES
UNIT - I FLEXIBILITY METHOD FOR INDETERMINATE FRAMES 1. What is meant by indeterminate structures? Structures that do not satisfy the conditions of equilibrium are called indeterminate structure. These
More informationA consistent dynamic finite element formulation for a pipe using Euler parameters
111 A consistent dynamic finite element formulation for a pipe using Euler parameters Ara Arabyan and Yaqun Jiang Department of Aerospace and Mechanical Engineering, University of Arizona, Tucson, AZ 85721,
More informationComputational Inelasticity FHLN05. Assignment A non-linear elasto-plastic problem
Computational Inelasticity FHLN05 Assignment 2017 A non-linear elasto-plastic problem General instructions A written report should be submitted to the Division of Solid Mechanics no later than October
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 informationModelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases
Modelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases Georg Haasemann Conrad Kloß 1 AIMCAL Conference 2016 MOTIVATION Wrinkles in web handling system Loss of
More informationNonlinear Static - 1D Plasticity - Isotropic and Kinematic Hardening Walla Walla University
Nonlinear Static - 1D Plasticity - Isotropic and Kinematic Hardening Walla Walla University Professor Louie L. Yaw c Draft date January 18, 2017 Copyright is reserved. Individual copies may be made for
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 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 informationSTRAIN ASSESSMENT USFOS
1 STRAIN ASSESSMENT IN USFOS 2 CONTENTS: 1 Introduction...3 2 Revised strain calculation model...3 3 Strain predictions for various characteristic cases...4 3.1 Beam with concentrated load at mid span...
More informationMulti Linear Elastic and Plastic Link in SAP2000
26/01/2016 Marco Donà Multi Linear Elastic and Plastic Link in SAP2000 1 General principles Link object connects two joints, i and j, separated by length L, such that specialized structural behaviour may
More informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More information202 Index. failure, 26 field equation, 122 force, 1
Index acceleration, 12, 161 admissible function, 155 admissible stress, 32 Airy's stress function, 122, 124 d'alembert's principle, 165, 167, 177 amplitude, 171 analogy, 76 anisotropic material, 20 aperiodic
More informationContinuation methods for non-linear analysis
Continuation methods for non-linear analysis FR : Méthodes de pilotage du chargement Code_Aster, Salome-Meca course material GNU FDL licence (http://www.gnu.org/copyleft/fdl.html) Outline Definition of
More informationGEOMETRIC NONLINEAR ANALYSIS
GEOMETRIC NONLINEAR ANALYSIS The approach for solving problems with geometric nonlinearity is presented. The ESAComp solution relies on Elmer open-source computational tool [1] for multiphysics problems.
More informationComparison of Models for Finite Plasticity
Comparison of Models for Finite Plasticity A numerical study Patrizio Neff and Christian Wieners California Institute of Technology (Universität Darmstadt) Universität Augsburg (Universität Heidelberg)
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 informationSoftware Verification
EXAMPLE 1-026 FRAME MOMENT AND SHEAR HINGES EXAMPLE DESCRIPTION This example uses a horizontal cantilever beam to test the moment and shear hinges in a static nonlinear analysis. The cantilever beam has
More informationLecture 4 Implementing material models: using usermat.f. Implementing User-Programmable Features (UPFs) in ANSYS ANSYS, Inc.
Lecture 4 Implementing material models: using usermat.f Implementing User-Programmable Features (UPFs) in ANSYS 1 Lecture overview What is usermat.f used for? Stress, strain and material Jacobian matrix
More informationComputational Inelasticity FHLN05. Assignment A non-linear elasto-plastic problem
Computational Inelasticity FHLN05 Assignment 2018 A non-linear elasto-plastic problem General instructions A written report should be submitted to the Division of Solid Mechanics no later than November
More informationNUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS
NUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS ABSTRACT : P Mata1, AH Barbat1, S Oller1, R Boroschek2 1 Technical University of Catalonia, Civil Engineering
More informationComputational Simulation of Dynamic Response of Vehicle Tatra T815 and the Ground
IOP Conference Series: Earth and Environmental Science PAPER OPEN ACCESS Computational Simulation of Dynamic Response of Vehicle Tatra T815 and the Ground To cite this article: Jozef Vlek and Veronika
More informationFinite Element Solutions for Geotechnical Engineering
Release Notes Release Date: July, 2015 Product Ver.: GTSNX 2015 (v2.1) Integrated Solver Optimized for the next generation 64-bit platform Finite Element Solutions for Geotechnical Engineering Enhancements
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 informationProject. First Saved Monday, June 27, 2011 Last Saved Wednesday, June 29, 2011 Product Version 13.0 Release
Project First Saved Monday, June 27, 2011 Last Saved Wednesday, June 29, 2011 Product Version 13.0 Release Contents Units Model (A4, B4) o Geometry! Solid Bodies! Parts! Parts! Body Groups! Parts! Parts
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 informationChapter 2 Examples of Optimization of Discrete Parameter Systems
Chapter Examples of Optimization of Discrete Parameter Systems The following chapter gives some examples of the general optimization problem (SO) introduced in the previous chapter. They all concern the
More informationSettlement and Bearing Capacity of a Strip Footing. Nonlinear Analyses
Settlement and Bearing Capacity of a Strip Footing Nonlinear Analyses Outline 1 Description 2 Nonlinear Drained Analysis 2.1 Overview 2.2 Properties 2.3 Loads 2.4 Analysis Commands 2.5 Results 3 Nonlinear
More informationPOST-BUCKLING BEHAVIOUR OF IMPERFECT SLENDER WEB
Engineering MECHANICS, Vol. 14, 007, No. 6, p. 43 49 43 POST-BUCKLING BEHAVIOUR OF IMPERFECT SLENDER WEB Martin Psotný, Ján Ravinger* The stability analysis of slender web loaded in compression is presented.
More informationLecture 15 Strain and stress in beams
Spring, 2019 ME 323 Mechanics of Materials Lecture 15 Strain and stress in beams Reading assignment: 6.1 6.2 News: Instructor: Prof. Marcial Gonzalez Last modified: 1/6/19 9:42:38 PM Beam theory (@ ME
More informationA NOTE ON RELATIONSHIP BETWEEN FIXED-POLE AND MOVING-POLE APPROACHES IN STATIC AND DYNAMIC ANALYSIS OF NON-LINEAR SPATIAL BEAM STRUCTURES
European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS 212) J. Eberhardsteiner et.al. (eds.) Vienna, Austria, September 1-14, 212 A NOTE ON RELATIONSHIP BETWEEN FIXED-POLE
More informationEngineering Sciences 241 Advanced Elasticity, Spring Distributed Thursday 8 February.
Engineering Sciences 241 Advanced Elasticity, Spring 2001 J. R. Rice Homework Problems / Class Notes Mechanics of finite deformation (list of references at end) Distributed Thursday 8 February. Problems
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 informationStructural Mechanics Module
Structural Mechanics Module User s Guide VERSION 4.3 Structural Mechanics Module User s Guide 1998 2012 COMSOL Protected by U.S. Patents 7,519,518; 7,596,474; and 7,623,991. Patents pending. This Documentation
More informationFinite Element Method
Finite Element Method Finite Element Method (ENGC 6321) Syllabus Objectives Understand the basic theory of the FEM Know the behaviour and usage of each type of elements covered in this course one dimensional
More informationProblem 1: Find the Equation of Motion from the static equilibrium position for the following systems: 1) Assumptions
Problem 1: Find the Equation of Motion from the static equilibrium position for the following systems: 1) Assumptions k 2 Wheels roll without friction k 1 Motion will not cause block to hit the supports
More informationSculptural Form finding with bending action
Sculptural Form finding with bending action Jens Olsson 1, Mats Ander 2, Al Fisher 3, Chris J K Williams 4 1 Chalmers University of Technology, Dept. of Architecture, 412 96 Göteborg, jens.gustav.olsson@gmail.com
More informationSoftware Verification
EXAMPLE 6-003 LINK GAP ELEMENT PROBLEM DESCRIPTION This example uses a single-bay, single-story rigid frame to test the gap link element. This link element carries compression loads only; it has zero stiffness
More informationNonlinear Analysis Of An EPDM Hydraulic Accumulator Bladder. Richard Kennison, Race-Tec
Nonlinear Analysis Of An EPDM Hydraulic Accumulator Bladder Richard Kennison, Race-Tec Agenda Race-Tec Overview Accumulator Experimental Testing Material Testing Numerical Analysis: 1. Linear Buckling
More information3. Overview of MSC/NASTRAN
3. Overview of MSC/NASTRAN MSC/NASTRAN is a general purpose finite element analysis program used in the field of static, dynamic, nonlinear, thermal, and optimization and is a FORTRAN program containing
More informationStretching of a Prismatic Bar by its Own Weight
1 APES documentation (revision date: 12.03.10) Stretching of a Prismatic Bar by its Own Weight. This sample analysis is provided in order to illustrate the correct specification of the gravitational acceleration
More informationComputational non-linear structural dynamics and energy-momentum integration schemes
icccbe 2010 Nottingham University Press Proceedings of the International Conference on Computing in Civil and Building Engineering W Tizani (Editor) Computational non-linear structural dynamics and energy-momentum
More informationDirect calculation of critical points in parameter sensitive systems
Direct calculation of critical points in parameter sensitive systems Behrang Moghaddasie a, Ilinca Stanciulescu b, a Department of Civil Engineering, Ferdowsi University of Mashhad, P.O. Box 91775-1111,
More informationDYNAMIC RESPONSE OF THIN-WALLED GIRDERS SUBJECTED TO COMBINED LOAD
DYNAMIC RESPONSE OF THIN-WALLED GIRDERS SUBJECTED TO COMBINED LOAD P. WŁUKA, M. URBANIAK, T. KUBIAK Department of Strength of Materials, Lodz University of Technology, Stefanowskiego 1/15, 90-924 Łódź,
More informationLINEAR AND NONLINEAR BUCKLING ANALYSIS OF STIFFENED CYLINDRICAL SUBMARINE HULL
LINEAR AND NONLINEAR BUCKLING ANALYSIS OF STIFFENED CYLINDRICAL SUBMARINE HULL SREELATHA P.R * M.Tech. Student, Computer Aided Structural Engineering, M A College of Engineering, Kothamangalam 686 666,
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 informationAs an example, the two-bar truss structure, shown in the figure below will be modelled, loaded and analyzed.
Program tr2d 1 FE program tr2d The Matlab program tr2d allows to model and analyze two-dimensional truss structures, where trusses are homogeneous and can behave nonlinear. Deformation and rotations can
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 informationGeometrically exact beam dynamics, with and without rotational degree of freedom
ICCM2014 28-30 th July, Cambridge, England Geometrically exact beam dynamics, with and without rotational degree of freedom *Tien Long Nguyen¹, Carlo Sansour 2, and Mohammed Hjiaj 1 1 Department of Civil
More informationLECTURE 14: DEVELOPING THE EQUATIONS OF MOTION FOR TWO-MASS VIBRATION EXAMPLES
LECTURE 14: DEVELOPING THE EQUATIONS OF MOTION FOR TWO-MASS VIBRATION EXAMPLES Figure 3.47 a. Two-mass, linear vibration system with spring connections. b. Free-body diagrams. c. Alternative free-body
More informationObtaining a Converged Solution with Abaqus
Obtaining a Converged Solution with Abaqus Day 1 Lecture 1 Workshop 1 Lecture 2 Lecture 3 Workshop 2 Lecture 4 Workshop 3 Lecture 5 Workshop 4 Workshop 5 Introduction to Nonlinear FEA Nonlinear Spring
More information2D Kirchhoff Thin Beam Elements
2D Kirchhoff Thin Beam Elements General Element Name Y,v X,u BM3 2 3 1 Element Group Element Subgroup Element Description Number Of Nodes 3 Freedoms Node Coordinates Geometric Properties Kirchhoff Parabolically
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 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 informationENGN2340 Final Project: Implementation of a Euler-Bernuolli Beam Element Michael Monn
ENGN234 Final Project: Implementation of a Euler-Bernuolli Beam Element Michael Monn 12/11/13 Problem Definition and Shape Functions Although there exist many analytical solutions to the Euler-Bernuolli
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 informationAdaptive Analysis of Bifurcation Points of Shell Structures
First published in: Adaptive Analysis of Bifurcation Points of Shell Structures E. Ewert and K. Schweizerhof Institut für Mechanik, Universität Karlsruhe (TH), Kaiserstraße 12, D-76131 Karlsruhe, Germany
More informationSolution of the Nonlinear Finite Element Equations in Static Analysis Part II
Topic 11 Solution of the Nonlinear Finite Element Equations in Static Analysis Part II Contents: Automatic load step incrementation for collapse and post-buckling analysis Constant arc-length and constant
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 informationLIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS
5 th GRACM International Congress on Computational Mechanics Limassol, 29 June 1 July, 2005 LIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS G.A. Drosopoulos I, G.E.
More informationBy drawing Mohr s circle, the stress transformation in 2-D can be done graphically. + σ x σ y. cos 2θ + τ xy sin 2θ, (1) sin 2θ + τ xy cos 2θ.
Mohr s Circle By drawing Mohr s circle, the stress transformation in -D can be done graphically. σ = σ x + σ y τ = σ x σ y + σ x σ y cos θ + τ xy sin θ, 1 sin θ + τ xy cos θ. Note that the angle of rotation,
More informationM.S Comprehensive Examination Analysis
UNIVERSITY OF CALIFORNIA, BERKELEY Spring Semester 2014 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials Name:......................................... M.S Comprehensive
More informationDynamic Model of a Badminton Stroke
ISEA 28 CONFERENCE Dynamic Model of a Badminton Stroke M. Kwan* and J. Rasmussen Department of Mechanical Engineering, Aalborg University, 922 Aalborg East, Denmark Phone: +45 994 9317 / Fax: +45 9815
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 information1 Slope Stability for a Cohesive and Frictional Soil
Slope Stability for a Cohesive and Frictional Soil 1-1 1 Slope Stability for a Cohesive and Frictional Soil 1.1 Problem Statement A common problem encountered in engineering soil mechanics is the stability
More informationGeneral Guidelines for Crash Analysis in LS-DYNA. Suri Bala Jim Day. Copyright Livermore Software Technology Corporation
General Guidelines for Crash Analysis in LS-DYNA Suri Bala Jim Day Copyright Livermore Software Technology Corporation Element Shapes Avoid use of triangular shells, tetrahedrons, pentahedrons whenever
More informationMechanics of Materials Primer
Mechanics of Materials rimer Notation: A = area (net = with holes, bearing = in contact, etc...) b = total width of material at a horizontal section d = diameter of a hole D = symbol for diameter E = modulus
More informationExample 37 - Analytical Beam
Example 37 - Analytical Beam Summary This example deals with the use of RADIOSS linear and nonlinear solvers. A beam submitted to a concentrated load on one extremity and fixed on the other hand is studied.
More information