Simcenter Testing Solutions Simcenter Testlab: Rigid Body Calculator

2022-05-07T22:46:12.000-0400
Simcenter Testing Solutions Marketing Simcenter SCADAS Simcenter Testlab Simcenter Testlab RT Teamcenter Share Simcenter Testlab Apps Simcenter Testxpress Simcenter Tecware Simcenter Sound Camera Other Hardware Simcenter Soundbrush Simcenter IQT

Summary


Details


Direct YouTube link: https://youtu.be/ISu4FcsV4fE


Simcenter Testlab: Rigid Body Property Calculator 

Need to determine the center of gravity and principal moments of inertia of an object? Simcenter Testlab can help! 

The inertias of the masses contained in a system are important in understanding how they will behave when set in motion. The inertia of an object determines how much it resists or continues to move when forces are applied. 

Some examples where inertial properties are important: 
  • Determining the amount of vibration imparted by a running motor mounted on elastomeric isolators. 
  • Controlling the motion of robot arm composed of an assembly of linkages (each linkage with its own inertia). A control algorithm for the robot arm benefits from knowing the inertia of the linkages so it operates as precisely as possible without overshooting its objective. 
  • Preventing a vehicle from rolling over during high-speed maneuvers (Figure 1). 
User-added image
Figure 1: As a vehicle maneuvers, the inertia of the vehicle body relative to the motion of suspension can cause it to rollover.

Using experimentally measured Frequency Response Functions (FRFs), the center of gravity and inertial properties of an object can be determined using the Rigid Body Property Calculator contained within Simcenter Testlab.  

This article describes the theoretical background, the FRF measurements required, and how to perform the inertial analysis in Simcenter Testlab software. 

Contents: 

1. Moment of Inertia Background 
2. Experimental Methods for Determination of Inertia 
   2.1 Trifilar Pendulum 
   2.2 Rigid Body Mode Mass Line Method 
3. Measurement Considerations 
   3.1 Test Object Suspension 
   3.2 FRF Measurements 
   3.3 Accurate Coordinates of Measurement Locations 
4. Simcenter Testlab Rigid Body Calculator Software 
   4.1 Getting Started 
   4.2 FRF Data and Modal Rigid Body Workbook 
   4.3 Band for Mass Line 
   4.4 Data Quality Check 
   4.5 Calculate Moment of Inertia 
   4.6 Simcenter Testlab Output Illustrated 
5. Conclusion 


1. Moment of Inertia Background 

The inertia of a point mass is given by the equation below (Equation 1): 
 
User-added image
Equation for inertia (I) of mass (m) about distance (r).

The inertia (I) is equal to the mass (m) of the object times the radius (r) squared.  The radius is the distance about which the object is rotating. Think of a weight being twirled on a string!  After being set in motion, the weight will want to keep going, even if it is suddenly accelerated in a different direction.
Inertia is often expressed in units of kilogram times meters squared (kg x m^2).

Many objects have complex mass distributions.  The inertia of a three-dimensional complex object requires more math than shown in Equation 1.  To determine the inertia of complex object, it is not as simple as multiplying a single mass value times the square of a specific distance.  Instead, the mass is distributed: the inertia varies depending on the axis of rotation.  For example, the inertia is different about a vertical rotation axis versus a lateral rotation axis, etc.

Fortunately, the inertia values of a three-dimensional complex object can be described by an ellipsoid shape as shown in Figure 2.
 
User-added image
Figure 2: Any physical object (left) has inertial properties that can be described by an ellipsoid shape (right). 

The ellipsoid describes how much inertia is present when the object is rotated in a specific direction about a point in space. Due to distribution of mass within the object (in this case a motor), the inertia is different in any given direction of rotation. 

The ellipsoid can be described from any frame of reference. However, the math needed to describe the ellipsoid is made simpler (has less terms) by judiciously selecting the origin and orientation of the axis system.  

