The document describes the implementation of geometrical nonlinearity in finite element analysis software FEASTSMT. It discusses total Lagrangian formulation and the Newton-Raphson method to solve nonlinear finite element equations arising from large deformations. Element formulations accounting for incremental strains, strain-displacement relationships, and stresses are developed. The implementation is validated on a solid element and beam bending problem by comparing results with MARC software and analytical solutions, showing good agreement.
Report
Share
Report
Share
1 of 6
Download to read offline
1
2
3
4
5
6
More Related Content
C012621520
1. IOSR Journal of Mechanical and Civil Engineering (IOSR-JMCE)
e-ISSN: 2278-1684,p-ISSN: 2320-334X, Volume 12, Issue 6 Ver. II (Nov. - Dec. 2015), PP 15-20
www.iosrjournals.org
DOI: 10.9790/1684-12621520 www.iosrjournals.org 15 | Page
Implementation Of Geometrical Nonlinearity in FEASTSMT
Kavya B S1
, Sarangapani G2
Ananya John3
1,3
( Department of Civil Engineering (Indira Gandhi Institute of Engineering & Technology for Women
Mahatma Gandhi University, Kottayam, India) 2
(Scientist / Engineer VSSC, Thiruvananthapuram)
Abstract:Analysis of the structures used in aerospace applications is done using finite element
method. These structures may face unexpected loads because of variable environmental situations.
These loads could lead to large deformation and inelastic manner. The aim of this research is to
formulate the finite elements considering the effect of large deformation and strain. Here total
Lagrangian method is used to consider the effect of large deformation. After deriving required
relations, implementation of formulated equation is done in FEASTSMT
(Finite Element Analysis of
Structures - Substructured and Multi-Threading). .Newton-Raphson method was utilized to solve
nonlinear finite element equations. The validation is carried out with the results obtained from the
Marc Software.
Keywordsâ Geometrical nonlinearity, Newton-Raphson scheme, Total Lagrangian formulation,
I. Introduction
Nonlinearity is natural in physical problems. In fact the linear assumptions made are only valid in
special circumstances and usually involve some measure of "smallness", for example, small strains, small
displacements, small rotations, small changes in temperature, and so on. The finite element equilibrium equation
derived for static analysis is ď ďďť ď˝ ďť ď˝FxK ď˝ .The fig 1 indicates the relation between force and displacement.
In linear analysis [K] is constant, that means if the force doubles, the displacement (and stresses) should be
double. This relation may not be linear in all cases. In these situations, nonlinear analysis may be required.
Fig.1: Force Displacement Diagram.
In nonlinear static analysis, the stiffness [K] is dependent on the displacement {x};
The resulting force-displacement curve may be nonlinear, as shown in fig 2. A nonlinear analysis is an iterative
solution because the relationship between load applied in the structure [F] and corresponding displacement (x) is
not known beforehand. Also no time dependent effects are considered.
Fig.2 : Force Displacement Diagram
Non-linear behaviour admits a wide variety of phenomena, possibly interacting with one another, and
each perhaps difficult to formulate. It is fortunate that the linear model provide satisfactory approximation for
many problems of practical interest such as collapse or buckling of strut due to sudden load overload,
progressive damage behaviour due to long lasting severe loads etc. Non-linear problem pose a difficulty of
describing phenomena by realistic mathematical and numerical models and difficulty of solving non-linear
equation.
2. Implementation Of Geometrical Nonlinearity in FEASTSMT
DOI: 10.9790/1684-12621520 www.iosrjournals.org 16 | Page
In structural mechanics, type of non-linearity include the following
ď Material non-linearity: Material nonlinearities occur when the stress-strain or force-displacement law is
not linear, or when material properties change with the applied loads.
ď Contact non-linearity: In which a gap between adjacent parts may open or closed, the contact area
between the parts changes as the contact force changes or there is a sliding contact with frictional
forces.
ď Geometrical non-linearity: Geometric nonlinearities involve nonlinearities in kinematic quantities such
as the strain-displacement relations in solids. Such nonlinearities can occur due to large displacements,
large strains, large rotations, and so on
The term large deflection[2] is somehow misleading since problems falling in this category need not
have actual deformations which are large. In fact, they can be as small as those in linear cases. However, for the
strain-displacement relationships to be accurate they must include the appropriate higher order nonlinear terms
which are taken to be small and negligible in the linear theory. Nonlinear solid mechanics theory[1] is certainly
more complex than the corresponding linear theory. Consequently, the application of the nonlinear theory to
physical problems can lead to complicated mathematical problems. In addition, analytic solutions are very
limited for the geometrical nonlinear problems in solid mechanics.
II. Total Lagrangian Method And Newton- Raphson Method
2.1 Total Lagrangian Method
There are different methods[4] available for non-linear analysis such as total lagrangian method, updated
lagrangian method depending on whether there exists large or small strain, deflection and rotation. Also the
method depends type non-linearity. In this report geometrical non-linearity is considered. Since the report deals
with implementation of a brick element with small strain and large deformation, total lagrangian method[1] is
used for the analysis and Newton-Raphson method for iteration.
1.Equation of motion
d0
V = ; where, = ;
= ( + + )
2.Incremental decompositions
(a) Stresses
= +
(b)Strains
= ; = +
= ( + + + ); ( + ;
3.Equation of motion with incremental decompositions
Noting that = the equation of motion is
d0
V + d0
V = - d0
V
= , =
4.Linearization of equation of motion
Using the approximations = , = , we obtain as appxoximate equation of motion
d0
V+ = - d0
V
2.2 Newton Raphson Method
Nonlinear solutions require several iterations. The actual relationship between load and displacement is
not known beforehand. Consequently , a series of linear approximations with correction is performed. The fig 3
is a simplified explanation of the Newton-Raphson method[2].
3. Implementation Of Geometrical Nonlinearity in FEASTSMT
DOI: 10.9790/1684-12621520 www.iosrjournals.org 17 | Page
Fig.3: force displacement diagram
ďˇ In the Newton- Raphson Method, the total load Fa is applied in iteration 1. The total load Fa is applied
in iteration 1.The result is x1. From the displacements, the internal forces F1 can be calculated. If Fa â
F1, then the system is not in equilibrium. Hence a new stiffness matrix (slope of red line) is calculated
based on the current conditions. The difference of Fa-F1 is the out-of-balance or residual forces. The
residual forces must be âsmallâ enough for the solution to converge.
ďˇ This process is repeated until Fa â Fi. In this example, after iteration 4, the system achieves equilibrium
and the solution is said to be converged.
III. Formulation Of Solid Elment
1.Incremental strains
( + )+ + + + ); i = 1,2,3; j = 1,2,3; k = 1,2,3; = (1)
2. Linear strain-displacement transformation matrix
Using = L ; where = ; (2)
= ; L = + (3)
= (4)
= ; = - ; = (5)
3.Nonlinear strain- displacement transformation matrix
= ; ; where = (6)
4. Second Piola-Kirchhoff stress matrix and vector
= ; = (7)
= ; where = (8)
IV. Implementation
ďˇ
Addition of incremental strains
ďˇ
Addition of Linear strain-displacement transformation matrix
ďˇ
Addition of nonlinear strain- displacement transformation matrix
4. Implementation Of Geometrical Nonlinearity in FEASTSMT
DOI: 10.9790/1684-12621520 www.iosrjournals.org 18 | Page
ďˇ Addition of Second Piola-Kirchhoff stress matrix and vector.
V. Validation Problems
5.1 Solid Element
Material properties:
Youngâs modulus, E = 1e8 N/mm2
,Poissonâs ratio = 0.3
Boundary and loading conditions:
ď Fully fixed condition (Ux = Uy = Uz = 0) on one face of the solid and (Uy = Uz = 0) on
opposite face as shown and a point Load (Fx) of 1.5e8N on four nodes as shown
Fig 4: Modelling in FEASTSMT
Fig 5: Modelling in MARC
Fig 6: Displacement contours of non-linear analysis from FEASTSMT
and MARC
LOAD DEFLECTION in mm
FEASTSMT
MARC
1 0.194302 0.193783
2 0.349164 0.349032
3 0.491028 0.481077
4 0.603111 0.597259
5 0.706461 0.701749
6 0.800946 0.797186
7 0.888453 0.885357
8 0.97014 0.967539
9 1.0469 1.04468
10 1.11943 1.11751
5. Implementation Of Geometrical Nonlinearity in FEASTSMT
DOI: 10.9790/1684-12621520 www.iosrjournals.org 19 | Page
Fig 7: Force-Displacement Curve
5.2 Simply Supported Beam
Material properties:
Youngâs modulus, E = 1e8 N/mm2
; Poissonâs ratio = 0.3
Boundary and loading conditions:
⢠Simply supported condition (Ux = Uy = Uz =0) on both ends of the beam
⢠A udl of 1000N/m on top face as shown
Fig 8: Modelling in FEASTSMT
Fig 9: Modelling in MARC
Fig 10: Displacement contours of non-linear analysis from FEASTSMT
and MARC
6. Implementation Of Geometrical Nonlinearity in FEASTSMT
DOI: 10.9790/1684-12621520 www.iosrjournals.org 20 | Page
Fig 11: Force-Displacement Curve
VI. Conclusion
The paper consists of implementation of geometrical non-linearity in FEASTSMT
. The formulation of
total lagrangian method and Newton-Raphson Method are implemented and the method is validated by using a
solid finite element for simply supported beam with udl and the analysis of a pressure vessel. The results were
compared with the analytical solution and the results obtained using the Marc software. The results are matching
well between FEASTSMT
and the Marc software as well as with the analytical solution.
References
[1]. Klaus-Jurgen Bathe and Said Bolourchit âLarge Displacement Analysis of Three-Dimensional Beam Structuresâ International
Journal For Numerical Methods In Engineering, Vol. 14,pp 961-986, 1979
[2]. A. Banerjee, B. B. Bhattacharya, A.K. Mallik,âLarge deflection of cantilever beams with geometric non-linearity:Analytical and
numerical approachesâ International Journal Of Non-Linear Mechanics June 2008
[3]. Dr. Seref Doguscan Akbas âThermal Deflections of a Simple Supported Beam with Temperature Dependet Physical Propertiesâ
International Journal of Solids and Structures ,2009
[4]. Izadpanah M HabibiA.R Zangeneh FarM, âNonlinear analysis of plane trusses using lagrangian descriptionsâ International Journal
Of Current Life Sciences - Vol.4, Issue, 5, pp. 644-648, May, 2015.
LOAD
DEFLECTION in mm
FEASTSMT
MARC
1 0.376326 0.3785
2 0.90927 0.9099
3 1.81396 1.752
4 2.89533 2.9443
5 3.52412 3.604
6 3.93534 3.964
7 4.24553 4.269
8 4.4977 4.519
9 4.71198 4.732
10 4.89943 4.919