WORKSHOP 6

Modal Frequency
Response Analysis

Objectives:
■ Define a frequency-varying excitation.

■ Produce a [Link] input file from a dynamic math

model created in Workshop 1.

■ Submit the file for analysis in [Link].

■ Compute nodal displacements for desired frequency

domain.

[Link] 102 Exercise Workbook6-1

6-2 [Link] 102 Exercise Workbook
WORKSHOP 6 Modal Frequency Response Analysis

Model Description:
Using the modal method, determine the frequency response of the flat
rectangular plate, created in Workshop 1, excited by a 0.1 psi pressure
load over the total surface of the plate and a 1.0 lb. force at a corner of
the tip lagging 45o. Use a modal damping of ξ = 0.03. Use a frequency
step of 20 hz between a range of 20 and 1000 hz; in addition, specify five
evenly spaced excitation frequencies between the half power points of
each resonant frequency between the range of 20-1000 hz.

Below is a finite element representation of the flat plate. It also contains

the loads and boundary constraints.

Figure 6.1-Loads and Boundary Conditions

0.1 psi over the total surface

1.0

[Link] 102 Exercise Workbook 6-3

Suggested Exercise Steps:
■ Reference a previously created dynamic math model, [Link], by
using the INCLUDE statement.

■ Specify modal damping as a tabular function of natural frequency

(TABDMP1).

■ Define the frequency-varying pressure loading (PLOAD2, LSEQ

and RLOAD2).

■ Define the frequency-varying tip load (DAREA and RLOAD2).

■ Define a set of frequencies to be used in the solution (FREQ1,

FREQ4).

■ Prepare the model for a modal frequency response analysis (SOL

111).

■ Define the dynamic load phase lead modal frequency response

(DPHASE).

■ Request response in terms of nodal displacement at Grids 11, 33,

and 55.

■ Generate an input file and submit it to the [Link] solver for

direct transient analysis.

■ Review the results, specifically the nodal displacements and phase

angles.

6-4 [Link] 102 Exercise Workbook

 WORKSHOP 6 Modal Frequency Response Analysis
ID SEMINAR,PROB6
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK
[Link] 102 Exercise Workbook 6-5
1 2 3 4 5 6 7 8 9 10

6-6 [Link] 102 Exercise Workbook

 WORKSHOP 6 Modal Frequency Response Analysis

1 2 3 4 5 6 7 8 9 10

ENDDATA

[Link] 102 Exercise Workbook 6-7

Exercise Procedure:
1. Users who are not utilizing [Link] for generating an input
file should go to Step 10, otherwise, proceed to step 2.

2. Create a new database called [Link].

File/New
New Database Name prob6
OK

In the New Model Preferences form, set the following:

Tolerance ◆ Default
Analysis Code: MSC/NASTRAN
Analysis Type: Structural
OK

3. Create the model by importing an existing [Link] input file,

([Link]).

◆ Analysis
Action: Read Input file
Object: Model Data
Method Translate
Select Input File...
Select File [Link]
OK
Apply
OK

4. Activate the entity labels by selecting the Show Labels icon on the
toolbar.

Show Labels

6-8 [Link] 102 Exercise Workbook

WORKSHOP 6 Modal Frequency Response Analysis

5. Create a frequency dependent load case for the frequency response.

◆ Load Cases
Action: Create
Load Case Name: frequency_Response
Load Case Type: Time Dependent
Assign/Prioritize Loads/BCs Displ_spc1.1
(Highlight the following:)

OK
Apply

6. Create the frequency dependent field for the frequency response.

◆ Fields
Action: Create
Object: Non Spatial
Method Tabular Input
Field Name frequency_dependent_load
Active Independent Variables ■ Frequency (f)
[Options... ]
Maximum Number of f 2
OK
Input Data...

Enter the Time/Frequency Scalar Table Data form with the data
below.

Freq (f) Value

1 10. 1.0
2 1000. 1.0
OK
Apply

[Link] 102 Exercise Workbook 6-9

 7. Create the frequency dependent unit force.

◆ Loads/BCs
Action: Create
Object: Pressure
Type: Element Uniform
New Set Name pressure
Target Element Type: 2D
Input Data...
Top Surf Pressure -0.1
* Time/Freq. Dependence f:frequency_dependent_load
(Select from the Time Dependent Fields box.)

OK
Select Application Region...
■ FEM
Select 2D Elements or Edge: Elem 1:40
Add
OK
Apply

8. Create forces.

