CAE Lab Manual
CAE Lab Manual
CAE Lab Manual
ANSYS
(Analysis Software)
1.0
FINITE ELEMENT CONCEPT
Finite element analysis simulates physical system and their loading conditions mathematically.
Analysis seeks to approximate the behavior of an arbitrary shaped structure under general
loading and constrain conditions. A continuum is divided into discrete number of small regions
called finite elements, whose behavior is easily understood. The entire system is then co-related
to such elements to study the integrated behavior.
1.2
ADVANTAGES OF FEM
Any complex structure can be analyzed
Different boundary conditions can be incorporated suitably
Complicated material properties such as anisotropy, non-linearity can be incorporated
The conventional method of analysis of beam, plates, shells etc are distinctly different from
one another, FEM on other hand adopts uniform approach for all type of structures
1.3
STEPS IN FEM
Discretization of continuum
Selection of displacement model
Derivation of element stiffness matrix
Assembly of element stiffness matrix & application of boundary
Solution for unknown displacements
Computation of element strains & stress from nodal displacement
P-Method (Polynomial Method): A different method to modify the subsequent finite element
analysis on the same problem is to increase the polynomial order in each element while
maintaining the original finite element size and mesh.
The increase of the interpolation order is internal, and the solution stops automatically once a
specified error tolerance is satisfied. This is known as the p-adaptivity method. The elements
associated with this capability are called the p-elements.
H-P Method: These two methods can be combined to modify the subsequent analysis on the
same model by simultaneously reducing the element size and increasing the interpolation order
in each element.
This combination is called mixed hp-adaptively.
All the
1.
2.
3.
4.
5.
1.5
PREPROCESSING
a) Create or import model geometry
b) Define material properties
c) Choose element type
d) Define geometric constants
e) Generate Finite element Mesh
2. SOLUTION
a) Apply boundary conditions
b) Apply load
c) Solve for unknowns
3. POSTPROCESSING
a) Review results like displacement, stresses, reactions etc.
b) Check validity of solution
Structural analysis is the most common application of the finite element method. The term
structural implies naval, aeronautical and mechanical & civil structures. Various types of
structural analyses are carried out using FEM.
Following are the various types of analysis
Structural Analysis
Thermal Analysis
Vibrations and Dynamics
Modal Analysis.
Buckling Analysis
Harmonic Analysis
Acoustics
Fluid flow simulations
Crash simulations
Mold flow simulations
The primary unknowns (nodal degrees of freedom) calculated in a structural analysis are
displacements. Other quantities, such as strains, stresses, and reaction forces, are then derived
from the nodal displacements.
1.6
Available FEM software Packages
ANSYS (General purpose, PC and workstations)
SDRC/I-DEAS (Complete CAD/CAM/CAE package)
NASTRAN (General purpose FEA on Mainframes)
LS-DYNA 3D (Crash/impact simulations)
ABAQUS (Nonlinear dynamic Analysis)
NISA (A General-purpose FEA tool)
PATRAN (Pre/post processor)
Expt.No: 1
Diagram:
Results:
1. The Deflection at each joint = -----------2. The stress in each member = ------------3. Reaction forces at the base = ------------
Date:
Practice Problems
1. Static Analysis of Truss Member
Aim: To perform static Analysis on Truss as shown in fig.
Discipline: Structural
Analysis Type: Static
Element Type: Link 180
Given data:
Cross-section area of Truss = 3x10-4 sq m
E= 2.07x1011 N /sq m
Diagram:
Results:
1. The Maximum stress = --------- Pa
2. Reaction forces at Node 1 and Node5.
3. Compare the above results with the theoretical values.
Diagram:
Results:
1. The Maximum stress = --------- Pa
2. Nodal Displacements
3. Reaction at the joints
Expt.No: 2
Date:
Results:
1. The Maximum stress = --------2. The Maximum displacement = ----------
Practice Problem
1. Static Analysis of 3D Space Truss
Expt.No: 3
Date:
Bending Stress Analysis on Beam
Results:
1. The Maximum stress = --------2. Deformation of the Beam
3. Maximum bending stress along the beam
4. Bending moment along the beam
Practice Problems
1. Static Analysis of Overhanging Beam
Results:
1. The Maximum stress = --------2. Deformation of the Beam
3. Reaction Forces
4. Bending moment along the beam
Results:
1. The Maximum stress = ---------
10
2. The Maximum deflections = ---------3. Compare the above results with the theoretical values.
4. Solve the problem by increasing the number of elements and compare the results
11
Expt.No: 4
Date:
Results:
1. The Maximum stress = --------2. The Maximum deflections = ----------
12
Expt.No: 5
Results:
1. The Maximum stress = --------2. The Maximum deflections = ----------
13
Date:
Expt.No:6
LOAD = 1 NEWTON
Results:
1. The fundamental frequency of the system is ------------2. Draw at least three modes of the system.
14
Date:
Expt.No: 7
Date:
Results:
1. The fundamental frequency of the system is ------------2. Draw at least three modes of frequency of the system.
3. Compare the above results with the theoretical values.
15
Expt.No: 8
Date:
Results:
1. The fundamental frequency of the system is ------------2. Draw at least three modes of frequency of the system.
3. Compare the above results with the theoretical values
Expt.No: 9
16
Date:
Test Case: This is a steady state transfer analysis of set of insulators, whose outer face
temperatures are defined. Conduction occurs in-between the insulators. The results of interest
are to plot nodal solutions and read the temperatures at the interfaces and also to obtain heat
flux value.
Results:
1. The interface temperatures = ---------------2
3. Plot the graphs between the Temperature Vs Distance & Heat flux Vs Distance.
4. Compare the above results with the theoretical results.
5. Consider Plane 2D quad-4-node element and check the results.
17
Expt.No: 10
Date:
Solidworks sketch
Results:
1. Maximum pressure = -----------2. Minimum pressure = -----------3. Maximum velocity = -------------
18
Expt.No: 11
Figure
Results:
1. Compute and plot the velocity distribution over the airfoil
2. Compute Lift and drag force
19
Date:
MATLAB Exercises
Expt.No: 12
Date:
Result:
Parameter
1.Nodal displacement
mm
2. Stress in N/m2
3. Reaction forces in N
in
MATLAB
q1
q2
q3
s1
s2
r1
r2
ANSYS
20
B)
Aim: To determine the nodal displacement, bending moments and bending stresses for a given
beam problem. Solve the problem in MATLAB and compare with ANSYS solution.
Test Case:
Result:
Parameter
1.Nodal Displacement in
MATLAB
ANSYS
mm
21
C)
Aim: To determine the displacement of node 1 and node 2 and the element stresses for plane
stress conditions using 2D CST elements. Solve the problem in MATLAB and compare with ANSYS
solution.
Test Case:
Results:
1.
Parameter
1.Displacement at node 2 in mm
2.Displacement at node 2 in mm
3.Element stress in element 1 in
mm
4.Element stress in element 2 in
mm
22
MATLAB
ANSYS
TUTORIALS
1. Nonlinear Analysis of a Cantilever Beam
Aim: To perform Nonlinear analysis of a cantilever beam
Element Type: Beam 188 and Beam 189
Discipline: Nonlinear
Analysis Type: Static
Problem Description:
Material Properties: Ex = 34000 MPa; Ey = 6530MPa; Ez= 6530MPa; vxy = 0.217;
Vyz = 0.366; vzx = 0.217; density = 2.6x10-6 kg/mm3
Diagram:
Results:
1. The Translation displacements = ---------- mm
2. The maximum Stress = ----------
23
Aim: To perform couple field analysis on a solid structure connected with a link
Element Types: Thermal mass (Link33)
Discipline: Structural-Thermal
Analysis Type: Static and Thermal
Problem Description:
Material Properties: E = 206.8 X 109 Pa; Thermal conductivity= 60.5W/m *K
Thermal Expansion Coefficient = 12e-06/K.
A steel link, with no internal stresses, is pinned between two solid structures at a reference
temperature of 0 C (273 K). One of the solid structures is heated to a temperature of 75 C (348
K). As heat is transferred from the solid structure into the link, the link will attempt to expand.
However, since it is pinned this cannot occur and as such, stress is created in the link.
Results:
1. The stress in the link = ---------------- MPa.
2. Compare the above result with theoretical result.
24
Problem Description:
A steel plate with a hole in the center is subjected to a tensile pressure. Due to the symmetric
nature of the geometry and loading, only of the structure modeled (as shown below).
Symmetry boundary conditions are to be applied, and local convergence criteria be specified.
Results:
1. The Maximum stress = --------- Pa
2. The Maximum displacement = ---------- mm
3. Stress concentration Factor deformations
4. Compare the above results with the theoretical values.
5. Solve the problem by increasing the number of elements and compare the results
25
Discipline: Thermal
Analysis Type: Steady state heat transfer (h method)
Element Type: Solid, Quad 4 node, Plane 55
Problem Description:
Material Properties: K=1.0 Btu/hr-ft-0F, hi=12 Btu/hr-ft2-0F,h0=3.0 Btu/hr-ft2-0F(1Btu=0.293
watt-hr)
Geometrical Properties: a=4 ft, b=1 ft (1 Ft=30.48 cm)
Loading: Tg=1000 F, Ta=0oF
Result:
5. Solidification of a Casting
Aim: To perform analysis during solidification of casting
26
Given data:
Material Properties for Sand: (1Btu=0.293 watt-hr), 1lb/in3=271 KN/m3, 1lb = 4.451N)
Conductivity (KXX) = 0.025Btu/ (hr-in-0 F);
Density (DENS) = 0.054 lb/in3
Specific heat (c) = 0.28 Btu/ (lb-0 F)
Conductivity for Steel (KXX):
At 00 F = 1.44Btu (hr-in-0 F) at 26430 F = 1.54 Btu (hr-in-0 F);
At 27500 F = 1.22 Btu (hr-in-0 F) at 28750 F = 1.22 Btu (hr-in-0 F)
Enthalpy for steel:
At 00 F = 0.0 Btu/in3; at 26430 F =128.1 Btu/in3; at 27500 F = 163.8 Btu/in3
At 28750 F =174.2 Btu/in3
Initial conditions:
Temperature of steel = 28750 F; Temperature of sand = 800 F
Convection Properties:
Film coefficient = 0.014 Btu (hr-in2 - 0 F)
Ambient temperature =800 F
Test Case: This is a transient heat transfer analysis of a casting process. The objective is to
track the temperature distribution in the steel casting and the mold during the solidification
process, which occurs over duration of 3 hours. The casting is made in an L-shaped sand mold
with 4 inch thick walls. Conduction occurs between the sand mold and the ambient air.
Results:
1. The maximum temperature is -----------------------2. Indicate the temperature distribution
27
Results:
1. Indicate the maximum stress locations in the component.
2. The Von-Misses Stress ------------------3. The maximum deflection in the component ---------------------
VIVA Questions
1. What are the different approximate solution methods?
2. What do you mean by continuum?
3. Define term node?
4. Define term element?
28
5. What is convergence?
6. What are the types of convergence?
7. What is p-convergence?
8. What is h convergence?
9. What is a higher order element?
10. Give example for higher order elements.
11. What do you mean by compatible elements?
12. What is geometric invariance?
13. Why do we use Pascals triangle in FEA?
14. What are the steps involved in FEA?
15. What is stiffness matrix?
16. How to obtain stiffness matrix?
17. What are the properties of stiffness matrix?
18. What is displacement function?
19. How to identify order of elements?
20. Mention different types of elements.
21. Mention some application of FEA.
22. What is connectivity?
23. What are the methods to improve problem solution?
24. Define symmetry in matrix.
25. What is plane stress?
26. What is plane strain?
27. Compare FEA with solid mechanics.
28. What are the packages available for FEA?
29. Define potential energy.
30. Define minimum potential energy.
31. Write potential energy equation for cantilever beam.
32. Mention two different methods to approach the model of physical system.
33. Difference between global coordinate and local coordinate?
34. What is local coordinate?
35. What is global coordinate?
36. What is shape function?
37. What are two general natural coordinate?
38. Mention the range of natural coordinate.
39. Number of shape function in CST
40. Number of shape function in quadrilateral.
41. Explain one point formula and Explain two point formula.
42. Why we are using polynomial equation in FEA?
43. Mention two schemes to represent band width?
44. What are forces involved in work potential?
45. What are anisotropic elements?
46. What are isotropic elements?
29
30