1.16.1 Analysis of a frame using substructures

Product: ABAQUS/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, Section 1.1.11, 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, Section 1.4.7, where the vibration frequencies are calculated for a prestressed structure.

The purpose of this example is to provide basic demonstration and verification of the substructuring capabilities in ABAQUS.

Problem description

The example is a frame consisting of two columns connected by a beam (Figure 1.16.1–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 the *BOUNDARY option. The n modes to be used in the substructure generation are then specified using the *RETAINED EIGENMODES option, 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 library files

By default, ABAQUS will write the substructures to a file whose name is defined by the job name. It is also possible to define the library name on the *SUBSTRUCTURE GENERATE and *ELEMENT options. The contents of a substructure library can be listed by using the *SUBSTRUCTURE DIRECTORY option. It is also possible to copy substructures from one file to another using the *SUBSTRUCTURE COPY option. substrframe_substrdirect.inp shows the use of these options.

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 on the *SUBSTRUCTURE PROPERTY option. The example illustrates this option with the three-dimensional frame structure shown in Figure 1.16.1–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 . 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.16.1–1 and Table 1.16.1–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 using the *RETAINED EIGENMODES option.

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_substrdirect.inp

Shows the use of *SUBSTRUCTURE DIRECTORY and the consolidation of the three substructure files from substrframe_2level.inp into two substructure files.

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 , 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 in substrframe_substrdirect.inp 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 -direction have been created by rotating substructure Z201 by 90° about the -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.16.1–1 Frequencies of the first three modes, planar model (all frequencies in cycles/time).

 Mode 1Mode 2Mode 3
Without substructures7.01767.243120.309
Three substructures7.260911.64033.291
Error3.5%61%64%
Three substructures, with two dynamic modes7.01737.244420.314
Error0.004%0.004%0.004%
Two levels of7.260911.64033.291
substructures   
Error3.5%61%64%

Table 1.16.1–2 Frequencies of the first three modes, three-dimensional model (all frequencies in cycles/time).

 Mode 1Mode 2Mode 3
Without substructures4.62774.62774.6320
Eight substructures5.04165.04165.4384
Error8.9%8.9%17.4%
Eight substructures, with two dynamic modes4.63014.63014.6321
Error0.05%0.05%0.05%


Figures

Figure 1.16.1–1 Frame model.

Figure 1.16.1–2 Three-dimensional model.