Ansys
Ansys
CHENNAI-600099
Date :
AIM:
To determine the displacement and bending stress of a given plate with a circular hole using Finite
Element Analysis bases ANSYS structure and view the displacement and bending stress plots.
PROCEDURE:
Fill in the window as shown above. This will create a rectangle where the bottom left corner has the
coordinates 0, 0, 0 and the top right corner has the coordinates 200, 100, 0.
Fill in the window as shown above. This will create a circle where the center has the
coordinates 100,50,0 (the center of the rectangle) and the radius of the circle is 20
mm.
a. Subtraction
Now we want to subtract the circle from the rectangle. Prior to this operation, your image should
resemble the following:
At this point a 'Subtract Areas' window will pop up and the ANSYS Input window will
display the following message: [ASBA] Pick or enter base areas from which to subtract
(as shown below)
Therefore, select the base area (the rectangle) by clicking on it. Note: The selected area will turn pink once it is
selected.
Ensure that the entire rectangular area is selected (otherwise click 'Next') and then click 'OK'. Click
Now you will be prompted to select the areas to be subtracted, select the circle by
clicking on it and then click 'OK'.
We are going to give the properties of Steel. Enter the following when prompted: EX
200000
PRXY 0.3
Mesh Size
To tell ANSYS how big the elements should be, Preprocessor > Meshing > Size Cntrls >
Manual Size > Areas > All Areas
o Select an element edge length of 25. We will return later to determine if this was adequate for
the problem.
Mesh
o In the 'Preprocessor' menu select Meshing > Mesh > Areas > Free and select the area when
prompted
o This location is fixed which means that all DOF's are constrained. Therefore, select 'All DOF'
by clicking on it and enter '0' in the Value field as shown above.
You will see some blue triangles in the graphics window indicating the displacement
contraints.
3. Apply Loads
o As shown in the diagram, there is a load of 20N/mm distributed on the right hand side of the
plate. To apply this load:
Solution > Define Loads > Apply > Structural > Pressure > On Lines
o When the window appears, select the line along the right hand edge of the plate and click
'OK'
o Calculate the pressure on the plate end by dividing the distributed load by the thickness of the
plate (1 MPa).
o Fill in the "Apply PRES on lines" window as shown below.
NOTE:
The pressure is uniform along the surface of the plate, therefore the last field is left
blank.
The pressure is acting away from the surface of the plate, and is therefore defined as
a negative pressure.
o The applied loads and constraints should now appear as shown below.
4. Solving the System
Solution > Solve > Current LS
5. Deformation
o General Postproc > Plot Results > Deformed Shape > Def + undeformd to view both the
deformed and the undeformed object.
o Observe the locations of deflection.
6. Deflection
o General Postproc > Plot Results > Nodal Solution... Then select DOF solution, USUM in
the window.
o Alternatively, obtain these results as a list. (General Postproc > List Results > Nodal
Solution...)
o Are these results what you expected? Note that all translational degrees of freedom were
constrained to zero at the left end of the plate.
7. Stresses
o General Postproc > Plot Results > Nodal Solution... Then select Stress, von Mises in the
window. You can list the von Mises stresses to verify the results at certain nodes
RESULT:
Ex. No. : 02 STRESS ANALYSIS OF SIMPLE BRACKET
Date :
AIM:
To determine the displacement and bending stress of a given simple bracket using Finite
Element Analysis bases ANSYS structure and view the displacement and bending stress plots.
PROCEDURE:
The 2D Plane Stress Bracket will introduce boolean operations, plane stress, and uniform pressure
loading. The problem to be modeled in this example is a simple bracket shown in the following figure. This
bracket is to be built from a 20 mm thick steel plate. A figure of the plate is shown below. This plate will be fixed
at the two small holes on the left and have a load applied to the larger hole on the right.
Preprocessing: Defining the Problem
1. Give the Bracket example a Title
Utility Menu > File > Change Title
2. Form Geometry
Again, Boolean operations will be used to create the basic geometry of the
Bracket.
a. Create the main rectangular shape
Preprocessor > Modeling > Create > Areas > Rectangle > By 2 Corners (X=0,Y=0, WIDTH= 80,
HEIGHT =100)
Preprocessor > Modeling > Create > Areas > Circle > Solid Circle
The center of the circle is located at (80,50) and has a radius of 50 mm The following code is used to
create a circular area:
c. Now create a second and third circle for the left hand side using the following
dimensions:
circle circle
parameter 2 3
XCENTER 0 0
YCENTER 20 80
RADIUS 20 20
d. Create a rectangle on the left hand end to fill the gap between the two small circles.
-
XCORNER
20
YCORNER 20
WIDTH 20
HEIGHT 60
Preprocessor > Modeling > Create > Areas > Circle > Solid Circle
Now select
Preprocessor > Modeling > Operate > Booleans > Subtract > Areas
Select the base areas from which to subract (the large plate that was created)
Next select the three circles that we just created. Click on the three circles that you
just created and click 'OK'.
Now you should have the following:
3. Apply Loads
As shown in the diagram, there is a single vertical load of 1000N, at the bottom of the large bolt hole.
Apply this force to the respective keypoint ( Solution > Define Loads > Apply > Structural >
Force/Moment > On Keypoints Select a force in the y direction of -1000)
The applied loads and constraints should now appear as shown below.
RESULT:
Ex. No. : 03. STRESS ANALYSIS OF RECTANGULAR L - BRACKET
Date :
AIM:
To determine the displacement and bending stress of a given L - bracket using Finite Element
Analysis bases ANSYS structure and view the displacement and bending stress plots.
PROCEDURE:
This is a simple, single load step, structural static
analysis of the (corner angle) L-bracket shown
below. The upper left-hand pin hole is constrained
(welded) around its entire circumference, and a
tapered pressure load is applied to the bottom of
the lower right-hand pin hole. The objective of the
problem is to demonstrate the typical ANSYS
analysis procedure.
Given
The bracket is made of A36 steel with a Young’s modulus of 30E6 psi and Poisson’s ratio of 0.27.
Approach and Assumptions
Assume plane stress for this analysis. Since the bracket is thin in the z direction (1/2 inch thickness) compared
to its x and y dimensions, and since the pressure load acts only in the x-y plane, this is a valid assumption.
Your approach is to use solid modeling to generate the 2-D model and automatically mesh it with nodes and
elements.
Build Geometry
1. Define rectangles.
2. Change plot controls and replot.
3. Change working plane to polar and create first circle.
4. Move working plane and create second circle.
5. Add areas.
6. Create line fillet.
7. Create fillet area.
8. Add areas together.
9. Create first pin hole.
10. Move working plane and create second pin hole.
11. Subtract pin holes from bracket.
12. Save the database as model.db.
Define Materials
13. Set Preferences.
14. Define Material Properties.
15. Define element types and options.
16. Define real constants.
Generate Mesh
17. Mesh the area.
Apply Loads
19. Apply displacement constraints.
20. Apply pressure load.
Obtain Solution
21. Solve.
Review Results
22. Enter the general postprocessor and read in the results.
23. Plot the deformed shape.
24. Plot the von Mises equivalent stress.
25. List the reaction solution.
Step 1: Define rectangles.
There are several ways to create the model geometry within ANSYS, some more convenient than others. The
first step is to recognize that you can construct the bracket easily with combinations of rectangles and circle
Primitives.
Decide where the origin will be located and then define the rectangle and circle primitives relative to that origin.
The location of the origin is arbitrary. Here, use the center of the upper left-hand hole. ANSYS does not need to know
where the origin is. Simply begin by defining a rectangle relative to that location. In ANSYS, this origin is called
the global origin.
1. Main Menu> Preprocessor> Modeling> Create> Areas> Rectangle> By Dimensions
2. Enter the following:
X1 = 0 (Note: Press the Tab key between entries)
X2 = 6
Y1 = -1
Y2 = 1
Before going to the next step, save the work you have done so far. ANSYS stores any input data in
memory to the ANSYS database. To save that database to a file, use the SAVE operation, available as a tool on
the Toolbar. ANSYS names the database file using the format jobname.db. If you started ANSYS using the
product launcher, you can specify a jobname at that point (the default jobname is file). You can check the
current jobname at any time by choosing Utility Menu> List> Status> Global Status. You can also save the
database at specific milestone points in the analysis (such as after the model is complete, or after the model is
meshed) by choosing Utility Menu> File> Save As and specifying different jobnames (model.db, or
mesh.db, etc.).
It is important to do an occasional save so that if you make a mistake, you can restore the model from the last
saved state. You restore the model using the RESUME operation, also available on the Toolbar. (You can also
find SAVE and RESUME on the Utility Menu, under File.)
4. Toolbar: SAVE_DB.
Step 3: Change working plane to polar and create first circle.
The next step in the model construction is to create the half circle at each end of the bracket. You will actually create
a full circle on each end and then combine the circles and rectangles with a Boolean "add" operation (discussed in
step 5.). To create the circles, you will use and display the working plane. You could have shown the working plane as
you created the rectangles but it was not necessary.
Before you begin however, first "zoom out" within the Graphics Window so you can see more of the circles as you
create them. You do this using the "Pan-Zoom-Rotate" dialog box, a convenient graphics control box you’ll use
often in any ANSYS session.
1. Utility Menu> PlotCtrls> Pan, Zoom, Rotate
6. Click on Polar.
Note
While you are positioning the cursor for picking, the "dynamic" WP X and Y values are displayed in the Solid
Circular Area dialog box. Also, as an alternative to picking, you can type these values along with the radius into
the dialog box.
Step 4: Move working plane and create second circle.
To create the circle at the other end of the bracket in the same manner, you need to first move the working
plane to the origin of the circle. The simplest way to do this without entering number offsets is to move the WP
to an average keypoint location by picking the keypoints at the bottom corners of the lower, right rectangle.
1. Utility Menu> WorkPlane> Offset WP to> Keypoints
2. Pick keypoint at lower left corner of rectangle.
3. Pick keypoint at lower right of rectangle.
4. OK to close picking menu.
3. Toolbar: SAVE_DB.
Step 9: Create first pin hole.
1. Utility Menu> WorkPlane> Display Working Plane (toggle on)
6. Material> Exit
Step 15: Define element types and options.
In any analysis, you need to select from a library of element types and define the appropriate ones for your
analysis.
For this analysis, you will use only one element type, PLANE82, which is a 2-D, quadratic, structural, higherorder
element. The choice of a higher-order element here allows you to have a coarser mesh than with lowerorder
elements while still maintaining solution accuracy. Also, ANSYS will generate some triangle shaped elements in
the mesh that would otherwise be inaccurate if you used lower-order elements (PLANE42). You will need to
specify plane stress with thickness as an option for PLANE82. (You will define the thickness as a real constant in the
next step.)
1. Main Menu> Preprocessor> Element Type> Add/Edit/Delete
3. OK for PLANE82.
Before clicking on the Help button in the next step, you should be aware that the help
information may appear in the same window as this tutorial, replacing the contents of the
tutorial. After reading the help information, click on the Back button to return to this
tutorial. If the help information appears in a separate window from the tutorial, minimize or
close the help window after you read the help information.
Generate Mesh
Step 17: Mesh the area.
One nice feature of the ANSYS program is that you can automatically mesh the model without specifying any
mesh size controls. This is using what is called a default mesh. If you’re not sure how to determine the mesh
density, let ANSYS try it first! Meshing this model with a default mesh however, generates more elements
than are allowed in the ANSYS ED program. Instead you will specify a global element size to control overall
mesh density.
1. Main Menu> Preprocessor> Meshing> Mesh Tool
3. Type in 0.5.
4. OK.
6. Click on Mesh.
7. Pick All for the area to be meshed (in picking menu). Close any warning messages that
appear.
Note
The mesh you see on your screen may vary slightly from the mesh shown here. As a result of this, you may see slightly
different results during postprocessing. For a discussion of results accuracy, see Planning Your
Approach in the ANSYS Modeling and Meshing Guide.
Step 18: Save the database as mesh.db.
Here again, you will save the database to a named file, this time mesh.db.
8. Toolbar: SAVE_DB.
Step 20: Apply pressure load.
Now apply the tapered pressure load to the bottom, right-hand pin hole. ("Tapered" here means varying
linearly.) Note that when a circle is created in ANSYS, four lines define the perimeter. Therefore, apply the
pressure to two lines making up the lower half of the circle. Since the pressure tapers from a maximum value (500
psi) at the bottom of the circle to a minimum value (50 psi) at the sides, apply pressure in two separate steps, with
reverse tapering values for each line.
The ANSYS convention for pressure loading is that a positive load value represents pressure into the surface
(compressive).
1. Main Menu> Solution> Define Loads> Apply> Structural> Pressure> On Lines
2. Pick line defining bottom left p art of the circle (line 6).
3. Apply.
6. Apply.
7. Pick line defining bottom right part of circle (line 7).
8. Apply.
Obtain Solution
Step 21: Solve.
3. OK.
6. OK.
7. Make choices in the Animation Controller (not shown), if necessary, then choose Close.
Step 24: Plot the von Mises equivalent stress.
1. Main Menu> General Postproc> Plot Results> Contour Plot> Nodal Solu
8. OK.
9. Make choices in the Animation Controller (not shown), if necessary, then choose Close.
Step 25: List reaction solution.
1. Main Menu> General Postproc> List Results> Reaction Solu
RESULT:
Ex. No. : 04 STRESS ANALYSIS OF AN AXISYMMETRIC COMPONENT
Date :
AIM:
To determine the displacement and bending stress of a given axisymmetric component using Finite
Element Analysis bases ANSYS structure and view the displacement and bending stress plots.
PROCEDURE:
The model will be that of a closed tube
made from steel. Point loads will be applied at the
center of the top and bottom plate to make an
analytical verification simple to calculate. A 3/4
cross section view of the tube is shown below.
As a warning, point loads will create
discontinuities in the your model near the point of
application. If you chose to use these types of
loads in your own modelling, be very careful and
be sure to understand the theory of how the FEA
package is appling the load and the assumption it
is making. In this case, we will only be concerned
about the stress distribution far from the point of
application, so the discontinuities will have a
negligable effect.
3. Create Areas
Preprocessor > Modeling > Create > Areas > Rectangle > By Dimensions
For an axisymmetric problem, ANSYS will rotate the area around the y-axis at x=0.
Therefore, to create the geometry mentioned above, we must define a U-shape.
We are going to define 3 overlapping rectangles as defined in the following table:
Rectangle X1 X2 Y1 Y2
1 0 20 0 5
2 15 20 0 100
3 0 20 95 100
4. Add Areas Together
Preprocessor > Modeling > Operate > Booleans > Add > Areas
6. Turn on Axisymmetry
While the Element Types window is still open, click the Options... button.
Under Element behavior K3 select Axisymmetric.
2. Apply Constraints
o Solution > Define Loads > Apply > Structural > Displacement > Symmetry B.C. > On Lines
Pick the two edges on the left, at x=0, as shown below. By using the symmetry B.C.
command, ANSYS automatically calculates which DOF's should be constrained for the line
of symmetry. Since the element we are using only has 2 DOF's per node, we could have
constrained the lines in the x-direction to create the symmetric boundary conditions.
o Utility Menu > Select > Entities
Select Nodes and By Location from the scroll down menus. Click Y coordinates and type 50 into
the input box as shown below, then click OK.
Solution > Define Loads > Apply > Structural > Displacement > On Nodes > Pick All
Constrain the nodes in the y-direction (UY). This is required to constrain the model in space,
otherwise it would be free to float up or down. The location to constrain the model in the y-
direction (y=50) was chosen because it is along a symmetry plane. Therefore, these nodes
won't move in the y-direction according to theory.
3. Utility Menu > Select > Entities
In the select entities window, click Sele All to reselect all nodes. It is important to always reselect all
entities once you've finished to ensure future commands are applied to the whole model and not just a few
entities. Once you've clicked Sele All, click on Cancel to close the window.
4. Apply Loads
o Solution > Define Loads > Apply > Structural > Force/Moment > On Keypoints
Pick the top left corner of the area and click OK. Apply a load of 100 in the FY direction.
o Solution > Define Loads > Apply > Structural > Force/Moment > On Keypoints
Pick the bottom left corner of the area and click OK. Apply a load of -100 in the FY direction.
o The applied loads and constraints should now appear as shown in the figure below.
o If you take the average of the stress in the y-direction over the thickness of the
tube, (0.18552 + 0.17866)/2, the stress in the tube is 0.182 MPa, matching the analytical solution.
The average is used because in the analytical case, it is assumed the stress is evenly distributed
across the thickness. This is only true when the location is far from any stress concentrators, such
as corners. Thus, to approximate the analytical solution, we must average the stress over the
thickness.
3. Plotting the Elements as Axisymmetric
Utility Menu > PlotCtrls > Style > Symmetry Expansion > 2-D Axi-symmetric
The following window will appear. By clicking on 3/4 expansion you can produce the figure
shown at the beginning of this tutorial.
RESULT:
Ex. No. : 05 STRESS ANALYSIS OF CANTILEVER BEAM
Date:
AIM:
To determine the displacement and bending stress of a given Cantilever Beam using Finite Element
Analysis bases ANSYS structure and view the displacement and bending stress plots.
PROCEDURE:
We are going to define 2 keypoints for the simplified structure as given in the following table
coordinate
keypoint
x y z
1 0 0 0
2 500 0 0
From the 'ANSYS Main Menu' select:
Preprocessor > Modeling > Create > Keypoints > In Active CS
3. Form Lines
The two keypoints must now be connected to form a bar using a straight line.
Select: Preprocessor > Modeling> Create > Lines > Lines > Straight Line.
Pick keypoint #1 (i.e. click on it). It will now be marked by a small yellow box. Now
pick keypoint #2. A permanent line will appear.
When you're done, click on 'OK' in the 'Create Straight Line' window.
4. Define the Type of Element
It is now necessary to create elements on this line.
From the Preprocessor Menu, select: Element Type > Add/Edit/Delete.
Click on the 'Add...' button. The following window will appear:
For this example, we will use the 3D elastic straight pipe element as selected in the above
figure. Select the element shown and click 'OK'. You should see 'Type 1 PIPE16' in the
'Element Types' window.
Click on the 'Options...' button in the 'Element Types' dialog box. The following window will
appear:
Click and hold the K6 button (second from the bottom), and select 'Include Output' and click
'OK'. This gives us extra force and moment output.
Click on 'Close' in the 'Element Types' dialog box and close the 'Element Type' menu.
The applied loads and constraints should now appear as shown below.
The maximum deflection occurs at the end of the rod and was found to be 6.2mm as shown above.
Stress:
The maximum stress occurs at the base of the rod and was found to be 64.9MPa as shown above (pure
bending stress).
Observe the value of the maximum deflection in the upper left hand corner (shown here
surrounded by a blue border for emphasis). This is identical to that obtained via hand
calculations.
Deflection :
For a more detailed version of the deflection of the beam,
From the 'General Postproc' menu select Plot results > Contour Plot > Nodal Solution.
Select 'DOF solution' and 'USUM'. Leave the other selections as the default
values. Click 'OK'.
o You may want to have a more useful scale, which can be accomplished by going to the Utility
Menu and selecting Plot Controls > Style > Contours > Uniform Contours
o The deflection can also be obtained as a list as shown below. General Postproc > List
Results > Nodal Solution ... select 'DOF Solution' and 'ALL DOFs' from the lists in the 'List
Nodal Solution' window and click 'OK'. This means that we want to see a listing of all
translational and rotational degrees of freedom from the solution. If we had only wanted to
see the displacements for example, we would have chosen 'ALL Us' instead of 'ALL DOFs'.
Are these results what you expected? Again, the maximum deflection occurs at node 2, the
right end of the rod. Also note that all the rotational and translational degrees of freedom were
constrained to zero at node 1.
If you wanted to save these results to a file, use the mouse to go to the 'File' menu (at the
upper left-hand corner of this list window) and select 'Save as'.
Stresses
For line elements (ie beams, spars, and pipes) you will need to use the Element Table to gain access to
derived data (ie stresses, strains).
From the General Postprocessor menu select Element Table > Define Table...
Click on 'Add...'
As shown above, in the 'Item,Comp' boxes in the above window, select 'Stress' and 'von
Mises SEQV'
Click on 'OK' and close the 'Element Table Data' window.
Plot the Stresses by selecting Plot Elem Table in the Element Table Menu
The following window will appear. Ensure that 'SEQV' is selected and click 'OK'
o If you changed the contour intervals for the Displacement plot to "User Specified" you may
need to switch this back to "Auto calculated" to obtain new values for VMIN/VMAX.
Utility Menu > PlotCtrls > Style > Contours > Uniform Contours...
o In the window,
A. Enter IMoment as the 'User label for item' - this will give a name to the data
B. Select 'By sequence num' in the Item box
C. Select 'SMISC' in the first Comp box
D. Enter SMISC,6 in the second Comp box
E. Click 'OK'
This will save all of the bending moment data at the left hand side (I side) of each element.
Now we need to find the bending moment data at the right hand side (J side) of each element.
Again, click 'Add...' in the 'Element Table Data' window.
Enter JMoment as the 'User label for item' - again, this will give a name to the data
A. Same as above
B. Same as above
C. For step D, enter SMISC,12 in the second Comp box
D. Click 'OK'
Click 'Close' in the 'Element Table Data' window and close the 'Element Table' Menu. Select Plot
Results > Contour Plot > Line Elem Res...
From the 'Plot Line-Element Results' window, select 'IMOMENT' from the pull down menu for LabI, and
'JMOMENT' from the pull down menu for LabJ. Click 'OK'. Note again that you can modify the intervals for the
contour plot.
Now, you can double check these solutions analytically. Note that the line between the I and J
point is a linear interpolation.
Before the explanation of the above steps, enter help pipe16 in the command line as shown below and then
hit enter.
Briefly read the ANSYS documentation which appears, pay particular attention to the Tables near the end
of the document (shown below).
Table 1. PIPE16 Item, Sequence Numbers, and Definitions for the ETABLE Commands
node I
Note that SMISC 6 (which we used to obtain the values at node I) correspond to MMOMZ - the Member
moment for node I. The value of 'e' varies with different Element Types, therefore you must check the ANSYS
Documentation files for each element to determine the appropriate SMISC corresponding to the plot you wish
to generate.
RESULT:
Ex. No. : 06 STRESS ANALYSIS OF SIMPLY SUPPORTED BEAM
Date:
AIM:
To determine the displacement and bending stress of a given simply supported beam using Finite
Element Analysis bases ANSYS structure and view the displacement and bending stress plots.
PROCEDURE:
Coordinates
Keypoint (x,y)
1 (0,0)
2 (1000,0)
4. Define Lines
Preprocessor > Modeling > Create > Lines > Lines > Straight Line
L,K#,K#
Create a line between Keypoint 1 and Keypoint 2.
5. Define Element Types
Preprocessor > Element Type > Add/Edit/Delete...
For this problem we will use the BEAM3 element. This element has 3 degrees of freedom
(translation along the X and Y axis's, and rotation about the Z axis).
6. Define Real Constants
Preprocessor > Real Constants... > Add...
In the 'Real Constants for BEAM3' window, enter the following geometric properties:
i. Cross-sectional area AREA: 100
ii. Area Moment of Inertia IZZ: 833.333
iii. Total beam height HEIGHT: 10
This defines an element with a solid rectangular cross section 10mm x 10mm.
7. Define Element Material Properties
Preprocessor > Material Props > Material Models > Structural > Linear > Elastic > Isotropic In
the window that appears, enter the following geometric properties for steel:
i. Young's modulus EX: 200000
ii. Poisson's Ratio PRXY: 0.3
8. Define Mesh Size
Preprocessor > Meshing > Size Cntrls > ManualSize > Lines > All Lines... For
this example we will use an element length of 100mm.
9. Mesh the frame
Preprocessor > Meshing > Mesh > Lines > click 'Pick All'
10. Plot Elements
Utility Menu > Plot > Elements
You may also wish to turn on element numbering and turn off keypoint numbering
Utility Menu > PlotCtrls > Numbering...
2. Apply Constraints
Solution > Define Loads > Apply > Structural > Displacement > On Keypoints
Pin Keypoint 1 (ie UX and UY constrained) and fix Keypoint 2 in the y direction (UY
constrained).
3. Apply Loads
We will apply a distributed load, of 1000 N/m or 1 N/mm, over the entire length of the beam.
o Select Solution > Define Loads > Apply > Structural > Pressure > On Beams
o Click 'Pick All' in the 'Apply F/M' window.
o As shown in the following figure, enter a value of 1 in the field 'VALI Pressure value at node
I' then click 'OK'.
The applied loads and constraints should now appear as shown in the
figure below.
Note: To have the constraints and loads appear each time you select
'Replot' you must change some settings. Select Utility Menu >
PlotCtrls > Symbols.... In the window that appears, select 'Pressures'
in the pull down menu of the 'Surface Load Symbols' section.
4. Solve the System
Solution > Solve > Current LS
As shown in the plot below, the maximum stress occurs in the middle of the beam with a value of 750
MPa.
RESULT:
Ex. No. : 07 STRESS ANALYSIS OF FIXED BEAM
Date:
AIM:
To determine the displacement and bending stress of a given FIXED BEAM using Finite Element
Analysis bases ANSYS structure and view the displacement and bending stress plots.
PROCEDURE:
coordinate
keypoint
x y z
1 0 0 0
2 500 0 0
From the 'ANSYS Main Menu' select:
Preprocessor > Modeling > Create > Keypoints > In Active CS
12. Form Lines
The two keypoints must now be connected to form a bar using a straight line.
Select: Preprocessor > Modeling> Create > Lines > Lines > Straight Line.
Pick keypoint #1 (i.e. click on it). It will now be marked by a small yellow box. Now
pick keypoint #2. A permanent line will appear.
When you're done, click on 'OK' in the 'Create Straight Line' window.
13. Define the Type of Element
It is now necessary to create elements on this line.
From the Preprocessor Menu, select: Element Type > Add/Edit/Delete.
Click on the 'Add...' button. The following window will appear:
For this example, we will use the 3D elastic straight pipe element as selected in the above
figure. Select the element shown and click 'OK'. You should see 'Type 1 PIPE16' in the
'Element Types' window.
Click on the 'Options...' button in the 'Element Types' dialog box. The following window will
appear:
Click and hold the K6 button (second from the bottom), and select 'Include Output' and click
'OK'. This gives us extra force and moment output.
Click on 'Close' in the 'Element Types' dialog box and close the 'Element Type' menu.
The applied loads and constraints should now appear as shown below.
The maximum deflection occurs at the end of the rod and was found to be 6.2mm as shown above.
Stress:
The maximum stress occurs at the base of the rod and was found to be 64.9MPa as shown above (pure
bending stress).
Observe the value of the maximum deflection in the upper left hand corner (shown here
surrounded by a blue border for emphasis). This is identical to that obtained via hand
calculations.
Deflection :
For a more detailed version of the deflection of the beam,
From the 'General Postproc' menu select Plot results > Contour Plot > Nodal Solution.
Select 'DOF solution' and 'USUM'. Leave the other selections as the default values. Click
'OK'.
o You may want to have a more useful scale, which can be accomplished by going to the Utility
Menu and selecting Plot Controls > Style > Contours > Uniform Contours
o The deflection can also be obtained as a list as shown below. General Postproc > List
Results > Nodal Solution ... select 'DOF Solution' and 'ALL DOFs' from the lists in the 'List
Nodal Solution' window and click 'OK'. This means that we want to see a listing of all
translational and rotational degrees of freedom from the solution. If we had only wanted to
see the displacements for example, we would have chosen 'ALL Us' instead of 'ALL DOFs'.
Are these results what you expected? Again, the maximum deflection occurs at node 2, the
right end of the rod. Also note that all the rotational and translational degrees of freedom were
constrained to zero at node 1.
If you wanted to save these results to a file, use the mouse to go to the 'File' menu (at the
upper left-hand corner of this list window) and select 'Save as'.
Stresses
For line elements (ie beams, spars, and pipes) you will need to use the Element Table to gain access to
derived data (ie stresses, strains).
From the General Postprocessor menu select Element Table > Define Table...
Click on 'Add...'
As shown above, in the 'Item,Comp' boxes in the above window, select 'Stress' and 'von
Mises SEQV'
Click on 'OK' and close the 'Element Table Data' window.
Plot the Stresses by selecting Plot Elem Table in the Element Table Menu
The following window will appear. Ensure that 'SEQV' is selected and click 'OK'
o If you changed the contour intervals for the Displacement plot to "User Specified" you may
need to switch this back to "Auto calculated" to obtain new values for VMIN/VMAX.
Utility Menu > PlotCtrls > Style > Contours > Uniform Contours ...
o In the window,
A. Enter IMoment as
the 'User label for item' - this will give a name to the data
B. Select 'By sequence num' in the Item box
C. Select 'SMISC' in the first Comp box
D. Enter SMISC,6 in the second Comp box
E. Click 'OK'
This will save all of the bending moment data at the left hand side (I side) of each element. Now we need to find
the bending moment data at the right hand side (J side) of each element.
Again, click 'Add...' in the 'Element Table Data' window.
Enter JMoment as the 'User label for item' - again, this will give a name to the data
A. Same as above
B. Same as above
C. For step D, enter SMISC,12 in the second Comp box
D. Click 'OK'
Click 'Close' in the 'Element Table Data' window and close the 'Element Table' Menu. Select Plot Results > Contour Plot
> Line Elem Res...
From the 'Plot Line-Element Results' window, select 'IMOMENT' from the
pull down menu for LabI, and 'JMOMENT' from the pull down menu for LabJ. Click 'OK'. Note again that you can modify the
intervals for the contour plot.
Now, you can double check these solutions analytically. Note that the line between the I and J point is a linear
interpolation.
Before the explanation of the above steps, enter help pipe16 in the command line as shown below and then hit enter.
Briefly read the ANSYS documentation which appears, pay particular attention to the Tables near the end of the
document (shown below).
Table 1. PIPE16 Item, Sequence Numbers, and Definitions for the ETABLE Commands
node I
Note that SMISC 6 (which we used to obtain the values at node I) correspond to MMOMZ - the Member
moment for node I. The value of 'e' varies with different Element Types, therefore you must check the ANSYS
Documentation files for each element to determine the appropriate SMISC corresponding to the plot you wish
to generate.
RESULT:
Ex. No.:08. MODE FREQUENCY ANALYSIS OF A 2D COMPONENT
Date:
AIM:
To determine the frequency and three mode shapes of a given 2D Component using Finite Element
Analysis based ANSYS software and view the three mode shapes.
PROCEDURE:
The purpose of this tutorial is to outline the steps
required to do a simple modal analysis of 2D
component.
3. Create Keypoints
We are going to define 2 keypoints (the beam vertices) for this structure as given in the
following table:
4. Define Lines
Preprocessor > Modeling > Create > Lines > Lines > Straight Line
For this problem we will use the BEAM3 (Beam 2D elastic) element. This element has 3
degrees of freedom (translation along the X and Y axis's, and rotation about the Z axis). With only
3 degrees of freedom, the BEAM3 element can only be used in 2D analysis.
In the 'Real Constants for BEAM3' window, enter the following geometric properties:
This defines an element with a solid rectangular cross section 0.01 m x 0.01 m.
Preprocessor > Material Props > Material Models > Structural > Linear > Elastic > Isotropic
In the window that appears, enter the following geometric properties for steel:
To enter the density of the material, double click on 'Linear' followed by 'Density' in the
'Define Material Model Behavior' Window
Note: For dynamic analysis, both the stiffness and the material density have to be specified.
Preprocessor > Meshing > Size Cntrls > ManualSize > Lines > All Lines...
For this example we will specify 10 element divisions along the line.
Preprocessor > Meshing > Mesh > Lines > click 'Pick All'
o As shown, select the Subspace method and enter 5 in the 'No. of modes to extract'
o Check the box beside 'Expand mode shapes' and enter 5 in the 'No. of modes to expand'
o Click 'OK'
Note that the default mode extraction method chosen is the Reduced Method. This is the
fastest method as it reduces the system matrices to only consider the Master Degrees of
Freedom (see below). The Subspace Method extracts modes for all DOF's. It is therefore
more exact but, it also takes longer to compute (especially when the complex geometries).
o The following window will then appear
The following table compares the mode frequencies in Hz predicted by theory and ANSYS.
Percent
Mode Theory ANSYS
Error
Note: To obtain accurate higher mode frequencies, this mesh would have to be refined even more
(i.e. instead of 10 elements, we would have to model the cantilever using 15 or more elements
depending upon the highest mode frequency of interest).
2. View Mode Shapes
o Select: General Postproc > Read Results > First Set
This selects the results for the first mode shape
o Select General Postproc > Plot Results > Deformed shape . Select 'Def + undef edge'
The first mode shape will now appear in the graphics window.
o To view the next mode shape, select General Postproc > Read Results > Next Set . As
above choose General Postproc > Plot Results > Deformed shape . Select 'Def + undef
edge'.
o The first four mode shapes should look like the following:
3. Animate Mode Shapes
o Select Utility Menu (Menu at the top) > Plot Ctrls > Animate > Mode Shape
The following window will appear
Mode 2
Mode 3
Mode 4
RESULT:
Ex.No.: 09. MODE FREQUENCY ANALYSIS OF CANTILEVER BEAM
Date :
AIM:
To determine the frequency and three mode shapes of a given CANTILEVER BEAM using Finite
Element Analysis based ANSYS software and view the three mode shapes.
PROCEDURE:
3. Create Keypoints
We are going to define 2 keypoints (the beam vertices) for this structure as given in the
following table:
4. Define Lines
Preprocessor > Modeling > Create > Lines > Lines > Straight Line
For this problem we will use the BEAM3 (Beam 2D elastic) element. This element has 3
degrees of freedom (translation along the X and Y axis's, and rotation about the Z axis). With only
3 degrees of freedom, the BEAM3 element can only be used in 2D analysis.
In the 'Real Constants for BEAM3' window, enter the following geometric properties:
This defines an element with a solid rectangular cross section 0.01 m x 0.01 m.
Preprocessor > Material Props > Material Models > Structural > Linear > Elastic > Isotropic
In the window that appears, enter the following geometric properties for steel:
To enter the density of the material, double click on 'Linear' followed by 'Density' in the
'Define Material Model Behavior' Window
Note: For dynamic analysis, both the stiffness and the material density have to be specified.
Preprocessor > Meshing > Size Cntrls > ManualSize > Lines > All Lines...
For this example we will specify 10 element divisions along the line.
Preprocessor > Meshing > Mesh > Lines > click 'Pick All'
o As shown, select the Subspace method and enter 5 in the 'No. of modes to extract'
o Check the box beside 'Expand mode shapes' and enter 5 in the 'No. of modes to expand'
o Click 'OK'
Note that the default mode extraction method chosen is the Reduced Method. This is the
fastest method as it reduces the system matrices to only consider the Master Degrees of
Freedom (see below). The Subspace Method extracts modes for all DOF's. It is therefore
more exact but, it also takes longer to compute (especially when the complex geometries).
o The following window will then appear
Note: To obtain accurate higher mode frequencies, this mesh would have to be refined even more
(i.e. instead of 10 elements, we would have to model the cantilever using 15 or more elements
depending upon the highest mode frequency of interest).
5. View Mode Shapes
o Select: General Postproc > Read Results > First Set
This selects the results for the first mode shape
o Select General Postproc > Plot Results > Deformed shape . Select 'Def + undef edge'
The first mode shape will now appear in the graphics window.
o To view the next mode shape, select General Postproc > Read Results > Next Set . As
above choose General Postproc > Plot Results > Deformed shape . Select 'Def + undef
edge'.
o The first four mode shapes should look like the following:
RESULT:
Ex. No. : 10 MODE FREQUENCY ANALYSIS OF SIMPLY SUPPORTED BEAM
Date:
AIM:
To determine the frequency and three mode shapes of a given SSB using Finite Element Analysis based
ANSYS software and view the three mode shapes.
PROCEDURE:
3. Create Keypoints
We are going to define 2 keypoints (the beam vertices) for this structure as given in the
following table:
4. Define Lines
Preprocessor > Modeling > Create > Lines > Lines > Straight Line
For this problem we will use the BEAM3 (Beam 2D elastic) element. This element has 3
degrees of freedom (translation along the X and Y axis's, and rotation about the Z axis). With only
3 degrees of freedom, the BEAM3 element can only be used in 2D analysis.
In the 'Real Constants for BEAM3' window, enter the following geometric properties:
This defines an element with a solid rectangular cross section 0.01 m x 0.01 m.
Preprocessor > Material Props > Material Models > Structural > Linear > Elastic > Isotropic
In the window that appears, enter the following geometric properties for steel:
To enter the density of the material, double click on 'Linear' followed by 'Density' in the
'Define Material Model Behavior' Window
Note: For dynamic analysis, both the stiffness and the material density have to be specified.
Preprocessor > Meshing > Size Cntrls > ManualSize > Lines > All Lines...
For this example we will specify 10 element divisions along the line.
Preprocessor > Meshing > Mesh > Lines > click 'Pick All'
o As shown, select the Subspace method and enter 5 in the 'No. of modes to extract'
o Check the box beside 'Expand mode shapes' and enter 5 in the 'No. of modes to expand'
o Click 'OK'
Note that the default mode extraction method chosen is the Reduced Method. This is the
fastest method as it reduces the system matrices to only consider the Master Degrees of
Freedom (see below). The Subspace Method extracts modes for all DOF's. It is therefore
more exact but, it also takes longer to compute (especially when the complex geometries).
o The following window will then appear
The following table compares the mode frequencies in Hz predicted by theory and ANSYS.
Percent
Mode Theory ANSYS
Error
o The first four mode shapes should look like the following:
9. Animate Mode Shapes
o Select Utility Menu (Menu at the top) > Plot Ctrls > Animate > Mode Shape
The following window will appear
Mode 3 and Mode 4
RESULT:
Ex. No. : 11. MODE FREQUENCY ANALYSIS OF FIXED BEAM
Date :
AIM:
To determine the frequency and three mode shapes of a given FIXED BEAM using Finite Element
Analysis based ANSYS software and view the three mode shapes.
3. Create Keypoints
We are going to define 2 keypoints (the beam vertices) for this structure as given in the
following table:
4. Define Lines
Preprocessor > Modeling > Create > Lines > Lines > Straight Line
For this problem we will use the BEAM3 (Beam 2D elastic) element. This element has 3
degrees of freedom (translation along the X and Y axis's, and rotation about the Z axis). With only
3 degrees of freedom, the BEAM3 element can only be used in 2D analysis.
In the 'Real Constants for BEAM3' window, enter the following geometric properties:
This defines an element with a solid rectangular cross section 0.01 m x 0.01 m.
Preprocessor > Material Props > Material Models > Structural > Linear > Elastic > Isotropic
In the window that appears, enter the following geometric properties for steel:
To enter the density of the material, double click on 'Linear' followed by 'Density' in the
'Define Material Model Behavior' Window
Note: For dynamic analysis, both the stiffness and the material density have to be specified.
Preprocessor > Meshing > Size Cntrls > ManualSize > Lines > All Lines...
For this example we will specify 10 element divisions along the line.
Preprocessor > Meshing > Mesh > Lines > click 'Pick All'
o As shown, select the Subspace method and enter 5 in the 'No. of modes to extract'
o Check the box beside 'Expand mode shapes' and enter 5 in the 'No. of modes to expand'
o Click 'OK'
Note that the default mode extraction method chosen is the Reduced Method. This is the
fastest method as it reduces the system matrices to only consider the Master Degrees of
Freedom (see below). The Subspace Method extracts modes for all DOF's. It is therefore
more exact but, it also takes longer to compute (especially when the complex geometries).
o The following window will then appear
The following table compares the mode frequencies in Hz predicted by theory and ANSYS.
Percent
Mode Theory ANSYS
Error
Mode 3 and Mode 4
RESULT:
Ex. No. : 12. HARMONIC ANALYSIS OF A 2D COMPONENT
Date:
AIM:
To determine the amplitude and frequency of a given 2D Component using Finite Element Analysis based
ANSYS software and view the three mode shapes.
PROCEDURE:
3. Create Keypoints
We are going to define 2 keypoints (the beam vertices) for this structure as given in the
following table:
Preprocessor > Modeling > Create > Lines > Lines > Straight Line
For this problem we will use the BEAM3 (Beam 2D elastic) element. This element has 3
degrees of freedom (translation along the X and Y axis's, and rotation about the Z axis). With only
3 degrees of freedom, the BEAM3 element can only be used in 2D analysis.
In the 'Real Constants for BEAM3' window, enter the following geometric properties:
This defines an element with a solid rectangular cross section 0.01 m x 0.01 m.
Preprocessor > Material Props > Material Models > Structural > Linear > Elastic > Isotropic
In the window that appears, enter the following geometric properties for steel:
To enter the density of the material, double click on 'Linear' followed by 'Density' in the
'Define Material Model Behavior' Window
Note: For dynamic analysis, both the stiffness and the material density have to be specified.
Preprocessor > Meshing > Size Cntrls > ManualSize > Lines > All Lines...
For this example we will specify 10 element divisions along the line.
o As shown, select the Full Solution method, the Real + imaginary DOF printout format and
do not use lumped mass approx.
o Click 'OK'
The following window will appear. Use the default settings (shown below).
3. Apply Constraints
o Select Solution > Define Loads > Apply > Structural > Displacement > On Nodes
The following window will appear once you select the node at x=0 (Note small changes in the
window compared to the static examples):
By doing this we will be subjecting the beam to loads at 1 Hz, 2 Hz, 3 Hz,........ 100 Hz. We
will specify a stepped boundary condition (KBC) as this will ensure that the same amplitude (100
N) will be applyed for each of the frequencies. The ramped option, on the other hand, would
ramp up the amplitude where at 1 Hz the amplitude would be 1 N and at 100 Hz the amplitude
would be 100 N.
You should now have the following in the ANSYS Graphics window
o Select Add (the green '+' sign in the upper left corner) from this window and the following
window should appear
o We are interested in the Nodal Solution > DOF Solution > Y-Component of displacement.
Click OK.
o Graphically select node 2 when prompted and click OK. The 'Time History Variables'
window should now look as follows
3. List Stored Variables
o In the 'Time History Variables' window click the 'List' button, 3 buttons to the left of 'Add'
The following window will appear listing the data:
Note that we get peaks at frequencies of approximately 8.3 and 51 Hz. This corresponds with the
predicted frequencies of 8.311 and 51.94Hz.
To get a better view of the response, view the log scale of UY.
o Select Utility Menu > PlotCtrls > Style > Graphs > Modify Axis
The following window will appear
o As marked by an 'A' in the above window, change the Y-axis scale to 'Logarithmic'
o Select Utility Menu > Plot > Replot
o You should now see the following
This is the response at node 2 for the cyclic load applied at this node from 0 - 100 Hz.
o For ANSYS version lower than 7.0, the 'Variable Viewer' window is not available. Use the
'Define Variables' and 'Store Data' functions under TimeHist Postpro. See the help file for
instructions.
RESULT:
Ex. No. : 13. THERMAL STRESS ANALYSIS OF A 2D COMPONENT - STATIC
Date :
AIM:
To determine the temperature distribution axial stress and tangential stress of a given rectangular
block using FEA based ANSYS software and view the temperature distribution axial stress and tangential
stress.
PROCEDURE:
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 attemp to expand. However, since it is pinned this
cannot occur and as such, stress is created in the link. A steady-state solution of the resulting stress will be found to
simplify the analysis.
Loads will not be applied to the link, only a temperature change of 75 degrees Celsius. The link is steel with a
modulus of elasticity of 200 GPa, a thermal conductivity of 60.5 W/m*K and a thermal expansion coefficient
of 12e-6 /K.
3. Define Keypoints
Preprocessor > Modeling > Create > Keypoints > In Active CS...
K,#,x,y,z
We are going to define 2 keypoints for this link as given in the following table:
Keypoint Coordinates (x,y,z)
1 (0,0)
2 (1,0)
4. Create Lines
Preprocessor > Modeling > Create > Lines > Lines > In Active Coord
Create a line joining Keypoints 1 and 2, representing a link 1 meter long.
5. Define the Type of Element
Preprocessor > Element Type > Add/Edit/Delete...
For this problem we will use the LINK33 (Thermal Mass Link 3D conduction) element. This
element is a uniaxial element with the ability to conduct heat between its nodes.
6. Define Real Constants
Preprocessor > Real Constants... > Add...
In the 'Real Constants for LINK33' window, enter the following geometric properties:
i. Cross-sectional area AREA: 4e-4
This defines a beam with a cross-sectional area of 2 cm X 2 cm.
7. Define Element Material Properties
Preprocessor > Material Props > Material Models > Thermal > Conductivity > Isotropic
In the window that appears, enter the following geometric properties for steel:
i. KXX: 60.5
8. Define Mesh Size
Preprocessor > Meshing > Size Cntrls > ManualSize > Lines > All Lines... For
this example we will use an element edge length of 0.1 meters.
9. Mesh the frame
Preprocessor > Meshing > Mesh > Lines > click 'Pick All'
10. Write Environment
The thermal environment (the geometry and thermal properties) is now fully described and
can be written to memory to be used at a later time.
Preprocessor > Physics > Environment > Write
In the window that appears, enter the TITLE Thermal and click OK.
If the Physics option is not available under Solution, click Unabridged Menu at the bottom of the
Solution menu. This should make it visible.
3. Apply Constraints
Solution > Define Loads > Apply > Thermal > Temperature > On Keypoints Set
the temperature of Keypoint 1, the left-most point, to 348 Kelvin.
4. Solve the System
Solution > Solve > Current LS
SOLVE
5. Close the Solution Menu
Main Menu > Finish
It is very important to click Finish as it closes that environment and allows a new one to be
opened without contamination. If this is not done, you will get error messages.
The thermal solution has now been obtained. If you plot the steady-state temperature on the link, you
will see it is a uniform 348 K, as expected. This information is saved in a file labelled Jobname.rth,
were .rth is the thermal results file. Since the jobname wasn't changed at the beginning of the analysis,
this data can be found as file.rth. We will use these results in determing the structural effects.
6. Read in the Structural Environment
Solution > Physics > Environment > Read
Choose struct and click OK.
7. Apply Constraints
Solution > Define Loads > Apply > Structural > Displacement > On Keypoints
Fix Keypoint 1 for all DOF's and Keypoint 2 in the UX direction.
8. Include Thermal Effects
Solution > Define Loads > Apply > Structural > Temperature > From Therm Analy
As shown below, enter the file name File.rth. This couples the results from the solution of the
thermal environment to the information prescribed in the structural environment and uses it during
the analysis.
The following list should appear. Note the stress in each element: -0.180e9 Pa, or 180 MPa in
compression as expected.
RESULT:
Ex. No. : 14. CONDUCTIVE HEAT TRANSFER ANALYSIS OF A 2D COMPONENT
Date:
AIM:
To determine the temperature distribution and heat flux of a given 2D component using FEA based
ANSYS software and view the temperature and heat flux.
PROCEDURE:
The Simple Conduction Example is constrained as
shown in the following figure. Thermal
conductivity (k) of the material is 10 W/m*C and
the block is assumed to be infinitely long.
3. Create geometry
Preprocessor > Modeling > Create > Areas > Rectangle > By 2 Corners > X=0, Y=0,
Width=1, Height=1
For this example, we will use PLANE55 (Thermal Solid, Quad 4node 55). This element has 4 nodes and a
single DOF (temperature) at each node. PLANE55 can only be used for 2 dimensional steadystate or
transient thermal analysis.
5. Element Material Properties
Preprocessor > Material Props > Material Models > Thermal > Conductivity > Isotropic >
KXX = 10 (Thermal conductivity)
6. Mesh Size
Preprocessor > Meshing > Size Cntrls > ManualSize > Areas > All Areas > 0.05
7. Mesh
Preprocessor > Meshing > Mesh > Areas > Free > Pick All
Solution Phase: Assigning Loads and Solving
1. Define Analysis Type
Solution > Analysis Type > New Analysis > Steady-State
2. Apply Constraints
For thermal problems, constraints can be in the form of Temperature, Heat Flow, Convection, Heat
Flux, Heat Generation, or Radiation. In this example, all 4 sides of the block have fixed temperatures.
o Solution > Define Loads > Apply
Note that all of the -Structural- options cannot be selected. This is due to the type of element
(PLANE55) selected.
o Thermal > Temperature > On Nodes
o Click the Box option (shown below) and draw a box around the nodes on the top line.
o Fill the window in as shown to constrain the side to a constant temperature of 500
o Using the same method, constrain the remaining 3 sides to a constant value of 100
Orange triangles in the graphics window indicate the temperature contraints.
3. Solve the System
Solution > Solve > Current LS
Note that due to the manner in which the boundary contitions were applied, the top corners are held at a
temperature of 100. Recall that the nodes on the top of the plate were constrained first, followed by the
side and bottom constraints. The top corner nodes were therefore first constrained at 500C, then
'overwritten' when the side constraints were applied. Decreasing the mesh size can minimize this
effect, however, one must be aware of the limitations in the results at the corners.
RESULT:
Ex. No. : 15. CONVECTIVE HEAT TRANSFER ANALYSIS OF A 2D COMPONENT
Date :
AIM:
To determine the temperature distribution and heat flux of a given 2D component using FEA based
ANSYS software and view the temperature and heat flux.
PEROCEDURE:
3. Create geometry
Preprocessor > Modeling > Create > Areas > Rectangle > By 2 Corners > X=0, Y=0,
Width=1, Height=1
As in the conduction example, we will use PLANE55 (Thermal Solid, Quad 4node 55). This element has 4
nodes and a single DOF (temperature) at each node. PLANE55 can only be used for 2
dimensional steady-state or transient thermal analysis.
5. Element Material Properties
Preprocessor > Material Props > Material Models > Thermal > Conductivity > Isotropic >
KXX = 10
7. Mesh
Preprocessor > Meshing > Mesh > Areas > Free > Pick All
o Fill in the window as shown. This will specify a convection of 10 W/m2*C and an ambient
temperature of 100 degrees Celcius. Note that VALJ and VAL2J have been left blank. This is
because we have uniform convection across the line.
4. Apply Insulated Boundary Conditions
o Solution > Define Loads > Apply > Thermal > Convection > On Lines
o Select the bottom of the block.
o Enter a constant Film coefficient (VALI) of 0. This will eliminate convection through the
side, thereby modeling an insulated wall. Note: you do not need to enter a Bulk (or ambient)
temperature
You should obtain the following:
5. Solve the System
RESULT: