Analysis of a frame using substructures

This example demonstrates and verifies the substructuring capabilities in Abaqus.

This page discusses:

ProductsAbaqus/Standard

Substructures provide significant savings in computer time in certain cases. They can be used in a linear analysis to divide a problem into several smaller analyses or to save parts of the structure that are not modified every time the analysis is done. Substructures can also be useful in a nonlinear analysis. As demonstrated in The Hertz contact problem, if any part of the structure remains linear during the analysis, the stiffness of that part can be calculated once and saved. Substructures can also be used in the analysis of linear perturbations around a state reached as a result of nonlinear response, as shown in Vibration of a rotating cantilever plate, where the vibration frequencies are calculated for a prestressed structure.

By default, Abaqus writes the substructures to the files whose name is defined by the job name and the step number. Alternatively, you can define the substructure name explicitly in the substructure generation analysis and when defining elements. The input files for this example show the use of these options.

Problem description

The example is a frame consisting of two columns connected by a beam (Figure 1). The columns and the beam are each modeled with 10 elements of type B21. The frame is loaded by a distributed load applied on the beam.

One-level substructure analysis

Once a substructure has been generated, it can be used just as any other element in an analysis. The substructure is connected to the rest of the model through specified degrees of freedom (the retained degrees of freedom). During substructure generation it is also possible to define load cases for the substructure. These load cases can then be scaled and applied during the analysis. In this example the loading on the beam is defined as a substructure load case and then applied at the global level.

Many structures include parts with the same geometry. In such cases the part can be generated as a single substructure that is then used several times. In this case the two columns are identical, so a substructure can be formed once for the column geometry and then used twice.

The analysis includes two steps. The first is a static load case in which the distributed load is applied to the structure. For this step the substructure analysis gives the same solution that is obtained by using the model without substructures, since the static reduction to a substructure stiffness matrix does not change that matrix.

In the second step a frequency analysis is performed to extract the first few vibrational modes of the structure. The technique of generating a substructure so that its dynamic response is modeled with reasonable accuracy is called “Component Mode Synthesis.” In this approach the dynamic response of the substructure is made up from a combination of its static modes of deformation and some of its natural modes of vibration obtained with its retained degrees of freedom fully constrained. In the case of a substructure with many internal degrees of freedom it may be expensive to extract these internal modes of vibration. In that case a sufficiently accurate dynamic representation of the substructure can often be obtained by “Guyan reduction.” This means that additional degrees of freedom that are not needed for the physical connection of the substructure to the rest of the mesh are retained simply to improve the dynamic modeling. The drawback to this method is that a judicious choice of these degrees of freedom is necessary to get the most benefit from the larger size of the substructure matrix. Component mode synthesis is, thus, a more reliable alternative but requires a frequency analysis of the substructure, which can be expensive for a large model. The modes are obtained by extracting some of the eigenfrequencies of the substructure in a frequency step, with all of the retained degrees of freedom restrained with a boundary condition. The n modes to be used in the substructure generation are specified with n not greater than the number of modes that have been extracted. These extra modes are then added when the substructure is used in the global analysis. In general, this technique improves the representation of the substructure's dynamic behavior significantly, even with just a few additional modes. In this example we have added two extra modes in each of the substructures. substrframe_1level_2mode_gen1.inp and substrframe_1level_2mode_gen2.inp show the input data when extra dynamic modes are added.

Multilevel substructure analysis

Substructures can be used within substructures. We illustrate this capability by dividing each column into two identical parts. A substructure is then generated for the half-column, and a higher-level substructure is then generated for the full column by combining two substructures on the lowest level. The same process is followed for the beam. In this way a total of nine substructures are used to model the structure, although the repetitive structure requires very little calculation. Abaqus places no limit on the number of levels of substructures that can be used. substrframe_2level_gen1.inp, substrframe_2level_gen2.inp, substrframe_2level_gen3.inp, and substrframe_2level_gen4.inp show the input data for this example with two levels of substructures.

Substructure rotation and mirroring

