3.19.1 Frequency extraction using the AMS eigensolver

Product: ABAQUS/Standard  

The tests in this section verify the frequency extraction procedure using the AMS eigensolver in ABAQUS/Standard by comparing the results with those obtained by the Lanczos eigensolver.

I. One-element tests

Elements tested

CPE4    C3D8   

Features tested

Eigenvalue extraction for a system with a symmetric stiffness matrix and multi-point constraints, selective modal recovery, full modal recover, and import.

Problem description

The two-dimensional model consists of a linear element of unit length. The nodes at one end (y = 0) are constrained, while the nodes at the other end are involved in a LINK MPC. The eigenvalue extraction is performed for the undeformed configuration. The three-dimensional model consists of a single linear element and is mainly used for testing the import feature.

Results and discussion

The eigenvalues obtained for both the AMS and Lanczos procedures are identical.

Input files

ams_1cpe4.inp

Eigenvalue extraction for a model with one element using the AMS eigensolver.

ams_import0.inp

Preloading of a single C3D8 element.

ams_import.inp

Frequency extraction of the import model using the AMS eigensolver.

II. Model with various lagrange-multiplier constraints (contact, connectors, distributing couplings)

Elements tested

C3D8I    C3D8R    C3D10M   

Features tested

Constraints with Lagrange multipliers and submodeling, mode-based steady-state dynamic restart, and selective modal recovery.

Problem description

The model consists of a semisphere pressed against a cube that is in contact with a rigid surface. The semisphere is also connected to the cube via four axial connectors.

In the preloading step the semisphere is pressed against the cube to establish contact. The load is applied at the reference node of the distributing coupling. In the second step the frequencies of the preloaded structure are extracted via the AMS procedure. Finally, the mode-based steady-state response is calculated in the third step using the results of the frequency extraction step. The results are compared with those obtained by the Lanczos eigensolver.

Results and discussion

In the following table the frequency extraction step results obtained by the Lanczos and AMS eigensolvers are compared.

ModeLanczosAMS
111.54711.551
211.91611.921
320.66420.690
425.79225.840
527.91627.963
628.80728.862
742.04842.110
842.37042.441

Input files

ams_conn_contact.inp

Full analysis using the AMS eigensolver.

ams_conn_contact_res.inp

Mode-based steady-state dynamic analysis restarted from the end of the frequency step.

ams_conn_contact_submodel.inp

Frequency extraction and mode-based steady-state dynamic analysis of a submodel driven entirely from the original model.

III. Model with coupled-temperature displacement

Elements tested

S8RT    B31H    B33H    B31    B33   

Features tested

Coupled temperature-displacement steps, hybrid Bernoulli and Timoshenko beams, full modal recovery, and mode-based steady-state dynamic analysis.

Problem description

The model consists of two rectangular parallel plates connected via beams at each corner. The structure is preloaded by applying a heat flux at the center of the top plate. The linear response is analyzed in a mode-based steady-state dynamic step preceded by a frequency extraction step using the AMS solver.

Results and discussion

In the following table the frequency extraction step results obtained by the Lanczos and AMS eigensolvers are compared.

ModeLanczosAMS
114.74314.745
214.74314.748
315.29615.301
417.15817.164
529.47629.505
638.68438.749
738.68438.773
852.77853.009
963.20163.545
1067.25367.621
1167.25367.641
1270.05570.555
1387.08087.166
1488.78989.594
1588.78989.720
1688.81890.292
1791.94692.735
1892.87793.825

Input file

ams_temp_plates.inp

Full analysis using the AMS eigensolver.

IV. Tire model with symmetric model generation and symmetric results transfer

Elements tested

CGAX3H    CGAX4H    SFMGAX1   

Features tested

Eigenvalue extraction for a tire model with hybrid and/or cylindrical elements, axisymmetric model followed by symmetric model generation with symmetric results transfer, and full modal recovery.

Problem description

The axisymmetric tire is inflated and then transferred to a full three-dimensional configuration. Subsequently, the rigid surface is brought in contact with the full tire, obtaining the footprint. Finally, the linear response is analyzed by performing a frequency extraction using the AMS eigensolver followed by a mode-based steady-state dynamic step.

