Chapter-1 1.1 Finite Element Analysis: Discretisation, and The Assembly of Elements Is Called A Mesh
Chapter-1 1.1 Finite Element Analysis: Discretisation, and The Assembly of Elements Is Called A Mesh
Chapter-1 1.1 Finite Element Analysis: Discretisation, and The Assembly of Elements Is Called A Mesh
INTRODUCTION
Elements can be of various shapes (as shown in Figure 1.2), in two dimensions,
quadrilateral or triangular, and in three-dimensions, brick-shaped (hexahedral), wedge-shaped
(pentahedral) or tetrahedral. This is, of course, not an exhaustive list.
To linear elastic analysis for the moment, then the quantity, that is, as a rule, first found
in the analysis is the displacement at series of points called nodes. The nodes are at the corners
of the elements and, depending on the element type, possibly at the midsides of the elements
or even within the element. Nodes on the boundaries of adjacent elements must belong to the
elements that meet there: examples of permitted and forbidden meshes are shown in Figure
1.3(A) and Figure 1.3(B) respectively. The analysis calculates the displacement at the nodes
for the particular loading applied to the FE model.
The displacement of each point within an element is fixed by the values of the
displacements of the nodes of the element, that is, it is a function of the nodal displacements.
In this way, the problem of finding the displacement of every point within the body is replaced
by the problem of finding the displacements of a finite number of points, namely, the nodes.
The displacement of each point is then defined in terms of the displacements of the nodes of
the element to which the point belongs. If we are considering a two dimensional model, then
the displacement of each node consists of two components, one parallel to a reference x axis
and a second parallel to the y axis: these are called degrees of freedom.
Pre-processing
Pre-processing is concerned with the creation of the model and the definition of the
way in which it is to be loaded. The pre-processor includes a graphics package that enables the
user to build up the model of the component to be analyzed and to display the model on the
computer screen. The successfulness of the entire analysis is largely determined at this stage
by the skill of the user in determining what simplifications (if any) are to be introduced into
the model as compared with the `real thing' and by the choice of the mesh and type of element
to be used. Appropriate mechanical or physical properties must be allocated to the material of
which the model is made and loads and possibly any restriction to the movement of certain
nodes (restraints) must be applied.
Analysis
The analysis part of the FE package takes in the input file carries out certain checks on
the information contained therein and then, if there are no errors in the input file, the analysis
is carried out and output files are produced. These contain an enormous amount of information
if the analysis is at all complex. These files can be examined and the relevant information
extracted but, as a rule, there is so much information that it needs to be presented to the user in
a more intelligible and user-friendly manner. This is the job of the post-processor. The pre and
post-processor are essentially the same software package.
Post-processing
The post-processor takes in the information from the output files and can present it to
the user in a range of different graphical and tabular forms. For example, depending on the
facilities available, colour may be used to indicate the value of some component of stress on
the surface of the component, or contour lines of equal stress may be drawn as in Figure 1.4,
or similar forms of display may be produced on sections through the model.
Fundamental equation is
σ = Eε
In this equation, the proportionality constant E between strain and stress is the “Modulus
of Elasticity” of the material.
In order to cover three-dimensional stress and strain states, in a first step we solve this
equation for ε and just look for the first principal strain
1
ε1 = .σ 1
E
Now, we add on the right side of this equation the missing terms from the two lateral
principal stresses σ 2 and σ 3 . Compared to σ 1 , these lateral stresses influence the first principal
strain ε1 much less: So, they are multiplied with “proportionality constant”, known as the
Poisson’s ratioν :
1
ε1 = {σ 1 −ν (σ 2 + σ 3 )}
E
A hyperelastic material is still an elastic material, which means it returns to its original
shape after the forces have been removed. Hyperelastic material also is Cauchy-elastic, which
means that the stress is determined by the current state of deformation, and not the path or
history of deformation.
The difference to linear elastic Material is that in hyperelastic material the stress-strain
relationship derives from a strain energy density function, and not a constant factor. This
definition says nothing about the Poisson's ratio or the amount of deformation that a material
will undergo under loading. However, often elastomers are modeled as hyperelastic.
Hyperelasticity may also be used to describe biological materials, like tissue
Elastomers are often modeled as hyperelastic
Elastomers (like rubber) typically have large strains (often some 100 %) at small loads
(means a very low modulus of elasticity, for example just 10MPa). The material is nearly
incompressible, so the Poisson’s ratio is very close to 0.5
Their loading and unloading stress-strain curve is not the same, depending on different
influence factors (time, static or dynamic loading, frequency, etc.). This viscous behavior is
ignored if the hyperelastic material model is used for description.
Elastomer material in comparison with metals and plastics
Energy-elasticity: Loading changes the distance of the atoms within the lattice of the metal
and so increases the internal energy. When unloading it, this energy is immediately set free,
the initial shape appears again.
Entropy-elasticity: Within an elastomer, its macromolecules are balled if unloaded. During
loading, a stretching and unballing appears. After unloading, more or less the unordered state
appears again. Viscous behavior: every loading leads to an even small remaining deformation
(creeping, relaxation).
1.3.1 MOONEY-RIVLIN
The Mooney-Rivlin model applies to current-technology shell, beam, solid, and plane
elements.
For NPTS = 2 (2 parameter Mooney-Rivlin option, which is also the default), the form of
the strain energy potential is:
1
= C10 ( I1 − 3) + C01 ( I 2 − 3) +
W
d ( I − 1) 2
Where,
C10 , C01 = material constants characterizing the deviatoric deformation of the material
LITERATURE SURVEY
Sedat Suslera & Halit Turkmen. S, In this study, the geometrically non-linear dynamic
behavior of simply supported tapered laminated composite plates subjected to the air blast
loading is investigated numerically. In-plane stiffness, inertia and the geometric non-linearity
effects are considered in the formulation of the problem. The equations of motion for the
tapered laminated plate are derived by the use of the virtual work principle. Approximate
solution functions are assumed for the space domain and substituted into the equations of
motion. Then, the Galerkin method is used to obtain the non-linear algebraic differential
equations in the time domain. The resulting equations are solved by using the finite difference
approximation over the time. The effects of the taper ratio, the stacking sequence and the fiber
orientation angle on the dynamic response are investigated. The displacement-time and strain-
time histories are obtained on certain points in the tapered direction. The results obtained by
using the present method are compared with the ones obtained by using a commercial finite
element software ANSYS. The results are found to be in an agreement. The method presented
here is able to determine the non-linear dynamic response of simply supported tapered
laminated plates to the air blast loading accurately.
Chang Chuan Xie & Jia Zhen Leng, A composite high-aspect-ratio wing of a high-altitude
long-endurance (HALE) aircraft was modeled with FEM by MSC/NASTRAN, and the non-
linear static equilibrium state is calculated under design load with follower force effect, but
without load redistribution. Assuming the little vibration amplitude of the wing around the
static equilibrium state, the system is linearized and the natural frequencies and mode shapes
of the deformed structure are obtained. Planar doublet lattice method is used to calculate
unsteady aerodynamics in frequency domain ignoring the bending effect of the deflected wing.
And then, the aeroelastic stability analysis of the system under a given load condition is
successively carried out. Comparing with the linear results, the non-linear displacement of the
wing tip is higher. The results indicate that the critical non-linear flutter is of the
flap/chordwise bending type because of the chord wise bending having quite a large torsion
component, with low critical speed and slowly growing damping, which does not appear in the
linear analysis. Furthermore, it is shown that the variation of the non-linear flutter speed
depends on the scale of the load and on the chord wise bending frequency. The research work
indicates that, for the very flexible HALE aircraft, the non-linear aeroelastic stability is very
important, and should be considered in the design progress. Using present FEM software as
the structure solver (e.g. MSC/NASTRAN), and the unsteady aerodynamic code, the non-
linear aeroelastic stability margin of a complex system other than a simple beam model can be
determined.
Majid Shahzad & Ali Kamran, The aim of research work is to characterize hyperelastic
material and to determine a suitable strain energy function (SEF) for an indigenously
developed rubber to be used in flexible joint use for thrust vectoring of solid rocket motor. In
order to evaluate appropriate SEF uniaxial and volumetric tests along with equi-biaxial and
planar shear tests were conducted. Digital image correlation (DIC) technique was utilized to
have strain measurements for biaxial and planar specimens to input stress-strain data in
Abaqus. Yeoh model seems to be right choice, among the available material models, because
of its ability to match experimental stress-strain data at small and large strain values. Quad lap
specimen test was performed to validate material model fitted from test data. FE simulations
were carried out to verify the behavior as predicted by Yeoh model and results are found to be
in good agreement with the experimental data.
Jadhav. A.N & Dr. S.R. Bahulikar, A Hyperelastic material is type of the ideally elastic
material for which elasticity shows non-linear behavior, because of that the stress strain
relationship for them expressed in terms of strain energy density. Hyperelastic models are used
to model the mechanical behavior of rubber-like materials ranging from elastomers, such as
natural rubber and silicon, to biologic materials, such as muscles and skin tissue. The
presented work is carried out to study the effect of the different Mooney-Rivlin hyperelastic
models used in commercial finite element software. All the models studied have under gone
same loading and boundary conditions and finally compared with testing. The final conclusion
based on their capturing hyper elasticity of material is stated
Nor Fazli Adull Manan & Siti Noor Azizzati Mohd Noor, Soft tissues, skin and rubber-like
material research has gained considerable attention over the last few years. Most related
research has adopted hyperelastic material models namely Mooney-Rivlin, Neo-Hookean and
Ogden models. The complex behavior (highly non-linear) makes it challenging to be analyzed.
Nevertheless, the fundamental understanding of a particular hyperelastic model could assist
researchers to have a better judgement on their findings. Therefore, performing a parametric
study is vital especially at the beginning of such numerical analysis. This study aims to
investigate numerically, the contribution of material parameters of two hyperelastic
constitutive models viz. Ogden and Mooney-Rivlin. This study is divided into three stages, i.e.
the derivation and simplification of the hyperelastic equations, the parametric analysis of the
hyperelastic models and the demonstration of the hyperelastic material parameters with
respect to stress and stretch. The second stage consists of the following investigations (1. The
influence of the first material parameter Constant 1, C1 and Ogden Exponent, α for Mooney-
Rivlin and Ogden Model respectively; 2. The influence of the second material parameter
Constant 2, C2 and Ogden Coefficient, µ for Mooney-Rivlin and Ogden Model
correspondingly and 3. The influence and sensitivity of the stress levels with different stretch,
on the three-dimensional stress-stretch curve). From this parametric study, the Mooney-Rivlin
Model indicates that an incompressible, isotropic Mooney-Rivlin Model is more sensitive to
C1 variations as compared to C2 variations. Nevertheless, the Ogden model shows results that
the higher Ogden exponent influence the stress level as well as the stress-stretch curve with an
accelerated stress increment at the beginning of deformation. It is also evident from the stress-
stretch curves illustrated in each investigation; the hyperelastic models are sensitive towards
different material coefficients.
Yihong Hong & Wenjuan Yao, Al-polymer laminated membranes are widely used in large
aerospace structures. When the laminated membranes are pressurized, wrinkles emerge, which
have an important effect on the performance of the structures during operation. This paper
describes the numerical simulation and experimental investigation of wrinkles in laminated
membranes. The non-linear post buckling analysis method, based on laminated thin-shell
elements, was used to simulate the onset, growth, and final configuration of wrinkles when
laminated membranes are subjected to external loads. The simulations are conducted with the
ANSYS finite element package. Changing regularities of number, wave length, and range for
the wrinkles during the onset and growth processes are investigated. The wrinkles of
laminated membranes with different design parameters such as material selection, ply number,
ply angle, and ply mode are predicted. Devices that can be used to clamp and load laminated
membranes in several load cases were designed and developed. A 3D photogrammetric system
was constructed to characterize wrinkling patterns of laminated membranes subjected to shear
displacement loads. By comparing the results of numerical analysis and experimental results,
the accuracy of the numerical analysis method was verified. This study work is expected to
inform wrinkling simulation and shape control of aerospace laminated membrane structures.
Marcel STERE & Daniela BARAN, The state – of - the art of aeronautical structures show
that parts are manufactured and subsequently assembled with the use of fasteners and/ or
bonding. Adhesive bonding is a key technology to low weight, high fatigue resistance,
robustness and an attractive design for cost structures. The paper results resolve significant
problems for two groups of end-users: 1) for the aerospace design office: a robust procedure
for the design of the hybrid joint structural components; 2) for the aeronautical repair centers:
a useful procedure for structural design and analysis with significant cost savings.
Sani. M.S.M & Abdullah. N.A.Z, Model updating is a process of making adjustment of
certain parameters of finite element model in order to reduce discrepancy between analytical
predictions of finite element (FE) and experimental results. Finite element model updating is
considered as an important field of study as practical application of finite element method
often shows discrepancy to the test result. The aim of this research is to perform model
updating procedure on a composite structure as well as trying improving the presumed
geometrical and material properties of tested composite structure in finite element prediction.
The composite structure concerned in this study is a plate of reinforced kenaf fiber with epoxy.
Modal properties (natural frequency, mode shapes, and damping ratio) of the kenaf fiber
structure will be determined using both experimental modal analysis (EMA) and finite element
analysis (FEA). In EMA, modal testing will be carried out using impact hammer test while
normal mode analysis using FEA will be carried out using MSC. Nastran/Patran software.
Correlation of the data will be carried out before optimizing the data from FEA. Several
parameters will be considered and selected for the model updating procedure.
CHAPTER-3
METHODOLOGY
The problems which have been selected for the analysis is been from ANSYS tutorial
problems. The ANSYS problem is chosen from Non-linear Hyperelastic Mooney-Rivlin. The
target material is made of polyethylene and the hyperelastic material (contact body) is made
up of Mooney-Rivlin.
1. Mooney-Rivlin 3 parameters
2. Mooney-Rivlin 5 parameters
3. Mooney-Rivlin 9 parameters
From this we have chosen Mooney-Rivlin 5 parameters for the analysis. The constant 5
parameters are solved from curve fitting. These all analysis data’s are selected from the
engineering data sources. This analysis is done by two-dimensional analysis.
The polyethylene edges are selected as target bodies and Elastomer material (Mooney-
Rivlin) are selected as contact body. These two bodies (contact and target) are initiated with
frictional contact. The constant frictional co-efficient is 0.14. The normal stiffness factor is
converted into manual. The normal stiffness factor in this analysis is 0.1. The stiffness factor
should be updated in each and every iteration to improve the accuracy of analysis. The mesh is
done by Refinement meshing for the increased mesh at the contact point. Due to the solution
steps will be highly increased. Now the Boundary conditions is been allotted. Polyethylene
material is been given fixed support and Elastomer material is been given displacement.
3.2.1 ANALYSIS
In this analysis total deformation, maximum principal stress and von misses strain are been
analyzed and the behavior of the Hyperelastic Mooney-Rivlin material is been identified. This
model is been used in the fishing rod and reels.
LITERATURE SURVEY
FINALISE THE
DESIGN
PARAMETERS
ANALYZING
DESIGN
MESH
Nodes: 2935
Elements: 904
FIXED SUPPORT
DISPLACEMENT
FIXED SUPPORT
Polyethylene
Contact body
Polyethylene
Part body 2
Part body 3
ANSYS16.0→ WORKBENCH
MATERIAL PROPERTIES
II. GEOMETRY
Tools→ Mid surface→ Click both surface of the material one by one →Generate the
DETAILS OF IMPORT
III. MODEL
Click manual contact region→ Click edge (ctrl+E) contact→ Select the contact edges
only
Click target→ Select edge→ Above the Mooney-Rivlin material→ Apply
Connections→ Contacts→ Right click→ Select duplicate→ Click target→ Select
the polyethylene material edges
Select→ Both Bonded part body to part body type→ Change bonded to frictional
• GEOMETRY
DEFINITION
STATISTICS
Bodies 3
Active bodies 3
Nodes 2935
Elements 904
IV. MESH
SCOPE
Table 4.5 Target and contact edges
Contact 5 Edges
Target 6 Edges
DEFINITION
Table 4.6 Conditions for Bonded surfaces
Type Frictional
Frictional coefficient 0.14
Behavior Symmetric
V. SETUP
Static structural→ insert→ Fixed support→ Top and Bottom edge of the polyethylene
material →Apply
Table 4.7 (A) & (B) Conditions for Analysis
ADVANCED
(A)
Formulation Pure penalty
Detection method Nodal-Normal to Target
Normal stiffness Manual
Normal stiffness factor 0.1
Update stiffness Each iteration
Pinball region Radius
Pinball radius 2mm
GEOMETRIC MODIFICATION
(B)
Interface Treatment Add Offset, No Ramping
VI. SOLUTION
• ANALYSIS SETTINGS
STEP CONTROLS
(A)
Auto time stepping ON
Initial sub steps 100
Minimum sub steps 50
Maximum sub steps 1000
SOLVER CONTROLS
(B)
Solver type Direct
Weak springs OFF
Large deflection ON
Inertial relief OFF
OUTPUT CONTROLS
(C)
Stress YES
Strain YES
Nodal force YES
RESULT
CONCLUSION
The literature survey has been carried out from the various Finite Element Analysis of
Non- linearity and Hyperelastic of structures. The non-linear analysis structure is designed
using CATIA. This CATIA designed structure is been imported into ANSYS. The structure
which is imported in ANSYS is been analyzed and interpenetration is been taken place. This
problem identifies the Interpenetration of targets and contact surfaces.
The problem is been solved at the present mesh but the result analysis is not accurate
due to the interpenetration at contact region. Generally interpenetration should be not taken
place because it is not physically possible. The interpenetration is occurred due to the pin ball
parameter and element mesh. In next simulation the meshing count will be increased and
interpenetration will be prevented.
REFERENCES
1. Chang Chuan Xie & Jia Zhen Leng, “Geometrical nonlinear aeroelastic stability
analysis of a composite high-aspect-ratio wing” Shock and Vibration 15 (2008) 325–
333 IOS Press
2. Houshyar. S & Shanks. R.A, “Modeling of polypropylene fiber-matrix composites
using finite element analysis” DOI: 10.3144/expresspolymlett.2009.2
3. Huei-Huang Lee, “Finite Element Simulations With ANSYS Workbench” ISBN:978-
1-58503-653-0
4. Jadhav. A.N & Dr. S.R.Bahulikar, “Comparative Study of Variation of Mooney Rivlin
Hyperelastic Material Models under Uniaxial Tensile Loading” Vol-2 Issue-4 2016
IJARIIE-ISSN(O)-2395-4396
5. Majid Shahzad & Ali Kamran, “Mechanical Characterization and FE Modeling of a
Hyperelastic Material” Materials Research. 2015; 18(5): 918-924 DOI:
http://dx.doi.org/10.1590/1516-1439.320414
6. Marcel STERE & Daniela BARAN, “Calculation of hybrid joints used in modern
aerospace structures” DOI: 10.13111/2066-8201.2011.3.4.1525
7. Nor Fazli Adull Manan & Siti Noor Azizzati Mohd Noor, “Numerical Investigation Of
Ogden and Mooney-Rivlin Material Parameters” VOL. 10, NO. 15, AUGUST 2015
ISSN 1819-6608
8. Sani. M.S.M & Abdullah. N.A.Z, “Finite element model updating of natural fiber
reinforced composite structure in structural dynamics” MATEC Web of Conferences
477, 83 CSNDD 2 016, DOI: 10.1051/matecconf/20168303007
9. Sani. M.S.M & Ouyang. H, “Identification of Nonlinearities in Joints of a Wing
Structure” MATEC Web of Conferences 477 , 83 CSNDD 2016
10. Sedat Susler & Halit Turkmen. S, “The nonlinear dynamic behavior of tapered
laminated plates subjected to blast loading” Shock and Vibration 19 (2012) 1235–1255
DOI 10.3233/SAV-2012-0667 IOS Press
11. Yihong Hong & Wenjuan Yao, “Numerical and Experimental Investigation of
Wrinkling Pattern for Aerospace Laminated Membrane Structures” Volume 2017,
Article ID 8476041, 11 pages https://doi.org/10.1155/2017/8476041