A tube model

This is an example of a tube model under nodal force power spectral density excitation. *DATABASE_FREQUENCY_BINARY_D3PSD and *DATABASE_FREQUENCY_BINARY_D3RMS are used for post-processing of the results. Example 6.6 from LS-Dyna training class NVH, Fatigue and Frequency Domain Analysis with LS-DYNA by Yun Huang.


*BOUNDARY_SPC_SET
*CONTROL_IMPLICIT_EIGENVALUE
*CONTROL_IMPLICIT_GENERAL
*CONTROL_IMPLICIT_NONLINEAR
*CONTROL_TERMINATION
*DATABASE_BINARY_D3PLOT
*DATABASE_BINARY_D3THDT
*DATABASE_FREQUENCY_BINARY_D3PSD
*DATABASE_FREQUENCY_BINARY_D3RMS
*DATABASE_GLSTAT
*DATABASE_MATSUM
*DEFINE_CURVE
*ELEMENT_SHELL
*END
*FREQUENCY_DOMAIN_RANDOM_VIBRATION
*KEYWORD
*MAT_ELASTIC
*NODE
*PART
*SECTION_SHELL
*SET_NODE_LIST
*TITLE

*KEYWORD
*TITLE
pipe simulation
*FREQUENCY_DOMAIN_RANDOM_VIBRATION
         1        20
      0.01
         8         1

      1537         0         3       100
*DATABASE_FREQUENCY_BINARY_D3PSD
         1
       10.      200.       191 
*DATABASE_FREQUENCY_BINARY_D3RMS
         1
*CONTROL_IMPLICIT_GENERAL
         1  1.000000         2         1         2         0         0         1
*CONTROL_IMPLICIT_EIGENVALUE
        20     0.000         0-1.000E+29         01.0000E+29         2     0.000
         0         0         0         0         1
*CONTROL_IMPLICIT_NONLINEAR
         1         0         0     0.000     0.000     0.000     0.000     0.000
         0         0         0         0
         0         0     0.000         1         2
*DEFINE_CURVE
       100         0     0.000   1.0e+03     0.000     0.000         0
           1.0000000           1.0000000
        3000.0000000           1.0000000
*CONTROL_TERMINATION
  1.000000         0     0.000     0.000     0.000
*DATABASE_GLSTAT
 1.0000E-3         0         0         1
*DATABASE_MATSUM
 1.0000E-3         0         0         1
*DATABASE_BINARY_D3PLOT
 1.0000E-3         0         0         0         0
         0
*DATABASE_BINARY_D3THDT
 100.00000         0         0         0         0
*BOUNDARY_SPC_SET
         1         0         1         1         1         1         1         1
*SET_NODE_LIST
         1     0.000     0.000     0.000     0.000MECH      
      1000      1001      1102      1153      1204      1255      1306      1357
      1408      1459      1510      1611      1662      1713      1764      1815
      1866      1917      1968      2019         0         0         0         0
*PART
pipe
         2         2         2         0         0         0         1         0
*SECTION_SHELL
         2        18   0.00000         3         0         0         0         0
  10.00000  10.00000  10.00000  10.00000     0.000     0.000     0.000         0
*MAT_ELASTIC
         2 2.7300E-9 32408.000  0.300000
*ELEMENT_SHELL
*NODE
*END