◆ Load/BCs
Action: Create
Object: Force
Type: Nodal
New Set Name force
Input Data...
Force <F1 F2 F3> <0, 0 , 1>
* Time/Freq. Dependence: f:frequency_dependent_load
(Select from the Time Dependent Fields box.)

OK

6-10 [Link] 102 Exercise Workbook

WORKSHOP 6 Modal Frequency Response Analysis

Select Application Region...

■ FEM
Select Nodes Node 11
Add
OK
Apply

To better visualize the model, hide the labels using the following
toolbar icon:

Hide Labels

Iso 3 View

Action: Plot Markers

Under Assigned Load/BC Sets, highlight:

Displ_spc1.1

Press_pressure

Force_force

Under Select Groups, highlight:

default_group

The model should be similar to Figure 6.2.

[Link] 102 Exercise Workbook 6-11

 Figure 6.2

12345
12345 45
12345 31 46
34 47
12345 21 35 32
23 33 48
22 36
12345 12 11 24 37 34 49
12 25 23 50
1 13 24 38 35
1 .1000 13 26 36 51
2 2 14 .1000 14 25 39
.1000 3 15 27 40 37 52
3 .1000 15 28 26 38 53
.1000 .1000 4 16 27 41
4 .1000 5 .1000 16 29 39 54
.1000 5 17 .1000 17 28 42
.1000 .1000 6 18 30 43 40 55
.1000 6 .1000 18 31 29
.1000 .1000 .1000 7 19 30 44
.1000 7 .1000 8 .1000 19 32
.1000 8 20 .1000 20
.1000 .1000 9 21 33
.1000 .1000 9 .1000
.1000 .1000 10 22 .1000
.1000 10 .1000
.1000 11
.1000 .1000
.1000 .1000
.1000
.1000
.1000
1.000

Y
X

In order to make the next step easier, turn the entity labels back on.

Show Labels

9. Now you are ready to generate an input file for analysis

Click on the Analysis radio button on the Top Menu Bar and complete
the entries as shown here.

◆ Analysis
Action: Analyze
Object: Entire Model
Method Analysis Deck
Job Name prob6
Translation Parameters...
Data Output: XDB and Print
OK
Solution Type...

6-12 [Link] 102 Exercise Workbook

WORKSHOP 6 Modal Frequency Response Analysis

Solution Type: ◆ FREQUENCY RESPONSE

Formulation: Modal
Solution Parameters...
Mass Calculation: Coupled
Wt.-Mass Conversion = 0.00259
Eigenvalue Extraction...
Frequency Range of Interest:
Lower = 10.
Upper = 2000.
OK
OK
OK

Subcase Create...
Available Subcases frequency_dependent
Subcase Parameters...
DEFINE FREQUENCIES...
Add Row

Incr. Type Start Freq. End Freq. No. Incr. Cluster/Spread

1 Linear 20 1000 49 Not Used
2 Linear 20 1000 5 Not Used

OK
Modal Damping Crit. Damp. (CRIT)
DEFINE MODAL DAMPING
Add Row

[Link] 102 Exercise Workbook 6-13

 Frequency Value
1 0 .03
2 10 .03

OK
OK
Output Requests...

Under Output Requests, highlight:

SPCFORCES(SORT1,Real)=All FEM

Delete
OK
Apply
Cancel
Subcase Select...
Subcases Selected: Default
(Click to de-select.)

Subcases for Solution Sequence: frequency_dependent

111
(Click to select.)

OK
Apply

NOTE: Not so Fast!! Before submitting the job, we still need to

define the phase lead.

An input file called [Link] will be generated. This process of

translating your model into an input file is called the Forward
Translation. The Forward Translation is complete when the Heartbeat
turns green.

10. However, since the phase lead term in the equation of the dynamic
loading function (DPHASE) is currently not supported by PATRAN,
you will need to manually edit the file, [Link], to insert the
appropriate phase for the point load.

Search for:
RLOAD1 8 9 1

6-14 [Link] 102 Exercise Workbook

WORKSHOP 6 Modal Frequency Response Analysis

Insert the identification number of the DPHASE entry in the 5th field.
The revised RLOAD1 card should look as follows;
RLOAD1 8 9 92 1

Also, insert the necessary DPHASE card;

DPHASE 92 11 3 -45.

(NOTE: The placement of the numbers must fit the within the alloted
8 character “cell” widths)

[Link] users should now proceed to Step 12.

[Link] 102 Exercise Workbook 6-15