Substructure stiffness and mass matrices are created based on the location and orientation of the elements when the substructure is generated. These matrices can be used anywhere in space as long as they are only translated. If the substructure is needed in a rotated configuration, the matrices must be transformed. Rotation and mirroring of a substructure is defined using a substructure property. The example illustrates this option with the three-dimensional frame structure shown in Figure 2. Element 82 on the usage level has been created by rotating substructure Z201 by 90°. Element 81 on the usage level has been created by mirroring substructure Z201 through the plane x=-y. Element 72 on the global level could have been created by a simple translation of substructure Z201 but has instead been created by both a rotation and a mirroring of the substructure. substrframe_b31_gen1.inp and substrframe_b31_gen2.inp show the input data for this analysis.

Results and discussion

The analysis without substructures provides “correct” results for this example (that is, the results are correct for the finite element discretization used). For the static loading case, substructure analysis does not imply any approximation; and the results with substructures are, therefore, identical. For the frequency analysis the substructure mass matrix is an approximate representation, so some accuracy is lost. Table 1 and Table 2 show comparisons of the lowest frequencies estimated by the planar and three-dimensional models, respectively. They show the significant improvement in results that is obtained by including extra modes.

Input files

substrframe_1level_2modes.inp

One level of substructures and two dynamic modes added for each substructure.

substrframe_1level_2mode_gen1.inp

Substructure generation referenced by the analysis substrframe_1level_2modes.inp.

substrframe_1level_2mode_gen2.inp

Substructure generation referenced by the analysis substrframe_1level_2modes.inp.

substrframe_2level.inp

Two levels of substructures.

substrframe_2level_gen1.inp

Substructure generation referenced by the analysis substrframe_2level.inp.

substrframe_2level_gen2.inp

Substructure generation referenced by the analysis substrframe_2level.inp.

substrframe_2level_gen3.inp

Substructure generation referenced by the analysis substrframe_2level.inp.

substrframe_2level_gen4.inp

Substructure generation referenced by the analysis substrframe_2level.inp.

substrframe_b31.inp

Analysis of the three-dimensional frame structure using substructures. In this case element 82 is created by a rotation, element 81 by mirroring in the plane x=-y, and element 72 by rotation followed by mirroring.

substrframe_b31_gen1.inp

Substructure generation referenced by the analysis substrframe_b31.inp.

substrframe_b31_gen2.inp

Substructure generation referenced by the analysis substrframe_b31.inp.

substrframe_nosubstr.inp

Analysis without substructures.

substrframe_1level.inp

Analysis with one level of substructures.

substrframe_1level_gen1.inp

Substructure generation referenced by the analysis substrframe_1level.inp.

substrframe_substrpath.inp

Shows the use of the two substructure files generated earlier in redoing the global analysis of the frame.

substrframe_3d_nosubstr.inp

Analysis of the three-dimensional frame structure without substructures.

substrframe_rotz201.inp

Analysis of the three-dimensional frame structure using substructures. The two structures in the y-direction have been created by rotating substructure Z201 by 90° about the z-axis.

substrframe_rotz201_gen1.inp

Substructure generation referenced by the analysis substrframe_rotz201.inp.

substrframe_rotz201_gen2.inp

Substructure generation referenced by the analysis substrframe_rotz201.inp.

substrframe_rotz201_2modes.inp

Same as substrframe_rotz201.inp with two extra dynamic modes added for each substructure.

substrframe_rotz201_2mode_gen1.inp

Substructure generation referenced by the analysis substrframe_rotz201_2modes.inp.

substrframe_rotz201_2mode_gen2.inp

Substructure generation referenced by the analysis substrframe_rotz201_2modes.inp.

Tables

Table 1. Frequencies of the first three modes, planar model (all frequencies in cycles/time).
  Mode 1 Mode 2 Mode 3
Without substructures 7.0176 7.2431 20.309
Three substructures 7.2609 11.640 33.291
Error 3.5% 61% 64%
Three substructures, with two dynamic modes 7.0173 7.2444 20.314
Error 0.004% 0.004% 0.004%
Two levels of 7.2609 11.640 33.291
substructures      
Error 3.5% 61% 64%
Table 2. Frequencies of the first three modes, three-dimensional model (all frequencies in cycles/time).
  Mode 1 Mode 2 Mode 3
Without substructures 4.6277 4.6277 4.6320
Eight substructures 5.0416 5.0416 5.4384
Error 8.9% 8.9% 17.4%
Eight substructures, with two dynamic modes 4.6301 4.6301 4.6321
Error 0.05% 0.05% 0.05%

Figures

Figure 1. Frame model.

Figure 2. Three-dimensional model.