To illustrate how the origin and axes can simplify the math to describe the inertia of an object, an example using a two-dimensional plane cut from the ellipsoid is shown in Figure 3.
 
User-added image
Figure 3: Using coordinate system (XY in left figure, red) with a specific origin and reference axes, the ellipsoid (blue) is difficult to explain mathematically.  Switching to a different coordinate system (12 in right figure, magenta), the ellipsoid is easier to describe mathematically.

If using a coordinate system with origin and reference axes (XY - red) as shown in the left figure above, the ellipsoid (blue) is more complex to explain mathematically.  This requires more mathematical terms (including something called the products of inertia) to describe. If the object rotates around this reference origin (red dot), it would not be in “balance”, meaning it would rotate and move away from the origin.

If the center of gravity of the object is used as the origin, the axes (12 - magenta) can be selected so there are no products of inertia required to describe the ellipsoid of inertia (figure on right above).  In this situation, if the object rotates around the center of gravity (black dot) origin, it is in “balance”.  This means that the object will continue to rotate around the origin.  In this situation, the ellipsoid can be described with three “principal moments of inertia”.

The ultimate output of Simcenter Testlab Rigid Body Calculator is the coordinates of the center of gravity and the three principal moments of inertia (and their orientation) for the test object.

2. Experimental Methods for Determination of Inertia

An overview of methods used to experimentally determine the inertial properties of an object are reviewed below.  The Simcenter Testlab software uses the mass line method.

2.1 Trifilar Pendulum

Traditionally, moments of inertia are measured using a pendulum to find the ellipsoid of inertia. The moment of inertia is a measurement of the resistance an object has to acceleration about a given axis. 

This is done by hanging a component on a pendulum and measuring the period it takes to complete one cycle as shown in Figure 4
 
User-added image
Figure 4: Test objects are hung from a rotating pendulum (left) in a trifilar pendulum inertia test, while in a mass line method inertia test (right) the test object is suspended on soft springs.

The trifilar pendulum equation can be found below (Equation 2): 
 
User-added image
Equation 2: Inertia from trifilar pendulum test.
 
Where:
   I = the moment of inertia of the bob (test component and pendulum)
   W = the weight of the pendulum fixture with test component
   R = the perpendicular distance from the axis of rotation to the center line of any wire
   T = the period of one full cycle
   L = length of wire from supporting frame of pendulum to bottom of pendulum

The object is hung so that its center of gravity is in the center of the pendulum. The time (T) is directly proportional to the inertial resistance of the test component.

While this is an accurate way of getting moment of inertia properties, this test can be highly time consuming and prone to calculation errors. Every axis of rotation must be measured separately, which can prove difficult in large, complex components. The support frames also need to be changed between measurements, and can lead to long, time-intensive experiments. As such, newer methods that can be quickly implemented have been designed to get around these problems. 

2.2 Rigid Body Mode and Mass Line Method

A mode shape is the motion that an object makes when it is excited at a natural frequency. There are two different types of modes: 
  • Rigid body modes: Where the entire object is moving in unison along a specific degree of freedom. The 6 degrees of freedom for an object include the 3 translational degrees of freedom, and 3 rotational degrees of freedom. 
  • Flexible modes: Object deforms or flexes itself.  Bending and torsion are common examples of flexible modes.
In a flexible, or deformation mode, the component is moving and twisting at different phases. Flexible modes are therefore not used to find rigid body properties. All the mass of the object participates in rigid body modes in all directions, making them suitable for the calculation of the ellipsoid of inertia.

The image below shows the rigid body modes of a disc brake (Figure 5).
 
User-added imageUser-added imageUser-added image
User-added imageUser-added imageUser-added image
Figure 5: The rigid body modes of a disk brake.  The top row is translational modes, the bottom row is rotational modes.

