5.5. Steps for Performing a Substructured Multibody Simulation

The methodology for performing a substructured multibody simulation assumes that you have generated the entire finite element model of the multibodies including the joints—using Ansys Workbench, for example—and want to take advantage of substructuring to reduce the solution time. This method is referred to as a top-down approach (as opposed to a bottom-up approach of defining the substructure first and then building the rest of the model around it).

Before proceeding, prepare the full multibody model (as described in Steps 1 through 4 in Overview of the Multibody Analysis Process). Verify that the bodies are connected to the joints as described in Connecting Bodies to Joints.

The multiple passes used in substructuring require that the files created and used in the process are handled appropriately. To aid in file management when performing a substructured multibody simulation, use the /FILNAME command to modify the current jobname as needed.

5.5.1. Step 1: Prepare the Full Model for a Substructured Multibody Analysis

Prepare the full model for a substructured multibody analysis, as follows:

StepActionCommentsCommand(s)
1.1Specify the full jobname.Example: /FILNAME,FULL /FILNAME
1.2Resume (or build) the full model.See Overview of the Multibody Analysis Process and Connecting Bodies to Joints. RESUME
1.3Make components of the flexible body. (Repeat for each flexible body.)Create an element component of the elements of the body, including any contact elements used to connect the body to a joint. (Do not include the joint elements.)
ESEL
CM,Ename,ELEM
1.4Select the entire model.--- ALLSEL
1.5Save the model.--- SAVE

5.5.2. Step 2: Create the Substructures (Generation Pass)

Perform the generation pass to create the CMS substructure (in the matrix .SUB file) characterizing the dynamic flexibility of the body.

You must decide how many modes to include in the CMS substructure. The number you determine depends on several factors including:

  • The driving frequency.

  • The frequencies to be excited (such as flexural, axial, torsional, etc.).

  • Whether displacements are of primary interest, or whether stresses/strains are of primary interest. (The latter require more modes to accurately capture their response.)

  • Whether impact (contact) is included. (Impact tends to excite higher frequencies.)

  • Whether acoustic frequencies are desired.

For most analyses, and particularly for rotating bodies, the fixed-interface method (CMSOPT,FIX) is sufficient. For analyses where higher frequencies are of interest (foe example, those involving acoustics or high-speed equipment), the residual-flexible free-interface method (CMSOPT,RFFB) provides more accuracy. For more information, see Supported CMS Methods in the Substructuring Analysis Guide.

For nonrotating bodies, you can apply constraints (D) in the generation pass to the degrees of freedom, but not the master degree of freedom. Set KEYOPT(4) = 1 for these superelements in the use pass; otherwise, your analysis will have convergence problems. For rotating bodies, do not apply constraints in the generation pass because the superelement must have six rigid body modes; you can, however, apply constraints to its master degree of freedom in the use pass.

Loading Considerations

When applying loads, be aware that:

  • The loads rotate with the rotating substructure by default. This behavior is valid for most load types (especially pressure loads). In the use pass, however, you can specify that the load vector not rotate with the substructure; disabling load rotation is useful in some cases, such as those involving nodal forces where you want to maintain their original direction.

  • When to apply gravity and other acceleration loads (such as those applied via ACEL and OMEGA commands) depends on whether the body is rotating or not. For a rotating body, apply the loads in the use pass. For a nonrotating body, you can apply the loads in this step and use it in the use pass; however, be careful not to specify it twice (for example, by issuing an ACEL command in the use pass). Issue the CMACEL command to apply the acceleration to the nonsubstructured elements only.

  • By applying a unit load in this step, you can easily scale it in the use pass and make use of tabular loads to apply a complex load-versus-time history in a single load step. Ansys, Inc. recommends this approach as it allows for straightforward creation of the full model results file.

Creating the Superelements

Follow these steps to create the superelements for a substructured multibody analysis:

StepActionCommentsCommand(s)
2.1Clear the database.Required only if performing this step in the same session as the prior step. /CLEAR
2.2Specify the generation pass jobname.Example: /FILNAME,BODY1 /FILNAME
2.3Resume the full model.Example: RESUME,FULL.db RESUME
2.4Define the analysis type.The analysis type is substructure.
/SOLU
ANTYPE,SUBSTR
2.5Define substructure options.Substructure name, and generate stiffness and mass, as in this example: SEOPT,BODY1SE,2SEOPT,Sename,2
CMS options, including the number of modes.CMSOPT,FIX,NMODE
2.6Select the substructure nodes and elements.Select the elements defined in Step 1.3.CMSEL,S,ELEM
Select the interface nodes defined in Step 1.3CMSEL,S,NODE
Create master degrees of freedom at all selected nodes.M,ALL,ALL
Select the nodes attached to the elements. NSLE
2.7Apply loads, if any.These are loads typically interior to the body (that is, not applied to a master degree of freedom).
F
SF
SFE
ACEL
...
2.8Create the substructure.Save the model. SAVE
Execute the creation. SOLVE

Repeat the steps above for each flexible body you wish to replace with CMS substructures. Use unique jobnames and substructure names for each flexible body.

Residual-Flexible Free-Interface CMS Method

If you are using the residual-flexible free interface method, use CMSOPT,RFFB,NMODE (rather than CMSOPT,FIX,NMODE) in Step 2.5. You must also define pseudo-constraints (D,,,SUPPORT).

For further information, see The CMS Generation Pass in the Substructuring Analysis Guide.

5.5.3. Step 3: Build the CMS-based Model (Use Pass)

Replace the flexible bodies with their corresponding CMS substructures.

