Cantilever beam II

This is a second example of a cantilever beam subjected to power spectral density. Gravity is used for acceleration. *DATABASE_FREQUENCY_BINARY_D3PSD and *DATABASE_FREQUENCY_BINARY_D3RMS are defined to get output. Example 6.3 from LS-Dyna training class NVH, Fatigue and Frequency Domain Analysis with LS-DYNA by Yun Huang.

*BOUNDARY_SPC_NODE
*CONTROL_IMPLICIT_EIGENVALUE
*CONTROL_IMPLICIT_GENERAL
*CONTROL_TERMINATION
*DATABASE_FREQUENCY_BINARY_D3PSD
*DATABASE_FREQUENCY_BINARY_D3RMS
*DATABASE_HISTORY_NODE
*DEFINE_COORDINATE_SYSTEM
*DEFINE_CURVE
*ELEMENT_BEAM_ORIENTATION
*END
*FREQUENCY_DOMAIN_RANDOM_VIBRATION
*KEYWORD  
*MAT_ELASTIC
*NODE
*PART
*SECTION_BEAM
*TITLE

*KEYWORD  
*TITLE
A beam example under PSD loading
*FREQUENCY_DOMAIN_RANDOM_VIBRATION
         0         0        0.     2000.
      0.05
         1         1         2   

                             2        12 
*DATABASE_FREQUENCY_BINARY_D3PSD
         1
       20.     2000.      1981
*DATABASE_FREQUENCY_BINARY_D3RMS
         1
*CONTROL_IMPLICIT_EIGENVALUE
       100     0.000         0     0.000         1   2000.00         2     0.000
         0                                       1
*CONTROL_IMPLICIT_GENERAL
         1  0.100000         0         0         0         0         0         0
*CONTROL_TERMINATION
  0.100000         0     0.000     0.000     0.000
*BOUNDARY_SPC_NODE
       101         0         1         1         1         1         1         1
*DATABASE_HISTORY_NODE
       101       149
*PART
Part number 101
       101       101       101         0         0         0         0         0
*SECTION_BEAM
       101        13  0.886525         0         0     0.000     0.000
  0.785070  0.049087  0.049087  0.098092     0.000
*MAT_ELASTIC
       101 2.5390E-4 9.9000E+6  0.330000     0.000     0.000         0
*DEFINE_COORDINATE_SYSTEM
         1     0.000     0.000     0.000  1.000000     0.000     0.000
     0.000  1.000000     0.000
*DEFINE_COORDINATE_SYSTEM
         2     0.000     0.000     0.000  1.000000     0.000     0.000
     0.000  1.000000     0.000
*DEFINE_CURVE
        12         0     0.000     0.010     0.000     0.000         0
          20.0000000           1.0000000
        2000.0000000           1.0000000
*ELEMENT_BEAM_ORIENTATION
*NODE
*END