FEA Lab PDF
FEA Lab PDF
1 D/F
Channabasaveshwara Institute of Technology
(Affiliated to VTU, Belagavi & Approved by AICTE, New Delhi)
(NAAC Accredited & ISO 9001:2015 Certified Institution)
NH 206 (B.H. Road), Gubbi, Tumakuru – 572 216. Karnataka.
LAB MANUAL
(2017-2018)
VI SEMESTER
NAME: __________________________________________
USN : __________________________________________
BATCH: _______________SECTION:_________________
Channabasaveshwara Institute of Technology
(Affiliated to VTU, Belagavi & Approved by AICTE, New Delhi)
(NAAC Accredited & ISO 9001:2015 Certified Institution)
NH 206 (B.H. Road), Gubbi, Tumakuru – 572 216. Karnataka.
Version 1.1
February 2018
Approved by:
HOD
Dept. of Mechanical Engineering
Channabasaveshwara Institute of Technology
(Affiliated to VTU, Belagavi & Approved by AICTE, New Delhi)
(NAAC Accredited & ISO 9001:2015 Certified Institution)
NH 206 (B.H. Road), Gubbi, Tumakuru – 572 216. Karnataka.
SYLLABUS
MODELLING AND ANALYSIS LAB (FEA)
PART – A
Study of a FEA package and modelling and stress analysis of:
1. Bars of constant cross section area, tapered cross section area and stepped bar.
2. Trusses – (Minimum 2 exercises of different types)
3. Beams – Simply supported, cantilever, beams with point load, UDL, beams with varying
load etc., (Minimum 6 exercises different nature)
4. Stress analysis of a rectangular plate with a circular hole.
PART - B
1. Thermal Analysis – 1D & 2D problem with conduction and convection boundary
conditions (Minimum 4 exercises of different types)
2. Dynamic Analysis to find
a) Fixed – fixed beam for natural frequency determination
b) Bar subjected to forcing function
c) Fixed – fixed beam subjected to forcing function
REFERENCE BOOKS:
1. A first course in the Finite element method, Daryl L Logan, Thomason, Third Edition
2. Fundaments of FEM, Hutton – McGraw Hill, 2004
3. Finite Element Analysis, George R. Buchanan, Schaum Series
Record Marks
Manual Marks
Dates
(Student)
Signature
Signature
(Faculty)
Sl.
No
Name of the Experiments
Submission of
Conduction Repetition
Record
Average
Note: If the student fails to attend the regular lab, the experiment
has to be completed in the same week. Then the
manual/observation and record will be evaluated for 50% of
maximum marks.
II
Channabasaveshwara Institute of Technology
(Affiliated to VTU, Belagavi & Approved by AICTE, New Delhi)
(NAAC Accredited & ISO 9001:2015 Certified Institution)
NH 206 (B.H. Road), Gubbi, Tumakuru – 572 216. Karnataka.
OBJECTIVES
OUTCOMES
2. General Steps
3. Bars of Constant Cross-section Area
4. Bars of Tapered Cross section Area
5. Stepped Bar
6. Trusses
7. Simply Supported Beam
8. Simply Supported Beam with Uniformly varying load
9. Simply Supported Beam with Uniformly distributed load
10. Beam with moment and overhung
11. Cantilever Beam
12. Beam with angular loads, one end hinged and at other end roller
support
13. Stress analysis of a rectangular plate with a circular hole
14. Corner angle bracket
15. Thermal analysis
Modal Analysis of Cantilever beam for natural frequency
16.
determination
17. Fixed- fixed beam subjected to forcing function
18. Bar subjected to forcing function and Part - C
19. Additional problems
20. Viva questions
V
MODELLING AND ANALYSIS LAB [15MEL68] VI SEM
PART-A
Problem 1.1: Bar of Constant Cross-section Area
Consider the bar shown in figure below. Young’s modulus is 2.1×105N/mm2 and Area is
500mm2. Determine the Nodal Displacement, Stress in each element, Reaction forces.
RESULT:
Analytical approach:
Calculation:
Displacement: ______________________
Stress: ____________________________
ANSYS Results:
ANSYS Theoretical
Deformation
Stress
Reaction
Solution: The tapered bar is modified into 2 elements as shown below with modified area of
cross section.
RESULT:
Analytical approach:
Calculation:
Displacement: ______________________
Stress: ____________________________
ANSYS Results:
ANSYS Theoretical
Deformation
Stress
Reaction
RESULT:
Analytical approach:
Calculation:
Displacement: ______________________
ANSYS Results:
ANSYS Theoretical
Deformation
Stress
Reaction
2. TRUSSES
Problem 2.1: 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.
RESULT:
Analytical approach:
Calculation:
Displacement: ______________________
Stress: ____________________________
ANSYS Results:
ANSYS Theoretical
Deformation
Stress
Reaction
Problem 2.2: Consider the two 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.
RESULT:
Analytical approach:
Calculation:
Displacement: ______________________
Stress: ____________________________
ANSYS Results:
ANSYS Theoretical
Deformation
Stress
Reaction
Problem 2.3: Determine the nodal deflections, reaction forces, and stress for the truss system
shown below (E = 200GPa, A = 3250mm2).
ANSYS Results:
ANSYS Theoretical
Deformation
Stress
Reaction
3. BEAMS
Problem 3.1: Simply Supported Beam
Compute the Shear force and bending moment diagrams for the beam shown and find the
maximum deflection. Assume rectangular c/s area of 100 mm * 100mm, Young’s modulus of
210 MPa, Poisson’s ratio 0.27.
14. Plot results – contour plot – Line Element Results – Elem table item at node I –
SMIS3 – Elem table item at node J – SMIS16 – ok (bending moment diagram will be
displayed).
15. Reaction forces: List Results – reaction solution – items to be listed – All items – ok
(reaction forces will be displayed with the node numbers).
NOTE: For Shear Force Diagram use the combination SMISC 6 & SMISC 19, for
Bending Moment Diagram use the combination SMISC 3 & SMISC 16.
16. Animation: PlotCtrls – Animate – Deformed results – DOF solution – USUM – ok.
RESULT:
Analytical approach:
Calculation:
Displacement: ______________________
Stress:_____________________
ANSYS Results:
ANSYS Theoretical
Deflection
Shear force
Bending moment
Stress
16. Reaction forces: List Results – reaction solution – items to be listed – All items – ok
(reaction forces will be displayed with the node numbers).
17. Animation: PlotCtrls – Animate – Deformed results – DOF solution – deformed +
undeformed – ok.
RESULT:
Analytical approach:
Calculation:
Deflection:______________________
Bending moment:___________________
Stress:_____________________________
ANSYS Results:
ANSYS Theoretical
Deflection
Shear force
Bending moment
Stress
RESULT:
Analytical approach:
Calculation:
Deflection:______________________
Stress:_____________________________
ANSYS results:
ANSYS Theoretical
Deflection
Shear force
Bending moment
Stress
RESULT:
Analytical approach:
Calculation:
Deflection: ______________________
Bending moment:____________________
Stress: ____________________
ANSYS Results:
ANSYS Theoretical
Deflection
SFD
BMD
Stress
RESULT:
Analytical approach:
Calculation:
Deflection: ______________________
Stress: _____________________
ANSYS Results:
ANSYS Theoretical
Deflection
SFD
BMD
Stress
15. Plot results – contour plot – Line Element Results – Elem table item at node I –
SMIS3 – Elem table item at node J – SMIS16 – ok (bending moment diagram will be
displayed).
RESULT:
Analytical approach:
Calculation:
Deflection:______________________
Stress:_____________________________
ANSYS Results:
ANSYS Theoretical
Deflection
Shear force
Bending moment
Stress
Problem 4.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.
RESULT:
Analytical approach:
Calculation:
ANSYS Results:
ANSYS Theoretical
Deformation
Stress
Problem 4.2: The corner angle bracket is shown below. The upper left hand pin-hole is
constrained around its entire circumference and a tapered pressure load is applied to the
bottom of lower right hand pin-hole. Compute Maximum displacement, Von-Mises stress.
500 – ok. Structural – Pressure – on Lines – Pick line defining bottom right part of the
circle – apply – load PRES value – 500 – optional PRES value – 50 – ok.
13. Solve – current LS – ok (Solution is done is displayed) – close.
14. Plot Results – Deformed Shape – def+undeformed – ok.
15. Plot results – contour plot – Element solu. – Stress – Von Mises Stress – ok (the stress
distribution diagram will be displayed).
16. PlotCtrls – Animate – Deformed shape – def+undeformed-ok.
RESULT:
PART B
THERMAL ANALYSIS
Problem 5.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).
RESULT:
Problem 5.2: A furnace wall is made up of silica brick (K=1.5W/moC) and outside magnesia
brick (K= 4.9 W/moC) each 10cm thick. The inner and outer surfaces are exposed to fluids at
temperatures of 820oC and 110oC respectively. The contact resistance is 0.001m2oC/W. The
heat transfer coefficient for inner and outer surfaces is equal to 35 W/m 2K. Find the heat flow
through the wall per unit area per unit time and temperature distribution across the wall.
Area= 1m2.
Given: K1=1.5W/moC, K2= 1 W/moC, K3= 4.9 W/moC, h1=h4=35 W/m2K, T1=820°C,
T6=110 0C, L1=L2= 10 cm & L= 1mm.
1. Preferences-thermal-h method-ok
2. Preprocessor-Element type-add/edit/delete-add-link, 3d conduction 33,element type
reference N0.=1-apply-link, convection 34 element type reference no.2=2-ok-close
3. Real constant- add/edit/delete-add-real constant set no=1-C/S area =1-ok-close.
4. Real constant- add/edit/delete-add-real constant set no=2-C/S area =1-ok-close.
5. Material properties-material model-thermal conductivity-isotropic-KXX=1.5-ok.
From the define material model behavior menu bar-material new model
Enter define material id=2-ok
Thermal-conductivity-isotropic-Kxx=1-ok
Define material id=3-ok
Thermal-conductivity-isotropic-Kxx=4.9-ok
Define material id=4-ok-convection or film coefficient HF= 35, close
Modeling-create-nodes-in active CS
Enter node no=1,x=0,y=0,z=0-apply
Enter node no=2, X=0.001, Y=0, Z=0-apply
Enter node no=3, X=0.101, Y=0, Z=0-apply
Enter node no=4, X=0.102, Y=0, Z=0-apply
Enter node no=5, X=0.202, Y=0, Z=0-apply
Enter node no=6, X=0.203, Y=0, Z=0-ok.
Modeling-create-element-element attributes
Enter element type no=2 LINK 34 (convection)
Material no=4 (convection or film coefficient)
Real constant set no=2 (convection)-ok
RESULT:
Analytical approach:
Calculation:
ANSYS Results:
ANSYS Theoretical
Nodal temperature
T1
T2
T3
T4
T5
T6
Heat flux
Problem 5.3: The exterior wall of a building is constructed of four materials, 12mm thick
gypsum board, 75mm thick fibre glass insulation, 20mm thick plywood and 20mm thick
hardboard. The inside and outside air temperatures are 20oC and -10oC respectively. The
convective heat transfer coefficients on the inner and outer surfaces of the wall are 6W/m2 oC
and 10 W/m2 oC respectively. Determine the heat flux and the temperature distribution. Take
K for gypsum=0.176W/moC, K for fibre glass=0.036W/moC, K for plywood=0.115 W/moC
and K for hardboard=0.215 W/moC. Area= 1m2.
1. Preferences-thermal-h method-ok
2. Element type-add/edit/delete-add-link, 3d conduction 33,element type reference
N0.=1-apply-link, convection 34 element type reference no.2=2-ok-close
3. Real constant- add/edit/delete-add-real constant set no=1-C/S area =1-ok-close.
4. Real constant- add/edit/delete-add-real constant set no=2-C/S area =1-ok-close.
5. Material properties-material model-convection film coefficient-enter HF=6-ok
From the define material model behaviour menu bar-material new model
Enter define material id=2-ok
Thermal-conductivity-isotropic-Kxx=0.176-ok
Define material id=3-ok
Thermal-conductivity-isotropic-Kxx=0.036-ok
Define material id=4-ok
Thermal-conductivity-isotropic-Kxx=0.115-ok
Define material id=5-ok
Thermal-conductivity-isotropic-Kxx=0.215-ok
Define material id=6-convection or film coefficient HF= 10, close
Modeling-create-nodes-in active CS
Enter node no=1,x=0,y=0,z=0-apply
Enter node no=2, X=0.001, Y=0, Z=0-apply
Enter node no=3, X=0.013, Y=0, Z=0-apply
Enter node no=4, X=0.088, Y=0, Z=0-apply
Enter node no=5, X=0.108, Y=0, Z=0-apply
Enter node no=6, X=0.128, Y=0, Z=0-apply
Enter node no=7, X=0.129, Y=0, Z=0-ok
Modeling-create-element-element attributes
Enter element type no=2 LINK 34 (convection)
Material no=1 (convection or film coefficient)
RESULT:
Analytical approach:
Calculation:
ANSYS Results:
ANSYS Theoretical
Nodal temperature
T1
T2
T3
T4
T5
T6
T7
Heat flux
Problem 5.4: A plane wall ‘X’ (K=75W/mK) is 60 mm thick and has volumetric heat
generation of 1.5×106W/m3. It is insulated on one side while the other side is in contact with
the surface of another wall ‘Y’ (K=150W/mK) which is 30mm thick and has no heat
generation. The free surface of wall ‘Y’ is exposed to a cooling fluid at 200C with a
convection coefficient of 950 W/m2K. Find steady state temperatures at salient points across
the composite wall. Area= 1m2.
1. Preferences-Thermal-h method-ok
2. Element type-add/edit/delete- add-solid, quad 4node 55-enter reference number=1-
select options-element behaviour,K3= plane thickness-ok
3. Real constants- add/edit/delete-add- enter real constant set no. 1, thickness=2-ok-close
4. Material properties-Material models-material number 1-thermal-conductivity-
isotropic-KXX=75-ok.
5. From the menu bar select material-new model-enter material no. ID 2=2-select
6. Material model no.2- thermal-conductivity-isotropic-KXX=150-ok.
7. Modeling-create-areas-rectangles-by dimensions-X1=0, X2=0.06, Y1=0,Y2=0.03-
apply- X1=0.06, X2=0.09, Y1=0,Y2=0.03-OK
Modeling-operate-Boolean-glue-areas-pick the material-ok
8. Meshing-size controls-manual size-picked lines-pick the first vertical line, middle line
and the last vertical line-ok-number of element divisions-2-apply
Meshing-size controls-manual size-picked lines-pick the first rectangle top and
bottom lines-number of element divisions=60-apply
Meshing-size controls-manual size-picked lines-pick the top and bottom lines of
second rectangle-number of element divisions=30-ok.
Meshing-mesh areas- free-pick all-ok.
9. Solution-Analysis type-new analysis-steady state-ok
define loads-apply-thermal-heat generated on areas-pick the first rectangle-ok
define loads-apply-thermal-heat generated on areas-apply Hgen on areas as constant
value- load Hgen value=1.5e6
define loads-apply-thermal-convection-on lines-pick the back corner line-ok-enter
film coefficient=950-bulk temperature=200C (Don’t change other attributes)-ok
RESULT:
Analytical approach:
Calculation:
ANSYS Results:
ANSYS Theoretical
Nodal temperature
T1
T2
T3
Problem 6.1: Modal Analysis of Cantilever beam for natural frequency determination.
Modulus of elasticity = 200GPa, Density = 7800 Kg/m3.
RESULT:
Analytical solution:
ANSYS Results:
RESULT:
The importing and exporting commands allow users to share 3D geometry with other
3D program. There are many 3D programs used in the world nowadays. Depending on the
purpose the software is chosen to meet the requirement. Autodesk 3ds Max is one of the
programs widely used in the media and entertainment industry all over the world. It is also
utilized in media education institutions such as, Metropolia US. It was also used in the project
that this thesis describes. There are several file formats that can be imported and exported
into 3ds Max. All acceptable files are automatically displayed in the file dialog box when
selecting the import option.
The available import formats include the following: • Autodesk (FBX) • 3D Studio
Mesh, Project and Shapes (3DS, PRJ; SHP) • Adobe Illustrator (AI) • Collada (DAE) • Land
XML/DEM/DDF • AutoCAD and Legacy AutoCAD (DWG, DXF) • Flight Studio Open
Flight (FLT) • Motion Analysis (HTR, TRC) • Initial Graphics Exchange Standard (IGE,
IGS, IGES) • Autodesk Inventor (IPT, WIRE, IAM) • Lights cape (LS, VW, LP) • OBJ
Material and Object (OBJ) • ACIS SAT (SAT) • Google Sketch Up (SKP) • Stereo Litho
(STL) • VIZ Material XML Import (XML) • STEP (STP, STEP) • Rhino (3DM).
Many available options can be imported into 3ds Max, but not all of them are
mentioned in this thesis. In the thesis, only those formats were utilized to import a 3D model
for the project of the case study are analyzed. STEP, IGES and 3DS are file formats used to
export the model. Though those format files are exported from the same CAD file, two of
them are not compatible. In the case study, the model designed by the CAD software is
utilized mostly in industries. The model cannot be used directly because it lacks textures and
materials. In order for it to work in video commercials, it needs to be transferred to other
software for more designs. Moreover, to give more options for using the model with fewer
errors and the best result, the model is exported into three different formats IGES, STEP and
3DS.
IGES File
IGES (Initial Graphics Exchange Specification) was the first specification for CAD
data exchange published in 1980 as a NBS (National Bureau of Standards) report in USA.
IGES was originally created for exchanging the drafting data like 2D/3D wireframe models,
text, dimensioning data, and a limited class of surfaces. Due to developing and requiring of
users, IGES has been ongoing improvement and now it can support more capabilities such as
entities, syntax, clarity and consistency.
The IGES specification describes the file format, language format, and the product
definition data. In the product definition, geometric, topological, and non-geometric data are
included. The geometric entities used to define the geometry are determined in the geometry
part. The topology part describes the entities that identify the relationship between the
geometric entities. There are three divisions such as annotation, definition, and organization
in the non-geometric part. The annotation category contains dimensions, drafting notations
and text. With the definition category, users can define specific properties of individuality or
collections of entities.
The organization category describes groupings of geometry, annotation, or property
elements. An IGES file has six sections, which are Flag, Start, Global, Directory Entry,
Parameter Data, and Terminate. A directory entry and parameter data entry are included in
each entity instance. The directory entry has an index and attributes to describe the data. The
parameter data defines the specific entity and is defined by fixed length records in accordance
with the corresponding entity.
The size of the IGES files and the processing time are practical problems. As IGES
files are composed of fixed format records and in both the directory entry section and the
parameter data section each entry has to have record, errors occurs in pre-and post-processor
implementations.
STEP File
STEP (Standard for the Exchange of Product model data) is a new international
standard (ISO 10303) for representing and exchanging product model information including
an object-flavored data specification language called EXPRESS. STEP also describes
implementation methods, for example, a physical transfer files, and offers different resources
such as geometric and topological representation.
STEP was developed in 1984 as a worldwide collaboration with the purpose of
defining a standard to cover all aspects of a product during its lifetime. It is a collection of
standards to represent and exchange product information. While the main parts of STEP are
already international standards, many parts remain to be under development. The
development is performed under the control of the International Standards Organization
(ISO), Technical Committee 184 (TC184, Industrial Automation Systems) and Subcommittee
4 (SC4, Industrial Data and Global Manufacturing Programming Languages).
STEP aims to offer system-independent mechanisms describe the product information
in computeraided systems throughout its lifetime. It separates the representation of product
information from the implementation methods that are used for data exchange. A basis for
archiving product information and a methodology for the conformance testing of
implementations are provided by STEP.
EXPRESS is used to specify the representation of product information. This facilitates
development of implementation and enables consistency of representation. STEP does not
only define the geometric sharpness of a product, but also includes topology, features,
tolerance specifications, material properties, so that the goal of the design, manufacturing,
testing, inspection and support of the product is completely defined. STEP covers the total
product life cycle in terms of sharing, storage and exchange. It is said that STEP is the most
important effort ever established in engineering and will take the position of current CAD
exchange standards.
STEP is well known and widely used as an exchange data form. It is supported by
many software tools such as ECAD or EDA (Altium Designer, Circuit Studio, Circuit Maker,
Cadence, etc.), MCAD (IDA step and express Engine), Dassault Systemes (Catia, Solid
Works and 3DVIA shape), PTC and Autodesk.
1. A =160 MM, B = 80 MM
2.
ADDITIONAL PROBLEMS
1. Calculate the stresses and displacement for the plate shown below. Let the load
be P = 100N applied at equal distance from both ends and E = 3e7 N/mm2.
3. Calculate the maximum value of Von-misses stresses in the stepped beam with
a rounded plate as shown in the figure. Where Young’s modulus, E=210Gpa,
Poisson’s ratio is 0.3 and the beam thickness is 10mm, the element size is 2mm
4. Loads will not be applied to the beam shown below in order to observe the
deflection caused by the weight of the beam itself. The beam is to be made of
steel with a Young’s modulus of elasticity of 200 Gpa.
5. A beam has a force of 1000N applied as shown below. The purpose of this
optimization problem is to minimize the weight of the beam without exceeding
the allowable stress. It is necessary to find the cross sectional dimensions of the
beam in order to minimize the weight of the beam. However, the width and
height of the beam cannot be smaller than 10mm. The maximum stress
anywhere in the beam cannot exceed 200 MPa. The beam is to be made of steel
with a modulus of elasticity of 200 GPa.