StepActionCommentsCommand(s)
3.1Clear the database.Required only if performing this step in the same session as the prior step. /CLEAR
3.2Specify the use pass jobname.Example: /FILNAME,USE /FILNAME
3.3Resume the full model.Example: RESUME,FULL.db RESUME
3.4Replace the flexible bodies.Deselect the flexible elements.
/PREP7
CMSEL,U,Ename
Define the substructure element type using an available type number (ITYPE).ET,ITYPE,50
If any loads were applied in Step 2 and you do not want them to rotate with the substructure, set the appropriate key option.KEYOPT, ITYPE,3,1
For nonrotating substructures that have constraints applied in the generation pass, set the appropriate key option. KEYOPT,ITYPE,4,1
Define the substructure.
TYPE,ITYPE
SE,Sename

Repeat the CMSEL,U and SE commands for each flexible body.


Caution:  Be careful not to select all elements (for example, via an ALLSEL command) before initiating the solution (SOLVE) in the next step. If you do so, Mechanical APDL solves for both sets of elements.


5.5.4. Step 4: Run the Multibody Analysis

Set up the multibody analysis and run it.

StepActionCommentsCommand(s)
4.1Specify the analysis type.Large deflection, transient analysis (multibody analysis).
/SOLU
ANTYPE,TRANS
4.2Specify the transient analysis options.HHT method with 0.1 numerical damping.
TRNOPT,FULL,,,,,HHT
TINTP,0.1
4.3Specify boundary conditions.Constraints on motion and initial conditions
D
DJ
IC
Applied loads, including applying loads from the generation pass Step 2.7 (SFE,,,SELV)
F
FJ
ACEL
SF
SFE
...
4.4Specify load step options and solve.Ending time and time step sizes.
Results file output controls. OUTRES
Run the analysis. SOLVE

To dampen out excessive solution noise, particularly in the velocities and accelerations, you typically use numerical damping. For more information, see Damping.

In Step 4.3, use tabular loads to specify complex load-versus-time histories. By default, loads are simply ramped (or step-applied (KBC)) over the time interval from one load step to the next. Tabular loads, however, allow a general load curve. To use multiple load steps to define the loading, repeat Steps 4.3 and 4.4 for each load configuration.

For more information about setting up and performing a multibody analysis, see Performing a Multibody Analysis.

5.5.5. Step 5: Expand all Solutions (Expansion Pass)

Using the solutions from the prior step (displacements at the master degrees of freedom at each time point), obtain the displacements and stresses (if desired) for all nodes and elements of the flexible bodies.

StepActionCommentsCommand(s)
5.1Clear the database.Required only if performing this step in the same session as the prior step. /CLEAR
5.2Specify the generation pass jobname from Step 2.Example: /FILNAME,BODY1 /FILNAME
5.3Resume that jobname's database.No file name required. RESUME
5.4Specify an expansion pass.---
5.5Specify the substructure to expand.Substructure name and the use pass jobname from Step 3. Example: SEEXP,BODY1SE,USESEEXP,Sename,Usefil
5.6Specify the solutions to expand, then expandExpand all time points, and indicate whether or not to compute stresses, strains, and forces.NUMEXP,ALL,,,Elcalc
Perform the expansion. SOLVE

Repeat all steps for each substructured body (including clearing the database (/CLEAR)).

5.5.6. Step 6: Create the Merged Results File

Merge all results files (one from the use pass and one from each of the expanded substructures) to create a results file with the full model data. After completing this part of the process, you can perform postprocessing as though you had run the full model in the multibody simulation.

StepActionCommentsCommand(s)
6.1Clear the database.Required only if performing this step in the same session as the prior step. /CLEAR
6.2Specify the full model jobname from Step 1.Example: /FILNAME,FULL /FILNAME
6.3Resume that jobname's database.No file name required. RESUME
6.4Delete the merged results file.If you fail to delete the merged results file, the program appends the results from this step to that file. /DELETE
6.5Merge the results for each time point.Loop through each time point (solution substeps).
/POST1
*DO,J,1,NSUBSTEPS
Bring in the use pass results.
FILE,USE
SET,1,J

Append the expanded substructure results.

Repeat both of these commands for each substructure.

FILE,BODY1
APPEND,1,J
Write the combined results and loop back for the next time point.

Understanding the example commands in this step:

  • NSUBSTEPS is the total number of substeps (time points) in the results files.

  • In the example commands, the jobname from the use pass (Step 3) is USE; therefore, its results file is named USE.rst. Likewise, the jobname from the expansion pass (Step 5) is BODY1; therefore, its results file is named BODY1.rst. Adjust the command arguments accordingly to accommodate your own jobnames.

  • As presented here, the analysis in the use pass is performed in one load step with NSUBSTEPS substeps. If such is not the case in your analysis, modify the *DO loop to use the appropriate SET command.

  • The expansion pass results files always have only one load step with all time points contained as NSUBSTEPS substeps, irrespective of the use pass load stepping and substepping.

5.5.7. Step 7: Postprocess the Results

Postprocess the full model as though you had run a nonsubstructured analysis.

Use the POST1 postprocessor (/POST1) to review the results over the entire model. Use the POST26 postprocessor (/POST26) to obtain time-history listings and plots. For more information, see Reviewing Multibody Analysis Results for specific multibody postprocessing.

StepActionCommentsCommand(s)
7.1Specify the full model jobname from Step 1.Example: /FILNAME,FULL /FILNAME
7.2Resume that jobname's database.No file name required. RESUME
7.3Review results at a specific point in time---
7.4Select the entire model.---

/POST26

Nodal velocity and acceleration nodal results are not available for the substructure interior nodes (non-master degrees of freedom).