Exavator
Exavator
Exavator
Edition Note
This document describes the release information of RecurDyn V9R1.
Table of Contents
Getting Started ....................................................................... 1
Objective ............................................................................................... 1
Approach ............................................................................................... 1
Audience ................................................................................................ 2
Prerequisites........................................................................................... 2
Procedures ............................................................................................. 3
Estimated Time to Complete ..................................................................... 3
Opening the Initial Model .......................................................... 4
Task Objective ........................................................................................ 4
Estimated Time to Complete ..................................................................... 4
Starting RecurDyn ................................................................................... 5
Running an Initial Simulation with the Rigid Boom ....................................... 7
Viewing the Results ................................................................................. 7
Swapping In the RFlex Body ..................................................... 8
Task Objective ........................................................................................ 8
Estimated Time to Complete ..................................................................... 8
Swapping In the RecurDyn RFlex Body ....................................................... 9
Viewing Stress Contour Plots .................................................................. 11
Plotting the Results ................................................................. 14
Task Objective ...................................................................................... 14
Estimated Time to Complete ................................................................... 14
Plotting the Dipper Stick’s Out-of-Plane Tilt .............................................. 15
RFlex Body Review and Tuning ................................................. 18
Task Objective ...................................................................................... 18
Estimated Time to Complete ................................................................... 18
Examining the RFlex Body ...................................................................... 19
Improving Simulation Performance .......................................................... 22
Appendix A : Creating the RecurDyn RFlex Input (RFI) File ........ 25
Task Objective ...................................................................................... 25
Estimated Time to Complete ................................................................... 25
Preparing the Nastran Bulk Data File........................................................ 26
Component Mode Reduction Method ........................................................ 26
Superelement Method ............................................................................ 29
Appendix B: Supported FE Elements ......................................... 34
Ansys Element Library ........................................................................... 34
MSC/NASTRAN Element Library ............................................................... 35
I-DEAS Element Library ......................................................................... 35
E X C A V A T O R T U T O R I A L ( R F L E X )
1
Chapter
Getting Started
Objective
In this tutorial, you will learn how to simulate a model which has a flexible body. You will
start with an existing model which has all rigid bodies, and replace one of them with a
RecurDyn RFlex body. RecurDyn RFlex uses a modal approach to represent flexible bodies, by
using the superposition of the body’s natural modal shapes and constraint modes at the
attachment points. This approach is ideal when the flexible body has only fixed connection
points to other bodies. If it has sliding or rolling contact, RecurDyn FFlex should be used,
which takes a nodal (or mesh-based) approach to represent the body.
Approach
You will start with an excavator model, shown below, that includes all of the mechanical
components. All the bodies in this model are rigid. The model is set up so that the excavator
goes through a digging and unloading motion, with the cab rotating around the vertical axis.
You will then import a RecurDyn RFlex Input (RFI) file, which represents a flexible boom. One
way to create an RFI file is by using NX Advance Simulation to mesh the boom geometry,
export the simulation to a Nastran bulk data file, and then use NX Nastran to read the bulk
data file and create the RFI file.
The RFI file is supplied with this tutorial. For information on how to create the RFI file from a
Nastran bulk data file, see the appendix.
1
E X C A V A T O R T U T O R I A L ( R F L E X )
Audience
This tutorial is intended for intermediate users of RecurDyn who previously learned how to
create geometry, joints, and force entities. All new tasks are explained carefully.
Prerequisites
You should first work through the 3D Crank-Slider and Engine with Propeller tutorials, or the
equivalent. We assume that you have a basic knowledge of physics.
You will need a license for the RFlex module of RecurDyn.
2
E X C A V A T O R T U T O R I A L ( R F L E X )
Procedures
The tutorial is comprised of the following procedures. The estimated time to complete each
procedure is shown in the table.
Total 45
3
E X C A V A T O R T U T O R I A L ( R F L E X )
2
Chapter
4
E X C A V A T O R T U T O R I A L ( R F L E X )
Starting RecurDyn
To start RecurDyn and open the initial
model:
5
E X C A V A T O R T U T O R I A L ( R F L E X )
6
E X C A V A T O R T U T O R I A L ( R F L E X )
From the Animation Control group in the Analysis tab, click Play/Pause.
You should see the excavator rotate about the vertical axis, and then goes through a digging
and then unloading motion. This is being driven by motion input to the revolute joint around
which the cab rotates, and translational joints in the hydraulic cylinders.
7
E X C A V A T O R T U T O R I A L ( R F L E X )
3
Chapter
8
E X C A V A T O R T U T O R I A L ( R F L E X )
1. From the RFlex group in the Flexible tab, click Import RFI.
2. Set the Creation Method to Body in Command Toolbar.
3. In the Working window, select
the Rigid_Boom body.
The RFlex boom body should now have replaced the rigid boom body, as shown below (may
be a different color than green). By checking in the Database window, you should notice
that the rigid boom body is no longer in the model.
9
E X C A V A T O R T U T O R I A L ( R F L E X )
Important: In the next step, it is important to save the model using a different filename,
as the results from the new file will be compared with the old.
10
E X C A V A T O R T U T O R I A L ( R F L E X )
6. Now adjust the maximum value so that the contour display will show stresses at a lower
range:
Change the Type to User Defined.
Type in 200 for the Max value.
11
E X C A V A T O R T U T O R I A L ( R F L E X )
7. Click OK
8. Play the animation.
You should see the contour plots on the flexible boom (animation frame 14 is shown below)
12
E X C A V A T O R T U T O R I A L ( R F L E X )
9. If the animation plays slowly on your computer, you may want to do one of the following:
Use the Fast Play button to display every fifth frame.
Record the animation to an avi file and view the animation with the Windows Media
Player.
13
E X C A V A T O R T U T O R I A L ( R F L E X )
4
Chapter
Task Objective
Learn how to plot and compare results from different models, and observe the influence of
modeling the boom as a flexible body.
14
E X C A V A T O R T U T O R I A L ( R F L E X )
1. From the Plot group in the Analysis tab, click Plot Result.
The current model’s results are automatically loaded into the Database window on the
right. You will now load the old rigid-body-only model’s results.
2. In the File Menu click Import file and select the file RD_Excavator_Start.rplt.
The old model’s results should now appear in the Database window, under
RD_Excavator_RFlex.
3. Plot the following results:
15
E X C A V A T O R T U T O R I A L ( R F L E X )
Note: The Request items appear in the Plot Database window because a request to create
plot data for an expression was created in the model. If you were to go back into the model,
you would see an expression called Ex_dipperStickTilt which has the following form:
AX(DipperStick.CM, Cab.CM)
This expression measures the rotation of the dipper stick’s center of mass marker around the
X-axis of the cab’s center of mass marker.
4. You should now see a plot similar to the one shown below.
The plot shows that the amplitude of tilt with the flexible boom is much higher than with the
rigid boom, and the system natural frequency (as reflected by the time interval between
some of the peaks) is lower as well.
To make more sense of this plot, you will now add the rotational acceleration of the cab --
the motion input which causes the dipper stick’s out-of-plane tilt. But first, you will create an
extra Y-axis so that both tilt and acceleration data can be displayed on the plot with similar
scales.
16
E X C A V A T O R T U T O R I A L ( R F L E X )
You can now see the motion input and the response of the flexible and rigid boom to that
input. From 0.5 to 2.0 sec, the flexible boom’s transient response takes much longer to die
out than the rigid boom’s response. The plot clearly shows the importance of including the
flexibility of the boom in this excavator model.
Note that there is some noise in the rotational acceleration curve of the cab. The noise is the
result of using the STEP function with a small number of integration steps. The noise will be
almost totally eliminated by decreasing the Maximum Time Step under the Parameter tab of
the Dynamic/Kinematic Analysis dialog box. However, with this change the simulation will
take 2-3 times as long to solve. Therefore, for purposes of this tutorial (fast simulations) the
default solver parameters will be used.
17
E X C A V A T O R T U T O R I A L ( R F L E X )
5
Chapter
18
E X C A V A T O R T U T O R I A L ( R F L E X )
19
E X C A V A T O R T U T O R I A L ( R F L E X )
Reviewing the low frequency modes will reveal which parts of the structure are the weakest
and most prone to vibration.
20
E X C A V A T O R T U T O R I A L ( R F L E X )
Notice that the areas of highest deformation are much localized and the structure is highly
distorted. Also, the frequency of this mode is quite high, at 1007.28 Hz. If this mode does
not appear realistic based on experimental results for similar parts, it would be reasonable to
throw this mode out of the solution to increase simulation performance.
Also notice that the Damping Ratio for this mode is 1, which is relatively high. The Damping
Ratio controls how much a particular mode contributes to the overall behavior of the flexible
body. If it is set to a low value, the mode will substantially contribute to the overall behavior.
If the damping is set high, the mode will be damped out quickly and will not substantially
contribute to the overall behavior. Starting at mode 7, the Damping Ratio is 0.01. For
mode 13 (107.46 Hz), the Damping Ratio increases to 0.1. Finally, at mode 68 (1007.28
Hz), the Damping Ratio increases to 1, meaning that as is, mode 68 will not contribute much
to the overall behavior of the structure.
By default, RecurDyn assigns Damping Ratio values based on modal frequency, as follows:
0 < f < 100 Hz: Damping Ratio = 0.01
100 ≤ f < 1000 Hz: Damping Ratio = 0.1
1000 Hz ≤ f: Damping Ratio = 1
This scheme works well with large structures such as those found in construction equipment
and automobiles. However, if you have a small application where high frequency modes are
important, or otherwise would like to assign the Damping Ratios yourself, the fields are
editable so that you can do that. You can also import a file containing the damping ratios for
every mode. For more information on how to do this, please see the RFlex User’s Guide.
21
E X C A V A T O R T U T O R I A L ( R F L E X )
The results from this simulation should now be compared to the original results. If there is
no significant change, then the current selection of modes can be used for further simulations.
22
E X C A V A T O R T U T O R I A L ( R F L E X )
This plot shows that the latest results for dipper stick tilt match the original results very well
(yellow line overlays blue line). In fact, if you were to compare the values for peak tilt at
roughly 2.85 sec, you would see that the difference is less than 1%. Therefore, removing
the modes above 1000 Hz successfully reduced the simulation time while retaining good
results.
23
E X C A V A T O R T U T O R I A L ( R F L E X )
24
E X C A V A T O R T U T O R I A L ( R F L E X )
A
Appendix
Task Objective
Learn what must be added to a Nastran bulk data file so that when processed with NX
Nastran, an RFI file complete with stress contour information will be created. Note that a
similar process is followed when creating a RFI file from Ansys output. This process is
described in the RFlex documentation.
25
E X C A V A T O R T U T O R I A L ( R F L E X )
26
E X C A V A T O R T U T O R I A L ( R F L E X )
$ - Output grid point stress and strain for all SURFACE and VOLUME commands.
GPSTRAIN=ALL
GPSTRESS=ALL
$ - Generate RecurDyn RFlex Input (RFI) file. Here, DMAP solution is turned
$ off, and grid point stress and strain are output to the RFI file.
$ - If your version of NX Nastran is 6.1 or later, use the following command:
MBDEXPORT RECURDYN FLEXBODY=YES,FLEXONLY=YES,OUTGSTRS=YES,OUTGSTRN=YES
$ - Otherwise, if your version of NX Nastran is earlier than 6.1, use this
$ command:
RECURDYNRFI FLEXBODY=YES,FLEXONLY=YES,OUTGSTRS=YES,OUTGSTRN=YES
$*
$ - Select the real eigenvalue extraction parameters for component mode
$ reduction.
RSMETHOD = 100
$ - Select the real eigenvalue extraction parameters.
METHOD = 101
$ - Output displacement of all points.
VECTOR(SORT1,REAL)=ALL
$ - Output all single-point forces of constrain
SPCFORCES(SORT1,REAL)=ALL
$ - Define Set 5 as the same as Set 2 defined above.
SET 5 = 1 THRU 4086
$ - Output stress and strain for elements defined in Set 5, above.
STRESS=5
STRAIN=5
$*
$ - Indicate beginning of surface or volume commands.
OUTPUT(POST)
$ - Define Set 6 as the same as Set 2 defined above.
SET 6 = 1 THRU 4086
$ - Set the volume for which strains and stresses are calculated. Here, direct
$ stresses and strains are requested.
VOLUME 1 SET 6,DIRECT,SYSTEM CORD 0
$ - NOTE: If shell elements are used in the mesh, the SURFACE command should
$ be used instead, as shown below:
$
$ SURFACE 1 SET 6,FIBRE ALL,SYSTEM CORD 0
$
$ If you will be displaying contour plots, though, it is recommended that
27
E X C A V A T O R T U T O R I A L ( R F L E X )
$ solid meshes be used to create RFlex bodies. This is because, for shell
$ elements, RecurDyn only displays contour plots for midplane stress, strain,
$ and displacement – data for the top or bottom surfaces cannot be displayed.
$*
$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$*
$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$*
$* BULK DATA
$*
$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$*
BEGIN BULK
$ - Define the units of measurement to be used.
DTI,UNITS,1,KG,MN,MM,S
$
$ - Select nodes as the connection points to the flexible body. Here, nodes
$ 9001 – 9008 are all master nodes of RBE2 elements in the mesh.
$
ASET1,123456,9001,THRU,9008
$*
$*
$* SOLUTION CARDS
$*
$ - First modal solution:
$ - Specify the frequency range or number of constrained normal
$ modes desired.
$ - Frequency range should be at least 2x the range of interest in the MBD
$ solution.
$ - In this solution, ASET DOF are constrained.
$ - This is selected by RSMETHOD in the Case Control section, above.
$
EIGRL 100 40 0 7 MASS
$
$ - Modal reduction DOFs:
$ - Number of SPOINTs requested (ns) should be as follows:
$ ns >= n + (6 + p)
$ where:
$ n = number of modes requested in first modal solution (in this case,
28
E X C A V A T O R T U T O R I A L ( R F L E X )
Superelement Method
The following code must be added to the Nastran file, before the GRID CARDS section, in
order to create an RFI file using the Superelement method. See comments in red for
descriptions of key commands:
29
E X C A V A T O R T U T O R I A L ( R F L E X )
$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$*
$* EXECUTIVE CONTROL
$*
$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$*
ID,NASTRAN,recurdyn_rfi_create_se
$ - Set the solution type to SEMODES, solving for the normal modes.
SOL 103
$ - Set the maximum CPU time to 999 sec.
TIME 999
CEND
$*
$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$*
$* CASE CONTROL
$*
$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$*
$ - Define a Set which contains all of the node IDs.
SET 1 = 1 THRU 9008
$ - Define a Set which contains all of the element IDs.
SET 2 = 1 THRU 4086
$ - Generate and assemble all superelements.
SEALL = ALL
$ - Assign the subcase to all superelements and loading conditions.
SUPER = ALL
$ - Turn printing of bulk data off.
ECHO = NONE
$ - Output grid point stress and strain for all SURFACE and VOLUME commands.
GPSTRAIN=ALL
GPSTRESS=ALL
$ - Generate RecurDyn RFlex Input (RFI) file. Here, DMAP solution is turned
$ off, and grid point stress and strain are output to the RFI file.
$ - If your version of NX Nastran is 6.1 or later, use the following command:
MBDEXPORT RECURDYN FLEXBODY=YES,FLEXONLY=YES,OUTGSTRS=YES,OUTGSTRN=YES
$ - Otherwise, if your version of NX Nastran is earlier than 6.1, use this
$ command:
RECURDYNRFI FLEXBODY=YES,FLEXONLY=YES,OUTGSTRS=YES,OUTGSTRN=YES
30
E X C A V A T O R T U T O R I A L ( R F L E X )
$*
$ - Select the real eigenvalue extraction parameters.
METHOD = 100
$ - Output displacement of all points.
VECTOR(SORT1,REAL)=ALL
$ - Output all single-point forces of constraint.
SPCFORCES(SORT1,REAL)=ALL
$ - Define Set 5 as the same as Set 2 defined above.
SET 5 = 1 THRU 4086
$ - Output stress and strain for elements defined in Set 5, above.
STRESS=5
STRAIN=5
$*
$ - Indicate beginning of surface or volume commands.
OUTPUT(POST)
$ - Define Set 6 as the same as Set 2 defined above.
SET 6 = 1 THRU 4086
$ - Set the volume for which strains and stresses are calculated. Here, direct
$ stresses and strains are requested.
VOLUME 1 SET 6,DIRECT,SYSTEM CORD 0
$ - NOTE: If shell elements are used in the mesh, the SURFACE command should
$ be used instead, as shown below:
$
$ SURFACE 1 SET 6,FIBRE ALL,SYSTEM CORD 0
$
$ If you will be displaying contour plots, though, it is recommended that
$ solid meshes be used to create RFlex bodies. This is because, for shell
$ elements, RecurDyn only displays contour plots for midplane stress, strain,
$ and displacement – data for the top or bottom surfaces cannot be displayed.
$*
$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$*
$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$*
$* BULK DATA
$*
$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$*
BEGIN BULK
31
E X C A V A T O R T U T O R I A L ( R F L E X )
32
E X C A V A T O R T U T O R I A L ( R F L E X )
PARAM GRDPNT 0
PARAM K6ROT 100.0
PARAM MAXRATIO 1.0+8
PARAM POST -2
PARAM POSTEXT YES
PARAM RESVEC YES
PARAM USETPRT 0
$*
The Nastran file is now ready to be processed by NX Nastran, which will produce the RFI file
which can then be imported into a RecurDyn model.
33
E X C A V A T O R T U T O R I A L ( R F L E X )
B
Appendix
Appendix B: Supported FE
Elements
In this chapter, you can see the supported FE elements in RecurDyn. To see more
information, refer to RFlex of RecurDyn Help
34
E X C A V A T O R T U T O R I A L ( R F L E X )
Rigid Element RBAR, RBE2, RBE3, RROD, CBUSH, CBUSH1D, CELAS1, CELAS2
2D Element Thin Shell Linear Triangle, Thin Shell Parabolic Triangle, Thin Shell
Linear Quadrilateral, Thin Shell Parabolic Quadrilateral,
Plane Stress Linear Triangle, Plane Stress Parabolic Triangle, Plane
Stress Linear Quadrilateral, Plane Stress Parabolic Quadrilateral,
Plane Strain Linear Triangle, Plane Strain Parabolic Triangle, Plane
Strain Linear Quadrilateral, Plane Strain Parabolic Quadrilateral
Rigid Element Rigid, Rigid Bar, Node To Node Translational Spring, Node To Node
Rotational Spring
35