Generating an input file for [Link] Users:
[Link] users can generate an input file using the data from pages
6-3 (general model description). The result should be similar to the
output below.

11. [Link] input file: [Link].

ID SEMINAR, PROB6
SOL 111
TIME 30
CEND
TITLE = FREQUENCY RESPONSE WITH PRESSURE AND POINT LOADS
SUBTITLE = USING THE MODAL METHOD WITH LANCZOS
ECHO = UNSORTED
SEALL = ALL
SPC = 1
SET 111 = 11, 33, 55
DISPLACEMENT(PHASE, PLOT) = 111
METHOD = 100
FREQUENCY = 100
SDAMPING = 100
SUBCASE 1
DLOAD = 100
LOADSET = 100
$
OUTPUT (XYPLOT)
$
XTGRID= YES
YTGRID= YES
XBGRID= YES
YBGRID= YES
YTLOG= YES
YBLOG= NO
XTITLE= FREQUENCY (HZ)
YTTITLE= DISPLACEMENT RESPONSE AT LOADED CORNER, MAGNITUDE
YBTITLE= DISPLACEMENT RESPONSE AT LOADED CORNER, PHASE
XYPLOT DISP RESPONSE / 11 (T3RM, T3IP)
YTTITLE= DISPLACEMENT RESPONSE AT TIP CENTER, MAGNITUDE
YBTITLE= DISPLACEMENT RESPONSE AT TIP CENTER, PHASE
XYPLOT DISP RESPONSE / 33 (T3RM, T3IP)
YTTITLE= DISPLACEMENT RESPONSE AT OPPOSITE CORNER, MAGNITUDE
YBTITLE= DISPLACEMENT RESPONSE AT OPPOSITE CORNER, PHASE

6-16 [Link] 102 Exercise Workbook

WORKSHOP 6 Modal Frequency Response Analysis

XYPLOT DISP RESPONSE / 55 (T3RM, T3IP)

$
BEGIN BULK
$
$ PARAMETERS FOR POST-PROCESSING
PARAM,COUPMASS,1
PARAM,WTMASS,0.00259
$
$ PLATE MODEL DESCRIBED IN NORMAL MODES EXAMPLE
$
INCLUDE ’[Link]’
$
$ EIGENVALUE EXTRACTION PARAMETERS
$
EIGRL, 100, 10., 2000.
$
$ SPECIFY MODAL DAMPING
$
TABDMP1, 100, CRIT,
+, 0., .03, 10., .03, ENDT
$
$ APPLY UNIT PRESSURE LOAD TO PLATE
$
LSEQ, 100, 300, 400
$
PLOAD2, 400, 1., 1, THRU, 40
$
$ APPLY PRESSURE LOAD
$
RLOAD2, 400, 300, , ,310
$
TABLED1, 310,
, 10., 1., 1000., 1., ENDT
$
$ POINT LOAD
$
$ IF 'DAREA' CARDS ARE REFERENCED, THEN
$ 'DPHASE' AND 'DELAY' CAN BE USED
$
RLOAD2, 500, 600, , 320, 310
$
DPHASE, 320, 11, 3, -45.
$

[Link] 102 Exercise Workbook 6-17

$
DAREA, 600, 11, 3, 1.0
$
$ COMBINE LOADS
$
DLOAD, 100, 1., .1, 400, 1.0, 500
$
$ SPECIFY FREQUENCY STEPS
$
FREQ1, 100, 20., 20., 49
FREQ4, 100, 20., 1000., .03, 5
$
ENDDATA

6-18 [Link] 102 Exercise Workbook

WORKSHOP 6 Modal Frequency Response Analysis

Submitting the input file for analysis:

12. Submit the input file to [Link] for analysis.

[Link] submit the [Link] .bdf file, find an available UNIX shell
window. At the command prompt enter nastran [Link]
scr=yes. Monitor the run using the UNIX ps command.

[Link] submit the [Link] .dat file, find an available UNIX

shell window and at the command prompt enter nastran prob6
scr=yes. Monitor the run using the UNIX ps command.

13. When the run is completed, use plotps utility to create a

postscript file, [Link], from the binary plot file
[Link]. The displacement response plots for Grids 11,
33 and 55 are shown in figures 6.2 to 6.4.

14. When the run is completed, edit the prob6.f06 file and
search for the word FATAL. If no matches exist, search for
the word WARNING. Determine whether existing
WARNING messages indicate modeling errors.

[Link] 102 Exercise Workbook 6-19