The first row shows the translational degrees of freedom, and the bottom shows the rotational. These rigid body modes occur at low natural frequencies and when found, can be used to find the center of gravity and moment of inertia for the test object. The idea is similar to the pendulum; the acceleration towards the center axis of the system can be measured because, like a pendulum, the entire body is moving in phase together. This method provides an accurate alternative to the pendulum test, without needing extra equipment than what is used in experimental modal testing. 

The main principle for obtaining inertia properties from rigid body modes is utilizing its frequency response function (FRF). FRFs are the frequency response ratio of a measured response (displacement, velocity or acceleration), over an input (force or other). In vibration, typically the acceleration is measured using an accelerometer and the corresponding frequency response function is acceleration over force, also known as accelerance. 

A full explanation of FRFs can be found in greater detail here: What is a Frequency Response Function (FRF)?

Frequency response curves show the expected response of single or multiple degree of freedom systems at a particular excitation frequency. The formula for a single degree of freedom system measuring acceleration over force can be written as shown in Equation 3:
 
User-added image
Equation 3: Accelerance equation 
 
Where A is acceleration, F is force, ω is the frequency in radians/second, k is the stiffness in N/m, c is the damping value in N/(m/s), and m is the mass in kg. 

Below shows an example of an FRF measuring a single degree of freedom system (Figure 6):
 
User-added image
Figure 6: The response of a dynamic system above its mode of vibration is proportional to the inverse of the mass.
 
The yellow zone shows the area of the frequency response curve that is controlled by the stiffness of the system, the orange zone shows the portion controlled by damping, and the blue zone shows the mass-controlled portion of the FRF. The mass-controlled zone is what is most important concerning rigid body modes. 

More about the mass line in a FRF in this knowledge article: Dynamic Stiffness, Compliance, Mobility, and more...

These mass lines correspond to 1/m, where m is the mass of the system in the measured direction. For example, if a SDOF had a weight of 10 kg, the value of the mass line would be 0.1. In the system above, the mass of the system was 25 kg, therefore, the flat portion of the mass line is at 0.04. In the mass line method, a user selects the frequency range of the FRF where the mass line is flat and unaffected by other modes (more guidance on selecting this region is offered further ahead in this document). 

Figure 7 below shows the highest frequency rigid body mode with a mass line that separates it from the first deformation mode. 
 
User-added image
Figure 7: With adequate frequency separation between the highest rigid body mode and lowest flexible mode, the FRF should contain a flat frequency band for determination of the mass from rigid body mode.

To use mass line methods to determine the ellipsoid of inertia of an object, the following is needed:
  • Geometry of the measurement locations on test object. Should “box in” the object and have accurately measured coordinates.
  • Total mass of the test object.
  • FRFs with flat mass lines at frequencies higher than the rigid body modes acquired at the measurement locations.
A good reference for how these calculations are performed is “Calculation of Rigid Body Properties from FRF Data: Practical Implementation and Cases” from the International Modal Analysis Conference (IMAC) proceedings of 1997.

Also see the document attached to this article.

3. Measurement Considerations

Some considerations when performing Frequency Response Function (FRF) measurements for moment of inertia calculations.

3.1 Test Object Suspension

The test object should be suspended in a manner that simulates free-free boundary conditions (or as much as possible). An example using soft bungee cords is shown in Figure 8.
 
User-added image
Figure 8: Powertrain test object suspended with soft bungee cords to simulate free-free boundary conditions.

In a free-free boundary condition, the test object is not affected by or attached to any other object. In other words, it is as if the test object is “floating in space”. In other words, it is suspended in a manner that does not restrict translational motion or rotational motion around any axis. In a dynamic simulation model, a free-free boundary condition results in six rigid body modes at zero Hertz.

Because a perfect free-free boundary condition is not possible in practice, there are rigid body modes in experimental conditions, but they are not located at zero Hertz. The rigid body modes are often located at frequencies just above zero Hertz.

