Space Frame Example: Problem Description
Space Frame Example: Problem Description
Space Frame Example: Problem Description
html
Modeling Tools
| Verification Example | | Preprocessing | | Solution | | Postprocessing | | Command Line |
Solid Modeling | Bicycle Example | | Preprocessing | | Solution | | Postprocessing | | Command Line |
Index Introduction
This tutorial was created using ANSYS 7.0 to solve a simple 3D space frame problem.
Contributions
Problem Description
Comments
MecE 563 The problem to be solved in this example is the analysis of a bicycle frame. The problem to be modeled in this example is
a simple bicycle frame shown in the following figure. The frame is to be built of hollow aluminum tubing having an
Mechanical Engineering
outside diameter of 25mm and a wall thickness of 2mm.
University of Alberta
ANSYS Inc.
Copyright © 2001
University of Alberta
Verification
The first step is to simplify the problem. Whenever you are trying out a new analysis type, you need something (ie
analytical solution or experimental data) to compare the results to. This way you can be sure that you've gotten the
correct analysis type, units, scale factors, etc.
The simplified version that will be used for this problem is that of a cantilever beam shown in the following figure:
1 of 18 5/2/2011 11:13 AM
U of A ANSYS Tutorials - Bicycle Space Frame http://www.mece.ualberta.ca/tutorials/ansys/BT/Bike/Bike.html
2. Enter Keypoints
For this simple example, these keypoints are the ends of the beam.
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
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.
When you're done, click on 'OK' in the 'Create Straight Line' window.
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.
2 of 18 5/2/2011 11:13 AM
U of A ANSYS Tutorials - Bicycle Space Frame http://www.mece.ualberta.ca/tutorials/ansys/BT/Bike/Bike.html
Click on 'Close' in the 'Element Types' dialog box and close the 'Element Type' menu.
Click Add... and select 'Type 1 PIPE16' (actually it is already selected). Click on 'OK'.
This defines an outside pipe diameter of 25mm and a wall thickness of 2mm.
Click on 'OK'.
'Set 1' now appears in the dialog box. Click on 'Close' in the 'Real Constants' window.
Double click Structural > Linear > Elastic and select 'Isotropic' (double click on it)
We are going to give the properties of Aluminum. Enter the following field:
EX 70000
PRXY 0.33
7. Mesh Size
In the Preprocessor menu select Meshing > Size Cntrls > ManualSize > Lines > All Lines
In the size 'SIZE' field, enter the desired element length. For this example we want an element length of 2cm,
therefore, enter '20' (i.e 20mm) and then click 'OK'. Note that we have not yet meshed the geometry, we
have simply defined the element sizes.
(Alternatively, we could enter the number of divisions we want in the line. For an element length of 2cm, we
would enter 25 [ie 25 divisions]).
NOTE
It is not necessary to mesh beam elements to obtain the correct solution. However, meshing is done in this case so
that we can obtain results (ie stress, displacement) at intermediate positions on the beam.
8. Mesh
In the 'Preprocessor' menu select Meshing > Mesh > Lines and click 'Pick All' in the 'Mesh Lines' Window
Utility Menu > File > Save as.... Select the name and location where you want to save your file.
From the Solution Menu, select 'Analysis Type > New Analysis'.
2. Apply Constraints
3 of 18 5/2/2011 11:13 AM
U of A ANSYS Tutorials - Bicycle Space Frame http://www.mece.ualberta.ca/tutorials/ansys/BT/Bike/Bike.html
In the Solution menu, select Define Loads > Apply > Structural > Displacement > On Keypoints
Select the left end of the rod (Keypoint 1) by clicking on it in the Graphics Window and click on 'OK' in the
'Apply U,ROT on KPs' window.
This location is fixed which means that all translational and rotational degrees of freedom (DOFs) are
constrained. Therefore, select 'All DOF' by clicking on it and enter '0' in the Value field and click 'OK'.
3. Apply Loads
As shown in the diagram, there is a vertically downward load of 100N at the end of the bar
Select the second Keypoint (right end of bar) and click 'OK' in the 'Apply F/M' window.
Enter a value of -100 in the 'Force/moment value' box and click 'OK'.
The applied loads and constraints should now appear as shown below.
Now, since the purpose of this exercise was to verify the results - we need to calculate what we should find.
Deflection:
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
4 of 18 5/2/2011 11:13 AM
U of A ANSYS Tutorials - Bicycle Space Frame http://www.mece.ualberta.ca/tutorials/ansys/BT/Bike/Bike.html
stress).
Deformation
from the Main Menu select General Postproc from the 'ANSYS Main Menu'. In this menu you will find a
variety of options, the two which we will deal with now are 'Plot Results' and 'List Results'
Select 'Def + undef edge' and click 'OK' to view both the deformed and the undeformed object.
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
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'.
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
5 of 18 5/2/2011 11:13 AM
U of A ANSYS Tutorials - Bicycle Space Frame http://www.mece.ualberta.ca/tutorials/ansys/BT/Bike/Bike.html
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'
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'
6 of 18 5/2/2011 11:13 AM
U of A ANSYS Tutorials - Bicycle Space Frame http://www.mece.ualberta.ca/tutorials/ansys/BT/Bike/Bike.html
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 ...
Note that a maximum stress of 64.914 MPa occurs at the fixed end of the beam as predicted analytically.
To further verify the simplified model, a bending moment diagram can be created. First, let's look at how ANSYS
defines each element. Pipe 16 has 2 nodes; I and J, as shown in the following image.
To obtain the bending moment for this element, the Element Table must be used. The Element Table contains most
of the data for the element including the bending moment data for each element at Node I and Node J. First, we
need to obtain obtain the bending moment data.
General Postproc > Element Table > Define Table... . Click 'Add...'.
7 of 18 5/2/2011 11:13 AM
U of A ANSYS Tutorials - Bicycle Space Frame http://www.mece.ualberta.ca/tutorials/ansys/BT/Bike/Bike.html
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.
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.
8 of 18 5/2/2011 11:13 AM
U of A ANSYS Tutorials - Bicycle Space Frame http://www.mece.ualberta.ca/tutorials/ansys/BT/Bike/Bike.html
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.
The above example was solved using the Graphical User Interface (or GUI) of ANSYS. This problem can also been
solved using the ANSYS command language interface. To see the benefits of the command line clear your current file:
From the Utility menu select: File > Clear and Start New
Ensure that 'Read File' is selected then click 'OK'
select 'yes' in the following window.
Copy the following code into the command line, then hit enter. Note that the text following the "!" are comments.
/PREP7 ! Preprocessor
K,1,0,0,0, ! Keypoint, 1, x, y, z
K,2,500,0,0, ! Keypoint, 2, x, y, z
L,1,2 ! Line from keypoint 1 to 2
!*
ET,1,PIPE16 ! Element Type = pipe 16
KEYOPT,1,6,1 ! This is the changed option to give the extra force and moment output
9 of 18 5/2/2011 11:13 AM
U of A ANSYS Tutorials - Bicycle Space Frame http://www.mece.ualberta.ca/tutorials/ansys/BT/Bike/Bike.html
!*
R,1,25,2, ! Real Constant, Material 1, Outside Diameter, Wall thickness
!*
MP,EX,1,70000 ! Material Properties, Young's Modulus, Material 1, 70000 MPa
MP,PRXY,1,0.33 ! Material Properties, Major Poisson's Ratio, Material 1, 0.33
!*
LESIZE,ALL,20 ! Element sizes, all of the lines, 20 mm
LMESH,1 ! Mesh the lines
FINISH ! Exit preprocessor
/SOLU ! Solution
ANTYPE,0 ! The type of analysis (static)
!*
DK,1, ,0, ,0,ALL ! Apply a Displacement to Keypoint 1 to all DOF
FK,2,FY,-100 ! Apply a Force to Keypoint 2 of -100 N in the y direction
/STATUS,SOLU
SOLVE ! Solve the problem
FINISH
Note that you have now finished Postprocessing and the Solution Phase with just these few lines of code. There are codes
to complete the Postprocessing but we will review these later.
Bicycle Example
Now we will return to the analysis of the bike frame. The steps which you completed in the verification example will not
be explained in great detail, therefore use the verification example as a reference as required. We will be combining the
use of the Graphic User Interface (GUI) with the use of command lines.
We are going to define the vertices of the frame using variables. These variables represent the various lengths of
the bicycle members. Notice that by using variables like this, it is very easy to set up a parametric description of
your model. This will enable us to quickly redefine the frame should changes be necessary. The quickest way to
enter these variables is via the 'ANSYS Input' window which was used above to input the command line codes for
the verification model. Type in each of the following lines followed by Enter.
x1 = 500
x2 = 825
y1 = 325
y2 = 400
z1 = 50
4. Enter Keypoints
For this space frame example, these keypoints are the frame vertices.
10 of 18 5/2/2011 11:13 AM
U of A ANSYS Tutorials - Bicycle Space Frame http://www.mece.ualberta.ca/tutorials/ansys/BT/Bike/Bike.html
We are going to define 6 keypoints for this structure as given in the following table (these keypoints are
depicted by the circled numbers in the above figure):
coordinate
keypoint
x y z
1 0 y1 0
2 0 y2 0
3 x1 y2 0
4 x1 0 0
5 x2 0 z1
6 x2 0 -z1
Now instead of using the GUI window we are going to enter code into the 'command line'. First, open the
'Preprocessor Menu' from the 'ANSYS Main Menu'. The preprocessor menu has to be open in order for the
preprocessor commands to be recognized. Alternatively, you can type /PREP7 into the command line. The
command line format required to enter a keypoint is as follows:
K, NPT, X, Y, Z
For a more detailed explanation, type help k into the command line
As with any programming language, you may need to add comments. The exclamation mark indicates that
anything following it is commented out. ie - for the second keypoint you might type:
K,2,0,y2,0 ! keypoint, #, x=0, y=y2, z=0
Enter the 4 remaining keypoints (listed in the table above) using the command line
Now you may want to check to ensure that you entered all of the keypoints correctly:
Utility Menu > List > Keypoints > Coordinates only
(Alternatively, type 'KLIST' into the command line)
If there are any keypoints which need to be re-entered, simply re-enter the code. A previously defined
keypoint of the same number will be redefined. However, if there is one that needs to be deleted simply enter
the following code:
KDELE,#
In this example, we defined the keypoints by making use of previously defined variables like y1 = 325. This was
simply used for convenience. To define keypoint #1, for example, we could have alternatively used the coordinates
x = 0, y = 325, z = 0.
To get a better view of our view of our model, we'll view it in an isometric view:
Select Utility menu bar > PlotCtrls > Pan, Zoom, Rotate...'
11 of 18 5/2/2011 11:13 AM
U of A ANSYS Tutorials - Bicycle Space Frame http://www.mece.ualberta.ca/tutorials/ansys/BT/Bike/Bike.html
6. Create Lines
keypoint
line
1st 2nd Again, we will use the command line to create the lines. The command format to create a
straight line looks like:
1 1 2
L, P1, P2
2 2 3 Line, Keypoint at the beginning of the line, Keypoint at the end of line
3 3 4 For example, to obtain the first line, I would write: ' L,1,2 '
4 1 4
Note: unlike 'Keypoints', 'Lines' will automatically assign themselves the next available
5 3 5 reference number.
6 4 5
7 3 6
8 4 6
Enter the remaining lines until you get a picture like that shown below.
Again, check to ensure that you entered all of the lines correctly: type ' LLIST ' into the command line
If there are any lines which need to be changed, delete the line by typing the following code: ' LDELE,# '
where # corresponds to the reference number of the line. (This can be obtained from the list of lines). And
then re-enter the line (note: a new reference number will be assigned)
12 of 18 5/2/2011 11:13 AM
U of A ANSYS Tutorials - Bicycle Space Frame http://www.mece.ualberta.ca/tutorials/ansys/BT/Bike/Bike.html
As in the verification model, define the type of element (pipe16). As in the verification model, don't forget to
change Option K6 'Include Output' to obtain extra force and moment output.
To set Young's Modulus and Poisson's ratio, we will again use the command line. (ensure that the preprocessor
menu is still open - if not open it by clicking Preprocessor in the Main Menu)
MP, LAB, MAT, C0
Material Property,Valid material property label, Material Reference Number, value
To enter the Elastic Modulus (LAB = EX) of 70000 MPa, type: ' MP,EX,1,70000 '
11. Mesh
In the 'Preprocessor' menu select 'Mesh' > 'Lines' and click 'Pick All' in the 'Mesh Lines' Window
Apply Constraints
13 of 18 5/2/2011 11:13 AM
U of A ANSYS Tutorials - Bicycle Space Frame http://www.mece.ualberta.ca/tutorials/ansys/BT/Bike/Bike.html
Once again, we will use the command line. We are going to pin (translational DOFs will be fixed) the first keypoint and
constrain the keypoints corresponding to the rear wheel attachment locations in both the y and z directions. The following
is the command line format to apply constraints at keypoints.
DK, KPOI, Lab, VALUE, VALUE2, KEXPND, Lab2, Lab3, Lab4, Lab5, Lab6
Displacement on K, K #, DOF label, value, value2, Expansion key, other DOF labels
Not all of the fields are required for this example, therefore when entering the code certain fields will be empty. For
example, to pin the first keypoint enter:
DK,1,UX,0,,,UY,UZ
The DOF labels for translation motion are: UX, UY, UZ. Note that the 5th and 6th fields are empty. These correspond to
'value2' and 'the Expansion key' which are not required for this constraint. Also note that all three of the translational
DOFs were constrained to 0. The DOFs can only be contrained in 1 command line if the value is the same.
Note that only UY and UZ are contrained to 0. UX is not constrained. Again, note that the 5th and 6th fields are empty
because they are not required.
Apply the constraints to the other rear wheel location (Keypoint 6 - UY and UZ).
Now list the constraints ('DKLIST') and verify them against the following:
If you need to delete any of the constraints use the following command: 'DKDELE, K, Lab' (ie 'DKDELE,1,UZ'
would delete the constraint in the 'z' direction for Keypoint 1)
Apply Loads
We will apply vertical downward loads of 600N at the seat post location (keypoint 3) and 200N at the pedal crank
location (keypoint 4). We will use the command line to define these loading conditions.
FK, KPOI, Lab, value, value2
Force loads at keypoints, K #, Force Label directions (FX, FY, FZ), value1, value2 (if req'd)
To apply a force of 600N downward at keypoint 3, the code should look like this: ' FK,3,FY,-600 '
Apply both the forces and list the forces to ensure they were inputted correctly (FKLIST).
If you need to delete one of the forces, the code looks like this: 'FKDELE, K, Lab' (ie 'FKDELE,3,FY' would delete the
force in the 'y' direction for Keypoint 3)
The applied loads and constraints should now appear as shown below.
14 of 18 5/2/2011 11:13 AM
U of A ANSYS Tutorials - Bicycle Space Frame http://www.mece.ualberta.ca/tutorials/ansys/BT/Bike/Bike.html
1. Deformation
Plot Results > Deformed Shape... 'Def + undef edge'
You may want to try plotting this from different angles to get a better idea what's going on by using the
'Pan-Zoom-Rotate' menu that was earlier outlined.
Try the 'Front' view button (Note that the views of 'Front', 'Left', 'Back', etc depend on how the object was
first defined).
Your screen should look like the plot below:
15 of 18 5/2/2011 11:13 AM
U of A ANSYS Tutorials - Bicycle Space Frame http://www.mece.ualberta.ca/tutorials/ansys/BT/Bike/Bike.html
2. Deflections
Now let's take a look at some actual deflections in the frame. The deflections have been calculated at the nodes of
the model, so the first thing we'll do is plot out the nodes and node numbers, so we know what node(s) we're after.
Go to Utility menu > PlotCtrls > Numbering... and turn on 'Node numbers'. Turn everything else off.
Note the node numbers of interest. Of particular interest are those nodes where the constraints were applied
to see if their displacements/rotations were indeed fixed to zero. Also note the node numbers of the seat and
crank locations.
List the Nodal Deflections (Main Menu > General Postproc > List Results > Nodal Solution...'). Are the
displacements and rotations as you expected?
3. Element Forces
We could also take a look at the forces in the elements in much the same way:
16 of 18 5/2/2011 11:13 AM
U of A ANSYS Tutorials - Bicycle Space Frame http://www.mece.ualberta.ca/tutorials/ansys/BT/Bike/Bike.html
Click on 'OK'.
For each element in the model, the force/moment values at each of the two nodes per element will be
displayed.
Close this list window when you are finished browsing.
Then close the 'List Results' menu.
4. Stresses
As shown in the cantilever beam example, use the Element Table to gain access to derived stresses.
Select Element Table > Define Table... to define the table (remember SMISC,6 and SMISC,12)
And, Plot Results > Line Elem Res... to plot the data from the Element Table
17 of 18 5/2/2011 11:13 AM
U of A ANSYS Tutorials - Bicycle Space Frame http://www.mece.ualberta.ca/tutorials/ansys/BT/Bike/Bike.html
Quitting ANSYS
To quit ANSYS, select 'QUIT' from the ANSYS Toolbar or select 'Utility Menu'/'File'/'Exit...'. In the dialog box that
appears, click on 'Save Everything' (assuming that you want to) and then click on 'OK'.
18 of 18 5/2/2011 11:13 AM