For [Link] users only. [Link] users
should skip to step 16.
15. While still editing prob6.f06, search for the word:

X Y - O U T P U T S U M M A R Y (spaces are necessary).

Displacement at Grid 11

Frequency (X) Displacement (Y)

140 = __________

440 = __________

Displacement at Grid 33

Frequency (X) Displacement (Y)

140 = __________

660 = __________

Displacement at Grid 55

Frequency (X) Displacement (Y)

140 = __________

1000 = __________

6-20 [Link] 102 Exercise Workbook

 WORKSHOP 6

results:

R E A L E I G E N V A L U E S
X Y - O U T P U T S U M M A R Y ( R E S P O N S E )
SUBCASE CURVE FRAME XMIN-FRAME/ XMAX-FRAME/ YMIN-FRAME/ X FOR YMAX-FRAME/ X FOR
ID TYPE NO. CURVE ID. ALL DATA ALL DATA ALL DATA YMIN ALL DATA YMAX
1 DISP 1 11( 5,--) 2.000000E+01 1.000000E+03 3.481836E-04 4.400000E+02 1.699121E-01 1.336996E+02
2.000000E+01 1.000000E+03 3.481836E-04 4.400000E+02 1.699121E-01 1.336996E+02
Comparison of Results

1 DISP 1 11(--, 11) 2.000000E+01 1.000000E+03 1.390213E+02 1.000000E+03 3.258276E+02 2.000000E+0

2.000000E+01 1.000000E+03 1.390213E+02 1.000000E+03 3.258276E+02 2.000000E+0
1 DISP 2 33( 5,--) 2.000000E+01 1.000000E+03 2.271459E-04 6.600000E+02 1.700317E-01 1.336996E+02
2.000000E+01 1.000000E+03 2.271459E-04 6.600000E+02 1.700317E-01 1.336996E+02
1 DISP 2 33(--, 11) 2.000000E+01 1.000000E+03 1.385571E+02 1.000000E+03 3.263339E+02 2.000000E+0
2.000000E+01 1.000000E+03 1.385571E+02 1.000000E+03 3.263339E+02 2.000000E+0
1 DISP 3 55( 5,--) 2.000000E+01 1.000000E+03 1.278678E-04 1.000000E+03 1.696787E-01 1.336996E+02
2.000000E+01 1.000000E+03 1.278678E-04 1.000000E+03 1.696787E-01 1.336996E+02
1 DISP 3 55(--, 11) 2.000000E+01 1.000000E+03 1.687413E+01 7.001384E+02 3.573561E+02 7.104853E+02
2.000000E+01 1.000000E+03 1.687413E+01 7.001384E+02 3.573561E+02 7.104853E+02

[Link] 102 Exercise Workbook 6-21

Modal Frequency Response Analysis

16. Compare the results obtained in the .f06 file with the following
 17. [Link] Users have finished this exercise. [Link] Users
should proceed to the next step.

18. Proceed with the Reverse Translation process, that is

attaching the [Link] results file into [Link]. To
do this, return to the Analysis form and proceed as
follows.

◆ Analysis
Action: Attach XDB
Object: Result Entities
Method Local
Select Results File...
Select File [Link]
OK
Apply

19. Plot the results in XY Plots.

The first plot is the Displacement versus Frequency plot at Node 11.

◆ Results
Action: Create
Object: Graph
Select Results Cases Frequency_dependent_0 of 54 subcases

Filter Method All

Filter
Apply
Close
Select y Result: Displacement, Translational
Quanity: Z Component
x: Global Variable
Variable: Frequency
Plot Options

6-22 [Link] 102 Exercise Workbook

WORKSHOP 6 Modal Frequency Response Analysis

Complex No. as: Magnitude

Display Attributes
y axis scale Log
Target Entities
Select Nodes: Node 11
Apply

Figure 6.3-Displacement Response at Loaded Corner

The second plot is the Displacement versus Frequency at Node 33..

Select Nodes: Node 33

Apply

[Link] 102 Exercise Workbook 6-23

 Figure 6.4-Displacement Response at Tip Center

6-24 [Link] 102 Exercise Workbook

WORKSHOP 6 Modal Frequency Response Analysis

The third plot is the Displacement versus Frequency at Node 55.

Select Nodes: Node 55

Apply

Figure 6.5-Displacement Response at Opposite Corner

Quit [Link] when you have completed this exercise.

[Link] 102 Exercise Workbook 6-25

6-26 [Link] 102 Exercise Workbook