For the rigid body calculation based on the FRFs to work properly, a flat mass line of the FRF is required. This means that there must be a large frequency separation between the flexible modes of the test object (torsion, bending, etc.) and the rigid body modes as shown in Figure 9.
 
User-added image
Figure 9: The first flexible mode of the test object should have a frequency that is ten times (10x) greater than the highest rigid body mode.
 
There are multiple options to approximate a free-free suspension experimentally (bungee cords, elastomeric mounts, foam, etc.). Whatever method is used for the suspension, the resulting FRF measurements need to be evaluated to ensure there is adequate separation between rigid and flexible modes in the structure and that the mass line of the rigid body modes of the FRF are flat.

3.2 FRF Measurements

To get the best results, FRFs need to be acquired “all around” the test object. It helps to think of creating a box around the test object as shown in Figure 10.  

User-added image
Figure 10: Minimal set of FRF measurement locations for a rigid body analysis.
 
While boxing in the object, it is also good practice to avoid the centerlines of any axis. 

In the above example, a set of FRFs will be measured by applying an excitation at all eight locations in multiple directions. If using a modal impact hammer, this would mean a possible 24 separate acquisitions could be acquired. Ideally, at minimum of nine FRF sets will be acquired with multiple directions being measured at some input nodes. Acquiring greater than nine will yield more possibilities during post analysis and is recommended.  The key is to excite all six rigid body modes (three translation, three rotation) with the inputs.

Some useful links on acquiring FRFs:

 3.3 Accurate Coordinates of Measurement Locations

The three-dimensional coordinates of the measurement's locations are used in the calculation of the moments of inertia. The orientations of the accelerometers to the axes system are also used. Because the locations and orientations are used in the calculations, they must be measured as accurately as possible.

A two-dimensional laser level (Figure 11) can be helpful in determining the exact coordinates of the measurement locations and orientations.
 
User-added image
Figure 11: Two dimensional laser level.
 
The laser site is placed next to the test object. It projects straight laser lines on the test object to make measuring distances as easy as possible in a precise manner.  This can also be used to make sure the accelerometers axes are oriented as expected. 

If a Computer Aided Drafting (CAD) model is available of the test object, coordinates could also be used from it.

4. Simcenter Testlab Rigid Body Calculator Software

To use the Simcenter Testlab Rigid Body Calculator, a Simcenter Testlab project containing the following is required (as described previously):
  • Set of Frequency Response Functions (FRFs) on a softly suspended object needs to have been measured as shown in Figure 12.  
 
User-added image
Figure 12: A set of Frequency Response Function (FRFs) for moment of inertia determination.
 
A minimum of 24 output locations and 9 input locations for the FRFs are recommended.
  • Geometry with precisely measured locations of the FRFs (Figure 13).
User-added image
Figure 13: Schematic (left) of test object used in this example. Coordinates are entered and stored in the Simcenter Testlab geometry (right) in the project file.
 
A good reference for the creation of geometry in Simcenter Testlab: Simcenter Testlab Geometry
  • Total mass of the test object.
The Simcenter Testlab Rigid Body Calculator is used to analyze the FRFs to determine the center of gravity and moments of inertia.

4.1 Getting Started

Start Simcenter Testlab Desktop and under “Tools -> Add-ins” turn on “Modal Analysis” and “Rigid Body Calculator” as shown in Figure 14.
 
User-added image
Figure 14: Calculating moments of inertia from FRFs requires two add-ins in Simcenter Testlab: Modal Analysis and Rigid Body Calculator.
 
If using token licensing, the two add-ins require 114 tokens total: 77 tokens for Modal Analysis and 37 tokens for Rigid Body Calculator.  

More about Simcenter Testlab token licensing in the knowledge article: Simcenter Testlab Tokens: What are they, and how do they work?

4.2 FRF Data and Modal Rigid Body Workbook

After turning on the appropriate add-ins, activate the section that contains the FRF measurements in Figure 15.
 
