Djs Fea Lab Manual
Djs Fea Lab Manual
FEA
LAB
MANUAL
. Page 1
Finite Element Analysis Lab
List of Experiments:
Students should use the commercial software or programmes form the text-books or self-developed programs, to
verify the results obtained by manual calculations. The input data and output results of the problem solved using
the computer programs should be included in the Journal. The proposed list is given below:
1. Any two problems using 1D element: Analysis of axially loaded stepped/Tapered bar and beam.
2. Any two problems using truss element: Analysis of plane trusses.
3. Any two problems using 2D element: Analysis of plate with hole/notch.
4. Any one problem on steady state heat conduction.
5. Any one problem using axisymmetric element: Analysis of pressure vessel.
6. Any one problem of free vibration analysis for dynamic system.
While performing the analysis the students should understand the concepts of selection of element type, meshing and
convergence of solution.
Course Outcomes:
1. Apply basic aspects of FEA to solve engineering problems.
2. Select appropriate element, suitable meshing, appropriate solver and perform Convergence test for given
problem.
3. Interpret and validate FEA result.
. Page 2
1.1 Introduction
The finite element method has been a powerful tool for the numerical solution of a wide
range of engineering problems. Application range from deformation and stress analysis of
automotive, aircraft, building, defense, missile, and bridge structures to the field analysis of
dynamics, stability, fracture mechanics, heat flux, fluid flow, magnetic flux, seepage and other
flow problems.
With the advances in computer technology and CAD systems Complex problems can be
modeled with relative ease. Several alternate configurations can be tried out on a computer before
the first prototype is built.
The basics in engineering field is must to idealize the given structure for the required
behavior. The proven knowledge in the typical problem area, modeling techniques, data transfer
and integration, computational aspects of the finite element method is essential.
Most often it is not possible to ascertain the behavior of complex continuous system
without some sort of approximation. For simple members like uniform beams, plates etc. Classical
solutions can be sought by forming differential and/or integral equations through structures like
machine tool frames, Pressure vessels, automobile bodies, ships, aircraft structures, domes etc,
need some approximate treatment to arrive at their behavior, in static deformation, dynamic
properties or heat conducting property.
Indeed these are continuous systems with their mass and elasticity being continuously
distributed.
The classical differential equation solution approach leads to intractability. To overcome
this engineers and mathematicians have form time to time proposed complex structure is defined
using a finite number of well-defined components. Such systems are then regarded as discrete
systems. The discretization method could be finite difference approximation, various residual
procedures etc.
Finite element method comes under this category of discretization methods. R.W. Clough
appears to be the first to use this term of finite element since 1960s there has been much progress
in the method. The method requires a large number of computations requiring a fast computer. In
fact digital computer advances have been responsible for the expanding usage of the finite element
method.
The finite element method was initially developed to solve structural problems. Its use of
late, has been rapidly extended to various fields. which still by no means can be claimed as
complete since fields of usage are being continuously diversified.
Page 3
1.2 Methods of Engineering Analysis
The general methods for solving the practical problems can be classified as shown in
figure:
Methods of analysis
Exact method Provides exact solution to the problems, but the limitation of this method is
that all practical problems cannot be solved and even if they can be solved, they may have
very complex solution.
Approximate analytical methods are alternative to the exact methods, in which certain
functions are assumed to satisfy the geometric boundary conditions, but not necessary the
governing equilibrium equations. These assumed functions, which are simpler, are then
solved by any conventional method available. The solutions obtained from these methods
Page 6
have limited range of values of variables for which the approximate solution is nearer
to the exact solution.
Numerical method has been developed to solve almost all types of practical problems.
The two common numerical methods used to solve the governing equations of practical
problems are numerical integration and finite element technique.
Numerical integration methods such as Runge-Kutta, Miline’s method etc. adopt
averaging of slopes of given function at the given initial values. These methods yield
better solutions over the entire field width but sometimes laborious.
Finite difference method the differential equations are approximated by finite
difference equation. Thus the given governing equation is converted to a set of
algebraic equations. These simultaneous equations can be solved by any simple
method such as Gauss Elimination, gauss-Seidal iteration method, Crout’s method
etc. the method of finite difference yield fairly good results and are relatively easy
to program. Hence, they are popular in solving heat transfer and fluid flow
problems. However, it is not suitable for problems with irregular geometry and not
suitable for problems of rapidly changing variables such as stress concentration
problems.
Finite Element Method (F.E.M) has emerged out to be powerful method for all kinds
of practical problems. In this method the solution region is considered to be built up of many
small-interconnected sub regions, called finite elements. These elements are applied with
an interpolation model, which is simplified version of substitute to the governing equation
of the material continuum property. The stiffness matrices obtained for these elements are
assembled together and the boundary conditions of the actual problems are satisfied to
obtain the solution all over the body or region FEM is well suited for computer
programming.
Boundary element method (B.E.M) like finite method is being used in all engineering
fields. In this approach, the governing differential equations are transformed into integral
identities applicable over the surface or boundary. These integral identities are integrated
over the boundary, which is divided into small boundary segments. As in finite element
method provided that the boundary conditions are satisfied, set of linear algebraic equations
emerges, for which a unique solution is obtained.
Page 8
element analysis. The general purpose mesh generation scheme can model a variety of complex
regions. The node numbering gives spars matrices and in many cases should given minimum
bandwidth proper block representation. Mesh plotting shows the element layout.
Finite Element Analysis is a way to simulate loading conditions on a design and determine
the design’s response to those conditions.
The design is modeled using discrete building blocks called elements
Each element has exact equations that describe how it responds to a
certain load.
The “sum” of the response of all elements in the model gives the total
response of the design.
The elements have a finite number of unknowns, hence the name finite
elements.
The finite element model, which has a finite number of unknowns, can only approximate
the response of the physical system, which has infinite unknowns so the question arises:
How good is the approximation?
Unfortunately, there is no easy answer to this question. It depends entirely on what you are
simulating and the tools you use for the simulation. We will, however, attempt to give you
guidelines throughout this training course.
Page 9
1.6 Degrees of Freedom
The point in space used for defining the geometry and connecting the elements. These points
Page
10
are useful to specify the displacement, force and boundary conditions of the structure
3. Displacement
Displacement defines the deformation or the movement caused by an external force. There
are six displacements for an element in space. These are three translations vectors u; v and w and
are three rotational vectors with respect to x, y and z. co-ordinate system. These nodal variables
are referred to as degrees of freedom.
4. Co-ordinate System
A co-ordinate system is the reference, which can specify the position of a point or an
element in a space. The co-ordinate systems used are.
Global co-ordinate system (Cartesian i.e. x. y. z.)
Local co-ordinate system.
Polar co-ordinate system(r. 0)
Cylindrical co-ordinate system(r. 0. z)
Spherical co-ordinate system(r. 0)
Can be Cartesian also (x. y. z)
5. Boundary Conditions/Constraints
Resistance (reaction) offered at grid points due to the motion of the structure being
constrained by the applied forces.
These reactions can be obtained with respect to the specified co-ordinate
System.
To prevent the rigid body motion there by reducing the continuous structure to discrete
system(avoiding the phenomena of mechanism)
These constraints can be zero (u. v. w. x. y. z = 0) or non-Zero and specified enforced
displacement at grid points.
The necessary conditions are imposed and the equations of equilibrium are solved for the
nodal displacements. There are two approaches of applying the boundary conditions.
a. Elimination technique
In this technique after imposing the boundary conditions, whichever displacement
is Zero, the corresponding row and column of the overall equilibrium is eliminated.
b. Penalty approach
A second approach for handling boundary condition is the penalty approach. This
approach easily implemented in a computer program and retains its simplicity even when
considering general boundary conditions where a large value of C is added to the highest
diagonal elements. On applying the boundary conditions, the overall equilibrium equation
gets modified.
c. Solution of the modified overall equilibrium equation by gauss elimination method
The value of nodal displacements will be obtained in this step. Thus knowing of all
the nodal displacements, stain of each element will be computed. From these values
elemental stresses are obtained.
Computation of Stains:
Page
12
Where [] = Strain Matrix
{b} = Strain Displacement Matrix where b is unique for
each element
{q} = Elemental Nodal Displacement Vector.
Computation of Stress:
Knowing the strain, the stress is calculated as follows
i.e., [] = {D} {}
Where [] = Stress Matrix
{D} = Constitutive Matrix for plane stress
{} = Strain Matrix
These are the main steps carried out in structural problems where the analysis is simple. It
becomes complicated thereafter for each and every problem, which includes dynamic analysis
where the forces and all the conditions vary with time. In thermal, it is transient analysis and in
fluid flow situations, it is turbulent flow analysis.
The finite element method was initially developed on a physical basis for the analysis of problems
in structural mechanics; however, it was soon recognized that the method can be applied equally
well to the solution of many other classes or problems
FEM can be easily interpreted in physical terms and it has a strong mathematical base.
Hence, finite element method can be easily applied to any problem with a proper
knowledge of the physical system under consideration and can be solved to a greater
accuracy by the application of proper mathematical tool.
Non-homogeneous continuum can also be dealt with by merely assigning different
properties to different elements. It is even possible to vary the properties within an element
according to the polynomial applied.
Finite element method accommodates complex geometry with ease and is capable of
handling non-linear and time dependent system also.
In finite element method, since boundary conditions are introduced in the assembled
equations, it requires only to specify the geometric boundary conditions without regarding
its effects on interior elements. Since the boundary conditions do not enter into the
individual finite element equations, the field variable models need not be changed, when
the boundary conditions change.
Finite element method considers the multidimensional continuity of the body. Hence it does not
Page
13
require separate interpolation process to extend the approximate solution to every point within the
continuum. Also it does not require the trial solutions that must all apply to the entire
multidimensional continuum.
FEM has been developed as a very high solution technique. However, the solution obtained
from FEM can be realistic if and only if the material properties are known precisely.
The major drawback of FEM is sensitivity of the solution on the geometry of the element
such as type, shape and orientation of elements used.
The compute programs yield a large amount of data as results it is very difficult to separate
out the required results from the pile of numbers.
1. Geometry is simplified.
Mistakes by Users
1. Elements are of the wrong type: e.g.) Shell elements are used where solid elements are
needed
Page
14
2. Distorted elements
Element Properties: The discretization of the structure or body into Finite Elements
Forms the basic first step in the analysis of a complicated structural system.
Mesh Generators:
The finite element models for practical analysis can contain tens of thousands or even
hundreds of thousands degrees of freedom. It is not possible to create such meshes manually. Mesh
generator is a software tool, which divides the solution domain into many Sub domains – finite
elements. Mesh generators can be of different types.
For two-dimensional problems, we want to mention two types: block mesh generators and
triangulates.
Block Mesh Generators:
Require some initial form of gross partitioning. The solution domain is partitioned in some
relatively small number of blocks. Each block should have some standard form The mesh inside
block is usually generated by mapping technique.
Triangulators:
Typically generate irregular mesh inside arbitrary domains. Voronoi polygons and
Delaunay triangulation are widely used to generate mesh. Later triangular mesh can be transformed
to the mesh consisting of quadrilateral elements. Delaney triangulation can be generalized for three
dimensional domains.
Page
15
Finite Element Model of Gear tooth
Page
16
Three-dimensional stress on an element
Page
17
2.1 Introduction to ANSYS
ANSYS/ ANSYS/
Structural ProFEA
ANSYS/
ED
ANSYS/
FLOTRAN
DesignSpace
ANSYS/
ANSYS/ ANSYS/
LS-DYNA University
PrepPost
Page 18
Other product lines:
o ANSYS/LS-DYNA - for highly nonlinear structural problems
o DesignSpace - an easy-to-use design and analysis tool meant for quick
analysis within the CAD environment
o ANSYS/ProFEA - for ANSYS analysis & design optimization within
Pro/ENGINEER
Page 20
Icon Toolbar Menu Utility Menu
Input Line Raise/Hidden Icon
Graphics Area
Main Menu
Output
Window
Utility Menu:-Contains functions which are available throughout the ANSYS session, such as
file controls, selecting, graphics controls, parameters, and exiting.
Toolbar Menu
Contains push buttons for executing commonly used ANSYS commands and
functions. Customized buttons can be created.
Graphics Area
Displays graphics created in ANSYS or imported into ANSYS.
Input Line
Displays program prompt messages and a text field for typing commands. All
previously typed commands appear for easy reference and access.
Main Menu
Contains the primary ANSYS functions, organized by processors (preprocessor,
solution, general postprocessor, etc.)
Output
Displays text output from the program. It is usually positioned behind the other
windows and can be raised to the front when necessary.
Resume:
This is opening a previously saved database. It is important to know that if you
Simply resume a database; it doesn’t change the job name. For example: You start ANSYS
Page 21
with a job name of “file”. Then you resume my model.db, do some work, and then save.
That save is done to file.db! Avoid this issue by always resuming using the icon on the
toolbar. If you open mymodel.db using this method, it resumes the model and automatically
changes the job name to my model.
Plotting:
Contrary to the name, this has nothing to do with sending an image to a plotter or
printer. Plotting in ANSYS refers to drawing something in the graphics window. Generally
you plot one type of entity (lines, elements, etc.) to the screen at a time. If you want to plot
more than one kind of entity use, “Plot → Multiplot”, which by default will plot everything
in your model at once.
Plot Controls:
This refers to how you want your “plot” to look on the screen (shaded, wire frame,
entity numbers on or off, etc). Other plot control functions include sending an image to a
graphics file or printer.
Creating Geometry:
Geometry in ANSYS is created from
“Main Menu → Preprocessor → Modeling → Create” and has the following
terminology,
KEYPOINTS: These are points, locations in 3D space.
LINES: This includes straight lines, curves, circles, spline curves, etc. Lines are typically
defined using existing keypoints.
AREAS: This is a surface. When you create an area, its associated lines and keypoints
are automatically created to border it.
VOLUMES: This is a solid. When you create a volume, its associated areas, lines and
keypoints are automatically created.
SOLID MODEL: In most packages this would refer to the volumes only, but in ANSYS
this refers to your geometry. Any geometry. A line is considered a “solid model”.
You can’t delete a child entity without deleting its parent, in other words you can’t delete
a line if it’s part of an area, can’t delete a key point if it’s the end point of a line, etc.
Boolean Operations:
Top Down style modeling can be a very convenient way to work. Instead of first
creating keypoints, then lines from those keypoints, then areas from the lines and so on
(bottom up modeling), start with volumes of basic shapes and use Boolean operations to
add them, subtract them, divide them etc. Even if you are creating a shell model, for
example a box, you could create the box as a volume (a single command) and then delete
the volume keeping the existing areas, lines and keypoints.
These kinds of operations are found under
“Main Menu → Preprocessor → Modeling → Operate →Booleans” with some common
ones being:
Add: Take two entities that overlap (or are at least touching) and make them one.
Page 22
Subtract: Subtract one entity from another. To make a hole in a plate, create the plate
(area of volume) then create a circular area or cylinder and subtract it from the plate.
Glue: Take two entities that are touching and make them contiguous or congruent so that
when meshed they will share common nodes. For example, using default mesh parameters,
Page 23
The Working Plane:
. Page 23
• You have many different areas at Z = 0 you want to constrain. You could select them all
one by one when applying the constraint, or select “By Location” beforehand, then say
“Pick All” in the picking dialog.
• You have a structure with many fastener holes that you want to constrain. Again, you could
select them all one by one when applying the constraint, or select lines “By
Length/Radius”, type in the radius of the holes to select all of them in one shot, then “Pick
All” in the picking dialog when applying the constraint.
After working with the selected set,
“Utility Menu → Select → Everything” to make the whole model active again.
Select Entities Dialog Box Terminology:
From Full: Select from the entire set of entities in the model.
Reselect: Select a subset from the currently selected entities.
Also Select: Select in addition to (from the whole model) the set you have currently
selected.
Unselect: Remove items from the selection set.
Select All: This is not the same as “Utility Menu → Select → Everything”. This selects
all of whatever entity you have specified at the top of the dialog.
Invert: Reverses the selected and unselected entities (just the entities specified at the top
of the dialog).
OK: This does the select operation (or brings up a picker dialog so that you can pick with
the mouse) and then dismisses the dialog.
Apply: This does the operation but keeps the dialog box. Typically use this so the dialog
stays active.
Replot: Replots whatever is active in the graphics window.
Plot: Plots only the entity specified at the top of the dialog.
Creating a Material:
Page 24
of Mechanical Engineering, MSEC, Bangalore.
Create the material properties for your model in
“Main Menu → Preprocessor → Material Props → Material Models”. This gives you this
dialog box where all materials can be created, Double click on items in the right hand pane
of this window to get to the type of material model you want to create. All properties can
be temperature dependant. Click OK to create the material and it will appear in the left
hand pane. Create as many different materials as you need for your analysis.
Selecting an Element Type:
ANSYS has a large library of element types. Why so many? Elements are organized
into groups of similar characteristics. These group names make up the first part of the
element name (BEAM, SOLID, SHELL, etc). The second part of the element name is a
number that is more or less (but not exactly) chronological. As elements have been created
over the past 30 years the element numbers have simply been incremented. The earliest
and simplest elements have the lowest numbers (LINK1, BEAM3, etc), the more recently
developed ones have higher numbers. The “18x” series of elements (SHELL181,
SOLID187, etc) are the newest and most modern in the ANSYS element library.
Tell ANSYS what elements you are going to use in your model using “Main Menu →
Element Type → Add/Edit/Delete”
Later, when meshing or creating elements manually you will need to tell ANSYS what
type of elements you want to create.
Creating Properties
A solid element (brick or tet) knows its thickness, length, volume, etc by virtue of
its geometry, since it is defined in 3D space. Shell, beam and link (truss) elements do not
know this information since they are a geometric idealization or engineering abstraction.
Properties in
Page 26
pick the element edge length depending on the sizes of
features in the geometry. It makes a
finer mesh around Smaller Features
in order to Capture them
adequately.
For example in the mesh below
no sizes were specified at all
except a Smart Size level of 4.
Applying Loads andBoundary Conditions:
Loads and boundary conditions can be applied in both the Preprocessor
(“Main Menu → Preprocessor → Loads → Define Loads → Apply”), and the Solution
processor (“Main Menu → Solution → Define Loads → Apply”).
1. Select the kind of constraint you want to apply.
2. Select the geometric entity where you want it applied.
3. Enter the value and direction for it.
There is no “modify” command for loads and B.C.’s. If you make a mistake simply apply
it again with a new value (the old one will be replaced if it’s on the same entity), or delete
it and reapply it.
Loads: Forces, pressures, moments, heat flows, heat fluxes, etc.
Constraints: Fixities, enforced displacements, symmetry and anti-symmetry conditions,
temperatures, convections, etc.
Although you can apply loads and boundary conditions to nodes or elements, it’s generally
better to apply all B.C.’s to your geometry. When the solve command is issued,
Page 27
they will be automatically transferred to the underlying nodes and elements. If B.C.’s are
put on the geometry, you can re-mesh that geometry without having to reapply them
Solving:
Solution is the term given to the actual simultaneous equation solving of the
mathematical model. The details of how this is done internally is beyond the scope of this
guideline but is addressed in a separate “ANSYS Tips” white paper. For the moment, it is
sufficient to say that the basic equation of the finite element method that we are solving is,
[K]{u}={F}
where [K] is the assembled stiffness matrix of the structure,
{u} is the vector of displacements at each node, and
{F} is the applied load vector. This is analogous to a simple spring and is the essence
of small deflection theory.
To submit your model to ANSYS for solving, go to “Main Menu → Solution → Solve
→ Current LS”. LS stands for load step. A load step is a loading “condition”.
This is a single set of defined loads and boundary conditions (And their associated solution
results. More on this in the next section). Within an interactive session the first solve you
do is load step 1, the next solution is load step 2, etc.
If you leave the solution processor after solving to do post-processing for example, the load
step counter gets set back to one. You can also define and solve multiple load steps all at
once.
There are several solvers in ANSYS that differ in the way that the system of equations is
solved for the unknown displacements. The two main solvers are the sparse solver and the
PCG solver. If the choice of solvers is left to “program chosen” then generally ANSYS will
use the sparse solver. The PCG (preconditioned conjugate gradient) solver works well for
models using all solid elements. From a practical perspective one thing to consider is that
the sparse solver doesn’t require a lot of RAM but swaps out to the disk a lot. Disk I/O is
very slow. If you have a solid model and lots of RAM the PCG solver could be significantly
faster since the solution runs mostly in core memory
Post processing:
The General Postprocessor is used to look at the results over the whole model at
one point in time.
This is the final objective of everything we have discussed so far; finding the
stresses, deflections, temperature distributions, pressures, etc. These results can then be
compared to some criteria to make an objective evaluation of the performance of your
design.
The solution results will be stored in the results file as result “sets”. For a linear
static analysis like we are talking about, the correlation between Load Step numbers and
Results Set numbers will be one to one as shown below.
Page 28
Only one set of results can be stored in the database at a time, so when you want to
look at a particular set, you have to read it in from the results file. Reading it in clears the
previous results set from active memory.
To read in a results set from the results file (not needed if you have run only a single
load step) use “Main Menu → General Postproc → Read Results → First Set, or By Pick”.
Most results are displayed as a contour plot as shown below. To generate a plot of stresses
use “Main Menu → General Postproc → Plot Results → Contour Plot → Nodal Solution”,
then pick the stresses you want to see
Building a Model
Building a finite element model requires more of an ANSYS user's time than any other
part of the analysis. First, you specify a jobname and analysis title. Then, you use the
PREP7 preprocessor to define the element types, element real constants, material
properties, and the
model geometry.
Specifying a Jobname and Analysis Title
This task is not required for an analysis, but is recommended.
Defining the Jobname
The jobname is a name that identifies the ANSYS job. When you define a jobname for
an analysis, the jobname becomes the first part of the name of all files the analysis creates.
(The extension or suffix for these files' names is a file identifier such as .DB.) By using a
jobname for each analysis, you insure that no files are overwritten. If you do not specify a
jobname, all files receive the name FILE or file, depending on the operating system.
Command(s): /FILNAME
GUI: Utility Menu>File>Change Jobname
Defining Element Types
The ANSYS element library contains more than 100 different element types. Each
element type has a unique number and a prefix that identifies the element category:
BEAM4, PLANE77,SOLID96, etc. The following element categories are available:
The element type determines, among other things:
A.1 .Consider the bar shown in figure below. Determine the Nodal Displacement,
Stress in each element, Reaction forces.
Procedure:
, Page 35
Step 7: General Post Processor
Element table – Define table – Add –‘ Results data item’ – By Sequence num –
LS – LS1 – ok.
Step 8: General Post Processor
Plot Results – Deformed Shape – def+undeformed – ok.
Plot results – contour plot – Line Element Results – Elem table item at node I –
LS1 – Elem table item at node J – LS1 – ok (Line Stress diagram will be displayed).
List Results – reaction solution – items to be listed – All items – ok (reaction forces
will be displayed with the node numbers).
List Results – Nodal loads – items to be listed – All items – ok (Nodal loads will
be displayed with the node numbers).
Step 9: PlotCtrls – Animate – Deformed shape – def+undeformed-ok
Ansys Theoretical
Particulars
Node 1 Node 2 Node 1 Node2
Deformation
Stress
Reaction
, Page 36
A.2.A steel rod subjected to compression is modeled by two bar elements, as shown in
Figure. Determine the nodal displacements and the axial stress in each Element.
E=207 GPa, A = 500 mm2
Procedure:
1. Utility Menu > File > Change Job name > Enter ‘Bar Analysis’
2. Utility Menu > File > Change Title > Enter Bar Analysis.
3. Preprocessor> Element type>Add/Edit/Delete>Add> Link1
4. Preprocessor> Real constant> Add/Edit/Delete>Add>Ok> Area>500
5. Preprocessor>Material Property> Young’s Modulus > Enter Ex = 2e5.
6. Preprocessor > Modeling > Create>node > In active CS > Create three nodes
at Location (0, 0, 0) (500, 0, 0) and (1000, 0, 0).
7. Preprocessor > Modeling > Create > Element > Auto numbered >
Thru node > (Pick nodes (1, 2), (2, 3) and create element).
8. Preprocessor > Load > Define load > Apply > Structural > Displacement>
Select Node number 1 and apply ‘ALL DOF’.
9. Preprocessor > Load > Define load > Apply> Structural > Force / moment on nodes >
Select node number 3 and Enter value Fx = -12000 N.
10. Solution > Solve > Current LS
11. General postprocessor > Plot Results > Deform shape
12. Postprocessor> Contour plot> Nodal solution> Displacement vector sum.
13. Postprocessor>Element Table>Define Table>Add>By sequence num>LS 1
14. Postprocessor>Element Table>List Elem table>LS 1>ok
, Page 37
Comparison between theoretical and Ansys values
Ansys Theoretical
Particulars
Node 1 Node 2 Node 3 Node 1 Node 2 Node 3
Deformation
Stress
Reaction
, Page 38
A.3. Load of P = 60 x 10 3 N is applied as shown. Determine the following, a) Nodal
Displacement, b) Stress in each member, c) Reaction Forces.
Given Data: E=20 x103N/m2. 0.3 (Poison's Ratio)
Procedure:
Ansys Theoretical
Particulars
Node 1 Node 2 Node 3 Node 1 Node 2 Node 3
Deformation
Stress
Reaction
, Page 41
A.4.For the tapered bar shown in the figure determine the displacement, stress and
reaction in the bar. Given A1 = 1000 mm2 and A2 = 500 mm2, E = 2x105 N/mm2
The Tapered bar is modified into 2 elements as shown below with modified area of cross section
(1000+500)/2 = 750 mm2
Areas of modified two elements:
A1 = (1000+750)/2 = 875 mm2
A2 = (750+500)/2 = 625 mm2
Procedure:
1. Utility Menu > File > Change Jobname
2. Utility Menu > File > Change Title
3. Preprocessor> Element type>Add/Edit/Delete>Add> Link1.
4. Preprocessor> Real constant> Add/Edit/Delete>Add>Ok> Area>
(Create 2 set of Real Constant) > A1= 875 mm2and A2 = 625 mm2.
5. Preprocessor>Material Property> Young’s Modulus > Enter Ex= 2x105
6. Preprocessor > modeling > Create > Node > In active CS> Create 3 Nodes > (0,
0, 0), (-187.5, 0, 0) and (-375, 0, 0).
7. Preprocessor>Modeling >Create>Element> Element Attributes
(You can set Real Constant number set to 1).
8. Preprocessor>Modeling >Create>Element >Auto numbered >Thru Node>
(Pick nodes and create element between nodes 1 and 2)
9. Repeat the procedures given in 7 and 8 set real constant number to 2 and create
Element between nodes 2 and 3.
10. Preprocessor>Load>Define load> Apply> Structural> Displacement> on nodes
and select node 1 and Apply ‘ALL DOF’ in the wizard.
11. Preprocessor>Load>Define load> Apply> Structural>Force / moments > on
Nodes > Select node number ‘3’ and Apply a force Fy = - 4000.
12. Solution> Solve> Current LS.
13. General postprocessor > Plot Result > Deform shape
14. Postprocessor > Contour plot > Nodal solution > Displacement vector sum.
Comparison between Ansys and Theoretical Results
Ansys Theoretical
Particulars
Node 1 Node 2 Node 3 Node 1 Node 2 Node 3
Deformation
Stress
Reaction
EXPERMENT TYPE: II
AII.Consider the four bar truss shown in figure. For the given data, find
Stress in each element, Reaction forces, Nodal displacement. E = 210 GPa, A
= 0.1 m2.
Procedure:
Page 59
Loads – Define loads – apply – Structural – Displacement – on Nodes – pick node
1 & 4 – apply – DOFs to be constrained – All DOF – ok – on Nodes – pick node 2
– apply – DOFs to be Constrained – UY – ok. Loads – Define loads – apply –
Structural – Force/Moment – on Nodes- pick node 2 – apply – direction of For/Mom
– FX – Force/Moment value – 2000 (+ve value) – ok – Structural Force/Moment –
on Nodes- pick node 3 – apply – direction of For/Mom – FY – Force/Moment value
– -2500 (-ve value) – ok.
Step 6: Solution
Solve – current LS – ok (Solution is done is displayed) – close.
Step 7: General Post Processor
Element table – Define table – Add – ‘Results data item’ – By Sequence num –
LS – LS1 – ok.
Step 8: General Post Processor
Plot Results – Deformed Shape – def+undeformed – ok.
Plot results – contour plot – Line Element Results – Elem table item at node I –
LS1 – Elem table item at node J – LS1 – ok (Line Stress diagram will be displayed).
Plot results – contour plot – Nodal solution – DOF solution – displacement vector
sum – ok. List Results – reaction solution – items to be listed
– All items – ok (reaction forces will be displayed with the node numbers).
List Results – Nodal loads – items to be listed – All items – ok (Nodal loads will
be displayed with the node numbers).
Step 91. PlotCtrls – Animate – Deformed shape – def+undeformed-ok
Page 60
1. Elemental Solution and Displacement
1. Reaction Solution
Page 61
AII.2Analyze the 2D truss as shown below. All the members have cross-sectional
Area of 5000 mm2 and are made of steel with Young’s modulus 210000 MPa. The
Settlement at support B is 10 mm. The roller at C is on a floor 45° from horizontal
direction.
(a) If the applied force P is 200 kN, determine the member forces and
Stresses.
(b) Determine the maximum value of P in which the maximum member
Force does not exceed 600 kN.
Procedure:
Click OK
Animate Deformation
PlotCtrls > Animate > Deformed Shape
List Member Forces & Stresses
Main Menu > General Postproc > Element Table > Define Element Table > Add > Select
By Sequence number in the left list box, and SMISC in the right list box. Type “1” after
the comma in the box at the bottom of the window.
Click Apply
For member stresses, choose By Sequence num> LS, 1
Main Menu > General Postproc > Element Table > List Element Table >
, Page 64
Select SMIS1 and LS1
Click OK
Now you can see the element forces and the stresses .
List the Deflections and Reaction Forces
Main Menu > General Postproc > List Results > Nodal Solution
Click DOF Solution and in the sub-list select Displacement vector sum
Click OK
(b)
From (a), the maximum member force (consider both from compression and tension) for P
equal to 200 kN is 1146.4 kN. Since this is an elastic problem, you can find the maximum
P that the member forces do not exceed 600 kN simply by
, Page 65
EXPERMENT TYPE: III
A.III.(1) In the plate with a hole under plane stress, find deformed shape of the hole
and Determine the maximum stress distribution along A-B (you may use t = 1 mm).
E = 210GPa, t = 1 mm, Poisson’s ratio = 0.3, Dia of the circle = 10 mm, Analysis
assumption – plane stress with thickness is used.
Procedure:
Page 88
Loads – Define loads – apply – Structural – Displacement – on Nodes – select box – drag
the left side of the area – apply – DOFs to be constrained – ALL DOF – ok.
Loads – Define loads – apply – Structural – Force/Moment – on Nodes – select box – drag
the right side of the area – apply – direction of For/Mom – FX – Force/Moment value –
2000 (+ve value) – ok.
Step 6: Solution
Solve – current LS – ok (Solution is done is displayed) – close.
Step 7: General Post Processor
Plot Results – Deformed Shape – def+undeformed – ok.
Plot results – contour plot – Element solu – Stress – Von Mises Stress – ok (the stress
distribution diagram will be displayed).
Step 8: PlotCtrls – Animate – Deformed shape – def+undeformed-ok
Page 89
EXPERMENT TYPE: IV
AIV.(1)A pipe of 100mm external dia. And 20mm thickness carries water at a
pressure of 20MPa. Determine the maximum and minimum intensities of hoop
stresses in the section of pipe. Also plot the variation of hoop and radial
stresses across the thickness of pipe.
Step1 : Preprocessor > Element Type > Add/Edit/Delete > Add > Solid > select Quad 4
Node 182 >
OK > Options (Element Behavior)> Axisymmetric > OK > Close
Step 2: Preprocessor > Material Props > Material Models > Structural > Linear > Elastic >
Isotropic >Enter Ex = 2.e5 and PRXY = 0.3 > OK > Close.
Note: In ANSYS the Y axis is always the axis of symmetry for axisymmetric problems.
Step 3: Preprocessor -> Create -> Areas -> Rectangle ->By 2 Corners
4. Preprocessor -> Meshing -> mesh tool -> (Click on smart size and reduce to fine mesh-3)
>pick the rectangular area>ok
Page 90
5. Solution -> Loads ->Apply ->Displacement -> On Lines (Pick the bottom line of the
rectangle)
uy = 0 along this line. This simply prevents rigid body motion in the Y direction. No other
displacement boundary conditions are required. The radial movement is prevented by the 'hoop'
tension in the cylinder.
6. Solutions -> Loads ->Apply ->Pressure -> On Lines (Pick the left hand line of the
rectangle).
Enter a pressure of 20
9. General Postprocessor -> Plot Results ->Contour plot> Element Solu >Stress>X
component
of stress>ok (Pick Sx, Sz and examine each).
Sx (Radial stresses)……………….
Page 91
The SX stress is the radial stress that is equal to the pressure (20MPa) on the interior of the
cylinder
and is zero on the exterior. We will need to examine the computed values more closely evaluate
this
further.
10. General Postprocessor -> Plot Results ->Contour plot> Element Solu >Stress>Z
component
of stress>ok (Pick Sx, Sz and examine each).
Sz (hoop stresses)……………….
SZ is the 'hoop' stress perpendicular to the plane of this rectangle and varies from
22.5613 to 42.2178 MPa in the legend above.
Page 92
11.Graph Stresses in Cross Section:
MAIN MENU: General Postproc > Path Operations > Define Path-By Nodes
DEFINE PATH: Pick nodes along the Radius (Horizontal ) from inner radius to outer radius.
Click OK.
In the Define Path Name box, type: Dist, click OK.
PATH OPERATIONS: Map Onto Path> Highlight Stress
MAP ONTO PATH: User Label Item >Hoop> StressHighlight Z-Direction SZ(Hoop) >OK
MAP ONTO PATH: User Label Item >Radial> StressHighlight X-Direction SX(Radial) >OK
PATH OPERATIONS: Plot Path Items -On Graph > Highlight SZ & SX > OK
A graph of the stresses along the chosen path will appear as shown in the figure below.
Page 93
At the inside of the cylinder SX(Radial Stresses) = 20MPa and SZ(Hoop Stresses) = 42.2178
At the outside of the cylinder SX = 0 MPa and SZ = 22.5613
Thus the ANSYS calculated results agree pretty well with the theory. Here, radial stress is
negative and Hoop stress is positive.
AIV.(2) A thick cylinder has inner and outer dia as 120mm and 180mm respectively, it is
subjected to an external pressure of 9MPa. Determine the maximum and minimum
intensities of Hoop stresses in the section of pipe also plot the variation of hoop and radial
stresses across the thickness of pipe.
Page 94
EXPERMENT TYPE: V
Procedure:
Page 96
EXPERMENT TYPE: VI
AVI.(1).Solve the 2-D heat conduction problem for the temperature distribution within the
rectangular plate. Thermal conductivity of the plate, KXX=401 W/(m-K).
Procedure:
Page 98
1. What is the total degree of freedom of Ansys commercial package.
2. What are the different menus in ANSYS.
3. What is Work space and Swap space.
4. What is default value of workspace and Swapspace?
5. What file format ANSYS can support
6. What is the current ANSYS version
7. What is optimization in ANSYS.
8. What is sub structuring and adaptive meshing?
9. Which package is better – ANSYS or NASTRAN.
10.What is Resume in ANSYS?
11.What is P-method and H-method
12.What are scalar parameters in ANSYS?
13.What are primary nodes and Secondary nodes?
14.What is Mirror (Reflection).
15.What are subroutines.
Page 99