Department of Mechanical Engineering: Laboratory Manual
Department of Mechanical Engineering: Laboratory Manual
Department of Mechanical Engineering: Laboratory Manual
MECHANICAL ENGINEERING
LABORATORY MANUAL
Name : _______________________________________
Section : _______________________________________
ME6711-SIMULATION AND ANALYSIS LABORATORY
CONTENTS
Analysis
Study of Basics in ANSYS
1 7
Stress analysis of a plate with a circular hole
2 13
Aim:
To study about the basic procedure to perform the analysis in ANSYS
CONTACt PLANE
HYPERelastic SOLID
INFINite SOURCe
INTERface SURFace
LINK TARGEt
MASS TRANSducer
MATRIX USER
VISCOelastic (or viscoplastic)
The element type determines, among other things:
The degree-of-freedom set (which in turn implies the discipline - structural, thermal, magnetic,
electric, quadrilateral, brick, etc.)
Whether the element lies in 2-D or 3-D space.
5. Defining Element Real Constants:
Element real constants are properties that depend on the element type, such as cross-
sectional properties of a beam element. For example, real constants for BEAM3, the 2-D beam
element, are area (AREA), moment of inertia (IZZ), height (HEIGHT), shear deflection constant
(SHEARZ), initial strain (ISTRN), and added mass per unit length (ADDMAS). Not all element types
require real constants, and different elements of the same type may have different real constant values.
6. Defining Material Properties:
Most element types require material properties. Depending on the application, material
properties can be linear (see linear material properties) or nonlinear (see nonlinear material properties).
As with element types and real constants, each set of material properties has a material
reference number. The table of material reference numbers versus material property sets is called the
material table. Within one analysis, you may have multiple material property sets (to correspond with
multiple materials used in the model). ANSYS identifies each set with a unique reference number.
7. Creating the model geometry:
Once you have defined material properties, the next step in an analysis is generating a finite
element model i.e., nodes and elements that adequately describes the model geometry. The graphic
below shows some sample finite element models.
There are two methods to create the finite element model; solid modeling and direct
generation. With solid modeling, you describe the geometric shape of your model, then instruct the
ANSYS program to automatically mesh the geometry with nodes and elements. You can control the
size and shape in the elements that the program creates. With direct generation, you "manually" define
the location of each node and the connectivity of each element.
In this step, you use the SOLUTION processor to define the analysis type and analysis
options, apply loads, specify load step options, and initiate the finite element solution. You also can
apply loads using the PREP7 preprocessor.
2. Applying Loads
The word loads as used in ANSYS documentation includes boundary conditions (constraints,
supports, or boundary field specifications) as well as other externally and internally applied loads.
Loads in the ANSYS program are divided into six categories:
DOF constraints
Forces
Surface loads
Body loads
Inertia loads
Coupled-field loads
You can apply most of these loads either on the solid model (keypoints, lines, and areas) or
the finite element model (nodes and elements).
Command(s): SOLVE
Main Menu> Solution>solve>current LS
Main Menu> Solution>solution method
When you issue this command, the ANSYS program takes model and loading information from the
database and calculates the results. Results are written to the results file (Jobname.RST,
Jobname.RTH, Jobname.RMG, or Jobname.RFL) and also to the database. The only difference is that
only one set of results can reside in the database at one time, while you can write all sets of results (for
all substeps) to the results file.
Structural analysis is probably the most the common application of the finite element
method such as piston, machine parts and tools.
Static Analysis:
It is the used to determine displacement, stress etc. under static loading conditions.
ANSYS can compute linear and non-linear types (e.g. the large strain hyper elasticity and
creep problems).
Buckling Analysis:
It is used to calculate buckling load and to determine the shape of the component after applying the
buckling load. Both linear buckling and non - linear buckling analysis are possible.
Thermal Analysis:
The steady state analysis of any solid under thermal boundary conditions calculates the
effect of steady thermal load on a system (or) component that includes the following.
a) Convection.
b) Radiation.
c) Heat flow rates.
d) Heat fluxes.
e) Heat generation rates.
f) Constant temperature boundaries.
Fluid Flow:
The ANSYS CFD offers comprehensive tools for analysis of two-dimensional and three
dimensional fluid flow fields.
Magnetic:
Magnetic analysis is done using ANSYS / Electromagnetic program. It can calculate the
magnetic field in device such as power generators, electric motor etc. Interest in
magnetic analysis is finding magnetic flux, magnetic density, power loss and magnetic
forces.
Acoustic / Vibrations:
Ansys is capable of modeling and analyzing vibration system. Acoustic is the study of the
generation, absorption and reflection of pressure waves in a fluid application.
A coupled field analysis is an analysis that takes into account the interaction
between two (or) more fields of engineering analysis. Pressure vessels, Induction
heating and micro electro mechanical systems are few examples.
Result:
Thus the basics of ANSYS are studied.
Expt. No.02 STRESS ANALYSIS OF A PLATE WITH A CIRCULAR
HOLE
Aim:
To conduct the stress analysis in a plate with a circular hole using ANSYS software
System configuration:
Ram : 8 GB
Processor : Core 2 Quad / Core 2 Duo
Operating system : Windows 7
Software : ANSYS (Version12.0/12.1)
Procedure:
The three main steps to be involved are
1. Pre Processing
2. Solution
3. Post Processing
Start - All Programs – ANSYS 12.0/12.1 - Mechanical APDL Product Launcher – Set the
Working Directory as E Drive, User - Job Name as Roll No., Ex. No. – Click Run
Preprocessing:
1. Preference - Structural- h-Method - ok
2. Preprocessor - Element type - Add/Edit/Delete – Add – Solid, 8 node 82 – ok – Option –
choose Plane stress w/thk - close
3. Real constants - Add/Edit/Delete – Add – ok – THK 0.5 – ok - close
4. Material props - Material Models – Structural – Linear – Elastic – Isotropic - EX 2e5, PRXY
0.3 - ok
5. Modeling – Create – Areas – Rectangle - by 2 corner - X=0, Y=0, Width=100, Height=50 -
ok- Circle - Solid circle - X=50, Y=25, Radius=10 - ok- operate –
Booleans – Subtract – Areas - Select the larger area (rectangle) – ok – Select Circle – Next
–ok
6. Meshing - Mesh Tool – Area – Set - Select the object – ok - Element edge length
2/3/4/5– ok - Mesh Tool -Select TRI or QUAD - Free/Mapped – Mesh - Select the
object –ok
Solution:
Post Processing:
9. General post proc - Plot Result - Contour plot - Nodal Solution – Stress –Von
mises stress - ok
10. Plot control – Animates - Mode Shape – Stress - Von mises - ok
11. Plot control – Animate - Save Animation - Select the proper location to save the file (E
drive-user) - ok
12. File – Report Generator – Choose Append – ok – Image Capture – ok - close
Result:
Thus the stress analysis of rectangular plate with a circular hole is done by using the ANSYS
Software.
Expt. No.03 STRESS ANALYSIS OF A RECTANGULAR
L BRACKET
Aim:
To conduct the stress analysis of a rectangular L section bracket using ANSYS software
System Configuration:
Ram : 8 GB
Processor : Core 2 Quad / Core 2 Duo
Operating system : Windows7
Software : ANSYS (Version12.0/12.1)
Procedure:
The three main steps to be involved are
1. Pre Processing
2. Solution
3. Post Processing
Start - All Programs – ANSYS 12.0/12.1 - Mechanical APDL Product Launcher – Set the
Working Directory as E Drive, User - Job Name as Roll No., Ex. No. – Click Run.
Preprocessing:
1. Preference - Structural- h-Method - ok
2. Preprocessor - Element type - Add/Edit/Delete – Add – Solid, 8 node 82 –ok– Option –
Choose Plane stress w/thk - close
3. Real constants - Add/Edit/Delete – Add – ok – THK 0.5 – ok - close
4. Material props - Material Models – Structural – Linear – Elastic – Isotropic - EX 2e5, PRXY
0.3 - ok
5. Modeling – Create – Key points - In active CS – enter the key point number and X, Y, Z
location for 6 key points to form the rectangular L-bracket.
Lines – lines - Straight line - Connect all key points to form as lines. Areas – Arbitrary - by
lines - Select all lines - ok. Lines - Line fillet - Select the two lines where the fillet is going to
be formed – ok – enter the Fillet radius=10- ok Areas – Arbitrary - through KPs - Select the
key points of the fillet - ok
Operate – Booleans – Add – Areas - Select the areas to be add (L Shape & fillet area) - ok.
Create – Areas – Circle - Solid circle - Enter the co-ordinates, radius of the circles at the
two ends(semicircles) -ok. Operate – Booleans – Add – Areas - Select the areas to be add
(L Shape & two circles) - ok
Create – Areas – Circle - Solid circle – Enter the coordinates, radius of the two circles
which are mentioned as holes - ok. Operate – Booleans – Subtract – Areas - Select the
area of rectangle – ok - Select the two circles - ok
6. Meshing - Mesh Tool – Area – Set - Select the object – ok - Element edge length 2/3/4/5 – ok
- Mesh Tool -Select TRI or QUAD - Free/Mapped – Mesh - Select the object - ok
Solution:
Post Processing:
9. General post proc - Plot Result - Contour plot - Nodal Solution – Stress –Von
mises stress - ok
12. Plot control – Animates - Mode Shape – Stress - Von mises - ok
13. Plot control – Animate - Save Animation - Select the proper location to save the file (E
drive-user) - ok
12. File – Report Generator – Choose Append – ok – Image Capture – ok - close
Young‟s Modulus : 200 GPa
Result:
Thus the stress analysis of rectangular L section bracket is done by using the ANSYS Software.
Expt. No.04(a) STRESS ANALYSIS OF BEAM
(CANTILEVER BEAM)
Aim:
To conduct the stress analysis in a cantilever beam using ANSYS software
System Configuration:
Ram : 8 GB
Processor : Core 2 Quad / Core 2 Duo
Operating system : Windows 7
Software : ANSYS (Version12.0/12.1)
Procedure:
The three main steps to be involved are
1. Pre Processing
2. Solution
3. Post Processing
Start - All Programs – ANSYS 12.0/12.1 - Mechanical APDL Product Launcher – Set the
Working Directory as E Drive, User - Job Name as Roll No., Ex. No. – Click Run.
Preprocessing:
1. Preference - Structural- h-Method - ok
2. Preprocessor - Element type - Add/Edit/Delete – Add – Beam, 2D elastic 3 – ok –
Options – ok - close
3. Sections – beam – Common sections – Select the correct section of the beam and input the
of “w1, w2,w3” and “t1, t2, t3” – Preview – Note down the values of area, Iyy
4. Real constants - Add/Edit/Delete – Add – ok – Enter the values of area=5500, Izz=0.133e8,
height=3 – ok -close
5. Material props - Material Models – Structural – Linear – Elastic – Isotropic - EX 2e5, PRXY 0.3
- ok
6. Modeling – Create – Key points – In active CS – Enter the values of CS of each key points –
Apply – ok. Lines – Lines – Straight line – Pick the all points – ok.
7. Meshing – Mesh attributes – All lines – ok. Meshing – Size cntrls – Manual size – Lines – All
lines – Enter the value of element edge length [or] Number of element divisions – ok. Mesh
tool – Mesh – Pick all
Solution:YARTHIPLUS.COM
8. Solution – Define Loads – Apply – Structural – Displacement - On key points – Select the 1st
key point – ALL DOF – ok. On key points – select the 2nd key point– UY – ok. Force/Moment
– On key points – Select the key point – ok – direction of force/moment FY, Value = -1,000 (-
sign indicates the direction of the force) – ok
9. Solve – Current LS – ok – Solution is done – close
Post Processing:
10. General post proc – Element table – Define table – Add – By sequence num – SMISC, 6 – ok
– SMISC,12 –ok – LS,2 – ok – LS,3 - ok – Close. Plot results – Contour plot – Nodal solution –
DOF solution – Y component of displacement – ok. Contour plot – Line Element Res – Node I
SMIS 6, Node J SMISC, 12 – ok. Contour plot– Line element Res – Node I LS 2, Node J LS 3
– ok
11. File–Report Generator–Choose Append–ok–Image Capture–ok - close
Result:
Thus the stress analysis of a cantilever beam is done by using the ANSYS Software.
Expt. No.04(b) STRESS ANALYSIS OF BEAM
(SIMPLY SUPPORTED)
Aim:
To conduct the stress analysis in a simply supported beam using ANSYS software
System Configuration:
Ram : 8 GB
Processor : Core 2 Quad / Core 2 Duo
Operating system : Windows 7
Software : ANSYS (Version12.0/12.1)
Procedure:
The three main steps to be involved are
1. Pre Processing
2. Solution
3. Post Processing
Start - All Programs – ANSYS 12.0/12.1 - Mechanical APDL Product Launcher – Set the
Working Directory as E Drive, User - Job Name as Roll No., Ex. No. – Click Run
Preprocessing:
1. Preference - Structural – h method - ok
2. Preprocessor - Element type - Add/Edit/Delete - Add - beam 2D elastic 3 –Options – ok– close
– real constant- Add/Edit/Delete- Add- area = 100, Izz = 833.33 & height =10- ok
3. Preprocessor – Material Properties – Material Model – Structural – Linear – Elastic – Isotropic –
EX 2e5, PRXY0.3 – ok
4. Preprocessor – Modeling – create – nodes – inactive CS
Node 1 X=0 Y=0
Node 2 X= 25 Y=0
Node 3 X= 50 Y=0
Node 4 X= 75 Y=0
Node 1 & 2
Node 2 & 3
Node 3 & 4
Node 4 & 5
Node 5 & 6 - ok
Solution:
7. Solution - define loads- apply- structural - displacement - on nodes - select node 1 & node 5 -
apply - UY - displacement = 0 -ok
Post Processing:
10. General post processor - plot result - deform shape - Deformed + Undeformed -ok
11. General post processor - element table - define table - add - user table for item
12. Plot result - line element result - Smax I- Smax J - first result -Evaluate table data –
Smax I, Smax J, Smin I, Smin J -ok
13. General postprocessor - list result - nodal solution - DOF solution - UY - displacement result (
Table 2)
Nodes UY
1
2
3
4
Result:
Thus the stress analysis of a simply supported beam is done by using the ANSYS Software.
Expt. No.04(c) STRESS ANALYSIS OF FIXED BEAM
Aim:
To conduct the stress analysis in a fixed beam using ANSYS software
System Configuration:
Ram : 8 GB
Processor : Core 2 Quad / Core 2 Duo
Operating system : Windows 7
Software : ANSYS (Version12.0/12.1)
Procedure:
The three main steps to be involved are
1. Pre Processing
2. Solution
3. Post Processing
Start - All Programs – ANSYS 12.0/12.1 - Mechanical APDL Product Launcher – Set the
Working Directory as E Drive, User - Job Name as Roll No., Ex. No. – Click Run.
Preprocessing:
Node 1 & 2
Node 2 & 3
Node 3 & 4
Node 4 & 5
Node 5 & 6 -ok
Solution:
7. Solution - define loads- apply- structural - displacement - on nodes - select node 1 & node 5 -
apply - UY - displacement = 0 -ok
Post Processing:
10. General post processor - plot result - deform shape - Deformed + Undeformed -ok
11. General post processor - element table - define table - add - user table for item
12. Plot result - line element result - Smax I- Smax J - first result -Evaluate table data –Smax I,
Smax J, Smin I, Smin J -ok
13. General postprocessor - list result - nodal solution - DOF solution - UY - displacement result (
Table 2)
Nodes UY
1
2
3
4
Result:
Thus the stress analysis of a fixed beam is done by using the ANSYS Software.
Expt. No.05 STRESS ANALYSIS OF AN AXI-SYMMETRIC
COMPONENT
Aim:
To obtain the stress distribution of an axisymmetric component
The model will be that of a closed tube made from steel. Point loads will be applied
at the centre of the top and bottom plate.
System Configuration:
Ram : 8 GB
Processor : Core 2 Quad / Core 2 Duo
Operating system : Windows 7
Software : ANSYS (Version12.0/12.1)
Procedure:
The three main steps to be involved are
1. Pre Processing
2. Solution
3. Post Processing
Preprocessing:
1. Utility Menu - Change Job Name - Enter Job Name. Utility Menu - File - Change
Title - Enter New Title
2. Preference - Structural –h method - ok
Rectangle X1 X2 Y1 Y2
1 0 20 0 5
2 15 20 0 100
3 0 20 95 100
Post Processing:
17. General postprocessor - List results - Nodal solution - stress - components SCOMP
Aim:
To conduct the thermal stress analysis of a 2D component by using ANSYS software
System Configuration:
Ram : 8 GB
Processor : Core 2 Quad / Core 2 Duo
Operating system : Windows 7
Software : ANSYS (Version12.0/12.1)
Procedure:
The three main steps to be involved are
1. Pre Processing
2. Solution
3. Post Processing
Start - All Programs – ANSYS 12.0/12.1 - Mechanical APDL Product Launcher – Set the
Working Directory as E Drive, User - Job Name as Roll No., Ex. No. – Click Run
Preprocessing:
1. Preference – Thermal - h-Method - ok
5.Modeling – Create – Areas - Rectangle – by 2 corners – Enter the coordinate values, height,
width - ok
6. Meshing – Mesh tool – Areas, set – select the object – ok – Element edge length 10 -
ok – Mesh tool- Tri, free - mesh – Select the object
Solution:
7. Solution – Define Loads – Apply – Structural – Displacement - On lines – Select the boundary
on the object –ok – Temperature – Uniform Temp – Enter the temp. Value 50 –ok.
8. Solve – Current LS – ok – Solution is done – close
Post Processing:
9. General post proc – Plot results – Contour plot – Nodal solution – Stress – 1st principal stress
– ok – Nodal solution – DOF Solution – Displacement vector sum - ok
10. File – Report Generator – Choose Append – ok – Image Capture – ok - close
Result:
Thus the thermal stress analysis of a 2D component is done by using the ANSYS
Software.
Expt. No.07 CONDUCTIVE HEAT TRANSFER ANALYSIS OF A 2D
COMPONENT
Aim:
To conduct the conductive heat transfer analysis of a 2D component by using ANSYS
software
System Configuration:
Ram : 8 GB
Processor : Core 2 Quad / Core 2 Duo
Operating system : Windows 7
Software : ANSYS (Version12.0/12.1)
Procedure:
The three main steps to be involved are
1. Pre Processing
2. Solution
3. Post Processing
Start - All Programs – ANSYS 12.0/12.1 - Mechanical APDL Product Launcher – Set the
Working Directory as E Drive, User - Job Name as Roll No., Ex. No. – Click Run
Preprocessing:
1. Preference – Thermal - h-Method - ok
5. Modeling – Create – Areas - Rectangle – by 2 corners – Enter the coordinate values, width - ok
6. Meshing – Mesh tool – Areas, set – select the object – ok – Element edge length 0.05 - ok –
Mesh tool- Tri, free - mesh – Select the object –ok
WWW.VIDYARTHIPLUS.COM
Solution:
7. Solution – Define Loads – Apply – Thermal – Temperature - On lines – Select the right and left
side of the object –ok – Temp. Value 100 – On lines – select the top and bottom of the object –
ok –Temp 500 – ok
Post Processing:
9. General post proc – Plot results – Contour plot – Nodal solution – DOF solution – Nodal
Temperature – ok
10. File–Report Generator–Choose Append–ok–Image Capture–ok - close
Result:
Thus the conductive heat transfer analysis of a 2D component is done by using ANSYS software.
Expt. No.08 CONVECTIVE HEAT TRANSFER ANALYSIS OF A
2D COMPONENT
Aim:
To conduct the convective heat transfer analysis of a 2D component by using ANSYS
software
System Configuration:
Ram : 8 GB
Processor : Core 2 Quad / Core 2 Duo
Operating system : Windows 7
Software : ANSYS (Version12.0/12.1)
Procedure:
The three main steps to be involved are
1. Pre Processing
2. Solution
3. Post Processing
Start - All Programs – ANSYS 12.0/12.1 - Mechanical APDL Product Launcher – Set the
Working Directory as E Drive, User - Job Name as Roll No., Ex. No. – Click Run
Preprocessing:
1. Preference – structural - h-Method - ok
5. Modeling – Create – Key points - In active CS – enter the key point number and X, Y, Z
location for 8 key points to form the shape as mentioned in the drawing. Lines – lines -
Straight line - Connect all the key points to form as lines. Areas – Arbitrary - by lines -
Select all lines - ok. .VIDYARTHIPLUS.COM
6. Meshing – Mesh tool – Areas, set – select the object – ok – Element edge length 0.05 - ok –
Mesh tool- Tri,free mesh – Select the object –ok
Solution:
7. Solution – Define Loads – Apply – Thermal – Temperature - On lines – Select the lines–ok
– Temp. Value 300 – ok – Convection – On lines – select the appropriate line – ok – Enter
the values of film coefficient 50, bulk temperature 40 –ok
Post Processing:
9. General post proc – List results – Nodal Solution – DOF Solution – Nodal temperature –
ok
10. Plot results – Contour plot – Nodal solution – DOF solution – Nodal Temperature – ok
11. File – Report Generator – Choose Append – ok – Image Capture – ok - Close
Result:
Thus the convective heat transfer analysis of a 2D component is done by using the ANSYS Software.
Expt. No.09(a) MODE FREQUENCY ANALYSIS OF A
CANTILEVER BEAM
Aim:
To conduct the Mode frequency analysis of a cantilever beam using ANSYS software
System Configuration:
Ram : 8 GB
Processor : Core 2 Quad / Core 2 Duo
Operating system : Windows 7
Software : ANSYS (Version12.0/12.1)
Procedure:
The three main steps to be involved are
1. Pre Processing
2. Solution
3. Post Processing
Start - All Programs – ANSYS 12.0/12.1 - Mechanical APDL Product Launcher – Set the
Working Directory as E Drive, User - Job Name as Roll No., Ex. No. – Click Run
Preprocessing:
Solution:
6. Solution – Define Loads – Apply – Structural – Displacement - On nodes – Select the node point
–ok – All DOF – ok- Analysis type – New analysis – Modal – ok
7. Analysis type – Analysis options – Block Lanczos – enter the value no of modes to extract as 3
or 4 or 5 – ok – End Frequency 10000 – ok.
8. Solve – Current LS – ok – Solution is done – close
Post Processing:www.VIDYARTHIPLUS.COM
9. General post proc – Read results – First set - Plot results – Deformed shape – Choose
Def+undeformed – ok. Read results – Next set - Plot results – Deformed shape – Choose
Def+undeformed – ok and so on
10. File–Report Generator–Choose Append–ok–Image Capture–ok – Close (Capture all images)
Result:
Thus the mode frequency analysis of a cantilever beam is done by using the ANSYS Software.
Expt. No.09(b) MODE FREQUENCY ANALYSIS OF A SIMPLY
SUPPORTED BEAM
Aim:
To conduct the Mode frequency analysis of a simply supported beam using ANSYS software
System Configuration:
Ram : 8 GB
Processor : Core 2 Quad / Core 2 Duo
Operating system : Windows 7
Software : ANSYS (Version12.0/12.1)
Procedure:
The three main steps to be involved are
1. Pre Processing
2. Solution
3. Post Processing
Start - All Programs – ANSYS 12.0/12.1 - Mechanical APDL Product Launcher – Set the
Working Directory as E Drive, User - Job Name as Roll No., Ex. No. – Click Run
Preprocessing:
Post Processing:www.Vidyarthiplus.Com
9. General post proc – Read results – First set - Plot results – Deformed shape – Choose
Def+undeformed – ok.Read results – Next set - Plot results – Deformed shape – Choose
Def+undeformed – ok and so on
10. File–Report Generator–Choose Append–ok–Image Capture–ok - Close.(Capture all images)
Result:
Thus the mode frequency analysis of a simply supported beam is done by using the ANSYS Software.
Expt. No.10 HARMONIC ANALYSIS OF A 2D COMPONENT
Aim:
To conduct the harmonic analysis of a 2D component by using ANSYS software
System Configuration:
Ram : 8 GB
Processor : Core 2 Quad / Core 2 Duo
Operating system : Windows 7
Software : ANSYS (Version12.0/12.1)
Procedure:
The three main steps to be involved are
1. Pre Processing
2. Solution
3. Post Processing
Start - All Programs – ANSYS 12.0/12.1 - Mechanical APDL Product Launcher – Set the
Working Directory as E Drive, User - Job Name as Roll No., Ex. No. – Click Run
Preprocessing:
1. Preprocessor - Element type - Add/Edit/Delete – Add – Beam, 2D elastic 3 – ok – Close
2. Real constants - Add/Edit/Delete – Add – ok – Area 0.1e-3, Izz 0.833e-9, Height 0.01 – ok – close
3. Material props - Material Models –Structural – Linear – Elastic - Isotropic – EX 206e9, PRXY
0.25 – ok –Density – DENS 7830 – ok
4. Modeling – Create – Key points – Inactive CS – Enter the coordinate values - ok. Lines – lines –
Straight Line – Join the two key points – ok
5. Meshing – Size Cntrls – manual size – lines – all lines – Enter the value of no of element
divisions 25 – ok. Mesh – Lines – Select the line – ok
Solution:
6. Solution - Analysis type – New analysis – Harmonic – ok. Analysis type – Analysis options –
Full, Real+ imaginary – ok– Use the default settings – ok
7. Solution – Define Loads – Apply – Structural – Displacement - On nodes – Select the node
point –ok – All DOF – ok. Force/Moment – On Nodes – select the node 2 – ok – Direction of
Viva-voce
Aim:
To study the capabilities of MATLAB software
Introduction:
Including graphical user interface building MATLAB is an interactive system whose basic data element is
an array that does not require dimensioning. It allows you to solve many technical computing problems,
especially those with matrix and vector formulations, in a fraction of the time it would take to write a
program in a scalar non-interactive language such as C or FORTRAN.
The name MATLAB stands for matrix laboratory. MATLAB was originally written to provide easy access to
matrix software developed by the LINPACK and EISPACK projects. Today, MATLAB engines incorporate
the LAPACK and BLAS libraries, embedding the state of the art in software for matrix computation.
Simulink Introduction:
Simulink is a graphical extension to MATLAB for modeling and simulation of systems. In Simulink,
systems are drawn on screen as block diagrams. Many elements of block diagrams are available, such as
transfer functions, summing junctions, etc., as well as virtual input and output devices such as function
generators and oscilloscopes. Simulink is integrated with MATLAB and data can be easily transferred
between the programs. In these tutorials, we will apply Simulink to the examples from the MATLAB
tutorials to model the systems, build controllers, and simulate the systems. Simulink is supported on Unix,
Macintosh, and Windows environments; and is included in the student version of MATLAB for personal
computers.
The idea behind these tutorials is that you can view them in one window while running Simulink in another
window. System model files can be downloaded from the tutorials and opened in Simulink. You will modify
and extend these systems while learning to use Simulink for system modeling, control, and simulation. Do
not confuse the windows, icons, and menus in the tutorialsfor your actual Simulink windows. Most images
in these tutorials are not live. They simply display what you should see in your own Simulink windows. All
Simulink operations should be done in your Simulink windows.
1. Starting Simulink
2. Model Files
3. Basic Elements
4. Running Simulations
5. Building Systems
Starting Simulink
Simulink is started from the MATLAB command prompt by entering the following command:
>> Simulink
When it starts, Simulink brings up the Simulink Library browser.
Open the modeling window with New then Model from the File menu on the Simulink Library
Browser as shown above.
You can open saved files in Simulink by entering the following command in the MATLAB command
window. (Alternatively, you can load a file using the Open option in the File menu in Simulink, or by hitting
Ctrl+O in Simulink)
>> filename
The following is an example model window.
A new model can be created by selecting New from the File menu in any Simulink window (or
by hitting Ctrl+N).
Basic Elements:
There are two major classes of items in Simulink: blocks and lines. Blocks are used to generate, modify,
combine, output, and display signals. Lines are used to transfer signals from one block to another.
Blocks:
There are several general classes of blocks:
Continuous
Discontinuous
Discrete
Look-up tables
Math operations
Model verification
Model-wide Utilities
Ports & subsystems
Signal attributes
Signal routing
Sinks: Used to output or display signals
Sources: Used to generate various signals
User-defined functions
Discrete: Linear, discrete-time system elements (transfer functions, state-space models, etc.)
Linear: Linear, continuous-time system elements and connections (summing junctions, gains, etc.)
Nonlinear: Nonlinear operators (arbitrary functions, saturation, delay, etc.)
Connections: Multiplex, Demultiplex, System Macros, etc.
Blocks have zero to several input terminals and zero to several output terminals. Unused input terminals
are indicated by a small open triangle. Unused output terminals are indicated by a small triangular point.
The block shown below has an unused input terminal on the left and an unused output terminal on the
right.
Lines:
Lines transmit signals in the direction indicated by the arrow. Lines must always transmit signals from the
output terminal of one block to the input terminal of another block. One exception to this is a line can tap
off of another line, splitting the signal to each of two destination blocks, as shown below.
Lines can never inject a signal into another line; lines must be combined through the use of a block such
as a summing junction.
A signal can be either a scalar signal or a vector signal. For Single-Input, Single-Output systems, scalar
signals are generally used. For Multi-Input, Multi-Output systems, vector signals are often used, consisting
of two or more scalar signals. The lines used to transmit scalar and vector signals are identical. The type of
signal carried by a line is determined by the blocks on either end of the line.
Running Simulations:
To run a simulation, we will work with the following model file:
simple2.mdl
Download and open this file in Simulink following the previous instructions for this file. You should see the
following model window.
Before running a simulation of this system, first open the scope window by double-clicking on the
scope block. Then, to start the simulation, either select Start from the Simulation menu (as shown below)
or hit Ctrl-T in the model window.
The simulation should run very quickly and the scope window will appear as shown below. If it doesn't, just
double click on the block labeled "scope."
Note that the simulation output (shown in yellow) is at a very low level relative to the axes of the scope. To
fix this, hit the autoscale button (binoculars), which will rescale the axes as shown below.
Note that the step response does not begin until t=1. This can be changed by double-clicking on
the "step" block. Now, we will change the parameters of the system and simulate the system again.
Double-click on the "Transfer Fcn" block in the model window and change the denominator to Re-run the
simulation (hit Ctrl-T) and you should see what appears as a flat line in the scope window. Hit the auto
scale button, and you should see the following in the scope window.
Notice that the autoscale button only changes the vertical axis. Since the new transfer function has
a very fast response, it compressed into a very narrow part of the scope window. This is not really a
problem with the scope, but with the simulation itself. Simulink simulated the system for a full ten seconds
even though the system had reached steady state shortly after one second.
To correct this, you need to change the parameters of the simulation itself. In the model window,
select Parameters from the Simulation menu. You will see the following dialog box
There are many simulation parameter options; we will only be concerned with the start and stop
times, which tell Simulink over what time period to perform the simulation. Change Start time from 0.0 to
0.8 (since the step doesn't occur until t=1.0. Change Stop time from 10.0 to 2.0, which should be only
shortly after the system settles. Close the dialog box and rerun the simulation.
After hitting the autoscale button, the scope window should provide a much better display of the
step response as shown below.
Result:
Thus the features of MATLAB are studied.
Expt. No.12 SIMULATION OF SPRING-MASS SYSTEM USING
MAT LAB
Aim:
To create a Simulink model for a mass attached to a spring with a linear damping force
System Configuration:
Ram : 8 GB
Processor : Core 2 Quad / Core 2 Duo
Operating system : Windows 7
Software : MATLAB
Procedure:
A mass on a spring with a velocity-dependent damping force and a time-dependent force acting upon it will
behave according to the following equation:
The model will be formed around this equation. In this equation, 'm' is the equivalent mass of the system; 'c' is
the damping constant; and 'k' is the constant for the stiffness of the spring. First we want to rearrange the
above equation so that it is in terms of acceleration; then we will integrate to get the expressions for velocity
and position. Rearranging the equation to accomplish this, we get:
To build the model, we start with a 'step' block and a 'gain' block. The gain block represents the mass, which
we will be equal to 5. We also know that we will need to integrate twice, that we will need to add these
equations together, and that there are two more constants to consider. The damping constant 'c' will act on
the velocity, that is, after the first integration, and the constant 'k' will act on the position, or after the second
integration. Let c = 0.35 and k = 0.5. Laying all these block out to get an idea of how to put them together, we
get:
By looking at the equation in terms of acceleration, it is clear that the damping term and spring term
are summed negatively, while the mass term is still positive. To add places and change signs of terms being
summed, double-click on the sum function block and edit the list of signs:
Once we have added places and corrected the signs for the sum block, we need only connect the
lines to their appropriate places. To be able to see what is happening with this spring system, we add a
'scope' block and add it as follows:
The values of 'm', 'c' and 'k' can be altered to test cases of under-damping, critical-damping and over-
damping. To accurately use the scope, right-click the graph and select "Autoscale".The mdl-file can now be
saved.
Result:
Then the simulation is verified for spring-mass system using MATLAB software, when the model is
run for 30 iterations.
Expt. No.13 SIMULATION OF CAM AND FOLLOWER
MECHANISM USING MAT LAB
Aim:
To simulate the cam and follower mechanism using MATLAB software
System Configuration:
Ram : 8 GB
Processor : Core 2 Quad / Core 2 Duo
Operating system : Windows 7
Software : MATLAB
Procedure:
A cam and follower system is system/mechanism that uses a cam and follower to create a specific
motion. The cam is in most cases merely a flat piece of metal that has had an unusual shape or profile
machined onto it. This cam is attached to a shaft which enables it to be turned by applying a turning action to
the shaft. As the cam rotates it is the profile or shape of the cam that causes the follower to move in a
particular way. The movement of the follower is then transmitted to another mechanism or another part of the
mechanism.
Examining the diagram shown above we can see that as some external turning force is applied to the
shaft (for example: by motor or by hand) the cam rotates with it. The follower is free to move in the Y plane
but is unable to move in the other two so as the lobe of the cam passes the edge of the follower it causes the
follower to move up. Then some external downward force (usually a spring and gravity) pushes the follower
down making it keep contact with the cam. This external force is needed to keep the follower in contact with
the cam profile.
Displacement Diagrams:
Displacement diagrams are merely a plot of two different displacements (distances). These two
dispalcements are:
Angle the cam has rotated through Distance moved by the follower
Start of the cycle 0º 0 mm
End of first half of the cycle 180º 30 mm
End of the full cycle 360º 0 mm
Uniform Velocity:
Uniform Velocity means travelling at a constant speed in a fixed direction and as long as the speed or
direction don't change then its uniform velocity. In relation to cam and follower systems, uniform velocity
refers to the motion of the follower.
Now let's consider a typical displacement diagram which is merely a plot of two different displacements
(distances). These two displacements are:
Let us consider the case of a cam imparting a uniform velocity on a follower over a displacement of 30mm for
the first half of its cycle.
We shall take the cycle in steps. Firstly if the cam has to impart a displacement of 30mm on follower over
half its cycle then it must impart a displacement of 30mm÷180º for every 1º turned by the cam i.e. it must
move the follower 0.167mm per degree turn. This distance is to much to small to draw on a displacement
diagram so we will consider the displacement of the follower at the start, at the end of the half cycle, the end
of the full cycle and at certain other intervals (these intervals or the length of these intervals will be decided
on later).
Once this is done then we can draw the displacement diagram as shown below. Note a straight line from the
displacement of the follower at the start of the motion to the displacement of the follower at the end of the
motion represents uniform velocity.