User-added image
Figure 15: Activate the section with the measured FRFs and then click on “Modal Data Selection” workbook. 

Then click on the “Modal Data Selection” workbook.

In the “Modal Data Selection” workbook allows quick viewing of the FRFs to be used in the analysis (Figure 16).
 
User-added image
Figure 16: View FRFs to make sure all desired are selected and then go to the Modal Rigid Body workbook.
 
On the middle left of the screen, green colored cells should appear. With the “Selected FRFs” checkbox on, FRFs can be interactively viewed by selecting cells.

When ready, click on the “Modal Rigid Body” workbook.

4.3 Band for Mass Line

In the “Modal Rigid Body” workbook, choose the “Change Modal Data Selection…” button in the upper left of the screen. Make sure “Active Section” is selected. Again, there should be green cells indicating FRFs are available for viewing.

Then, in the upper right of the screen, adjust the cursors to select a flat portion of the FRF curve. If necessary, zoom in to select the area between the rigid body modes and flexible modes as shown in Figure 17.

User-added image
Figure 17: After selecting some FRFs to be shown (upper left), move the cursors to select a band where the FRFs are flat.

Ideally, the mass line area should be flat (as described elsewhere in this document). The number of data points in the band (i.e., spectral lines) is also shown. The more lines in the flat area, the better the calculation.

By default, the “Mass line method” setting (middle left in software interface) is set to “Unchanged” which assumes the FRFs indeed have a flat mass line. This is always the preferred method.

However, there are situations where a flat mass line can be difficult to realize: modern aircraft, for example, have the first flexible mode of the wings well below one Hertz.  Designing a suspension system which would allow an adequate separation between rigid and flexible modes is expensive, if feasible at all. In this situation, dedicated methods have been developed to compensate for the lack of a flat mass line.  These correction methods require a modal analysis with lower and upper residual calculations:
  • Corrected FRFs: When insufficient bandwidth exists between the rigid body modes and the flexible modes, the measured FRF’s must be corrected.  First the flexible modes are estimated from measured FRF’s using a modal curvefitter. Then FRF’s are synthesized from these modes (without lower residual) in the selected frequency band and subtracted from the measured FRF.   
  • Lower residual: In the situation where accurate measured FRF’s are not available in the frequency range between the rigid body and the flexible modes, lower residual terms can be used. Lower residual terms are estimated from the measured FRF’s using a modal curvefitter.   The lower residual represents the influence of the rigid body modes and is used to approximate the mass line.
Next some automated checks on data quality can be run.

4.4 Data Quality Check

Press the “Check” button in lower left to perform an automated assessment of the data quality to be used for the moment of inertia calculation (Figure 18).
 
User-added image
Figure 18: Press “Check” in the lower left to initialize the automated data check. Then press “Logging” to see the results of the check.
 
Checks on the data include:
  • The recommended minimum of inputs and outputs for the FRFs are met.
  • Data quality expressed in correlation percentage (100% being the best).
In this example, the data check is “orange” indicating that an improvement in the FRF data quality (69.02%) is desirable. 
The FRFs in question and their correlation can be identified by pressing the “List” button and scrolling left to right to find the low percentage FRF as shown in Figure 19 below.
 
User-added image
Figure 19: In “Modal Data Selection” right click and exclude the previously identified low percentage FRFs.
 
The excluded FRFs will turn red to indicate they are removed from the analysis.

Returning to the “Modal Rigid Body” workbook, the data quality check can be run again as shown in Figure 20.
 
User-added image
Figure 20: Performing the automated data quality check should ideally result in all criteria being marked green. 

After excluding the problem data, the status of the quality check should turn green, indicating that all criteria were met. It is helpful to acquire more than the minimum nine inputs in case data needs exclusion.

Now the Moments of Inertia can be calculated.

4.5 Calculate Moment of Inertia