Results and discussion

The following table shows the comparison of eigenfrequencies obtained by the Lanczos and AMS eigensolvers.

ModeLanczosAMS
147.55247.590
248.99249.042
354.39154.445
456.74956.795
577.58277.743
682.15382.265
785.12385.268
885.55385.694
998.55498.802
10103.73104.06
11112.37112.77
12116.90117.47
13118.64119.08
14119.71120.04
15124.68125.18
16130.75131.43
17132.16132.60
18136.05136.61
19137.41138.03
20138.30139.02
21140.35140.97
22140.58141.23
23143.88144.66
24144.98145.75
25148.05148.99
26152.60153.74

Input files

ams_tire_axisymm.inp

Axisymmetric tire model.

ams_tire_full3d.inp

Three-dimensional tire model.

V. Model with map solution

Element tested

CPS3   

Features tested

Solution mapping and selective modal recovery.

Problem description

The first model is subject to a static preload. The solution is mapped onto a second mode with different elements, and the structure is further loaded statically. Finally, the eigenvalues of the loaded structure are extracted via the AMS eigensolver.

Results and discussion

The following table shows the comparison of eigenfrequencies obtained by the Lanczos and AMS eigensolvers.

ModeLanczosAMS
114.92514.926
243.61443.617
348.56648.571
495.49095.540

Input files

ams_mapsolution_1.inp

Original model preloaded statically.

ams_mapsolution_2.inp

Solution-mapped model with further preloading and frequency extraction using the AMS eigensolver.

VI. Models with material orientations, nodal transformations, and initial conditions

Elements tested

C3D8    SFM3D4R    S4    S8R   

Features tested

Material orientations, nodal transformations, initial conditions, selective modal recovery, and full modal recovery.

Problem description

Relatively small problems with simple topologies constructed for testing the features mentioned above.

Results and discussion

The eigenfrequencies obtained by the AMS and Lanczos eigensolvers are identical for the model with material orientations and initial conditions. The model with nodal transformations exhibits differences smaller than 1%.

Input files

ams_material_ori.inp

Model with material orientations and initial conditions.

ams_nodal_transf.inp

Model with nodal transformations.

VII. Models with residual modes

Elements tested

CPE4R    C3D20R   

Features tested

Residual modes, selective modal recovery, and full modal recovery.

Problem description

Models of simple topology to test the accuracy of residual modes using the AMS eigensolver.

Results and discussion

The following table compares the eigenmodes obtained using the Lanczos and AMS eigensolvers.

ResidualModeLanczosAMS
no14992.34993.1
no25430.05430.9
no37340.87344.6
no410875.10877.
no513716.13724.
yes625445.24359.
yes736445.34198.
yes840925.35377.

The maximum displacement in the steady-state dynamic step at 13kHz is 1.949 units with the Lanczos procedure, versus 1.848 units with the AMS eigensolver.

Input files

ams_resmod_c3d20r.inp

Three-dimensional model with residual modes, AMS, and full modal recovery.

lanczos_resmod_c3d20r.inp

Three-dimensional model with residual modes and Lanczos.

ams_resmod_cpe4r.inp

Two-dimensional model with residual modes, AMS, and selective modal recovery.

VIII. Miscellaneous models

Elements tested

SAXA12    M3D4   

Features tested

Motion of material through the mesh and section distributions.

Problem description

Models with simple topology to test the features mentioned above.

Results and discussion

The results are identical using both the Lanczos and AMS eigensolvers for the model with material motion. For the model with section distributions and SAXA12 elements the results differ slightly in the fifth eigenvalue, as shown in the table below.

ModeLanczosAMS
1531.33531.33
2771.31771.31
31017.61017.6
41129.41129.4
51217.01217.1
61639.21639.2
71754.11754.1
82275.62275.6
93382.03382.0
103490.53490.5
113556.73556.7
123994.93994.9

Input files

ams_motion.inp

Model with material motion.

lanczos_resmod_c3d20r.inp

Model with section distributions and SAXA12 elements.