Simulation Lab Manual New
Simulation Lab Manual New
TIRUNELVELI – 627003.
Name : .....................................................................................
Reg No : ......................................................................................
Branch : MECHANICAL ENGINEERING
Semester : III
Subject : ME 6711 / SIMULATION AND ANALYSIS
LABORATORY
College vision and mission:
Vision:
1
Our Vision is "To create innovative and vibrant young leaders and entrepreneurs in
Engineering and Technology for building India as a super knowledge power and blossom into a
University of excellence recognized globally".
Mission:
To provide education in Engineering with excellence and ethics and to reach the
unreached.
QUALITY POLICY
We aim at continuous pursuit for excellence through
Quality education tapped from National and International Resources 14
Modular approach to channelize knowledge and programmed evaluation of
knowledge accumulated
Continuous reviewable and renewal of quality systems leading to quality output
Producing Engineers with strong ethical and moral background
2
DEPARTMENT OF MECHANICAL ENGINEERING
The Department of Mechanical Engineering was started in the year 2005 and offers a full time 4
year B.E Degree Program. The Department has conducted much national level Technical Symposium and
Workshops. The Department is strongly built with excellent lab facilities and enriched by senior faculty
members. The students are given theoretical and practical input.
3
Programme Outcomes (PO’s)
S.No Programme Outcomes
Pa . Engineering Knowledge: Apply knowledge of mathematics, science, engineering
fundamentals and an engineering specialization for building engineering models.
Pb Problem Analysis: Identify and solve engineering problems reaching conclusions
using mathematics and engineering sciences.
Pc Design/Development of Solutions: Design and develop solutions for engineering
problems that meet specified needs.
Pd Conduct Investigations of Complex Problems: Conduct investigations of complex
problems including design of experiments and analysis to provide valid solutions.
Pe Modern Tool Usage: Create and apply appropriate techniques, resources, and
modern engineering tools for executing engineering activities.
Pf The Engineer and Society: Apply reasoning of the societal, safety issues and the
consequent responsibilities relevant to engineering practice.
Pg Environment and Sustainability: Understand the impact of engineering solutions in
the environment and exhibit the knowledge for sustainable development.
Ph Ethics: Apply ethical principles and commit to professional ethics, responsibilities
and norms of engineering practice.
Pi Individual and Team Work: Function effectively as an individual, and as a member
or leader in diverse teams in multi-disciplinary settings.
Pj Communication: Communicate effectively to the engineering community and the
outside world and also to write effective reports.
Pk Project Management and Finance: Understand engineering and management
principles and apply them to handle projects in multi disciplinary environments.
Pl Life-Long Learning: Recognize the need for life-long learning and apply in the
context of technological change.
4
INTRODUCTION
Finite element analysis (FEA) is the modeling of products and systems in a virtual
environment, for the purpose of finding and solving potential (or existing) structural or
performance issues. FEA is the practical application of the finite element method (FEM), which
is used by engineers and scientist to mathematically model and numerically solve very complex
structural, fluid, and multiphysics problems. FEA software can be utilized in a wide range of
industries, but is most commonly used in the aeronautical, biomechanical and automotive
industries.
Femap is an advanced engineering simulation software program that creates finite
element analysis models of complex engineering products and systems, and displays solution
results. Femap can virtually model components, assemblies or systems and determine the
behavioral response for a given operating environment.
Using Femap’s digital simulation capabilities you can:
OVERVIEW
The purpose of this manual is to help participants to build the knowledge in software tools.
This course involves use of simulation and analysis tools. This manual provides different types
of analysis and simulation listed below
Linear statics: linear analysis with applied loads and constraints that are static
Nonlinear statics and dynamics: effects due to contact (where one part of the model
comes into contact with another), nonlinear material definitions (plasticity, elasticity,
etc.) and large displacement (strains that exceed small displacement theory that limits a
linear analysis approach)
Normal modes: natural frequencies of vibration
5
Dynamic response: loads or motions that vary with time and frequency
Buckling: critical loads at which a structure becomes unstable
PREFACE
is proud to be in its 8th year of existence since its establishment in
Simulation and analysis laboratory
July 2009. It had a series of events by Starting with the analysis like stress, mode frequency and
temperature. The matlab is used to simulate spring mass system and hydraulic cylinder.
Instruction to Teachers:
Instruction to students:
Student should maintain discipline & silence inside the lab.
Follow the instructions given by the teacher
Students should bring their textbook, observation and record to the lab class and no other
materials are allowed.
HEADINGS and DETAILS should be neatly written
6
i. Aim of the experiment
ii. Apparatus / Tools / Instruments required
iii. Procedure / Theory / Algorithm / Program
iv. Neat Diagram
v. Result / discussions.
Students should maintain the cleanliness of the lab.
Be PATIENT, STEADY, SYSTEMATIC AND REGULAR.
Lab Code of Conduct:
Students should enter the lab with proper Dress code (i.e. with I.D card, shoes & lab coat)
Breakdown of apparatus /Equipment should be immediately intimated to the staff In-
charge and proper entry should be made in the breakdown Register
Students should arrange the chairs before leaving the lab.
Major Lab Equipments with Specifications:
SMPROFESSIONALI
DEVELOPMENT
IMPARTING
PEO
SKILL
S.
KNOWLEDGE
No.
PO
7
DISCIPLINEMULTI -
SMPROFESSIONALI
DEVELOPMENT
IMPARTING
PEO
SKILL
S.
KNOWLEDGE
No.
PO
8
ME6711 SIMULATION AND ANALYSIS LABORATORY
LIST OF EXPERIMENTS
A. SIMULATION
1. MATLAB basics, Dealing with matrices, Graphing-Functions of one variable and two variables
2. Use of Matlab to solve simple problems in vibration
3. Mechanism Simulation using Multibody Dynamic software
B. ANALYSIS
1. Force and Stress analysis using link elements in Trusses, cables etc.
2. Stress and deflection analysis in beams with different support conditions.
3. Stress analysis of flat plates and simple shells.
4. Stress analysis of axi – symmetric components.
5. Thermal stress and heat transfer analysis of plates.
6. Thermal stress analysis of cylindrical shells.
7. Vibration analysis of spring-mass systems.
8. Model analysis of Beams.
9. Harmonic, transient and spectrum analysis of simple systems.
TOTAL: 45 PERIODS
COURSE OUTCOMES:
CO 1: upon Completion of this course, the students can model, analyse and simulate experiments to
meet real world system and evaluate the performance.
CO 2: to give exposure to software tools needed to analyze engineering problems.
CO 3: to conduct thermal stress analysis of cylindrical shells.
CO 4: ability to use matlab to solve simple problems in vibrations.
CONTENT BEYOND SYLLABUS:
1.Simulation of hydraulic cylinder
2.Conductive heat transfer analysis in 2D component
3.Convective heat transfer analysis in 2D component
CO/PO Mapping
H/M/L indicates strength of correlation) H-High, M-Medium, L-Low
COs/POs POs
PO PO PO PO PO PO PO PO PO PO PO PO
COs
a b c d e f g h i j k l
CO 1 M M H H H H M M
CO2 M M H H H H M M
9
CO 3 M M H H H H M M
CO 4 M M H H H H M M
Direct Indirect
1. Model Examination 1. Course end survey
2. End Semester Exam
10
LIST OF EXPERIMENTS
Index
MARKS STAFF
EX.NO DATE EXERCISES PAGE NO
AWARDED SIGNATURE
11
Aim:
To determine the stress acting on a rectangular plate with a circular hole due
to the applied external load
Hardware required:
1. Pentium 4 processor.
2. 512 MB ram.
3. VGA colour monitor.
4. 2 GB hardisk free space.
5. Colour printer.
Software required:
1. Windows xpo.s.
2. FEMAP v 9.3.1 with NX NASTRAN.
Procedure:
3. To draw a Circle within a rectangle geometry -> curve -> circle ->
radius -> give x,y,z as 50,25,0 -> Ok -> 50,35,0 -> Ok.
4. To make a rectangle with a circular hole as a surface,
Geometry -> Boundary surface -> from curve -> select all -> Ok.
5. Material selection
Model -> Material -> give title -> load -> select any one available
material -> Ok -> cancel.
6. Assigning Property
Model -> property -> give title -> select the material -> give thickness
->Ok.
12
7. Meshing
Mesh -> geometry -> surface -> select all -> Ok -> Property (Take
title) -> Ok.
9. Constraints
Model -> constraints -> nodal -> give title -> Ok -> method -> on-
curve -> select any one side of the plate -> Ok -> Fixed -> Ok.
10.Loading
Model -> load -> Nodal ->give title -> Ok -> Method -> on-curve ->
Select any one side of the plate ->Ok -> select force -> give values on
Fz only -> Ok.
12.Analysis
Model -> Analysis -> new -> static -> Ok -> Analyze -> load results
-> continue.
15.Animation
Press F5 -> Animate -> Ok.
13
14
15
Result:
16
Ex. No: 2 Date:
Aim:
To determine the stress acting on a Truss due to the applied external load
Hardware required:
1. Pentium 4 processor.
2. 512 MB ram.
3. VGA colour monitor.
4. 2 GB hardisk free space.
5. Colour printer.
Software required:
1. Windows xpo.s.
2. FEMAP v 9.3.1 with NX NASTRAN.
Procedure:
4. Material selection
Model -> Material -> give title -> load -> select any one available
material -> Ok -> cancel.
5. Assigning Property
Model -> property -> give title -> select the material ->Change
Element property-> select bar element -> Define vector as feed 1 in Y
direction (Tip) ->Ok.
17
6. Meshing
Mesh -> geometry -> surface -> select all -> Ok -> Property (Take
title) -> Ok.
8. Constraints
Model -> constraints -> nodal -> give title -> Ok -> method -> on-
curve -> select top sided nodes -> Ok -> Fixed -> Ok.
9. Loading
Model -> load -> Nodal ->give title -> Ok -> Method -> on-curve ->
Select bottom side nodes ->Ok -> select force -> give values on Fy
only -> Ok.
11.Analysis
Model -> Analysis -> new -> static -> Ok -> Analyze -> load results
-> continue.
13.Animation
Press F5 -> Animate -> Ok.
18
19
Result:
20
Thus the performance of the stress analysis of a Truss was analyzed and
animated.
21
Ex. No: 3 Date:
Aim:
Hardware required:
6. Pentium 4 processor.
7. 512 MB ram.
8. VGA colour monitor.
9. 2 GB hardisk free space.
10.Colour printer.
Software required:
3. Windows xpo.s.
4. FEMAP v 9.3.1 with NX NASTRAN.
Procedure:
Geometry -> Boundary surface -> from curve -> select all -> Ok.
3. Material selection
Model -> Material -> give title -> load -> select any one available
material -> Ok -> cancel.
4. Assigning Property
22
Model -> property -> give title -> Elem/property type ->Axisymmetric->
Ok ->OK to show Axisymmetric Axis in all Views: Yes -> select the
material -> Ok.
5. Meshing
Meshes -> geometry ->Surface->Pick the boundary. -> Ok -> select the
property -> all triangles -> Ok.
6. Constraints
Model -> constraints -> nodal -> give title -> Ok ->Pick the four
nodes along the bottom of the part -> Ok ->Create Nodal
Constraints/DOFdialog box:TY-> Ok.
7. Loading
Model -> load -> Nodal ->give title -> Ok -> Method -> on-curve ->
Select inner side of the shell ->Ok -> select pressure load-> give
values -> Ok.
9. Analysis
Model -> Analysis -> new ->Static-> Ok -> Analyze -> yes (verify
that it's OK to flip the model) -> yes (force all element normals to lie
along the correct global axis) ->Scale Factor for Axisym Forces dialog
box:OK->Yes (To permanently convert the elements to triangular
elements with midside nodes) ->Yes again (to save the model so you
can view the results in the correct model).
Press F5 -> deform -> contour -> deform and contour data
->Deformation: 1. Total Translation ->Contour: 6029: Axisym Azimuth
Stress-> Ok.
23
12.Animation
Press F5 -> Animate -> Ok.
24
Result:
25
Ex. No: 4 Date:
Aim:
To determine the stress acting on abeam due to the applied external load
Hardware required:
1. Pentium 4 processor.
2. 512 MB ram.
3. VGA colour monitor.
4. 2 GB hardisk free space.
5. Colour printer.
Software required:
1. Windows xpo.s.
2. FEMAP v 9.3.1 with NX NASTRAN.
Procedure:
3. Material selection
Model -> Material -> give title -> load -> select any one available
material -> Ok -> cancel.
4. Assigning Property
26
Model -> property -> give title -> Elem/property type -> beam -> Ok
-> shape ->select rectangular bar from shape edit box -> give height
and width -> select material from the material box -> Ok.
5. Meshing
Mesh -> mesh control -> size along curve ->select the beam -> Ok ->
enter the number of elements -> Ok.
Meshes -> geometry -> curve -> select the beam -> Ok -> Property
(Take title) -> Ok -> give 0,0,0 for Base and then give 0,1,0 for Tip ->
Ok.
6. Constraints
Model -> constraint -> nodal -> give title -> Ok -> select the ends of
the beam -> Ok -> fixed -> Ok.
7. Loading
Model -> load -> Nodal ->give title -> Ok -> Method -> on-curve ->
Select the beam ->Ok -> select force -> give values on Fy only -> Ok.
9. Analysis
Model -> Analysis -> new -> static -> Ok -> Analyze -> load results
-> continue.
12.Animation
Press F5 -> Animate -> Ok.
27
28
Fixed Beam
Continuous Beam
29
Simply Supported Beam
30
Result:
Thus the performance of the stress analysis of a beam was analyzed and
animated.
31
Ex. No: 5 Date:
Aim:
Hardware required:
1. Pentium 4 processor.
2. 512 MB ram.
3. VGA colour monitor.
4. 2 GB hardisk free space.
5. Colour printer.
Software required:
1. Windows xpo.s.
2. FEMAP v 9.3.1 with NX NASTRAN.
Procedure:
Geometry -> Boundary surface -> from curve -> select all -> Ok.
3. Material selection
Model -> Material -> give title -> load -> select any one available
material -> Ok -> cancel.
4. Assigning Property
32
Model -> property -> give title -> Elem/property type -> plate -> Ok
-> Give thickness -> select the material -> Ok.
5. Meshing
Meshes -> geometry -> surface-> select all -> Ok -> select the
property -> all triangles -> Ok.
6. Constraints
Model -> constraints -> nodal -> give title -> Ok -> method -> on-
curve -> select any one side of the plate -> Ok -> Fixed -> Ok.
7. Loading
Model -> load -> Nodal ->give title -> Ok -> Method -> on-curve ->
Select any one side of the plate ->Ok -> select force -> give values on
Fz only -> Ok.
9. Analysis
Model -> Analysis -> new -> normal modes/Eigen value -> Ok ->
Analyze -> load results -> continue.
12.Animation
Press F5 -> Animate -> Ok.
33
34
Result:
35
Ex. No: 6 Date:
Aim:
Hardware required:
1. Pentium 4 processor.
2. 512 MB ram.
3. VGA colour monitor.
4. 2 GB hardisk free space.
5. Colour printer.
Software required:
1. Windows xpo.s.
2. FEMAP v 9.3.1 with NX NASTRAN.
Procedure:
3. Material selection
Model -> Material -> give title -> load -> select any one available
material -> Ok -> cancel.
4. Assigning Property
36
Model -> property -> give title -> Elem/property type -> beam -> Ok
-> shape ->select rectangular bar from shape edit box -> give height
and width -> select material from the material box -> Ok.
5. Meshing
Meshes -> mesh control -> size along curve ->select the beam -> Ok
-> enter the number of elements -> Ok.
Meshes -> geometry -> curve -> select the beam -> Ok -> Property
(Take title) -> Ok -> give 0, 0, 0 for Base and then give 0, 1, 0 for Tip
-> Ok.
6. Constraints
Model -> constraint -> nodal -> give title -> Ok -> select the ends of
the beam -> Ok -> fixed -> Ok.
7. Loading
Model -> load -> Nodal ->give title -> Ok -> Method -> on-curve ->
Select the beam ->Ok -> select force -> give values on Fy only -> Ok.
9. Analysis
Model -> Analysis -> new ->Normal modes and eigen values -> Ok
-> Analyze -> load results -> continue.
12.Animation
Press F5 -> Animate -> Ok.
37
Fixed Beam
38
Continuous Beam
Result:
Thus the performance of the mode frequencies on a beam was analyzed and
animated.
39
40
Ex. No: 7 Date:
Aim:
Hardware required:
1. Pentium 4 processor.
2. 512 MB ram.
3. VGA colour monitor.
4. 2 GB hardisk free space.
5. Colour printer.
Software required:
1. Windows xpo.s.
2. FEMAP v 9.3.1 with NX NASTRAN.
Procedure:
Geometry -> Boundary surface -> from curve -> select all -> Ok.
3. Material selection
Model -> Material -> give title -> load -> select any one available
material -> Ok -> cancel.
4. Assigning Property
41
Model -> property -> give title -> Elem/property type -> plate -> Ok
-> Give thickness -> select the material -> Ok.
5. Meshing
Meshes -> geometry -> curve -> select all -> Ok -> select the property
-> all triangles -> Ok.
6. Constraints
Model -> constraints -> nodal -> give title -> Ok -> method -> on-
curve -> select any one side of the plate -> Ok -> Fixed -> Ok.
7. Loading
Model -> load -> Nodal ->give title -> Ok -> Method -> on-curve ->
Select any one side of the plate ->Ok -> select temperature -> give the
temperature value -> Ok. Method -> on-curve -> Select another edge
->Ok -> select temperature -> give the temperature value -> -> Ok
9. Analysis
Model -> Analysis -> new -> Steady state heat transfer -> Ok ->
Analyze -> load results -> continue.
10.Loading II
Model -> Load ->From output -> Select Temperature under Nodal ->
Select Temperature for X vector -> OK -> OK.
11.Analysis II
Model -> Analysis -> new -> Static -> Ok -> Analyze -> load results
-> continue.
42
14.Animation
Press F5 -> Animate -> Ok.
43
44
Result:
45
Ex. No: 8 Date:
Aim:
Hardware required:
1. Pentium 4 processor.
2. 512 MB ram.
3. VGA colour monitor.
4. 2 GB hardisk free space.
5. Colour printer.
Software required:
1. Windows xpo.s.
2. FEMAP v 9.3.1 with NX NASTRAN.
Procedure:
Geometry -> Boundary surface -> from curve -> select all -> Ok.
3. Material selection
Model -> Material -> give title -> load -> select any one available
material -> Ok -> cancel.
4. Assigning Property
46
Model -> property -> give title -> Elem/property type ->Axisymmetric->
Ok ->OK to show Axisymmetric Axis in all Views: Yes -> select the
material -> Ok.
5. Meshing
Meshes -> geometry ->Surface->Pick the boundary. -> Ok -> select the
property -> all triangles -> Ok.
6. Constraints
Model -> constraints -> nodal -> give title -> Ok ->Pick the four
nodes along the bottom of the part -> Ok ->Create Nodal
Constraints/DOFdialog box:TY-> Ok.
7. Loading
Model -> load ->Elemental->Select All-> Ok ->Pressure:Value:100->Ok ->
select temperature -> give the temperature value ->Adjacent Faces-> Ok
->Method ->Pick an element face on the inside of the part ->Tolerance:
85 -> Ok
9. Analysis
Model -> Analysis -> new ->Static-> Ok -> Analyze -> yes (verify
that it's OK to flip the model) -> yes (force all element normals to lie
along the correct global axis) ->Scale Factor for Axisym Forces dialog
box:OK->Yes (To permanently convert the elements to triangular
elements with midside nodes) ->Yes again (to save the model so you
can view the results in the correct model).
Press F5 -> deform -> contour -> deform and contour data
->Deformation: 1. Total Translation ->Contour: 6029: Axisym Azimuth
Stress-> Ok.
47
Press F6 -> Tools -> view style ->Postporcessing->Unreformed model ->
Ok.
12.Animation
Press F5 -> Animate -> Ok.
48
Result:
49
Ex. No: 9 Date:
Aim:
Hardware required:
1. Pentium 4 processor.
2. 512 MB ram.
3. VGA colour monitor.
4. 2 GB hardisk free space.
5. Colour printer.
Software required:
1. Windows xpo.s.
2. FEMAP v 9.3.1 with NX NASTRAN.
Procedure:
Geometry -> Boundary surface -> from curve -> select all -> Ok.
3. Material selection
Model -> Material -> give title -> load -> select any one available
material -> Ok -> cancel.
4. Assigning Property
50
Model -> property -> give title -> Elem/property type -> plate -> Ok
-> Give thickness -> select the material -> Ok.
5. Meshing
Meshes -> geometry -> curve -> select all -> Ok -> select the property
-> all triangles -> Ok.
6. Constraints
Model -> constraints -> nodal -> give title -> Ok -> method -> on-
curve -> select any one side of the plate -> Ok -> Fixed -> Ok.
7. Loading
Model -> load -> Nodal ->give title -> Ok -> Method -> on-curve ->
Select any one side of the plate ->Ok -> select temperature -> give the
temperature value -> Ok. Method -> on-curve -> Select another edge
->Ok -> select temperature -> give the temperature value -> -> Ok
9. Analysis
Model -> Analysis -> new -> Steady state heat transfer -> Ok ->
Analyze -> load results -> continue.
12.Animation
Press F5 -> Animate -> Ok.
51
52
Result:
53
54
Ex. No: 10 Date:
Aim:
Hardware required:
1. Pentium 4 processor.
2. 512 MB ram.
3. VGA colour monitor.
4. 2 GB hardisk free space.
5. Colour printer.
Software required:
1. Windows xpo.s.
2. FEMAP v 9.3.1 with NX NASTRAN.
Procedure:
3. Assigning Property
Model -> property -> give title -> Elem/property type -> plate -> Ok
-> Give thickness -> select the material -> Ok.
4. Meshing
Meshes -> geometry -> curve -> select all -> Ok -> select the property
-> all triangles -> Ok.
5. Constraints
55
Model -> constraints -> nodal -> give title -> Ok -> method -> on-
curve -> select any one side of the plate -> Ok -> Fixed -> Ok.
6. Loading
Model -> load -> elemental ->select all -> ok -> in create loads on
elements dialog box -> highlight convection from the selection list ->
Enter "4.1E-5" into the Coefficient field, then -> Enter "25" in the
Temperature field -> Ok -> Choose Adjacent Faces radio button in
Method section -> Click in Face field to activate -> Select the top
Element face of any Element on the top surface of the circuit board
->Ok.
8. Analysis
Model -> Analysis -> new -> Steady state heat transfer -> Ok ->
Analyze -> load results -> continue.
11.Animation
Press F5 -> Animate -> Ok.
56
Result:
57
58
Ex. No: 11 Date:
Aim:
Hardware required:
1. Pentium 4 processor.
2. 512 MB ram.
3. VGA colour monitor.
4. 2 GB hardisk free space.
5. Colour printer.
Software required:
1. Windows xpo.s.
2. FEMAP v 9.3.1 with NX NASTRAN.
Procedure:
FEMAP93/Examples/Dynamics/hinge.neu
Locate hinge.neu
Click Open
59
Click OK
3. Constraints
Click OK
Select the nodes around the edge of the hole either one node at a time
or using alternative picking methods discussed in previous examples.
Click OK
Click Cancel
4. Analysis
Model-> Analysis
60
Select "36.NX Nastran" from the Analysis Program drop down menu,
then...
Click Next
Click OK
Model, Function
X = 1, Y = 1
Click Cancel
61
Click OK
Click OK
Click Cancel
Model, Function
Select "7.Critical Damping vs. Freq" from Type drop down menu.
X = 1, Y = 0.1
62
Click Cancel
Select "1.Mode 1, 89.81622 Hz" from the First Freq drop down menu,
then...
Select "3.Mode 3, 568.1597 Hz" from the Last Freq drop down menu
Click OK
Click OK
Press F5.
63
Choose XY of Function radio button
Click OK
Click OK
Model, Analysis
Select "36..NX Nastran" from the Analysis Program drop down menu,
then...
Click Next
Click the Browse button ("..." button) and locate the "*.MASTER"
file that was created earlier in the tutorial. Remember, this should be
located in the same directory as your Normal Modes Analysis Results
File. Then.
64
Click Open
Click OK
Press F5
Enter "44" in the Node field located in the Output Location section
Enter "11" in the from field in the Show Output Sets section
65
Result:
Thus the Harmonic response of the 2D component has been analyzed and
animated.
66
Ex. No: 12 Date:
Aim:
Hardware Required
Software Required
Problem Description
Solve the matrix problem for the AX=B. The elements for matrix [A] and
[B] have been given. Find the values of [X] matrix
Procedure
1. Open Matlab
67
A=[1 2 3; 0 1 -1; 4 2 0]
1 2 3
0 1 -1
4 2 0
B=[20;10;40]
20
10
40
inv(A)
ans=
AX=B
X=inv(A)*B
X=
5.5556
8.8889
-1.1111
Result
68
Ex. No: 13 Date:
Aim:
Hardware Required
Software Required
Problem Description
Procedure
1. Open Matlab
69
ONE VARIABLE
>> x=[-3:1:3]
x= -3 -2 -1 0 1 2 3
>>y=x.^2
y= 9 4 1 0 1 4 9
>> plot(x,y)
TWO VARIABLES
>> x= linespace(-3:3);
>> y= linespace(-3:3);
>>[x,y]=meshgrid(x,y);
>>z=x.^2+y.^2;
>>mesh(x,y,z)
ONE VARIABLE
70
TWO VARIABLES
Result
71
Ex. No: 14 Date:
Aim:
Hardware Required
Software Required
Problem Description
Procedure
1. Open Matlab
3. Go to Library Browser
72
73
Result
Thus the Simple Cam follower and belt system is simulated using MATLAB
Simulink.
74
Ex. No: 15 Date:
Aim:
Hardware Required
Software Required
Problem Description
Procedure
1. Open Matlab
3. Go to Library Browser
75
76
77
Result
78
VIVA QUESTIONS
A node is a specific point in the finite element at which the value of the field
variable is to be determined.
Nodes are the selected finite points at which basic unknowns (displacements
in elasticity problems) are to be determined in the finite element analysis
5. What is convergence?
Convergence refers to how close the FEM solution is to the exact solution
7. What is p-convergence?
79
Large elements and complex shape functions are used in p-method
problems. In order to increase the accuracy of the solution, the complexity of the
shape function must be increased. The mesh does not need to be changed when
using the p-method.
8. What is h convergence?
Simple shape functions and many small elements are used in h-method
problems. In order to increase the accuracy of the solution, more elements must be
added. This means creating a finer mesh.
When using an h-method finite element program (such as ANSYS), the user
must run two or more solutions to ensure that the solution has converged. The user
runs the solution with one mesh and then changes the mesh and reruns the solution.
80
10. Give example for higher order elements.
For an element, Stiffness matrix is a matrix such that { f } = [K] {Q}, [K]
relates nodal displacements to nodal force of a single element.
81
The displacement function, uniquely defines strain within an element in
terms of nodal displacements.
Plane stress is defined to be a state of stress in which the normal stress and
shear stress directed perpendicular to the plane are assumed to be zero e.g. thin
plate with hole
Plane strain is defined to be a state of strain in which normal strain and shear
strain normal to the XY plane are assumed to be zero.
82
27. Compare FEA with solid mechanics.
Finitie element analysis can be applied to any continous matter where you
can divide the situation into small elements (usually triangular) and apply a set of
edge constraints and then use a computer to solve for the area of concern for
whatever the value under investigation is e.g. temperarture, flow rate, stress, shear,
bending moment etc.
So Solid mechanics is the study of things as shear, stress, etc. and they use
FEA as a tool but FEA can be applied to many other fields e.g fluid mechanics
thermodynamics, etc.
83
For the convenience of deriving element properties, in FEM many times for
each element a separate coordinate system is used known as local coordinate
system
The coordinate system used to define the points in the entire structure is
called global coordinate system.
For the convenience of deriving element properties, in FEM many times for
each element a separate coordinate system is used known as local coordinate
system
The coordinate system used to define the points in the entire structure is
called global coordinate system.
Function which relates the field variable at any point within the element to
the field variables of nodal points is called shape function.
-1 to +1
41. Explain one point formula and Explain two point formula.
84
w1f(ξ1)+w2f(ξ2), w1= 1,ξ1 = 1/√3, w2= 1,
2 point formula ∫ ( ) ξ2 = -1/√3
42. Why we are using polynomial equation in FEA?
The property of the material is not same along all the directions; such
materials are called anisotropic elements.
The property of the material is same along all the directions; such materials
are called isotropic elements.
Shape function at a specified point is unity and other than the specified point
it is zero.
85
48. Define truss.
The principle of virtual work (PVW) states that the stress, body force and
traction are in equilibrium if and only if the IVW equals the EVW for every virtual
displacement field.
In classical methods exact equations are formed and exact solutions are
obtained where as in finite element analysis exact equations are formed but
approximate solutions are obtained.
Solutions have been obtained for few standard cases by classical methods,
where as solutions can be obtained for all problems by finite element analysis.
Whenever the following complexities are faced, classical method makes the
drastic assumptions’ and looks for the solutions: Shape, Boundary conditions,
Loading
86
To get the solution in the above cases, rectangular shapes, same boundary
condition along a side and regular equivalent loads are to be assumed. In FEM no
such assumptions are made. The problem is treated as it is.
When material property is not isotropic, solutions for the problems become
very difficult in classicalmethod. Only few simple cases have been tried
successfully by researchers. FEM can handle structures with anisotropic properties
also without any difficulty.
Within the limits of elasticity, the ratio of the linear stress to the linear strain
is termed the modulus of elasticity or Young's Modulus and may be written
Young's Modulus, or E =(Stress/Strain) It is this property that determines how
much a bar will sag under its own weight or under a loading when used as a beam
within its limit of proportionality. For steel, Young's Modulus is of the order of
205000 N/mm2.
a material, around 0.28 for ordinary steels. Named after its discoverer, the French
mathematician Siméon-Davis Poisson (1781-1840).
87
(ii)Modulus of Rigidity or Shear Modulus (G)
Modulus of Rigidity is the ratio of shear stress to shear strain within elastic
limit. It is denoted by N,C or G. if q is the shear stress within elastic limit and f the
corresponding shear strain, then G = q / f.
Bulk Modulus is the ratio of volumetric stress to volumetric strain within the
elastic limit. If pv is the volumetric stress within elastic limit and ev the
corresponding volumetric strain, we have K = pv / ev.
A bar subjected to forces and couples that lie in a plane containing its
longitudinal axis is called a beam
Convection is the mode of heat transfer b/w a solid surface and the adjacent
fluid that is in motion and it involves the combined effects of conduction and fluid
motion.
88
Radiation is the energy emitted by matter in the form of electromagnetic
waves as a result of the changes in the electronic configurations of the atoms or
molecules.
Heat flux is defined as the rate of heat transfer per unit area.\
Adiabatic surfaces are surfaces which do not allow the flow of heat either
into the body or out the body.
For a fluid confined in a vessel, the rate of flow of heat out of the fluid, per
unit area of vessel wall divided by the difference between the temperature in the
interior of the fluid and the temperature at the surface of the wall. Also known as
convection coefficient.
89
Types are modal analysis, harmonic analysis etc
90