To calculate the moments of inertia from the FRF data, click on the “Calculate” minor sheet (located near the top) of the Modal Rigid Body workbook (Figure 21).
 
User-added image
Figure 21: With the data quality status green and having entered the mass of the test object, the moments of inertia can be calculated.

In the “Calculate” minor workbook, the “Calculate” button (middle left) can be pressed. The moments of inertia will be shown in the upper right section of the screen, and the center of gravity location will be shown on the geometry in the lower right.

The output includes:
  • Coordinates of center of gravity (cog)
  • Moments and products of inertia about the center of gravity and user defined reference point
  • Principal moments of inertia and their direction
  • Synthesis of six scaled rigid body modes with user defined frequency and damping for use in simulation models
The calculated results can be copied to Microsoft Excel. Highlight the cells, right click, choose “Copy” as shown in Figure 22.
 
User-added image
Figure 22: Rigid Body Calculator results can be copied to Excel by highlighting all cells (upper left corner) and selecting “Copy”.
 
Open Excel and then choose “Paste”.

In the “Advanced...” button of the workbook, it is possible to calculate the inertia with respect to any origin.

4.6 Simcenter Testlab Output Illustrated

There are several items output in the Rigid Body Calculator results.  Some of the results are oriented to the original measurement coordinates, while others are oriented to the principal axes.

The origin used to describe the original measurement locations is listed (also called the “reference”), and the location of the calculated center of gravity (Figure 23).
 
User-added image
 Figure 23: Measurement reference origin (blue dot) and center of gravity (black dot).

The reference coordinate system are the coordinates used to do the measurements.  Because the location of the center of gravity and the orientation of the principal axes were not known when beginning the measurements, there is a good chance that the two coordinate systems (reference and principal) will not align.

The moment of inertias about the reference origin are also listed as shown in Figure 24:
 
User-added image
Figure 24: Moments of inertia around the reference origin and reference axes (light blue).

Because these moments (Ixx, Iyy, Izz) are not calculated around the center of gravity and in the principal axes, there are products of inertia (Ixy, Ixz, Iyz). As mentioned at the beginning of the document, this means that if the object were sent into a rotation around the reference origin, it would move away from the reference origin due to imbalance. In practice, it is even possible that the reference origin could be located outside the ellipsoid of inertia.

The moments of inertia are also calculated around the center of gravity (cog) and in the original reference axes (Figure 25):
User-added image
Figure 25: Moments of inertia around the center of gravity (cog) and reference axes (green).

Just like the moments of inertia (Ixx, Iyy, Izz) calculated around the reference origin, these also have products of inertia (Ixy, Ixz, Iyz).

Finally, the principal moment of inertia around the center of gravity and in the principal axis system is also calculated and displayed (Figure 26):
 
User-added image
Figure 26: Principal moments of inertia (red) around the center of gravity (cog) and principal axes.

The three principal moments (I11, I22, I33) are listed in numerical order, from highest to lowest.  There are no products of inertia because they have a value of zero when using the principal axis system.

Lastly, the rotation between the reference axis system orientation and the principal axis orientation is listed (Figure 27):
 
User-added image
 Figure 27: Translations and rotations between reference and principal axes systems (magenta).

The rotation is expressed two different ways: angles (to be performed in order from left to right) and as a rotation matrix.

5. Conclusion

Experimental determination of the moments of inertia can be useful in verifying a simulation model.  To get the proper motion (whether vibration or other) in a model, the mass properties of each component must be correct.  This can be particularly useful when the Computer Aided Drawing (CAD) model of some components do not exist, but physical parts are available.  These components can be represented by their experimentally verified mass inertial properties.

Hope this background and tutorial helps!  Questions? Email steff.nelson@siemens.com

 
Related Links

KB Article ID# KB000049612_EN_US

Contents

SummaryDetails

Associated Components

SCADAS RS SCADAS System Simcenter Testing Solutions Simcenter Testlab Simcenter Testlab Neo