Computer Aided Engineering Lab Manual
Computer Aided Engineering Lab Manual
LAB MANUAL
OF
Third Year
Bachelor of Engineering
(Choice Based Credit System)
Mechanical Engineering
(2019 Course)
Board of Studies – Mechanical and Automobile
Engineering
(With Effect from Academic Year 2021-22)
2
Department of Mechanical Engineering, ICEM, Pune
Laboratory Plan
Prac
tical Title of Experiment Session
No.
Note:
Practical 1, 2, 3, 4 and 10 are compulsory.
Perform any three practical among 5, 6, 7, 8 and 9.
The lab report shall consist of completion of Practical’s and Presentations.
Practical examination shall be based on the practical undertaken during
the semester.
Practical examination shall cover 20 marks among practical 1 to 4 and 30
marks among practical 5 to 9.
Use of Open-Source Analysis Software is highly recommended.
3
Department of Mechanical Engineering, ICEM, Pune
PRACTICAL 1
1D Bar Element – Structural Linear Analysis
4
Department of Mechanical Engineering, ICEM, Pune
Title: 1D Bar Element – Structural Linear Analysis
Description:
Find deformation, stress, shear force and bending moment diagram for the 1D bar element using
following specifications.
Density =2700 kg/m3; Young Modulus, E=7E10 N/m2; Cross section area, A=0.5m*0.5m; Length of
Element, Le=5m
Solution:
5
Department of Mechanical Engineering, ICEM, Pune
Click on project on upper left corner to go back.
v. Click on concept on menu bar – lines from sketch – select sketch – apply – generate
vii. Click on part – line body – select cross section Rect1 – generate
6
Department of Mechanical Engineering, ICEM, Pune
ix. Close design modeler environment
7
Department of Mechanical Engineering, ICEM, Pune
Check number of nodes and element in left bottom corner. This is 1d bar element
with rectangular cross section have 1 element only.
8
Department of Mechanical Engineering, ICEM, Pune
i. On upper side select – Deformation
9
Department of Mechanical Engineering, ICEM, Pune
iii. Click on shear-moment diagram in solution – select – path – click on solve.
10
Department of Mechanical Engineering, ICEM, Pune
More Examples-
1. For the figure shown below find displacement of each node and reaction at support. E=210Gpa.
All dimensions are in mm.
2. A horizontal step bar consists of two steps and made of steel material. An axial load P = 5 KN
is applied as shown in figure. Find:
a) Displacement at nodal point
b) Stresses in each element.
c) Reactions at support.
11
Department of Mechanical Engineering, ICEM, Pune
PRACTICAL 2
Truss Analysis using 1D Element
12
Department of Mechanical Engineering, ICEM, Pune
Title: Computer program for stress analysis 2-D truss subjected to plane forces
Description:
Determine the displacement, stresses at each element and reaction at support for the
following 2-D truss subjected to plane forces.
Solution:
13
Department of Mechanical Engineering, ICEM, Pune
1. Engineering Data: Double click on engineering data to select suitable material
from library or add material and its material properties. Default material
selected is structural steel. For stress analysis, check or add isotropic elasticity
i.e., young’s modulus and poisons ratio. Default material is structural steel.
iii. Click on Draw – use line to draw given truss model in workspace of ansys
v. Click on concept on menu bar – lines from sketch – select sketch – apply – generate
vi. Click on concept – cross-section – select circular corss-section – give radius 0.01 m
vii. Click on part – line body – select cross section Rect1 – generate
14
Department of Mechanical Engineering, ICEM, Pune
ix. Close design modeler environment
15
Department of Mechanical Engineering, ICEM, Pune
To change number of elements for truss structure give element size 3000mm –
generate mesh.
16
Department of Mechanical Engineering, ICEM, Pune
iv. Select top left, middle and right point – Click on force – Apply force 50 KN,
100 KN and 50 KN respectively.
ii. Right click on solution – insert – beam tool – solve for direct stress,
minimum and maximum combine stress.
Direct Stress
17
Department of Mechanical Engineering, ICEM, Pune
Maximum Combine Stress
18
Department of Mechanical Engineering, ICEM, Pune
Again Right click on solution – insert – probe – Force reaction – select boundary
condition displacement
19
Department of Mechanical Engineering, ICEM, Pune
More Examples-
1. For a truss shown in figure, determine the displacement of nodes, stresses in members and
reactions at support. Cross-sectional area of all members is 250 mm2 and E = 201 GPa.
2. Using FEA software, list and plot the axial stresses and reactions in each member of
given truss, also find out the maximum displacement. Truss details are given below:
20
Department of Mechanical Engineering, ICEM, Pune
PRACTICAL 3
Plate/Shell Element – Structural Linear Analysis
21
Department of Mechanical Engineering, ICEM, Pune
Title: Plate with hole – Structural Linear Analysis.
Problem:
22
Department of Mechanical Engineering, ICEM, Pune
1. Engineering Data: Double click on engineering data to select suitable material
from library or add material and its material properties. Default material
selected is structural steel.
v. Click on dimensions – general – select upper horizontal and put dimension – select
left vertical line and put dimension
23
Department of Mechanical Engineering, ICEM, Pune
viii. Edit dimension in left dimension box to prepare rectangle as per given
dimension
ix. Click on Draw – circle – draw circle at centre, left and right side under rectangle
xi. Click on concept – surface from sketches- click on XY Plane – sketch – apply
24
Department of Mechanical Engineering, ICEM, Pune
xiv. Close design modeler environment
Check number of nodes and element in left bottom corner. Meshing can vary by
click on sizing- course/medium/fine.
25
Department of Mechanical Engineering, ICEM, Pune
4. Static Structural: To apply boundary condition.
i. Select edge in upper toolbox
ii. Click on static structural – select left edge – click on support – fixed support
26
Department of Mechanical Engineering, ICEM, Pune
Equivalent or Von-mises Stress = 48.17 MPa
27
Department of Mechanical Engineering, ICEM, Pune
More Examples-
1. Determine the stress distribution over 1-hole plate shown in fig. Assume
dimensions.
3. A bicycle wrench shown below is made of stainless steel with a Young’s modulus of
193 GPa and a Poisson’s ratio of 0.27. The wrench is 2 mm thick. Using ANSYS
Workbench, determine the location and magnitude of the maximum deformation
and the maximum von Mises stress under the given load and boundary conditions.
28
Department of Mechanical Engineering, ICEM, Pune
PRACTICAL 4
Beam Element – Non-Linear Buckling Analysis
29
Department of Mechanical Engineering, ICEM, Pune
Title: Beam Element – Non-Linear Buckling Analysis
Description:
A bar of length 4 meters when used as a simply supported beam and subjected to a
uniformly distributed load of 30 kN/m over the whole span deflects 15 mm at the
Centre. Determine the crippling load when it is used as a column with the following
conditions. External Diameter is 137 mm and internal diameter is 103 mm.
Solution:
⚫ W =30 kN/m = 30 N/mm, l = 4 m = 4000 mm, = 15 mm
= EI = =
One end fixed, other end free
30
Department of Mechanical Engineering, ICEM, Pune
Click on project on upper left corner to go back.
iv. Click on Draw – Circle – draw two concentric circles in workspace of ansys
v. Click on dimensions – diameter – select inner circle and put dimension – select
outer circle and put dimension
vi. Edit dimension in left dimension box to prepare circle as per given
dimension
31
Department of Mechanical Engineering, ICEM, Pune
vii. Click extrude – select geometry – apply – edit depth 4000mm
Check number of nodes and element in left bottom corner. Meshing can vary by
click on sizing- course/medium/fine.
32
Department of Mechanical Engineering, ICEM, Pune
ii. Click on static structural – select left face – click on support – fixed support
33
Department of Mechanical Engineering, ICEM, Pune
Normal Stess in Z-direction = 229.2 MPa
34
Department of Mechanical Engineering, ICEM, Pune
ii. Right click on setup – click on edit
iii. Pre-stressed static structural data will be imported automatically to buckling analysis.
v. Click on solution – take five total deformations for 4 modes – Click on Solve
Deformation at Mode 1
35
Department of Mechanical Engineering, ICEM, Pune
Deformation at Mode 2
Deformation at Mode 3
Deformation at Mode 4
36
Department of Mechanical Engineering, ICEM, Pune
8. Again, select Static Structural Analysis over Eigenvalue Buckling
37
Department of Mechanical Engineering, ICEM, Pune
ii. Click on static structural – select left face – click on support – fixed support
38
Department of Mechanical Engineering, ICEM, Pune
Force Convergence
Displacement Convergence
39
Department of Mechanical Engineering, ICEM, Pune
More Examples-
1. Non-linear analysis of L-shape frame with rectangular cross-section
2. A bar of length 5 meters when used as a simply supported beam and subjected
to a uniformly distributed load of 50 kN/m over the whole span deflects 20 mm
at the Centre. Determine the crippling load when it is used as a column with
the following conditions. Consider square cross-section of column with
30mm×30mm dimesion.
40
Department of Mechanical Engineering, ICEM, Pune
PRACTICAL 5
Thermal Analysis – Static/Transient Analysis
41
Department of Mechanical Engineering, ICEM, Pune
Title: Thermal Analysis – Static/Transient Analysis.
For the rectangular cross-session of the element B = 1m and H = 0.5m and Length of Element is 5m.
Assign thermal Properties, Thermal Conductivity = 60.5 W/m∘C; Apply a constant temperature of
100∘C to the left side of the block (face 1) and a constant temperature of 300∘C to the right side of
the block (face 3). All other faces are insulated by default. Find temperature distribution in the
element.
Solution:
42
Department of Mechanical Engineering, ICEM, Pune
Click on project on upper left corner to go back.
v. Click on concept on menu bar – lines from sketch – select sketch – apply – generate
43
Department of Mechanical Engineering, ICEM, Pune
vi. Click on concept – cross-section – select rectangular corss-section – give
dimensions
vii. Click on part – line body – select cross section Rect1 – generate
44
Department of Mechanical Engineering, ICEM, Pune
Check number of nodes and element in left bottom corner. This is 1d bar element
with rectangular cross section have 1 element only.
ii. Click on steady state thermal – select left edge – click on temperature –
3000C.
iii. Click on steady state thermal – select left edge – click on temperature –
1000C.
45
Department of Mechanical Engineering, ICEM, Pune
iv. Boundary condition complete
46
Department of Mechanical Engineering, ICEM, Pune
More Examples-
1. For the square cross-session of the element 2m × 2m and Length of Element is 5m. Assign
thermal Properties, Thermal Conductivity = 75 W/m∘C; Apply a constant temperature of 200∘C
to the left side of the block (face 1) and a constant temperature of 400∘C to the right side of the
block (face 3). All other faces are insulated by default. Find temperature distribution in the
element.
2. Check temperature distribution of the cylindrical block having external diameter 50mm, internal
diameter 30 mm and Length 5m. Assign thermal Properties, Thermal Conductivity = 75 W/m∘C;
Apply a constant temperature of 300∘C to the left side of the block (face 1).
3. Check thermal performance of piston. Prepare CAD model of piston in any modelling software
and import to the analysis software.
47
Department of Mechanical Engineering, ICEM, Pune
PRACTICAL 6
Coupled Analysis- (Structural + Thermal)
48
Department of Mechanical Engineering, ICEM, Pune
Title: Coupled Thermal-Structural Analysis using FEA software.
Problem:
Material properties
Solution:
49
Department of Mechanical Engineering, ICEM, Pune
1. Engineering Data: Double click on engineering data to select suitable material
from library or add material and its material properties. Default material
selected is structural steel.
50
Department of Mechanical Engineering, ICEM, Pune
iii. Click on sketching on left middle
vi. Edit dimension in left dimension box to prepare circle as per given
dimension
vii. Click on extrude – select sketch – add material – give depth as length of bar.
51
Department of Mechanical Engineering, ICEM, Pune
x. Click on Draw – circle – draw small circle at centre
xii. Edit dimension in left dimension box to prepare circle as per given
dimension
xiii. Click on extrude – select sketch – add material – give depth as length of bar.
52
Department of Mechanical Engineering, ICEM, Pune
Right click on mesh – insert method – Sweep – algorithm axisymmetric –
Element order linear – O-Grid division – 10.
Check number of nodes and element in left bottom corner. Meshing can vary by
click on sizing- course/medium/fine.
ii. Click on steady state thermal – select left edge – click on temperature –
3000C.
iii. Click on steady state thermal – select left edge – click on temperature –
1000C.
53
Department of Mechanical Engineering, ICEM, Pune
iv. Boundary condition complete
54
Department of Mechanical Engineering, ICEM, Pune
ii. Select right face – Apply force
Equivalent stress
Total Deformation
55
Department of Mechanical Engineering, ICEM, Pune
More Examples-
1. Thermo-mechanical analysis of disc used in disc-brake. Assume suitable
dimension and data.
56
Department of Mechanical Engineering, ICEM, Pune
PRACTICAL 7
Analysis of Machine Component using 3D Elements
57
Department of Mechanical Engineering, ICEM, Pune
Title: Stress and deflection analysis of any machine component using 3-D elements.
Diameters:
D1 = 50 mm
D2 = 25 mm
D3 = 30 mm
D4 = 15 mm
E = 200 Gpa
Solution:
58
Department of Mechanical Engineering, ICEM, Pune
1. Engineering Data: Double click on engineering data to select suitable material
from library or add material and its material properties. Default material
selected is structural steel.
59
Department of Mechanical Engineering, ICEM, Pune
iv. Click on Draw – circle – draw in two circle in left and two on right side of y-axis
with center on x-axis in workspace of ansys.
vii. Click on Dimension – horizontal – select y-axis and center of circle on both side
60
Department of Mechanical Engineering, ICEM, Pune
xi. Close design modeler environment
Right click on mesh - generate mesh – Model mesh with hexahedral element
Right click on mesh - generate mesh – Model mesh with tetrahedral element
61
Department of Mechanical Engineering, ICEM, Pune
Check number of nodes and element in left bottom corner. Meshing can vary by
click on sizing- course/medium/fine.
62
Department of Mechanical Engineering, ICEM, Pune
Figure: Equivalent stress using hexahedral element
63
Department of Mechanical Engineering, ICEM, Pune
Figure: Total deformation using hexahedral element
64
Department of Mechanical Engineering, ICEM, Pune
More Examples-
1. The following figure shows a simple bracket made of structural steel with Young’s modulus
E = 200 GPa and Poisson’s ratio ν = 0.3. The bracket has a cross-section profile shown below
and an extrusion depth of 12 cm. A through-hole of diameter 1.5 cm is centered on the
bracket’s flat ends, one on each side, as shown in the figure. Suppose that the bracket is fixed
on the bottom, and lifted by a pressure load of 1000 N/cm2 from the inside surface of the top
flat portion. Compute the deformation and stresses in the bracket using the FEA.
2. The following figure made of structural steel with young’s modulus E = 200 GPa and
Poisson’s ratio ν = 0.3. Compute the deformation and stresses in the bracket using the FEA.
The dimensions are given figure. Total length is 2.3 cm.
65
Department of Mechanical Engineering, ICEM, Pune
PRACTICAL 08
Non-Linear Analysis of Assembly using Contact Elements
66
Department of Mechanical Engineering, ICEM, Pune
Title: Non-Linear Analysis of Assembly using Contact Elements
Problem: Non-linear analysis of L-Shape cylindrical arm welded with pipe.
67
Department of Mechanical Engineering, ICEM, Pune
Click on project on upper left corner to go back.
68
Department of Mechanical Engineering, ICEM, Pune
4. Model (Meshing): Double click on model to open analysis environment
ii. Right click on mesh – insert – method – tetrahedral for L-Shape bend due to
hole slot – increase resolution in mesh setting – generate mesh
Check number of nodes and element in left bottom corner. Meshing can vary by
click on sizing- course/medium/fine.
69
Department of Mechanical Engineering, ICEM, Pune
ii. Click on static structural – select hole internal face – click on support – fixed
support
iii. Click on left face of short pipe – Load – force – Vector – component
iv. In Y-direction – (1500 N)
v. Boundary condition complete
vi. Click on analysis setting – Auto time stepping on – Define by substep (Can
also use Time) – Specify supsteps as per requirement (For Ex. As follows)
vii. Weak spring – on
viii. Large deflection – on (For non-linear analysis)
70
Department of Mechanical Engineering, ICEM, Pune
7. Results
Total Deformation
71
Department of Mechanical Engineering, ICEM, Pune
Equivalent Stress
Force Convergence
72
Department of Mechanical Engineering, ICEM, Pune
Displacement Convergence
73
Department of Mechanical Engineering, ICEM, Pune
More Example-
1. One end of a steel pin is rigidly fixed to the solid plate while its other end is
force fitted to the steel arm. The configuration is shown in Figure.
74
Department of Mechanical Engineering, ICEM, Pune
PRACTICAL 9
Modal Analysis of a Cantilever Beam
75
Department of Mechanical Engineering, ICEM, Pune
Title: Modal Analysis of a Cantilever Beam
Problem:
Modulus of Elasticity=200GPa
Poisson’s Ratio=0.3
Section of Beam=20x20mm
Density=7850kg/m3
Solution:
76
Department of Mechanical Engineering, ICEM, Pune
Click on project on upper left corner to go back.
77
Department of Mechanical Engineering, ICEM, Pune
viii. Edit dimension in left dimension box to prepare rectangle as per given
dimension
Check number of nodes and element in left bottom corner. Meshing can vary by
click on sizing- course/medium/fine.
78
Department of Mechanical Engineering, ICEM, Pune
4. Modal: To apply boundary condition.
ii. Click on Modal – select left face – click on support – fixed support
iii. Click on Analysis setting under Modal – Select Max modes to find in Details
of Analysis setting – Default 6 – Change to 5
i. Deformation – Total for First mode – Again Total for second mode – again
total deformation for all five modes
79
Department of Mechanical Engineering, ICEM, Pune
6. Click on solve on upper toolbar
Results:
TABLE 1
Frequency at each mode
Mode Frequency [Hz]
1.
1597.7
2.
3. 7215.2
4.
8615.1
5.
80
Department of Mechanical Engineering, ICEM, Pune
Figure 2: Total deformation in First mode
81
Department of Mechanical Engineering, ICEM, Pune
More Examples-
1.
2.
82
Department of Mechanical Engineering, ICEM, Pune
EXPT. NO. 10
Presentation on advanced applications of FEA, NVH, CFD, Crash,
Fatigue, Manufacturing, etc.
83
Department of Mechanical Engineering, ICEM, Pune
Select application for the below any type of analysis and prepare presentation.
1. Crash Analysis
2. Fatigue analysis
9. Multibody dynamics
84
Department of Mechanical Engineering, ICEM, Pune