0% found this document useful (0 votes)
427 views363 pages

Total NAS102

1) The document describes a modal analysis workshop where students will analyze the natural frequencies and mode shapes of a flat rectangular plate using the Lanczos method in MSC NASTRAN. 2) The plate is fixed on one edge and free on the other edges, and students are asked to generate a finite element model of the plate including grids, elements, material properties, and boundary conditions. 3) Students will then create an input file and submit it to the MSC NASTRAN solver to obtain the first five natural frequencies and review the results.
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
427 views363 pages

Total NAS102

1) The document describes a modal analysis workshop where students will analyze the natural frequencies and mode shapes of a flat rectangular plate using the Lanczos method in MSC NASTRAN. 2) The plate is fixed on one edge and free on the other edges, and students are asked to generate a finite element model of the plate including grids, elements, material properties, and boundary conditions. 3) Students will then create an input file and submit it to the MSC NASTRAN solver to obtain the first five natural frequencies and review the results.
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 363

1-2

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 1

Modal Analysis of a Flat Plate

Model Description:
For this example, use Lanczos method to find the first five natural
frequencies and mode shapes of a flat rectangular plate. One of the edges
is fixed, (See Figure 1.2.). Below is a finite element representation of the
rectangular plate. It also contains the geometric dimensions and the loads
and boundary constraints. Table 1.1 contains the necessary parameters to
construct the input file.

Figure 1.1-Grid Coordinates and Element Connectivities

MSC/NASTRAN 102 Exercise Workbook

1-3

Figure 1.2-Loads and Boundary Conditions

1-4

Table 1.1
Length (a)

5 in

Height (b)

2 in

Thickness

0.100 in

Weight Density

0.282 lbs/in3

Mass/Weight Factor

2.59E-3 sec2/in

Elastic Modulus

30.0E6 lbs/in2

Poissons Ratio

0.3

MSC/NASTRAN 102 Exercise Workbook

Modal Analysis of a Flat Plate

WORKSHOP 1

Natural Frequency: Hertz

2
12
3
ij
Eh
f ij = ------------ ----------------------------2
2
2a 12 ( 1 )
where i= 1,2,3, ...
j= 1,2,3, ...
Description: Clamped-Free-Free-Free

a
a = length of plate
b = width of plate
h = thickness of plate
i = number of half-waves in mode shape along horizontal axis
j = number of half-waves in mode shape along vertical axis
C = clamped edge
E = modulus of elasticity
F = free edge
S = simply supported edge

= mass per unit area of plate ( h for a plate material with density )
= Poisson ratio
MSC/NASTRAN 102 Exercise Workbook

1-5

ij2 and (ij)

Mode Sequence
a/b

0.40

3.511

4.786

8.115

13.88

21.64

23.73

(11)

(12)

(13)

(14)

(21)

(22)

3.502

6.406

14.54

22.04

26.07

31.62

(11)

(12)

(13)

(21)

(22)

(14)

3.492

8.525

21.43

27.33

31.11

54.44

(11)

(12)

(21)

(13)

(22)

(23)

3.477

11.68

21.62

39.49

53.88

61.99

(11)

(12)

(21)

(22)

(13)

(31)

3.456

17.99

21.56

57.46

60.58

106.5

(11)

(12)

(21)

(22)

(31)

(32)

2/3

1.0

1.5

2.5

= 0.3

1-6

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 1

Modal Analysis of a Flat Plate

Suggested Exercise Steps


Explicitly generate a finite element representation of the plate
structure. (i.e., The grids (GRID) and element connectivities
(CQUAD4) should be defined manually.)
Define material (MAT1) and element (PSHELL) properties.
Apply the fixed boundary constraints (SPC1).
Prepare the model for a normal modes analysis (SOL 103 and
PARAMs).
PARAM, WTMASS, 0.00259
PARAM, COUPMASS, 1
Generate an input file and submit it to the MSC/NASTRAN
solver for normal modes analysis.

Review the results, specifically the eigenvalues.

MSC/NASTRAN 102 Exercise Workbook

1-7

ID SEMINAR,PROB1
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________

CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK

1-8

MSC/NASTRAN 102 Exercise Workbook

Modal Analysis of a Flat Plate

WORKSHOP 1

MSC/NASTRAN 102 Exercise Workbook

10

1-9

ENDDATA

1-10

MSC/NASTRAN 102 Exercise Workbook

10

WORKSHOP 1

Modal Analysis of a Flat Plate

Exercise Procedure:
1. Users who are not utilizing MSC/PATRAN for generating an input file
should go to Step 11, otherwise, proceed to step 2.
2. Create a new database named prob1.db.
File/New Database
New Database Name

prob1

OK
In the New Model Preference form set the following:
Tolerance

Default

Analysis Code:

MSC/NASTRAN

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

4. Create a surface.
Geometry
Action:

Create

Object:

Surface

Method

XYZ

Vector Coordinates List

<5, 2, 0>

Origin Coordinates List

[ 0, 0, 0]

Apply

MSC/NASTRAN 102 Exercise Workbook

1-11

Figure 1.3-The surface should resemble the output below.

Surface 1.2

Surface 1.3
Y

5. Create the finite element model and mesh the surface.


Finite Elements
Action:

Create

Object:

Mesh Seed

Type:

Uniform

Number of Elements
Number =

10

Curve List

Surface 1.2

(see Figure 1.3)

Apply
5a. Change the number of mesh seeds to 4 and select the right edge.
Number =

Curve List

Surface 1.3

(see Figure 1.3)

1-12

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 1

Modal Analysis of a Flat Plate

Apply
5b. Mesh the surface.
Action:

Create

Object:

Mesh

Type:

Surface

Surface List

Surface 1

Apply
Figure 1.4-The model should appear as below.

6. Create a set of material properties for the plate.


Materials
Action:

Create

MSC/NASTRAN 102 Exercise Workbook

1-13

Object:

Isotropic

Method:

Manual Input

Material Name

mat_1

Input Properties...
Elastic Modulus =

30.0E6

Poisson Ratio =

.3

Density =

.282

Apply
Cancel
7. Define the plate thickness.
Properties
Action:

Create

Dimension:

2D

Type:

Shell

Property Set Name

plate

Input Properties...
Material Name

m:mat_1

(Select from Material Property Sets box.)

Thickness

0.100

OK
Select Members

Surface 1

Add
Apply
8. Apply constraints to the model.
8a. Constrain the left edge from moving through all degrees of freedom.
Load/BCs
Action:
1-14

MSC/NASTRAN 102 Exercise Workbook

Create

WORKSHOP 1

Modal Analysis of a Flat Plate

Object:

Displacement

Type:

Nodal
fixed

New Set Name


Input Data...
Translations <T1 T2 T3>

<0, 0, 0>

Rotations <R1 R2 R3>

<0, 0, >

Analysis Coordinate Frame

Coord 0

OK
Select Application Region...
Select the curve or edge icon.
Curve

Select Geometry Entities

Surface 1.1

(see Figure 1.5)

Add
OK
Apply

MSC/NASTRAN 102 Exercise Workbook

1-15

Figure 1.5

Surface 1.1

9.

Run the analysis.

Before the complete input deck is generated for this analysis, a file that contains
only the model data needs to be created. This file is to be used in later workshops.
Analysis
Action:

Analyze

Object:

Entire Model

Method

Model Only

Job Name

plate

Apply
10.

Now, you will generate the input file for analysis.

Analysis

1-16

Action:

Analyze

Object:

Entire Model

Method

Analysis Deck

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 1

Modal Analysis of a Flat Plate


prob1

Job Name
Solution Type...

NORMAL MODES

Solution Type:
Solution Parameters...
Mass Calculation:

Coupled

Data Deck Echo:

Unsorted

Wt. -Mass Conversion =

.00259

OK
OK
Subcase Create...
Available Subcases

Default

Subcase Parameters...
Number of Desired Roots =

OK
Output Requests...
Under Output Requests, highlight:
SPCFORCES(SORT1,Real)=All FEM
Delete
OK
Apply
Cancel
Apply
An MSC/NASTRAN input file called prob1.bdf will be generated. The
process of translating your model into an input file is called Forward
Translation. The Forward Translation is complete when the Heartbeat
turns green. MSC/PATRAN Users should proceed to step 12.

MSC/NASTRAN 102 Exercise Workbook

1-17

Generating an input file for MSC/NASTRAN Users:


MSC/NASTRAN users can generate an input file using the data from
Table 1.1. The result should be similar to the output below.
11. MSC/NASTRAN Input File: prob1.dat
ID SEMINAR, PROB1
SOL 103
TIME 600
CEND
TITLE = NORMAL MODES EXAMPLE
ECHO = UNSORTED
SUBCASE 1
SUBTITLE= USING LANCZOS
METHOD = 1
SPC = 1
VECTOR=ALL
BEGIN BULK
PARAM
COUPMASS 1
PARAM
WTMASS .00259
EIGRL
1
PSHELL
1
1
.1
CQUAD4
1
1
1
=,*1,=,*1,*1,*1,*1
=8
CQUAD4
11
1
12
=,*1,=,*1,*1,*1,*1
=8
CQUAD4
21
1
23
=,*1,=,*1,*1,*1,*1
=8
CQUAD4
31
1
34
=,*1,=,*1,*1,*1,*1
=8
MAT1
1
3.+7
GRID
1
0.
=,*1,=,*0.5,==
=9
GRID
12
0.
=,*1,=,*0.5,==
=9
GRID
23
0.
=,*1,=,*0.5,==
=9
1-18

5
1
2

13

1
12

13

24

23

24

35

34

35

46

45

.3
0.

.282
0.

.5

0.

1.

0.

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 1

GRID
34
=,*1,=,*0.5,==
=9
GRID
45
=,*1,=,*0.5,==
=9
SPC1
1
ENDDATA

Modal Analysis of a Flat Plate

12345

1.5

0.

0.

2.

0.

0.

12

23

34

45

MSC/NASTRAN 102 Exercise Workbook

1-19

11a.

We will also create an input file plate.bdf, which contains all


the relevant model data. This file is to be used in later
workshops.

GRID
1
=,*1,=,*0.5,==
=9
GRID
12
=,*1,=,*0.5,==
=9
GRID
23
=,*1,=,*0.5,==
=9
GRID
34
=,*1,=,*0.5,==
=9
GRID
45
=,*1,=,*0.5,==
=9
PSHELL
1
1
CQUAD4
1
1
=,*1,=,*1,*1,*1,*1
=8
CQUAD4
11
1
=,*1,=,*1,*1,*1,*1
=8
CQUAD4
21
1
=,*1,=,*1,*1,*1,*1
=8
CQUAD4
31
1
=,*1,=,*1,*1,*1,*1
=8
MAT1
1
3.+7
SPC1
1
12345

1-20

0.

0.

0.

0.

.5

0.

0.

1.

0.

0.

1.5

0.

0.

2.

0.

.1
1

1
2

13

1
12

12

13

24

23

23

24

35

34

34

35

46

45

.3
12

.282
23

34

MSC/NASTRAN 102 Exercise Workbook

45

WORKSHOP 1

Modal Analysis of a Flat Plate

Submitting the input file for analysis:


12.

Submit the input file to MSC/NASTRAN for analysis.


12a. To submit the MSC/PATRAN .bdf file for analysis, find an
available UNIX shell window. At the command prompt enter:
nastran prob1.bdf scr=yes. Monitor the run using the UNIX
ps command.
12b. To submit the MSC/NASTRAN .dat file for analysis, find an
available UNIX shell window. At the command prompt enter:
nastran prob1 scr=yes. Monitor the run using the UNIX ps
command.

13. When the run is completed, edit the prob1.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.
14. While still editing prob1.f06, search for the word:
R E A L (spaces are necessary)
1st = __________Hz
2nd = __________Hz
3rd = __________Hz
4th = __________Hz
5th = __________Hz

MSC/NASTRAN 102 Exercise Workbook

1-21

Comparison of Results
15.

1-22

Compare the results obtained in the .f06 file with the results on
the following page:

MSC/NASTRAN 102 Exercise Workbook

EIGENVALUE
7.056994E+05
1.878432E+07
2.811467E+07
1.931709E+08
2.234434E+08

R E A L
RADIANS
8.400591E+02
4.334088E+03
5.302327E+03
1.389859E+04
1.494802E+04

E I G E N V A L U E S
CYCLES
GENERALIZED
MASS
1.336996E+02
1.000000E+00
6.897916E+02
1.000000E+00
8.438915E+02
1.000000E+00
2.212030E+03
1.000000E+00
2.379052E+03
1.000000E+00

GENERALIZED
STIFFNESS
7.056994E+05
1.878432E+07
2.811467E+07
1.931709E+08
2.234434E+08

MSC/NASTRAN 102 Exercise Workbook

16. MSC/NASTRAN Users have finished this exercise. MSC/PATRAN Users should
proceed to the next step.
17. Proceed with the Reverse Translation process, that is importing the prob1.op2 results file into MSC/
PATRAN. To do this, return to the Analysis form and proceed as follows:
Analysis
Action:

Read Output2

Object:

Result Entities

Method

Translate

Select Results File...


Select Results File
1-23

OK

prob1.op2

Modal Analysis of a Flat Plate

EXTRACTION
ORDER
1
2
3
4
5

WORKSHOP 1

MODE
NO.
1
2
3
4
5

Apply
To simplify the view, turn off the entity labels using the toolbar.
Hide Labels
In addition, switch to a 3 view isometric view point.
Iso 3 View
When the translation is complete bring up the Results form.
Results
Form Type:

Basic

Select Results Cases

1.1-Default, Mode 1

Select Deformation Result

1.1 Eigenvectors, Translational

Apply
The results should resemble Figure 1.6.

1-24

MSC/NASTRAN 102 Exercise Workbook

Modal Analysis of a Flat Plate

WORKSHOP 1

Figure 1.6

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

55

Z
Y
X

To reset the graphics, click on this icon:


Reset Graphics
Repeat the procedure to view the other mode shapes.
Quit MSC/PATRAN when you are finished with this exercise.

MSC/NASTRAN 102 Exercise Workbook

1-25

1-26

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP PROBLEM 2

Modal Analysis of A Flat


Plate using Static Reduction

Objectives
Reduce the dynamic math model, created in Workshop 1, to
one with fewer degrees of freedom.
Produce a MSC/NASTRAN input file.
Submit the file for analysis in MSC/NASTRAN.
Find the first five natural frequencies and mode shapes of the
flat plate.

MSC/NASTRAN 102 Exercise Workbook

2-1

2-2

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 2

Modal Analysis of a Flat Plate using Static Reduction

Model Description:
For this example, reduce the dynamic math model created in Workshop 1,
using static reduction. Then find the first five natural frequencies and
mode shapes using the Automatic Givens method. Use the points
indicated in Figure 2.2 for the A-set.
Figure 2.1-Grid Coordinates and Element Connectivities

MSC/NASTRAN 102 Exercise Workbook

2-3

Figure 2.2-Loads and Boundary Conditions

2-4

Table 2.1
Length (a)

5 in

Height (b)

2 in

Thickness

0.100 in

Weight Density

0.282 lbs/in3

Mass/Weight Factor

2.59E-3 sec2/in

Elastic Modulus

30.0E6 lbs/in2

Poissons Ratio

0.3

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 2

Modal Analysis of a Flat Plate using Static Reduction

Suggested Exercise Steps


Reference a previously created dynamic math model, plate.bdf,
by using the INCLUDE statement.
Prepare the model for a normal modes analysis (SOL 103 and
PARAMs).
PARAM, WTMASS, 0.00259
PARAM, COUPMASS, 1
Define degrees of freedom in the analysis set (ASET) for grids
indicated in Figure 2.2.
Generate an input file and submit it to the MSC/NASTRAN
solver for normal modes analysis.

Review the results, specifically the eigenvalues.

MSC/NASTRAN 102 Exercise Workbook

2-5

ID SEMINAR,PROB2
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________

CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK

2-6

MSC/NASTRAN 102 Exercise Workbook

Modal Analysis of a Flat Plate using Static Reduction

WORKSHOP 2

MSC/NASTRAN 102 Exercise Workbook

10

2-7

ENDDATA

2-8

MSC/NASTRAN 102 Exercise Workbook

10

WORKSHOP 2

Modal Analysis of a Flat Plate using Static Reduction

Exercise Procedure:
1. Users who are not utilizing MSC/PATRAN for generating an input file
should go to Step 7, otherwise, proceed to step 2.
2. Create a new database named prob2.db.
File/New Database
New Database Name

prob2

OK
In the New Model Preference form set the following:
Tolerance

Default

Analysis Code:

MSC/NASTRAN

OK
3. Create the model by importing an existing MSC/NASTRAN input file,
(plate.bdf).
Analysis
Action:

Read Input File

Object:

Model Data

Method

Translate

Select Input File...


Select File

plate.bdf

OK
Apply
OK

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

MSC/NASTRAN 102 Exercise Workbook

2-9

5. Add the pre-defined constraints into the Default load case.


Load Cases
Action:
Select Load Case to Modify

Modify
Default

(Highlight the following:)

Assign/Prioritize Loads/BCs
Select Load/BCs to Add to
Spreadsheet

Displ_spc1.1

(Highlight the following:)

OK
Apply
6. Create the new analysis deck.
Analysis
Action:

Analyze

Object:

Entire Model

Method

Analysis Deck

Solution Type...
Solution Type:

NORMAL MODES

Solution Parameters...
Mass Calculation:

Coupled

Data Deck Echo:

Unsorted

Wt. -Mass Conversion =

.00259

OK
OK
Direct Text Input...
In the Bulk Data Section, type in the following:
ASET1, 345, 3, 5, 7, 9, 11
ASET1, 345, 25, 27, 29, 31, 33
2-10

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 2

Modal Analysis of a Flat Plate using Static Reduction

ASET1, 345, 47, 49, 51, 53, 55

OK
Subcase Create...
Available Subcases

Default

Subcase Parameters...
Extraction Method:

Automatic Givens

Number of Desired Roots =

OK
Output Requests...
Under Output Requests, highlight:
SPCFORCES(SORT1,Real)=All FEM
Delete
OK
Apply
Cancel
Apply
An MSC/NASTRAN input file called prob2.bdf 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. MSC/PATRAN Users should proceed to step 8.

MSC/NASTRAN 102 Exercise Workbook

2-11

Generating an input file for MSC/NASTRAN Users:


MSC/NASTRAN users can generate an input file using the data from
Table 2.1. The result should be similar to the output below.
7.

MSC/NASTRAN input file: prob2.dat

ID SEMINAR, PROB2
SOL 103
TIME 10
CEND
TITLE = REDUCTION PROCEDURES, NORMAL MODES EXAMPLE
SUBTITLE = USING STATIC REDUCTION
ECHO = UNSORTED
SUBCASE 1
SUBTITLE=USING LANCZOS
METHOD = 1
SPC = 1
VECTOR=ALL
BEGIN BULK
EIGR,1,AGIV,,,,5
PARAM, COUPMASS, 1
PARAM, WTMASS, 0.00259
INCLUDE plate.bdf
$
$ SELECT A-SET, STATIC REDUCTION IS DONE AUTOMATICALLY
$
ASET1,345,3,5,7,9,11
ASET1,345,25,27,29,31,33
ASET1,345,47,49,51,53,55
ENDDATA

2-12

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 2

Modal Analysis of a Flat Plate using Static Reduction

Submitting the input file for analysis:


8.

9.

Submit the input file to MSC/NASTRAN for analysis.


8a.

To submit the MSC/PATRAN .bdf file for analysis, find an


available UNIX shell window. At the command prompt
enter: nastran prob2.bdf scr=yes. Monitor the run using
the UNIX ps command.

8b.

To submit the MSC/NASTRAN .dat file for analysis, find an


available UNIX shell window. At the command prompt
enter: nastran prob2 scr=yes. Monitor the run using the
UNIX ps command.

When the run is completed, edit the prob2.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.

10. While still editing prob2.f06, search for the word:


R E A L (spaces are necessary)
1st = __________Hz
2nd = __________Hz
3rd = __________Hz
4th = __________Hz
5th = __________Hz

Comparison of Results
11.

Compare the results obtained in the .f06 file with the results on
the following page:

MSC/NASTRAN 102 Exercise Workbook

2-13

2-14
MSC/NASTRAN 102 Exercise Workbook

MODE
NO.
1
2
3
4
5

EXTRACTION
ORDER
43
45
44
42
41

EIGENVALUE
7.057452E+05
1.880877E+07
2.818009E+07
1.956108E+08
2.367820E+08

R E A L
E I G E N V A L U E S
RADIANS
CYCLES
GENERALIZED
MASS
8.400864E+02
1.337039E+02
1.000000E+00
4.336908E+03
6.902404E+02
1.000000E+00
5.308492E+03
8.448727E+02
1.000000E+00
1.398609E+04
2.225956E+03
1.000000E+00
1.538772E+04
2.449032E+03
1.000000E+00

GENERALIZED
STIFFNESS
7.057452E+05
1.880877E+07
2.818009E+07
1.956108E+08
2.367820E+08

WORKSHOP 2

Modal Analysis of a Flat Plate using Static Reduction

12. MSC/NASTRAN Users have finished this exercise. MSC/


PATRAN Users should proceed to the next step.
13. Proceed with the Reverse Translation process, that is importing
the prob2.op2 results file into MSC/PATRAN. To do this,
return to the Analysis form and proceed as follows:
Analysis
Action:

Read Output2

Object:

Result Entities

Method

Translate

Select Results File...


Select Results File

prob2.op2

OK
Apply
To simplify the view, turn off the entity labels using the toolbar.
Hide Labels
In addition, switch to a 3 view isometric view point.
Iso 3 View

MSC/NASTRAN 102 Exercise Workbook

2-15

When the translation is complete bring up the Results form.


Results
Form Type:

Basic

Select Results Case

1.1-Default, Mode1

Select Deformation Result

1.1 Eigenvectors, Translational

Apply
Reset the graphics by clicking on this icon:
Reset Graphics
You can go back and select any Results Case, Fringe Results or
Deformation Results you are interested in.
Quit MSC/PATRAN when you are finished with this exercise.

2-16

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP PROBLEM 3

Direct Transient Response


Analysis

Objectives
Define time-varying excitation.
Produce a MSC/NASTRAN input file from dynamic math
model created in Workshop 1.
Submit the file for analysis in MSC/NASTRAN.
Compute nodal displacements for desired time domain.

MSC/NASTRAN 102 Exercise Workbook

3-1

3-2

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 3

Direct Transient Response Analysis

Model Description:
Using the direct method, determine the transient response of the flat
rectangular plate, created in Workshop 1, under time-varying excitation.
This example structure shall be excited by 1 psi pressure load over the
total surface of the plate varying at 250Hz. In addition, a 50 lb force is
applied at a corner of the tip also varying at 250Hz but out-of-phase with
the pressure load. Both time dependent dynamic loads are applied for the
duration of 0.008 seconds only. Use structural damping of g=0.06 and
convert this damping to equivalent viscous damping at 250Hz. Carry the
analysis for 0.04 seconds.
Below is a finite element representation of the flat plate. It also contains
the loads and boundary constraints.
Figure 3.1-Loads and Boundary Conditions
1 psi over the total surface

50.00

MSC/NASTRAN 102 Exercise Workbook

3-3

Suggested Exercise Steps


Reference previously created dynamic math model, plate.bdf,
by using the INCLUDE statement.
Define the time-varying pressure loading (PLOAD2, LSEQ and
TLOAD2). (Hint, be certain to specify phase angle since the
applied loads are out-of-phase).
Define the time-varying tip load (DAREA and TLOAD2).
(Again, be certain to specify the phase angle).
Combine the time-varying loads (DLOAD).
Specify integration time steps (TSTEP).
Prepare the model for a direct transient analysis (SOL 109).
Specify the structural damping and convert this damping to
equivalent viscous damping.
PARAM, G, 0.06
PARAM, W3, 1571.0
Request response in terms of nodal displacement at grid points
11, 33 and 55.
Generate an input file and submit it to the MSC/NASTRAN
solver for direct transient analysis.
Review the results, specifically the nodal displacements and
xy-plot output.

3-4

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 3

Direct Transient Response Analysis

ID SEMINAR,PROB3
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________

CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK

MSC/NASTRAN 102 Exercise Workbook

3-5

3-6

MSC/NASTRAN 102 Exercise Workbook

10

Direct Transient Response Analysis

WORKSHOP 3

10

ENDDATA

MSC/NASTRAN 102 Exercise Workbook

3-7

Exercise Procedure:
1. Users who are not utilizing MSC/PATRAN for generating an input file
should go to Step 13, otherwise, proceed to step 2.
2. Open a new database named prob3.db.
File/New Database
prob3

New Database Name


OK

In the New Model Preference form set the following:


Tolerance

Default

Analysis Code:

MSC/NASTRAN

OK
3. Create the model by importing an existing MSC/NASTRAN input file,
(plate.bdf).
Analysis
Action:

Read Input File

Object:

Model Data

Method

Translate

Select Input File


plate.bdf

Select File
OK
Apply
OK

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

3-8

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 3

Direct Transient Response Analysis

5. Add the pre-defined constraints into the default load case.


Load Cases
Action:

Create

Load Case Name

transient_response

Load Case Type:

Time Dependent

Assign/Prioritize Loads/BCs
Select Load/BCs to Add to
Spreadsheet

Displ_spc1.1

(Select from menu.)

OK
Apply
6. Create a time-dependent field for the transient response of the pressure
loading.
Fields
Action:

Create

Object:

Non Spatial

Method

Tabular Input
time_dependent_pressure

Field Name
[Options ...]
Maximum Number of t

21

OK
Input Data ...
Map Function to Table...
PCL Expression f(t):

sind(90000.*t)

Start Time

0.0

End Time

0.008

Number of Points

20

Apply
MSC/NASTRAN 102 Exercise Workbook

3-9

Cancel
In the Time/Frequency Scalar Table Data window, add the following to
Row 21:
Time(t)
21

0.04

Value
0.0

OK
Apply
7. Create another time-dependent field for the transient response of the
nodal force.
Fields
Action:

Create

Object:

Non Spatial

Method

Tabular Input

Field Name

time_dependent_force

[Options ...]
Maximum Number of t

21

OK
Input Data ...
Map Function to Table...
PCL Expression f(t)

-sind(90000.*t)

Start Time

0.0

End Time

0.008

Number of Points

20

Apply
Cancel

3-10

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 3

Direct Transient Response Analysis

In the Time/Frequency Scalar Table Data window, add the following to


Row 21:
Time(t)
21

0.04

Value
0.0

OK
Apply
8. Create the time dependent pressure.
Loads/BCs
Action:

Create

Object:

Pressure

Type:

Element Uniform

New Set Name

pressure

Target Element Type:

2D

Input Data...
-1

Top Surf Pressure

Note: The default direction of pressure in MSC/PATRAN is opposite


from default MSC/NASTRAN assumption.
Time Dependence:

f:time_dependent_pressure

(Select from the Time Dependent Fields


box)

OK
Select Application Region...
FEM
Select 2D Elements or Edge

Elem 1:40

(Select all elements)

Add
OK
Apply

MSC/NASTRAN 102 Exercise Workbook

3-11

9. Create the time-dependent nodal force.


Loads/BCs
Action:

Create

Object:

Force

Type:

Nodal
force

New Set Name


Input Data...
Spatial Dependence
Force <F1 F2 F3>

<0 0 50>

Time Dependence:

f:time_dependent_force

(Select from the Time Dependent Fields


box)

OK
Select Application Region...
FEM
Node 11

Select Nodes
Add
OK
Apply

To simplify the view, turn off the entity labels using the toolbar.
Hide Labels
In addition, switch to a 3 view isometric view point.
Iso 3 View
The result should be similar to Figure 3.2.

3-12

MSC/NASTRAN 102 Exercise Workbook

Direct Transient Response Analysis

WORKSHOP 3

Figure 3.2-The model with loads and boundary conditions applied.

12345
12345
12345
12345
12345
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
50.00
Z
Y
X

10. Create the analysis.


Analysis
Action:

Analyze

Object:

Entire Model

Method:

Analysis Deck
prob3

Job Name
Solution Type...

TRANSIENT RESPONSE

Solution Type:
Solution Parameters...
Formulation:

Direct

Mass Calculation:

Coupled

Wt.-Mass Conversion =

.00259

Struct. Damping Coeff. =

0.06

W3, Damping Factor =

1571

OK
OK

MSC/NASTRAN 102 Exercise Workbook

3-13

Direct Text Input...


Clear
OK
Subcase Create...
Available Subcases

transient_response

(Select from menu.)

Subcase Parameters...
Ending Time =

.04

Number of Time Steps =

100

OK
Output Requests...
Form Type:

Advanced

Under Output Requests, highlight:


SPCFORCES(SORT1,Real)=All FEM
Delete
Output Requests:

select DISPLACEMENT(...

Sorting:

By Freq/Time

Modify
OK
Apply
Cancel
Subcase Select...
Subcases Selected:

Default

(Click to de-select.)

Subcases for Solution


Sequence: 109
(Click to select.)

OK
3-14

MSC/NASTRAN 102 Exercise Workbook

transient_response

WORKSHOP 3

Direct Transient Response Analysis

Apply
An MSC/NASTRAN input file called prob3.bdf 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. MSC/PATRAN Users should proceed to step 14.

MSC/NASTRAN 102 Exercise Workbook

3-15

Generating an input file for MSC/NASTRAN Users:


MSC/NASTRAN users can generate an input file using the data
previously stated. The result should be similar to the output below.
11. MSC/NASTRAN input file: prob3.dat
ID SEMINAR, PROB3
SOL 109
TIME 30
CEND
TITLE= TRANSIENT RESPONSE WITH TIME DEPENDENT PRESSURE AND POINT LOADS
SUBTITLE= USE THE DIRECT METHOD
ECHO= PUNCH
SPC= 1
SET 1= 11, 33, 55
DISPLACEMENT= 1
SUBCASE 1
DLOAD= 700 $ SELECT TEMPORAL COMPONENT OF TRANSIENT LOADING
LOADSET= 100 $ SELECT SPACIAL DISTRIBUTION OF TRANSIENT LOADING
TSTEP= 100 $ SELECT INTEGRATION TIME STEPS
$
OUTPUT (XYPLOT)
XGRID=YES
YGRID=YES
XTITLE= TIME (SEC)
YTITLE= DISPLACEMENT RESPONSE AT LOADED CORNER
XYPLOT DISP RESPONSE / 11 (T3)
YTITLE= DISPLACEMENT RESPONSE AT CENTER TIP
XYPLOT DISP RESPONSE / 33 (T3)
YTITLE= DISPLACEMENT RESPONSE AT OPPOSITE CORNER
XYPLOT DISP RESPONSE / 55 (T3)
$
BEGIN BULK
PARAM, COUPMASS, 1
PARAM, WTMASS, 0.00259
$
$ PLATE MODEL DESCRIBED IN NORMAL MODES EXAMPLE
$
INCLUDE plate.bdf
$
$ SPECIFY STRUCTURAL DAMPING
$ 3 PERCENT AT 250 HZ. = 1571 RAD/SEC.
$
PARAM, G, 0.06
3-16

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 3

Direct Transient Response Analysis

PARAM, W3, 1571.


$
$ APPLY UNIT PRESSURE LOAD TO PLATE
$
LSEQ, 100, 300, 400
$
PLOAD2, 400, 1., 1, THRU, 40
$
$ VARY PRESSURE LOAD (250 HZ)
$
TLOAD2, 200, 300, , 0, 0., 8.E-3, 250., -90.
$
$ APPLY POINT LOAD OUT OF PHASE WITH PRESSURE LOAD
$
TLOAD2, 500, 600, , 0, 0., 8.E-3, 250., 90.
$
DAREA, 600, 11, 3, 1.
$
$ COMBINE LOADS
$
DLOAD, 700, 1., 1., 200, 50., 500
$
$ SPECIFY INTERGRATION TIME STEPS
$
TSTEP, 100, 100, 4.0E-4, 1
$
ENDDATA

MSC/NASTRAN 102 Exercise Workbook

3-17

Submitting the input file for analysis:


12.

Submit the input file to MSC/NASTRAN for analysis.


12a. To submit the MSC/PATRAN .bdf file for analysis, find an
available UNIX shell window. At the command prompt
enter: nastran prob3.bdf scr=yes. Monitor the run using
the UNIX ps command.
12b. To submit the MSC/NASTRAN .dat file for analysis, find
an available UNIX shell window. At the command prompt
enter: nastran prob3 scr=yes. Monitor the run using the
UNIX ps command.

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


postscript file, prob3.ps, from the binary plot file prob3.plt.
The displacement response plots for Grids 11, 33 and 55 are
shown in figures 3.2, 3.3 and 3.4.
14. Edit the prob3.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.
15. While still editing prob3.f06, search for the word:
D I S P L (spaces are necessary)
Displacement at Grid 11
Time

T3

.0024 = __________
.0052 = __________
.02

= __________

Displacement at Grid 33
Time

T3

.0024 = __________
.0052 = __________
.02

3-18

= __________

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 3

Direct Transient Response Analysis

Displacement at Grid 55
Time

T3

.0024 = __________
.0052 = __________
.02

= __________

MSC/NASTRAN 102 Exercise Workbook

3-19

Comparison of Results
16.

3-20

Compare the results obtained in the .f06 file with the results
on the following page:

MSC/NASTRAN 102 Exercise Workbook

POINT-ID =

11
D I S P L A C E M E N T
TYPE
G
G
G
G

T1

POINT-ID =

0.0
0.0
0.0
0.0

T2

T3
0.0
-2.173625E-02
-7.204904E-02
-1.433462E-01

R1
0.0
1.104167E-02
2.847414E-02
4.082027E-02

R2
0.0
1.050818E-02
2.852519E-02
4.915178E-02

0.0
0.0
0.0
0.0

0.0
0.0
0.0
0.0

0.0
0.0
0.0
0.0

-3.726422E-02
-2.122380E-02
-2.998187E-03
1.535974E-02

-6.629907E-05
-1.431050E-05
-7.089762E-06
5.380207E-06

1.039267E-02
5.916678E-03
8.371174E-04
-4.281030E-03

0.0
0.0
0.0
0.0

33
D I S P L A C E M E N T

MSC/NASTRAN 102 Exercise Workbook

TIME
0.0
4.000000E-04
8.000000E-04
1.200000E-03
.
.
.
3.879996E-02
3.919996E-02
3.959996E-02
3.999996E-02

TYPE
G
G
G
G

T1

G
G
G
G

POINT-ID =

T2

V E C T O R

0.0
0.0
0.0
0.0

0.0
0.0
0.0
0.0

T3
0.0
-1.122398E-02
-4.424753E-02
-1.030773E-01

R1
0.0
9.220218E-03
2.576699E-02
3.819036E-02

R2
0.0
6.138594E-03
2.014980E-02
3.922388E-02

0.0
0.0
0.0
0.0

0.0
0.0
0.0
0.0

0.0
0.0
0.0
0.0

-3.729695E-02
-2.121863E-02
-3.002583E-03
1.535096E-02

1.898676E-05
3.488550E-05
-2.228106E-07
-3.032754E-05

1.037927E-02
5.907703E-03
8.361273E-04
-4.274252E-03

0.0
0.0
0.0
0.0

R3

55
D I S P L A C E M E N T

3-21

TIME
0.0
4.000000E-04
8.000000E-04
1.200000E-03
.
.
.
3.879996E-02
3.919996E-02
3.959996E-02
3.999996E-02

R3

TYPE
G
G
G
G

G
G
G
G

T1

T2

V E C T O R

0.0
0.0
0.0
0.0

0.0
0.0
0.0
0.0

T3
0.0
-2.849185E-03
-1.992890E-02
-6.643156E-02

R1
0.0
7.791447E-03
2.322436E-02
3.540079E-02

R2
0.0
4.611430E-03
1.681028E-02
3.501805E-02

0.0
0.0
0.0
0.0

R3

0.0
0.0
0.0
0.0

0.0
0.0
0.0
0.0

-3.722652E-02
-2.115454E-02
-2.998628E-03
1.529953E-02

1.035188E-04
8.268487E-05
6.654292E-06
-6.482315E-05

1.039059E-02
5.912832E-03
8.371378E-04
-4.277684E-03

0.0
0.0
0.0
0.0

Direct Transient Response Analysis

G
G
G
G

0.0
0.0
0.0
0.0

WORKSHOP 3

TIME
0.0
4.000000E-04
8.000000E-04
1.200000E-03
.
.
.
3.879996E-02
3.919996E-02
3.959996E-02
3.999996E-02

V E C T O R

17. MSC/NASTRAN Users have finished this exercise. MSC/


PATRAN Users should proceed to the next step.
18. Proceed with the Reverse Translation process, that is importing
the prob3.op2 results file into MSC/PATRAN. To do this,
return to the Analysis form and proceed as follows:
Analysis
Action:

Read Output2

Object:

Result Entities

Method

Translate

Select Results File...


Select File

prob3.op2

OK
Apply
When the translation is complete bring up the Results form.
Results
Form Type:

Advanced

Select Results Cases


(Highlight all.)

Get Results
Select Result

1.1 Displacements, Translational

Plot Type

XY Plot

Plot Type Options...


Global Var...
Global Variables

1-Time

Apply
Result(Y)...
Results
Vector Component

1.1-Displacements, Translational
X

OK
Node IDs
3-22

MSC/NASTRAN 102 Exercise Workbook

Node 11

WORKSHOP 3

Direct Transient Response Analysis

Apply
New Title or Title Filter

Displacement Response At
Loaded Corner

Rename
Apply
You may reset the graphics by clicking on this icon:
Reset Graphics

Figure 3.3-Displacement Response at Node 11

MSC/NASTRAN 102 Exercise Workbook

3-23

Repeat the above steps for plotting the xy plots of Node 11, Node 33 and
Node 55. Return to the Results Display form. If the Curves for XY Plot
form and the Result XY Plot Options form are still open, close them by
pushing the Cancel button.
Plot Type Options...
Global Var...
Global Variables

1-Time

Apply
Result(Y)...
Results
Vector Component

1.1-Displacements, Translational
X

OK
Node IDs

Node 33

Apply
New Title or Title Filter

Displacement Response at Tip Center

Rename
Apply
Figure 3.4-Displacement Response at Node 33

3-24

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 3

Direct Transient Response Analysis

Return to the Results Display form. If the Curves for XY Plot form and
the Result XY Plot Options form are still open, close them by pushing the
Cancel button.
Plot Type Option...
Global Var...
Global Variables

1. Time

Apply
Result(Y)...
Results

1.1-Displacements, Translational
X

Vector Component

OK
Node 55

Node IDs
Apply
New Title or Title Filter

Displacement Response at
Opposite Corner

Rename
Apply
Figure 3.5-Displacement Response at Node 55

Quit MSC/PATRAN when you are finished with this exercise.


MSC/NASTRAN 102 Exercise Workbook

3-25

3-26

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP PROBLEM 4

Modal Transient Response


Analysis

Z
Y
X

Objectives
Define time-varying excitation.
Produce a MSC/NASTRAN input file from a dynamic math
model, created in Workshop 1.
Submit the file for analysis in MSC/NASTRAN.
Compute nodal displacements for desired time domain.

MSC/NASTRAN 102 Exercise Workbook

4-1

4-2

MSC/NASTRAN 102 Exercise Workbook

Modal Transient Response Analysis

WORKSHOP 4

Model Description:
Using the Modal Method, determine the transient response of the flat
rectangular plate, created in Workshop 1, under time-varying excitation.
This example structure shall be excited by a 1 psi pressure load over the
total surface of the plate varying at 250Hz. In addition, a 25 lb force is
applied at a corner of the tip also varying at 250Hz but starting 0.004
seconds after the pressure load begins. Both time-dependent dynamics
loads are applied only for the duration of 0.008 seconds only. Use a modal
damping of = 0.03 for all nodes. Carry out the analysis for 0.04 seconds.
Below is a finite element representation of the flat plate. It also contains
the loads and boundary constraints.
Figure 4.1-Loads and Boundary Conditions

12345
12345
12345
12345
12345

1.0 psi over the total surface

25.0
Z
Y
X

MSC/NASTRAN 102 Exercise Workbook

4-3

Suggested Exercise Steps


Reference previously created dynamic math model, plate.bdf, by
using the INCLUDE statement.
Specify modal damping as a tabular function of natural frequency
(TABDMP1).
Define the time-varying pressure loading (PLOAD2, LSEQ and
TLOAD2).
Define the time-varying tip load (DAREA and TLOAD2).
Define the time delay term in the equations of the dynamic loading
function (DELAY).
Combine the time-varying loads (DLOAD).
Specify integration time steps (TSTEP).
Prepare the model for a modal transient analysis (SOL 112).
Request response in terms of nodal displacement at grid 11, 33, and
55.
Generate an input file and submit it to the MSC/NASTRAN solver
for normal modes analysis.

4-4

Review the results, specifically the nodal displacements.

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 4

Modal Transient Response Analysis

ID SEMINAR,PROB4
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________

CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK

MSC/NASTRAN 102 Exercise Workbook

4-5

4-6

MSC/NASTRAN 102 Exercise Workbook

10

Modal Transient Response Analysis

WORKSHOP 4

10

ENDDATA

MSC/NASTRAN 102 Exercise Workbook

4-7

Exercise Procedure:
1. Users who are not utilizing MSC/PATRAN for generating an input file
should go to Step 11, otherwise, proceed to step 2.
2. Create a new database and named prob4.db.
File/New Database
prob4

New Database Name


OK

In the New Model Preference form set the following:


Tolerance

Default

Analysis code:

MSC/NASTRAN

OK
3. Create the model by importing an existing MSC/NASTRAN input file,
(plate.bdf).
Analysis
Action:

Read Input File

Object:

Model Data

Method:

Translate

Select Input File


plate.bdf

Select Input File


OK
Apply
OK

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

4-8

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 4

Modal Transient Response Analysis

5. Add the pre-defined constraints into a newly defined load case.


Load Cases
Action:

Create

Load Case Name:

transient_response

Load Case Type:

Time Dependent

Assign/Prioritize Loads/BCs
Select Load/BCs to Add to
Spreadsheet

Displ_spc1.1

(Select from menu.)

OK
Apply
6. Create a time-dependent field for the pressure loading.
Fields
Action:

Create

Object:

Non Spatial

Method:

Tabular Input
time_dependent_pressure

Field Name:
Options ...
Maximum Number of t:

21

OK
Input Data ...
Map Function to Table...
PCL Expression f(t):

sind(90000.*t)

Start time:

0.0

End time:

0.008

Number of Points:

20

Apply
Cancel
MSC/NASTRAN 102 Exercise Workbook

4-9

Go back to the Time/Frequency Scalar Table Data window, go down to


row 21, and add the following:
Time(t)
21

Value

0.04

0.0

OK
Apply
7. Create a time-dependent field for the nodal force.
7a.First, define the PCL function manually.
7b.The text below defines a PCL function called nodal_force. Using a
text editor, input the text into a file called prob4.pcl.
Function nodal_force(t)
real t
if (t < 0.004 || t > .012) then
return 0.0
else
return sind(90000.*t)
end if
End Function

7c.To compile PCL function, go into the command line and type:
!!input prob4
Fields
Action:

Create

Object:

Non Spatial

Method:

Tabular Input

Field Name:

time_dependent_force

Options...
Maximum Number of t:
OK
Input Data...
4-10

MSC/NASTRAN 102 Exercise Workbook

32

WORKSHOP 4

Modal Transient Response Analysis

Map Function to Table...


PCL Expression f(t):

nodal_force(t)

Start time:

0.000

End time:

0.012

Number of Points:

31

Apply
Cancel
Go back to the Time/Frequency Scalar Table Data window, go down to
row 32, and add the following:
Time(t)
32

0.04

Value
0.0

OK
Apply
8. Create the time-dependent pressure.
Loads/BCs
Action:

Create

Object:

Pressure

Type:

Element Uniform

New Set Name:

pressure

Target Element Type:

2D

Input Data...
Top Surf Pressure

-1

Time Dependence

f:time_dependent_pressure

(Select from the Time Dependent Fields


box.)

OK
Select Application Region ...
FEM
MSC/NASTRAN 102 Exercise Workbook

4-11

Elm 1:40

Select 2D Elements or Edge


(Select all elements.)

Add
OK
Apply
9. Create the time-dependent nodal force.
Loads/BCs
Action:

Create

Object:

Force

Type:

Nodal
force

New Set Name:


Input Data ...
Force <F1 F2 F3>

<0 0 25>

Time Dependence

f:time_dependent_force

(Select from the Time Dependent Fields


box.)

OK
Select Application Region ...
FEM
Node 11

Select Nodes
Add
OK
Apply

To simplify the view, turn off the entity labels using the toolbar.
Hide Labels
In addition, switch to a 3 view isometric view point.
Iso 3 View
4-12

MSC/NASTRAN 102 Exercise Workbook

Modal Transient Response Analysis

WORKSHOP 4

The result should be similar to Figure 4.2.


Figure 4.2

12345
12345
12345
12345
12345
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
1.000
25.00
Z
Y
X

10. Generate the input file.


Analysis
Action:

Analyze

Object:

Entire Model

Method:

Analysis Deck
prob4

Jobname:
Solution Type...

TRANSIENT RESPONSE

Solution Type:
Solution Parameters ...
Formulation

Modal

Mass Calculation

Coupled

Wt.-Mass Conversion

.00259

Eigenvalue Extraction...
Number of Desired Roots

OK

MSC/NASTRAN 102 Exercise Workbook

4-13

OK
OK

Direct Text Input...


Case Control Section

SDAMPING = 100

Bulk Data Section

TABDMP1, 100, CRIT,


+, 0., .03, 10., .03, ENDT

(Note that these are two separate lines.)

OK
Subcase Create...
Available Subcases

transient_response

(Select from menu.)

Subcase Parameters...
Ending Time =

.04

Number of Time Steps =

100

OK
Output Requests...
Form Type:

Advanced

under Output Requests, highlight:


SPCFORCES(SORT1,Real)=ALL FEM
Delete
Output Requests:

select DISPLACEMENT(...

Sorting:

By Freq/Time

Modify
OK
Apply
Cancel
4-14

MSC/NASTRAN 102 Exercise Workbook

Modal Transient Response Analysis

WORKSHOP 4

Subcase Select ...


Subcases Selected:

default

(Click to deselect.)

Subcases for Solution


Sequence: 112
( Click

transient_response

to select.)

OK
Apply
An MSC/NASTRAN input file called prob4.bdf 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. MSC/PATRAN Users should proceed to step 12.

MSC/NASTRAN 102 Exercise Workbook

4-15

Generating an input file for MSC/NASTRAN Users:


MSC/NASTRAN users can generate an input file using the data
previously stated. The result should be similar to the output below.
11. MSC/NASTRAN input File: prob4.dat
ID SEMINAR, PROB4
SOL 112
TIME 30
CEND
TITLE = TRANSIENT RESPONSE WITH TIME DEPENDENT PRESSURE AND POINT LOADS
SUBTITLE = USE THE MODAL METHOD
ECHO = UNSORTED
SPC = 1
SET 111 = 11, 33, 55
DISPLACEMENT(SORT2) = 111
SDAMPING = 100
SUBCASE 1
METHOD = 100
DLOAD = 700
LOADSET = 100
TSTEP = 100
$
OUTPUT (XYPLOT)
XGRID=YES
YGRID=YES
XTITLE= TIME (SEC)
YTITLE= DISPLACEMENT RESPONSE AT LOADED CORNER
XYPLOT DISP RESPONSE / 11 (T3)
YTITLE= DISPLACEMENT RESPONSE AT TIP CENTER
XYPLOT DISP RESPONSE / 33 (T3)
YTITLE= DISPLACEMENT RESPONSE AT OPPOSITE CORNER
XYPLOT DISP RESPONSE / 55 (T3)
$
BEGIN BULK
PARAM, COUPMASS, 1
PARAM, WTMASS, 0.00259
$
$ PLATE MODEL DESCRIBED IN NORMAL MODES EXAMPLE PROBLEM
$
INCLUDE plate.bdf
$
$ EIGENVALUE EXTRACTION PARAMETERS
$
4-16

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 4

Modal Transient Response Analysis

EIGRL, 100, , ,5
$
$ 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
$
$ VARY PRESSURE LOAD (250 HZ)
$
TLOAD2, 200, 300, , 0, 0., 8.E-3, 250., -90.
$
$ APPLY POINT LOAD (250 HZ)
$
TLOAD2, 500, 600,610, 0, 0.0, 8.E-3, 250., -90.
$
DAREA, 600, 11, 3, 1.
DELAY, 610, 11, 3, 0.004
$
$ COMBINE LOADS
$
DLOAD, 700, 1., 1., 200, 25., 500
$
$ SPECIFY INTERGRATION TIME STEPS
$
TSTEP, 100, 100, 4.0E-4, 1
$
ENDDATA

MSC/NASTRAN 102 Exercise Workbook

4-17

Submitting the input file for analysis:


12. Submit the input file to MSC/NASTRAN for analysis.
12a.To submit the MSC/PATRAN .bdf file for analysis, find an
available UNIX shell window. At the command prompt enter:
nastran prob4.bdf scr=yes. Monitor the run using the UNIX ps
command.
12b.To submit the MSC/NASTRAN .dat file for analysis, find an
available UNIX shell window. At the command prompt enter:
nastran prob4 scr=yes. Monitor the run using the UNIX ps
command.
13. When the run is completed, use plotps utility to create a
postscript file, prob4.ps, from the binary plot file prob4.plt.
The displacement response plots for Grids 11, 33 and 55 are
shown in figures 4.3, 4.4, and 4.5.
14. When the run is completed, edit the prob4.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.
15. While still editing prob4.f06, search for the word:
D I S P L (spaces are necessary)
Displacement at Grid 11.
Time

T3

.0064= __________
.0092= __________
.02 = __________
Displacement at Grid 33.
Time

T3

.0068= __________
.0092= __________
.02 = __________

4-18

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 4

Modal Transient Response Analysis

Displacement at Grid 55.


Time

T3

.0068= __________
.0092= __________
.02 = __________

MSC/NASTRAN 102 Exercise Workbook

4-19

Comparison of Results
16. Compare the results obtained in the .f06 file with the results on the following page:

4-20

MSC/NASTRAN 102 Exercise Workbook

POINT-ID =

11
D I S P L A C E M E N T

POINT-ID =

T1
T2
T3
R1
R2
.0
.0
.0
.0
.0
3.138503E-15 5.333171E-15 1.873720E-04 -6.340404E-06 5.161942E-05
9.901832E-13 2.054362E-12 1.177721E-01 3.009433E-03 -3.418436E-02
1.301660E-12 2.610037E-12 1.495051E-01 6.088505E-03 -4.068905E-02
1.361787E-12 2.701773E-12 1.469920E-01 1.272909E-02 -3.955902E-02
-1.399555E-12 -2.778663E-12 -1.534481E-01 -1.183270E-02 4.150834E-02
-1.565439E-12 -3.143368E-12 -1.692225E-01 -1.568289E-02 4.710494E-02
-1.396290E-12 -2.830759E-12 -1.533197E-01 -1.109114E-02 4.241555E-02
1.762308E-13 3.609815E-13 2.043042E-02 1.705799E-04 -5.402198E-03
5.328810E-14 1.099402E-13 6.485358E-03 -1.067494E-05 -1.817145E-03
33
D I S P L A C E M E N T

MSC/NASTRAN 102 Exercise Workbook

TIME
TYPE
.0
G
4.000000E-04
G
5.999999E-03
G
6.399998E-03
G
6.799998E-03
G
7.199998E-03
G
8.799998E-03
G
9.199998E-03
G
2.000000E-02
G
3.959996E-02
G
POINT-ID =

4-21

TIME
.0
4.000000E-04
6.399998E-03
6.799998E-03
7.199998E-03
8.799998E-03
9.199998E-03
9.599999E-03
2.000000E-02
3.959996E-02

R3
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0

TYPE
G
G
G
G
G
G
G
G
G
G

V E C T O R

T1
T2
T3
R1
.0
.0
.0
.0
-1.482405E-15 5.013036E-15 1.835858E-04 -3.682543E-14
-7.504724E-13 1.795108E-12 1.207052E-01 2.766146E-03
-9.796515E-13 2.269019E-12 1.556774E-01 6.135463E-03
-1.045427E-12 2.321332E-12 1.599741E-01 1.308315E-02
-8.759517E-13 1.829195E-12 1.278863E-01 1.804037E-02
1.073612E-12 -2.391129E-12 -1.655026E-01 -1.212454E-02
1.217997E-12 -2.696051E-12 -1.851369E-01 -1.596580E-02
-1.280910E-13 3.184533E-13 2.059413E-02 1.496853E-04
-3.946867E-14 9.686225E-14 6.468208E-03 -2.682333E-05
55
D I S P L A C E M E N T
V E C T O

R2
.0
5.086755E-05
-3.436569E-02
-4.140196E-02
-4.083382E-02
-3.306541E-02
4.270243E-02
4.844051E-02
-5.411018E-03
-1.813416E-03

R3
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0

T1
T2
T3
R1
R2
.0
.0
.0
.0
.0
-4.974573E-15 5.435887E-15 1.873720E-04 6.340404E-06 5.161942E-05
-2.151553E-12 2.189475E-12 1.615714E-01 5.634375E-03 -4.218743E-02
-2.229502E-12 2.212659E-12 1.728057E-01 1.248924E-02 -4.215589E-02
-1.757434E-12 1.697970E-12 1.456460E-01 1.736620E-02 -3.459810E-02
2.297826E-12 -2.283533E-12 -1.773652E-01 -1.152392E-02 4.395273E-02
2.561482E-12 -2.555566E-12 -2.007833E-01 -1.524836E-02 4.984966E-02
2.290908E-12 -2.331874E-12 -1.754856E-01 -1.067918E-02 4.413952E-02
-2.917673E-13 3.108817E-13 2.072625E-02 1.192453E-04 -5.428383E-03
-8.886282E-14 9.431532E-14 6.432103E-03 -4.190110E-05 -1.814187E-03

R3
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0

Modal Transient Response Analysis

TYPE
G
G
G
G
G
G
G
G
G
G

WORKSHOP 4

TIME
.0
4.000000E-04
5.999999E-03
6.399998E-03
6.799998E-03
8.799998E-03
9.199998E-03
9.599999E-03
2.000000E-02
3.959996E-02

V E C T O R

17. MSC/NASTRAN Users have finished this


exercise. MSC/PATRAN Users should proceed to
the next step.
18.

Proceed with the Reverse Translation process, that is


importing the prob4.op2 results file into MSC/PATRAN. To
do this, return to the Analysis form and proceed as follows:

Analysis
Action:

Read Output2

Object:

Result Entities

Method:

Translate

Select Results File...


Select Results File

prob4.op2

OK
Apply
When the translation is complete bring up the Results form.
Results
Form Type:

Advanced

Select Results Case


(Select all.)

Get Results
Select Result

1.1 Displacements, Translational

Plot Type

XY Plot

Plot Type Options...


Global Variable
Global Variables:

1. Time

Apply
Result (Y)

4-22

Results:

1.1 Displacements, Translational

Vector Component

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 4

Modal Transient Response Analysis

OK
Node 11

Node IDs
Apply

Displacement Response At
Loaded Corner

New Title or Title Filter:


Rename
Apply

The output should look similar to Figure 4.3.


Figure 4.3-Displacement Response at Loaded Corner
.180

.120

.0600

0.

-.0600

-.120

-.180
0.

.00700

.0140

.0210

.0280

.0350

.0420

19. Repeat the procedure to find the nodal displacement for Node 33.
Results
Form Type:

Advanced

Select Results Case


(Select all.)

Get Results
Select Result

1.1 Displacements, Translational


MSC/NASTRAN 102 Exercise Workbook

4-23

Plot Type

XY Plot

Plot Type Options...


Global Variable
Global Variables:

1. Time

Apply
Result(Y)...
Results:

1.1 Displacements, Translational

Result(Y)...

OK
Node 33

Node IDs
Apply

Displacement Response at Tip


Center

New Title or Title Filter:


Rename
Apply

The output should look similar to Figure 4.4.


Figure 4.4-Displacement Response at Tip Center

.180

.120

.0600

0.

-.0600

-.120

-.180
0.

4-24

.00700

.0140

.0210

MSC/NASTRAN 102 Exercise Workbook

.0280

.0350

.0420

WORKSHOP 4

Modal Transient Response Analysis

20. Repeat the procedure to find the nodal displacement for Node 55.
Results
Form Type:

Advanced

Select Results Case


(Highlight all.)

Get Results
Select Result

1.1 Displacements, Translational

Plot Type

XY Plot

Plot Type Options...


Global Variable
Global Variables:

1. Time

Apply
Result(Y)...
Results:

1.1 Displacement, Translational

Result(Y)

OK
Node 55

Node IDs
Apply
New Title or Title Filter:

Displacement Response at
Opposite Corner

Rename
Apply
The output should look similar to Figure 4.5.

MSC/NASTRAN 102 Exercise Workbook

4-25

Figure 4.5-Displacement Response at Opposite Corner


.210

.140

.0700

0.

-.0700

-.140

-.210
0.

.00700

.0140

.0210

.0280

.0350

.0420

Quit MSC/PATRAN when you are finished with this exercise.

4-26

MSC/NASTRAN 102 Exercise Workbook

5-2

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 5

Direct Frequency Response Analysis

Model Description:
Using the direct method, determine the frequency response of the flat
rectangular plate, created in Workshop 1, under frequency-varying
excitation. This example structure shall be excited by a unit load at a
corner of the tip. Use a frequency step of 20 Hz between a range of 20
and 1000 Hz. Use structural damping of g=0.06.
Below is a finite element representation of the flat plate. It also contains
the loads and boundary constraints.
Figure 5.1-Loads and Boundary Conditions

1.00

MSC/NASTRAN 102 Exercise Workbook

5-3

Suggested Exercise Steps:


Reference previously created dynamic math model, plate.bdf, by
using the INCLUDE statement
Define the frequency-varying tip load (DAREA and RLOAD2).
Define a set of frequencies to be used in the solution (FREQ1).
Prepare the model for a direct frequency response analysis (SOL
108).
Specify the structural damping.
PARAM, G, 0.06
Request response in terms of nodal displacement at Grids 11, 33
and 55.
Generate an input file and submit it to the MSC/NASTRAN solver
for direct transient analysis.

5-4

Review the results, specifically the nodal displacements and phase


angles.

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 5

Direct Frequency Response Analysis

ID SEMINAR,PROB5
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________

CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK
MSC/NASTRAN 102 Exercise Workbook

5-5

5-6

MSC/NASTRAN 102 Exercise Workbook

10

Direct Frequency Response Analysis

WORKSHOP 5

ENDDATA
MSC/NASTRAN 102 Exercise Workbook

5-7

10

Exercise Procedure:
1. Users who are not utilizing MSC/PATRAN for generating an input
file should go to Step 9, otherwise, proceed to step 2.
2. Create a new database called prob5.db.
File/New Database
prob5

New Database Name:


OK

In the New Model Preference form set the following:


Tolerance:

Default

Analysis Code:

MSC/NASTRAN

Analysis Type:

Structural

OK
3. Create the model by importing an existing MSC/NASTRAN input
file, (plate.bdf).

Analysis
Action:

Read Input file

Object:

Model Data

Method:

Translate

Select Input File ...


plate.bdf
OK
Apply
OK
4. Activate the entity labels by selecting the Show Labels icon on the
toolbar.
Show Labels

5-8

MSC/NASTRAN 102 Exercise Workbook

Direct Frequency Response Analysis

WORKSHOP 5

5. Create a time dependent load case for the transient 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 a frequency dependent field for the frequency dependent


load.

Fields
Action:

Create

Object:

Non Spatial

Method:

Tabular Input
frequency_dependent_load

Field Name

Frequency (f)
[Options ... ]
Maximum Number of f:

OK
Input Data...
Using the data in the table below, enter the values describing the time
dependent force into the Time/Frequency Scalar Table Data form.
Freq (f)

Value

1.0

1000

1.0
MSC/NASTRAN 102 Exercise Workbook

5-9

OK
Apply
7. Create the frequency dependent unit force.

Load/BCs
Action:

Create

Object:

Force

Type:

Nodal
unit_force

New Set Name:


Input Data...
Spatial Dependence/Force:

<0 0 1>

Time Dependence:

f:frequency_dependent_load

(Select from the Time Dependent Fields box)

OK
Select Application Region...

FEM
Node 11

Select Nodes:
Add
OK
Apply

To better visualize the model, hide the entity labels and switch to an
isometric view using the icons below:
Hide Labels

Iso 3 View

The model should be similar to Figure 5.2.

5-10

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 5

Direct Frequency Response Analysis


Figure 5.2

8. 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
prob5

Job Name
Solution Type...

FREQUENCY RESPONSE

Solution Type:
Solution Parameters...
Formulation:

Direct

Mass Calculation:

Coupled

Wt.-Mass Conversion=

0.00259

MSC/NASTRAN 102 Exercise Workbook

5-11

Struct. Damping Coeff. =

0.06

OK
OK
Subcase Create ...
Available Subcases

frequency_response

Subcase Parameters...
Starting Frequency =

20

Ending Frequency =

1000

# of Freq. Increments =

49

OK
Output Requests...
Form Type:

Advanced

under Output Request highlight: SPCFORCES(SORT1,Real)=All FEM


Delete
Output Requests:

select DISPLACEMENT(...

Sorting:

By Freq/Time

Modify
OK
Apply
Cancel
Subcase Select ...
Subcases Selected:

Default

(Click to de-select.)

Subcases for Solution


Sequence: 108
( Click

frequency_response

to select.)

OK
Apply
5-12

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 5

Direct Frequency Response Analysis


An input file called prob5.bdf 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. MSC/PATRAN users should now proceed to Step 10.

MSC/NASTRAN 102 Exercise Workbook

5-13

Generating an input file for MSC/NASTRAN Users:


MSC/NASTRAN users can generate an input file using the data from
page 5-3 (Model Description). The result should be similar to the
output below.
9. MSC/NASTRAN input file: prob5.dat
ID SEMINAR, PROB5
SOL 108
TIME 30
CEND
TITLE = FREQUENCY RESPONSE DUE TO UNIT FORCE AT TIP
ECHO = UNSORTED
SPC = 1
SET 111 = 11, 33, 55
DISPLACEMENT(SORT2, PHASE) = 111
SUBCASE 1
DLOAD = 500
FREQUENCY = 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
XYPLOT DISP RESPONSE / 55 (T3RM, T3IP)
$
BEGIN BULK
PARAM, COUPMASS, 1
PARAM, WTMASS, 0.00259
$
$ PLATE MODEL DESCRIBED IN NORMAL MODES EXAMPLE
$
INCLUDE plate.bdf
5-14

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 5

Direct Frequency Response Analysis

$
$ SPECIFY STRUCTURAL DAMPING
$
PARAM, G, 0.06
$
$ APPLY UNIT FORCE AT TIP POINT
$
RLOAD2, 500, 600, , ,310
$
DAREA, 600, 11, 3, 1.0
$
TABLED1, 310,
, 0., 1., 1000., 1., ENDT
$
$ SPECIFY FREQUENCY STEPS
$
FREQ1, 100, 20., 20., 49
$
ENDDATA

MSC/NASTRAN 102 Exercise Workbook

5-15

Submitting the input file for analysis:


10. Submit the input file to MSC/NASTRAN for analysis.
10a.To submit the MSC/PATRAN .bdf file, find an available UNIX
shell window. At the command prompt enter nastran prob5.bdf
scr=yes. Monitor the run using the UNIX ps command.
10b.To submit the MSC/NASTRAN .dat file, find an available UNIX
shell window and at the command prompt enter nastran prob5
scr=yes. Monitor the run using the UNIX ps command.
11. When the run is completed, use plotps utility to create a postscript
file, prob5.ps, from the binary plot file, prob5.plt. The displacement response plots for Grids 11, 33 and 55 are shown in figures 52 to 5-7.
12. When the run is completed, edit the prob5.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.

5-16

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 5

Direct Frequency Response Analysis

For MSC/NASTRAN users only. MSC/PATRAN


users should skip to step 16.
13. While still editing prob5.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

__________

380

__________

Displacement at Grid 33
Frequency (X)

Displacement (Y)

140

__________

600

__________

Displacement at Grid 55
Frequency (X)

Displacement (Y)

140

__________

1000

__________

MSC/NASTRAN 102 Exercise Workbook

5-17

Comparison of Results
14. Compare the results obtained in the .f06 file with the results on the
following page:

5-18

MSC/NASTRAN 102 Exercise Workbook

11
C O M P L E X

TYPE
G

4.000000E+01

9.799999E+02

1.000000E+03

POINT-ID =

T1
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0

T2
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0

V E C T O R

T3
R1
R2
8.817999E-03
6.435859E-04
2.632016E-03
356.4954
176.5664
176.5000
9.404316E-03
6.434991E-04
2.795561E-03
356.2596
176.5677
176.2785
9.965085E-04
2.691742E-04
4.097779E-04
187.6832
7.8008
15.1581
8.803170E-04
2.354656E-04
3.317750E-04
186.9299
8.2146
14.6645

R3
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0

33
C O M P L E X

MSC/NASTRAN 102 Exercise Workbook

FREQUENCY
2.000000E+01

TYPE
G

4.000000E+01

9.799999E+02

1.000000E+03

POINT-ID =

T1
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0

T2
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0

D I S P L A C E M E N T
(MAGNITUDE/PHASE)

V E C T O R

T3
R1
R2
8.183126E-03
5.993295E-04
2.443290E-03
356.4899
176.5639
176.4950
8.768992E-03
6.006200E-04
2.606561E-03
356.2376
176.5565
176.2581
6.867234E-04
3.836353E-04
5.393046E-04
188.0180
5.5597
10.0794
6.062436E-04
3.454144E-04
4.648783E-04
186.8358
5.4959
8.8514

R3
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0

55
C O M P L E X

5-19

FREQUENCY
2.000000E+01

TYPE
G

4.000000E+01

9.799999E+02

1.000000E+03

T1
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0

T2
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0

D I S P L A C E M E N T
(MAGNITUDE/PHASE)

V E C T O R

T3
R1
R2
7.606255E-03
5.587703E-04
2.371172E-03
356.4844
176.5612
176.4928
8.190030E-03
5.613805E-04
2.534562E-03
356.2155
176.5442
176.2492
2.558788E-04
4.612964E-04
5.702980E-04
193.1958
4.6290
9.0143
2.144666E-04
4.204372E-04
4.981144E-04
190.6200
4.3746
7.6762

R3
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0

Direct Frequency Response Analysis

FREQUENCY
2.000000E+01

D I S P L A C E M E N T
(MAGNITUDE/PHASE)

WORKSHOP 5

POINT-ID =

15. MSC/NASTRAN Users have finished this exercise.


MSC/PATRAN Users should proceed to the next step.
16. Proceed with the Reverse Translation process, that is importing the
prob5.op2 results file into MSC/PATRAN. To do this, return to the
Analysis form and proceed as follows.

Analysis
Action:

Read Output2

Object:

Result Entities

Method:

Translate

Select Results File...


Select Available Files

prob5.op2

OK
Apply
17. Plot the results in XY plots.
The first plot is to make the Displacement versus Frequency plot at
Node 11.

Results
Form Type:

Advanced

Select Result Cases


(Highlight all cases.)

Get Results
Select Result

1.1-Displacements, Translational

Plot Type:

XY Plot

Plot Type Options...


Result XY Plot Types

Results Versus Global Variables

Global Var...
Global Variable:
Apply
Result (Y)...
5-20

MSC/NASTRAN 102 Exercise Workbook

1-Frequency

Direct Frequency Response Analysis

WORKSHOP 5

Results

1.1-Displacements, Translational

Vector Component

Numerical Form for Complex


Results

Mag.

OK
Node 11

Node IDs
Apply...
New Title or Title Filter

Displacement vs Frequency
at Node 11

Rename
Apply
Figure 5.3-Displacement Response at Node 11

MSC/NASTRAN 102 Exercise Workbook

5-21

The next step is to make the plot of Phase versus Frequency. Return
to the Results Display form. If the Curves for XY Plot form and the
Result XY Plot Options form are still open, close them by pushing the
Cancel button.
Plot Type Options...
Result (Y)...
Result

1.1-Displacements, Translational

Vector Component

Numerical Form for Complex


Results

Phase

OK
Node IDs

Node 11

Apply...
Result XY Window Name:

XYWindow2

New Title or Title Filter

Phase vs Frequency at Node 11

Rename
Apply
Figure 5.4-Phase Angle at Node 11

5-22

MSC/NASTRAN 102 Exercise Workbook

Direct Frequency Response Analysis

WORKSHOP 5

Repeat the above steps of plotting the XY plots of Node 11 for Node
33 and 55. Once again, push Cancel to remove any miscellaneous
forms until the Results Display form.
Plot Type Options...
Result (Y)...
Numerical Form for Complex
Results

Mag.

OK
Node 33

Node IDs
Apply...
Result XY Window Name:

XYWindow3

New Title or Title Filter

Displacement vs Frequency
at Node 33

Rename
Apply

MSC/NASTRAN 102 Exercise Workbook

5-23

Figure 5.5-Displacement Response at Node 33

Plot Type Options...


Result (Y)...
Numerical Form for Complex
Results

Phase

OK
Node IDs

Node 33

Apply
Result XY Window Name:

XYWindow4

New Title or Title Filter

Phase vs Frequency
at Node 33

Rename
Apply

5-24

MSC/NASTRAN 102 Exercise Workbook

Direct Frequency Response Analysis

WORKSHOP 5

Figure 5.6-Phase Angle at Node 33

Plot Type Options...


Result (Y)...
Numerical Form for Complex
Results

Mag.

OK
Node 55

Node IDs
Apply...
Result XY Window Name:

XYWindow5

New Title or Title Filter

Displacement vs Frequency
at Node 55

Rename
Apply

MSC/NASTRAN 102 Exercise Workbook

5-25

Figure 5.7-Displacement Response at Node 55

Plot Type Options...


Result (Y)...
Numerical Form for Complex
Results

Phase

OK
Node IDs

Node 55

Apply...
Result XY Window Name:

XYWindow6

New Title or Title Filter

Phase vs Frequency
at Node 55

Rename
Apply

5-26

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 5

Direct Frequency Response Analysis


Figure 5.8-Phase Angle at Node 55

Quit MSC/PATRAN when you have completed this exercise.

MSC/NASTRAN 102 Exercise Workbook

5-27

5-28

MSC/NASTRAN 102 Exercise Workbook

6-2

MSC/NASTRAN 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

MSC/NASTRAN 102 Exercise Workbook

6-3

Suggested Exercise Steps:


Reference a previously created dynamic math model, plate.bdf, 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 direct transient 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 MSC/NASTRAN solver
for direct transient analysis.

6-4

Review the results, specifically the nodal displacements and phase


angles.

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 6

Modal Frequency Response Analysis

ID SEMINAR,PROB6
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________

CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK
MSC/NASTRAN 102 Exercise Workbook

6-5

6-6

MSC/NASTRAN 102 Exercise Workbook

10

Modal Frequency Response Analysis

WORKSHOP 6

10

ENDDATA
MSC/NASTRAN 102 Exercise Workbook

6-7

Exercise Procedure:
1.

Users who are not utilizing MSC/PATRAN for


generating an input file should go to Step 10,
otherwise, proceed to step 2.

2.

Create a new database called prob6.db.

File/New Database
prob6

New Database Name


OK

In the New Model Preference form set the following:


Tolerance

Default

Analysis Code:

MSC/NASTRAN

Analysis Type:

Structural

OK
3.

Create the model by importing an existing MSC/


NASTRAN input file, (plate.bdf).

Analysis
Action:

Read Input file

Object:

Model Data

Method

Translate

Select Input File...


plate.bdf

Select File
OK
Apply
OK
4.

Activate the entity labels by selecting the Show Labels


icon on the toolbar.
Show Labels

6-8

MSC/NASTRAN 102 Exercise Workbook

Modal Frequency Response Analysis

WORKSHOP 6

5.

Create a time dependent load case for the transient


response.

Load Cases
Action:

Create

Load Case Name:

frequency_dependent

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 transient


response.

Fields
Action:

Create

Object:

Non Spatial

Method

Tabular Input
frequency_dependent_load

Field Name

Frequency (f)
[Options... ]
Maximum Number of f

OK
Input Data...
Enter the Time/Frequency Scalar Table Data form with the data
below.
Freq (f)

Value

10.

1.0

1000.

1.0

MSC/NASTRAN 102 Exercise Workbook

6-9

OK
Apply
7.

Create the frequency dependent unit force.

Loads/BCs
Action:

Create

Object:

Pressure

Type:

Element Uniform
pressure

New Set Name


Target Element Type:

2D

Input Data...
Top Surf Pressure

-0.1

Time 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
To better visualize the model, hide the labels using the following
toolbar icon:
Hide Labels

Iso 3 View
The model should be similar to Figure 6.2.

6-10

MSC/NASTRAN 102 Exercise Workbook

Modal Frequency Response Analysis

WORKSHOP 6

Figure 6.2

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

55

Z
Y
X

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

Create forces.

Load/BCs
Action:

Create

Object:

Force

Type:

Nodal
force

New Set Name


Input Data...
Force <F1 F2 F3>

<0, 0 , 1>

Time Dependence:

f:frequency_dependent_load

(Select from the Time Dependent Fields box.)

OK
MSC/NASTRAN 102 Exercise Workbook

6-11

Select Application Region...

FEM
Select Nodes

Node 11

Add
OK
Apply
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

Solution Type...
Solution Type:

FREQUENCY RESPONSE

Solution Parameters...
Formulation:

Modal

Mass Calculation:

Coupled

Wt.-Mass Conversion =

0.00259

Eigenvalue Extraction...
Frequency Range of Interest:
Lower =

10.

Upper =

2000.

OK
OK
OK
Direct Text Input...

6-12

MSC/NASTRAN 102 Exercise Workbook

Modal Frequency Response Analysis

WORKSHOP 6

Case Control Section

SDAMPING = 100

Bulk Data Section

TABDMP1, 100, CRIT,


+, 0., .03, 10., .03, ENDT
FREQ4, 2, 20., 1000., .03, 5

(Each line in the box is a separate line to


input!)

OK
Subcase Create...
frequency_dependent

Available Subcases
Subcase Parameters...
Starting Frequency =

20

Ending Frequency =

1000

# of Freq. Increments =

49

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:


111

frequency_dependent

(Click to select.)

OK
Apply

MSC/NASTRAN 102 Exercise Workbook

6-13

An input file called prob6.bdf 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 to insert the appropriate phase for the point load.

Search for:
RLOAD1

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

92

Also, insert the necessary DPHASE card;


DPHASE

92

11

-45.

(NOTE: The placement of the numbers must fit the within the alloted
8 character cell widths)
MSC/PATRAN users should now proceed to Step 12.

6-14

MSC/NASTRAN 102 Exercise Workbook

Modal Frequency Response Analysis

WORKSHOP 6

Generating an input file for MSC/NASTRAN Users:


MSC/NASTRAN 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.

MSC/NASTRAN input file: prob6.dat.

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
MSC/NASTRAN 102 Exercise Workbook

6-15

YBTITLE= DISPLACEMENT RESPONSE AT OPPOSITE CORNER, PHASE


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 plate.bdf
$
$ 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.
$
6-16

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 6

Modal Frequency Response Analysis

$
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

MSC/NASTRAN 102 Exercise Workbook

6-17

Submitting the input file for analysis:


12.

Submit the input file to MSC/NASTRAN for analysis.


12a. To submit the MSC/PATRAN .bdf file, find an available
UNIX shell window. At the command prompt enter nastran
prob6.bdf scr=yes. Monitor the run using the UNIX ps
command.
12b. To submit the MSC/NASTRAN .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.

6-18

13.

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


postscript file, prob6.ps, from the binary plot file
prob6.plt. 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.

MSC/NASTRAN 102 Exercise Workbook

Modal Frequency Response Analysis

WORKSHOP 6

For MSC/NASTRAN users only. MSC/PATRAN


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

__________

MSC/NASTRAN 102 Exercise Workbook

6-19

Comparison of Results
16.

6-20

Compare the results obtained in the .f06 file with the


results on the following page:

MSC/NASTRAN 102 Exercise Workbook

E )
YMAX-FRAME/
ALL DATA
1.699121E-01
1.699121E-01
3.258276E+02
3.258276E+02
1.700317E-01
1.700317E-01
3.263339E+02
3.263339E+02
1.696787E-01
1.696787E-01
3.573561E+02
3.573561E+02

X FOR
YMAX
1.336996E+02
1.336996E+02
2.000000E+01
2.000000E+01
1.336996E+02
1.336996E+02
2.000000E+01
2.000000E+01
1.336996E+02
1.336996E+02
7.104853E+02
7.104853E+02

Modal Frequency Response Analysis

X FOR
YMIN
4.400000E+02
4.400000E+02
1.000000E+03
1.000000E+03
6.600000E+02
6.600000E+02
1.000000E+03
1.000000E+03
1.000000E+03
1.000000E+03
7.001384E+02
7.001384E+02

WORKSHOP 6

MSC/NASTRAN 102 Exercise Workbook

X Y - O U T P U T S U M M A R Y ( R E S P O N S
SUBCASE CURVE FRAME
XMIN-FRAME/ XMAX-FRAME/ YMIN-FRAME/
ID
TYPE NO. CURVE ID.
ALL DATA
ALL DATA
ALL DATA
1 DISP
1
11( 5,--) 2.000000E+01 1.000000E+03 3.481836E-04
2.000000E+01 1.000000E+03 3.481836E-04
1 DISP
1
11(--, 11) 2.000000E+01 1.000000E+03 1.390213E+02
2.000000E+01 1.000000E+03 1.390213E+02
1 DISP
2
33( 5,--) 2.000000E+01 1.000000E+03 2.271459E-04
2.000000E+01 1.000000E+03 2.271459E-04
1 DISP
2
33(--, 11) 2.000000E+01 1.000000E+03 1.385571E+02
2.000000E+01 1.000000E+03 1.385571E+02
1 DISP
3
55( 5,--) 2.000000E+01 1.000000E+03 1.278678E-04
2.000000E+01 1.000000E+03 1.278678E-04
1 DISP
3
55(--, 11) 2.000000E+01 1.000000E+03 1.687413E+01
2.000000E+01 1.000000E+03 1.687413E+01

6-21

17.

MSC/NASTRAN Users have finished this exercise.


MSC/PATRAN Users should proceed to the next step.

18.

Proceed with the Reverse Translation process, that is


importing the prob6.op2 results file into MSC/PATRAN.
To do this, return to the Analysis form and proceed as
follows.

Analysis
Action:

Read Output2

Object:

Result Entities

Method

Translate

Select Results File...


Select File

prob6.op2

OK
Apply
19.

Plot the results in XY Plots.

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

Results
Form Type:

Advanced

Select Result Cases


(Highlight all cases.)

Get Results
Select Result

1.1-Displacements, Translational

Plot Type:

XY Plot

Plot Type Options...


Result XY Plot Types

Results Versus Global Variables

Global Var...
Global Variable:

1-Frequency

Apply
Result (Y)...
Results:
6-22

MSC/NASTRAN 102 Exercise Workbook

1.1-Displacements, Translational

Modal Frequency Response Analysis

WORKSHOP 6

Vector Component

Numerical Form for Complex


Results

Mag.

OK
Node 11

Node IDs
Apply...
New Title or Title Filter

Displacement Response
at Loaded Corner

Rename
Apply

XY Plot
Action:

Modify

Object:

Axis

Active Axis:

Scale...
Scale:

Logarithmic

Apply
Cancel

MSC/NASTRAN 102 Exercise Workbook

6-23

Figure 6.3-Displacement Response at Loaded Corner

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


Repeat the above steps of plotting the XY plots of Node 11 for Nodes
33 and 55. Push Cancel to remove any miscellaneous forms until the
Results Display form.

Results
Plot Type Options...
Result (Y)...
Numerical Form for Complex
Results

Mag.

OK
Node IDs

Node 33

Apply...
Result XY Window Name:

XYWindow2

New Title or Title Filter:

Displacement Response
at Tip Center

Rename

6-24

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 6

Modal Frequency Response Analysis


Apply

XY Plot
Action:

Modify

Object:

Axis

Active Axis:

Scale...
Scale:

Logarithmic

Apply
Cancel

Figure 6.4-Displacement Response at Tip Center

MSC/NASTRAN 102 Exercise Workbook

6-25

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

Results
Plot Type Options...
Result (Y)...
Numerical Form for Complex
Results

Mag.

OK
Node IDs

Node 55

Apply...
Result XY Window Name:

XYWindow3

New Title or Title Filter

Displacement Response
at Opposite Corner

Rename
Apply

XY Plot
Action:

Modify

Object:

Axis

Active Axis:

Scale...
Scale:
Apply
Cancel

6-26

MSC/NASTRAN 102 Exercise Workbook

Logarithmic

WORKSHOP 6

Modal Frequency Response Analysis


Figure 6.5-Displacement Response at Opposite Corner

Quit MSC/PATRAN when you have completed this exercise.

MSC/NASTRAN 102 Exercise Workbook

6-27

6-28

MSC/NASTRAN 102 Exercise Workbook

7-2

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 7

Direct Transient Response Analysis

Model Description:
Using the direct method, determine the transient response to a unit
acceleration sine pulse of 250 Hz applied at the base in the z direction. A
large mass of 1000 lb is applied to the base. Use a structural damping
coefficient of g = 0.06 and convert this damping to equivalent viscous
damping at 250 Hz.
Below is a finite element representation of the flat plate. It also contains
the loads and boundary constraints.
Figure 7.1-Loads and Boundary Conditions
Concentrated Mass
1000 lbs
1245
1245
1245
1245
1245

Z
Y
X

MSC/NASTRAN 102 Exercise Workbook

7-3

Suggested Exercise Steps


Reference previously created dynamic math model, plate.bdf,
by using the INCLUDE statement.
Modify base constraints and release displacements in the
Z-direction.
Define the time-varying unit acceleration (TLOAD2 and
DAREA).
Create the large mass at the base (CMASS2 and RBE2).
Specify the structural damping and convert this damping to
equivalent viscous damping.
PARAM, G, 0.06
PARAM, W3, 1571
Specify integration time steps (TSTEP).
Prepare the model for a direct transient analysis (SOL 109).
Request response in terms of nodal displacement, velocity and
acceleration at Grids 11, 33 and 55.
Generate an input file and submit it to the MSC/NASTRAN
solver for normal modes analysis.

7-4

Review the results, specifically the nodal displacements,


velocities, and acceleration.

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 7

Direct Transient Response Analysis

ID SEMINAR,PROB7
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________

CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK

MSC/NASTRAN 102 Exercise Workbook

7-5

7-6

MSC/NASTRAN 102 Exercise Workbook

10

WORKSHOP 7

Direct Transient Response Analysis

10

ENDDATA

MSC/NASTRAN 102 Exercise Workbook

7-7

Exercise Procedure:
1. Users who are not utilizing MSC/PATRAN for generating an input file
should go to Step 11, otherwise, proceed to step 2.
2. Create a new database and named prob7.db
File/New Database
prob7

New Database Name


OK

In the New Model Preference form set the following:


Tolerance

Default

Analysis code:

MSC/NASTRAN

OK
3. Create the model by importing an existing MSC/NASTRAN input file,
(plate.bdf).
Analysis
Action:

Read Input File

Object:

Model Data

Method:

Translate

Select Input File


plate.bdf

Select Input File


OK
Apply
OK

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

7-8

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 7

Direct Transient Response Analysis

5. Create a time dependent load case for the transient response..


Load Cases
Action:

Create

Load Case Name:

transient_response

Load Case Type:

Time Dependent

Assign/Prioritize loads/BCs

Displ_spc1.1

(Highlight the following:)

OK
Apply
6. Place the large mass at the base. To do this, it will be necessary to create
a point element at Node 23.

Finite Elements
Action:

Create

Object:

Element

Method:

Edit

Shape:

Point

Node 1 =

Node 23

Apply
7. Next, define the scalar mass.

Properties
Action:

Create

Dimension:

0D

Type:

Mass

Property Set Name:

scalar_mass

Option(s)

Grounded
MSC/NASTRAN 102 Exercise Workbook

7-9

Input Properties...
Mass

1000

DOF at Node 1

UZ

(Select String (UZ) from the button on the


right.)

OK
Select Members:

Elm 41

(Change the select menu from point to point


element, then select the point element which
was previously created.)

Point element
Add
Apply
8. Create the RBE Mass, which will be connected to the remaining base
points.

Finite Elements
Action:

Create

Object:

MPC

Type:

RBE2

Define Terms ...


(Turn off Auto Execute.)

Auto Execute

Create Independent
Node List:

Node 23

Apply

Create Dependent
Node List:
7-10

MSC/NASTRAN 102 Exercise Workbook

Node 1, 12, 34, 45

WORKSHOP 7

Direct Transient Response Analysis


UZ

DOFs
Apply
Cancel
Apply

9. Create the time-dependent field for the unit acceleration.


Fields
Action:

Create

Object:

Non Spatial

Method:

Tabular Input

Field Name:

time_dependent_acceleration

Options ...
Maximum Number of t:

21

OK
Input Data ...
Map Function to Table...
PCL Expression f(t):

sind(90000.*t)

Start time:

0.0

End time:

0.004

Number of Points:

20

Apply
Cancel
Go back to the Time/Frequency Scalar Table Data window.
following to row 21:
Time(t)
21

0.04

Add the

Value
0.0

OK
Apply
MSC/NASTRAN 102 Exercise Workbook

7-11

9a. Now define the unit acceleration.


Loads/BCs
Action:

Create

Object:

Force

New Set Name:

unit_acceleration

Input Data...
Spatial Dependence:

<0, 0, 2.588>

Time Dependence:

f:time_dependent_acceleration

OK
Select Application Region ...
FEM
Node 23
Add
OK
Apply
9b. Finally, modify the existing constrants by releasing DOF3.

Load/BCs
Action:

Modify

Object:

Displacement

Type:

Nodal

Select Set to Modify

spc1.1

Modify Data...
Translations < T1 T2 T3 >

< 0, 0, >

Rotations < R1 R2 R3 >

< 0, 0, >

OK
Apply
7-12

MSC/NASTRAN 102 Exercise Workbook

Direct Transient Response Analysis

WORKSHOP 7

9c. The result should look like Figure 7.2


Figure 7.2

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

55

Z
Y
X

10. Generate the input file.


Analysis
Action:

Analyze

Object:

Entire Model

Method:

Analysis Deck
prob7

Job Name:
Solution Type...

TRANSIENT RESPONSE

Solution Type:
Solution Parameters ...
Formulation

Direct

Mass Calculation

Coupled

Wt.-Mass Conversion

.00259

MSC/NASTRAN 102 Exercise Workbook

7-13

Struct. Dampening Coeff. =

.06

W3, Damping Factor =

1571

OK
OK
Subcase Create ...
Available Subcases

transient_response

Subcase Parameters...
Ending Time =

.04

Number of Time Steps =

200

OK
Output Request
under Select Result Type,
highlight:

Velocities
Accelerations

under Output Requests, highlight:


SPCFORCES(SORT1,Real)=ALL FEM
Delete
OK
Apply
Cancel
Subcase Select ...
Subcases Selected:

Default

(Click to deselect.)

Subcases for Solution


Sequence: 109

transient_response

(Click to select.)

OK
Apply

7-14

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 7

Direct Transient Response Analysis

An MSC/NASTRAN input file called prob7.bdf 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. MSC/PATRAN Users should proceed to step 12.

MSC/NASTRAN 102 Exercise Workbook

7-15

Generating an input file for MSC/NASTRAN Users:


MSC/NASTRAN users can generate an input file using the data
previously stated. The result should be similar to the output below.
11. MSC/NASTRAN Input File: prob7.dat
ID SEMINAR, PROB7
SOL 109
TIME 30
CEND
TITLE = TRANSIENT RESPONSE WITH BASE EXCITATION
SUBTITLE = USING DIRECT TRANSIENT METHOD, NO REDUCTION
ECHO = UNSORTED
SPC = 200
SET 111 = 23, 33
DISPLACEMENT (SORT2) = 111
VELOCITY (SORT2) = 111
ACCELERATION (SORT2) = 111
SUBCASE 1
DLOAD = 500
TSTEP = 100
$
OUTPUT (XYPLOT)
XGRID=YES
YGRID=YES
XTITLE= TIME (SEC)
YTITLE= BASE ACCELERATION
XYPLOT ACCELERATION RESPONSE / 23 (T3)
YTITLE= BASE DISPLACEMENT
XYPLOT DISP RESPONSE / 23 (T3)
YTITLE= TIP CENTER DISPLACEMENT RESPONSE
XYPLOT DISP RESPONSE / 33 (T3)
$
BEGIN BULK
$
$ PLATE MODEL DESCRIBED IN NORMAL MODES EXAMPLE
$
INCLUDE plate.bdf
PARAM, COUPMASS, 1
PARAM, WTMASS, 0.00259
$
$ SPECIFY STRUCTURAL DAMPING
$
PARAM, G, 0.06
7-16

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 7

Direct Transient Response Analysis

PARAM, W3, 1571.


$
$ APPLY EDGE CONSTRAINTS
$
SPC1, 200, 12456, 1, 12, 23, 34, 45
$
$ PLACE BIG FOUNDATION MASS (BFM) AT BASE
$
CMASS2, 100, 1000., 23, 3
$
$ RBE MASS TO REMAINING BASE POINTS
$
RBE2, 101, 23, 3, 1, 12, 34, 45
$
$ APPLY LOADING TO FOUNDATION MASS
$
TLOAD2, 500, 600, , 0, 0.0, 0.004, 250., -90.
$
DAREA, 600, 23, 3, 2.588
$
$ SPECIFY INTEGRATION TIME STEPS
$
TSTEP, 100, 200, 2.0E-4, 1
$
ENDDATA

MSC/NASTRAN 102 Exercise Workbook

7-17

Submitting the input file for analysis:


12. Submit the input file to MSC/NASTRAN for analysis.
12a.To submit the MSC/PATRAN .bdf file for analysis, find an
available UNIX shell window. At the command prompt enter:
nastran prob7.bdf scr=yes. Monitor the run using the UNIX ps
command.
12b.To submit the MSC/NASTRAN .dat file for analysis, find an
available UNIX shell window. At the command prompt enter:
nastran prob7 scr=yes. Monitor the run using the UNIX ps
command.
13. When the run is completed, use plotps utility to create a
postscript file, prob7.ps, from the binary plot file prob7.plt.
The displacement, velocity, and acceleration response plots for
Grids 11, 33 and 55 are shown in figures 7.3 to 7.8.
14. When the run is completed, edit the prob7.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.

7-18

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 7

Direct Transient Response Analysis

15. While still editing prob7.f06, search for the word:


D I S P L (spaces are necessary)
Displacement at Grid 23
Time
.0

T3

= __________

.02 = __________
.04 = __________
Displacement at Grid 33
Time
.0

T3

= __________

.02 = __________
.04 = __________
V E L O C (spaces are necessary)
Velocity at Grid 23
Time
.0

T3

= __________

.02 = __________
.04 = __________
Velocity at Grid 33
Time
.0

T3

= __________

.02 = __________
.04 = __________

MSC/NASTRAN 102 Exercise Workbook

7-19

A C C E L (spaces are necessary)


Acceleration at Grid 23
Time
.0

T3

= __________

.02 = __________
.04 = __________
Acceleration at Grid 33
Time
.0

T3

= __________

.02 = __________
.04 = __________

7-20

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 7

Direct Transient Response Analysis

Comparison of Results
16. Compare the results obtained in the .f06 file with the results on the following page:

MSC/NASTRAN 102 Exercise Workbook

7-21

POINT-ID =

23
D I S P L A C E M E N T

7-22

TIME
.0
1.999998E-02
3.999996E-02

TYPE
G
G
G

T1

POINT-ID =

.0
.0
.0

T2
.0
.0
.0

T3
.0
2.511800E-06
2.523053E-06

MSC/NASTRAN 102 Exercise Workbook

TYPE
G
G
G

T1
.0
.0
.0

POINT-ID =

T2
.0
.0
.0

T3
.0
4.559710E-06
7.363514E-07

TYPE
G
G
G

T1
.0
.0
.0

POINT-ID =

T2
.0
.0
.0

T3
1.029922E-05
4.035838E-07
4.712748E-07

R1
.0
7.870200E-20
-8.475326E-19

R2
.0
-5.699044E-07
4.972196E-07

.0
.0
.0

R1
.0
.0
.0

R2
.0
.0
.0

.0
.0
.0

TYPE
G
G
G

T1
.0
.0
.0

T2
.0
.0
.0

V E C T O R

T3
-4.379307E-07
1.220225E-03
6.578242E-04

R1
1.492096E-20
-2.594927E-15
3.797481E-15

R2
-1.481305E-07
-3.394969E-04
-1.829342E-04

.0
.0
.0

23
A C C E L E R A T I O N
TYPE
G
G
G

T1
.0
.0
.0

POINT-ID =

T2
.0
.0
.0

T3
1.029922E-01
1.615566E-04
-1.353748E-04

V E C T O R
R1
.0
.0
.0

R2
.0
.0
.0

.0
.0
.0

33
A C C E L E R A T I O N

TIME
.0
1.999998E-02
3.999996E-02

V E C T O R

33

POINT-ID =

TIME
.0
1.999998E-02
3.999996E-02

.0
.0
.0

V E C T O R

V E L O C I T Y
TIME
.0
1.999998E-02
3.999996E-02

R2
.0
.0
.0

23
V E L O C I T Y

TIME
.0
1.999998E-02
3.999996E-02

R1
.0
.0
.0

33
D I S P L A C E M E N T

TIME
.0
1.999998E-02
3.999996E-02

V E C T O R

TYPE
G
G
G

T1
.0
.0
.0

T2
.0
.0
.0

T3
-4.379307E-03
-1.464230E+00
1.231636E+00

V E C T O R
R1
1.492096E-16
-1.464654E-11
1.381082E-11

R2
-1.481305E-03
4.074106E-01
-3.427888E-01

.0
.0
.0

WORKSHOP 7

Direct Transient Response Analysis

17. MSC/NASTRAN Users have finished this exercise. MSC/PATRAN


Users should proceed to the next step.
18. Proceed with the Reverse Translation process, that is importing the
prob7.op2 results file into MSC/PATRAN. To do this, return to the
Analysis form and proceed as follows:
Analysis
Action:

Read Output2

Object:

Result Entities

Method:

Translate

Select Results File...


prob7.op2

Select Results File


OK
Apply

When the translation is complete bring up the form.


Results
Form Type:

Advanced

Select Result Cases


(Highlight all cases.)

Get Results
Select Result

1.1 Displacements, Translational

Plot Type

XY Plot

Plot Type Options ...


Global Var...
Global Variables:

1. Time

Apply
Result (Y)...
Results

1.1 Displacements, Translational

Vector Component

OK
MSC/NASTRAN 102 Exercise Workbook

7-23

Node IDs:

Node 23

Apply...
New Title or Title Filter:

Base Displacement At
Node 23

Rename
Apply
The output should look similar to Figure 7.3
Figure 7.3-Base Displacement at Node 23

19. Repeat the procedure to find the nodal displacement for Node 33.
Results
Form Type:

Advanced

Select Result Cases


(Highlight all cases.)

Get Results

7-24

Select Result

1.1 Displacements, Translational

Plot Type

XY Plot

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 7

Direct Transient Response Analysis

Plot Type Options ...


Global Var...
Global Variables:

1. Time

Apply
Result (Y)...
Results

1.1 Displacements, Translational

Vector Component

OK
Node 33

Node IDs:
Apply...
Result XY Window Name:

XYWindow2

New Title or Title Filter:

Tip Displacement at Node 33

Rename
Apply
The output should look similar to Figure 7.4

MSC/NASTRAN 102 Exercise Workbook

7-25

Figure 7.4-Tip Displacement at Node 33

20. Repeat the procedure to find the nodal velocity for Node 23.
Results
Form Type:

Advanced

Select Result Cases


(Highlight all cases.)

Get Results
Select Result

1.1 Velocity, Translational

Plot Type

XY Plot

Plot Type Options ...


Global Var...
Global Variables:

1. Time

Apply
Result (Y)...

7-26

Results

1.1 Velocity, Translational

Vector Component

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 7

Direct Transient Response Analysis

OK
Node 23

Node IDs:
Apply...
Result XY Window Name:

XYWindow3

New Title or Title Filter:

Base Velocity at Node 23

Rename
Apply
The output should look similar to Figure 7.5
Figure 7.5-Base Velocity at Node 23

21. Repeat the procedure to find the velocity at Node 33.


Results
Form Type:

Advanced

Select Result Cases


(Highlight all cases.)

Get Results
MSC/NASTRAN 102 Exercise Workbook

7-27

Select Result

1.1 Velocity, Translational

Plot Type

XY Plot

Plot Type Options ...


Global Var...
Global Variables:

1. Time

Apply
Result (Y)...
Results

1.1 Velocity, Translational

Vector Component

OK
Node IDs:

Node 33

Apply...
Result XY Window Name:

XYWindow4

New Title or Title Filter:

Tip Velocity at Node 33

Rename
Apply
The output should look similar to Figure 7.6

7-28

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 7

Direct Transient Response Analysis

Figure 7.6-Tip Velocity at Node 33

22. Repeat the procedure to find the nodal acceleration for Node 23.
Results
Form Type:

Advanced

Select Result Cases


(Highlight all cases.)

Get Results
Select Result

1.1 Acceleration, Translational

Plot Type

XY Plot

Plot Type Options ...


Global Var...
Global Variables:

1. Time

Apply...
Result (Y)...
Results

1.1 Acceleration, Translational

Vector Component

MSC/NASTRAN 102 Exercise Workbook

7-29

OK
Node IDs:

Node 23

Apply...
Result XY Window Name:

XYWindow5

New Title or Title Filter:

Base Acceleration at Node 23

Rename
Apply
The output should be similar to Figure 7.7
Figure 7.7-Base Acceleration at Node 23

23. Repeat the procedure to find the acceleration at Node 33.


Results
Form Type:

Advanced

Select Result Cases


(Highlight all cases.)

Get Results
Select Result
7-30

MSC/NASTRAN 102 Exercise Workbook

1.1 Acceleration, Translational

WORKSHOP 7

Direct Transient Response Analysis

Plot Type

XY plot

Plot Type Option ...


Global Var...
Global Variables:

1. Time

Apply
Result (Y)...
Results:

1.1 Acceleration, Translational

Vector Component

OK
Node IDs:

Node 33

Apply...
Result XY Window Name:

XYWindow6

New Title or Title Filter:

Tip Acceleration at Node 33

Rename
Apply
The output should look similar to Figure 7.8

MSC/NASTRAN 102 Exercise Workbook

7-31

Figure 7.8-Tip Acceleration at Node 33

Quit MSC/PATRAN when you are finished with this exercise.

7-32

MSC/NASTRAN 102 Exercise Workbook

8-2

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 8

Enforced Motion with Direct Frequency Response

Model Description:
Using the direct method, determine the frequency response of the flat
rectangular plate, created in Workshop 1, under a 0.1 displacement at a
corner of the tip. Use a frequency step of 20 Hz in the range of 20 to 1000
Hz. Use a structural damping of g = 0.06.
Below is a finite element representation of the flate plate. It also contains
the loads and boundary constraints.
Figure 8.1-Loads and Boundary Conditions

12345
12345
12345
12345
12345

Concentrated Mass
(100,000 lbs)
Z
Y
X

MSC/NASTRAN 102 Exercise Workbook

8-3

Suggested Exercise Steps


Reference previously created dynamic math model, plate.bdf,
by using the INCLUDE statement.
Create the large mass at a corner of the tip (CMASS2).
Define the frequency-varying tip displacement (RLOAD2,
TABLED4, DAREA).
Define a set of frequencies to be used in the solution (FREQ1).
Prepare the model for a direct frequency analysis (SOL 108).
Sprecify the structural damping.
PARAM, G, 0.06
Request response in terms of nodal displacement and grid points
11, 33, and 55.
Generate an input file, and submit it to the MSC/NASTRAN
solver for direct transient analysis.
Review the results, specifically the grid displacements.

8-4

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 8

Enforced Motion with Direct Frequency Response

ID SEMINAR,PROB8
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________

CEND
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
BEGIN BULK
MSC/NASTRAN 102 Exercise Workbook

8-5

8-6

MSC/NASTRAN 102 Exercise Workbook

10

WORKSHOP 8

Enforced Motion with Direct Frequency Response

10

ENDDATA
MSC/NASTRAN 102 Exercise Workbook

8-7

Exercise Procedure:
1. Users who are not utilizing MSC/PATRAN for generating an input file
should go to Step 9, otherwise, proceed to step 2.
2. Create a new database and named prob8.db.
File/New Database
prob8

New Database Name


OK

In the New Model Preference form set the following:


Tolerance

Default

Analysis code:

MSC/NASTRAN

OK
3. Create the model by importing an existing MSC/NASTRAN input file,
(plate.bdf).
Analysis
Action:

Read Input File

Object:

Model Data

Method:

Translate

Select Input File


plate.bdf

Select Input File


OK
Apply
OK

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

8-8

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 8

Enforced Motion with Direct Frequency Response

5. Create the frequency dependent load case.


Load Cases
Action:

Create

Load Case Name:

frequency response

Load Case Type:

time_dependent

Assign/Prioritize loads/BCs

Displ_spc1.1

(Highlight the following:)

Apply
6. Place a large mass at a corner of the tip (Node 11). However, a point
element must be created first.
Finite Element
Action:

Create

Object:

Element

Method:

Edit

Shape:

Point

Node 1 =

Node 11

Apply
7. Then define the scalar mass.
Properties
Action:

Create

Dimension:

0D

Type:

Mass

Property Set Name:

scalar_mass

Option(s):

Grounded

Input Properties ...


Mass:

1.0E+5

Dof at Node 1 (Value Type)

UZ

MSC/NASTRAN 102 Exercise Workbook

8-9

OK
Select Members:

Elm 41

(Click on point element in select menu.


Then select Elm 41.)

Point Element
Add
Apply
8. Start the analysis.
Analysis
Action:

Analyze

Object:

Entire Model

Method:

Analysis Deck

Jobname:

prob8

Solution Type...
Solution Type:

FREQUENCY RESPONSE

Solution Parameters ...


Formulation:

Direct

Mass Calculation

Coupled

Wt.-Mass Conversion

.00259

Structure Damping Coeff:

0.06

OK
OK
Direct Text Input ...
Bulk Data Section

8-10

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 8

Enforced Motion with Direct Frequency Response

Bulk Data Section:


(Each line in the box is a separate line to
input!)

RLOAD2,500,600, , ,310
TABLED4,310,0.,1.,0.,10000.,
+,0.,0.,-39.4784,ENDT
DAREA,600,11,3,25.8799

Case Control Section


Case Control Section:

DLOAD=500

OK
Subcase Create...
Available Subcases

frequency_response

Subcase Parameters...
Starting Frequency

20

Ending Frequency

1000

# of Freq. Increments

49

OK
Output Requests...
under Output Request highlight:
SPCFORCES(SORT1,Real)=All FEM
Delete
OK
Apply
Cancel
Subcase Select ...
Subcases Selected: click on

Default

Subcases for Solution


Sequence: 108 click on

frequency_response

OK
Apply

MSC/NASTRAN 102 Exercise Workbook

8-11

An MSC/NASTRAN input file called prob8.bdf 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. MSC/PATRAN Users should proceed to step 10.

8-12

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 8

Enforced Motion with Direct Frequency Response

Generating an input file for MSC/NASTRAN Users:


MSC/NASTRAN users can generate an input file using the data
previously stated. The result should be similar to the output below.
9.

MSC/NASTRAN input file: prob8.dat

ID SEMINAR, PROB8
SOL 108
TIME 30
CEND
TITLE= FREQUENCY RESPONSE DUE TO .1 DISPLACEMENT AT TIP
SUBTITLE= DIRECT METHOD
ECHO= UNSORTED
SPC= 1
SET 111= 11, 33, 55
DISPLACEMENT(PHASE, SORT2)= 111
SDISP(PHASE, SORT2)= ALL
set 222 = 11
OLOAD= 222
SUBCASE 1
DLOAD= 500
FREQUENCY= 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
XYPLOT DISP RESPONSE / 55 (T3RM, T3IP)
$
BEGIN BULK
$
MSC/NASTRAN 102 Exercise Workbook

8-13

$ PLATE MODEL DESCRIBED IN NORMAL MODES EXAMPLE


$
INCLUDE plate.bdf
PARAM, COUPMASS, 1
PARAM, WTMASS, 0.00259
$
$ SPECIFY STRUCTURAL DAMPING
$
PARAM, G, 0.06
$
$ APPLY UNIT DISPLACEMENT AT TIP POINT
$
CMASS2, 5000, 1.0E+5, 11, 3
$
RLOAD2, 500, 600, , ,310
$
TABLED4, 310, 0., 1., 0., 10000.,
,0., 0., -39.4784, ENDT
$
DAREA, 600, 11, 3, 25.8799
$
$ SPECIFY FREQUENCY STEPS
$
FREQ1, 100, 20., 20., 49
$
ENDDATA

8-14

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 8

Enforced Motion with Direct Frequency Response

Submitting the input file for analysis:


10.

Submit the input file to MSC/NASTRAN for analysis.


10a. To submit the MSC/PATRAN .bdf file for analysis, find an
available UNIX shell window. At the command prompt
enter: nastran prob8.bdf scr=yes. Monitor the run using
the UNIX ps command.
10b. To submit the MSC/NASTRAN .dat file for analysis, find an
available UNIX shell window. At the command prompt
enter: nastran prob8 scr=yes. Monitor the run using the
UNIX ps command.

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


postscript file, prob8.ps, from the binary plot file prob8.plt.
The displacement response plots for Grids 11, 33 and 55 are
shown in figures 8.2 to 8.7.
12. Edit the prob8.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.

MSC/NASTRAN 102 Exercise Workbook

8-15

13.

While still editing prob8.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

__________

380

__________

Displacement at Grid 33
Frequency (X)

Displacement (Y)

140

__________

600

__________

Displacement at Grid 55
Frequency (X)

8-16

Displacement (Y)

140

__________

1000

__________

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 8

Enforced Motion with Direct Frequency Response

Comparison of Results
Compare the results obtained in the .f06 file with the results on the
following page:

MSC/NASTRAN 102 Exercise Workbook

8-17

8-18
MSC/NASTRAN 102 Exercise Workbook

X Y - O U T P U T S U M M A R Y ( R E S P O
SUBCASE CURVE FRAME
XMIN-FRAME/
XMAX-FRAME/
YMIN-FRAME/
ID
TYPE
NO.
CURVE ID.
ALL DATA
ALL DATA
ALL DATA
1 DISP
1
11( 5,--) 2.000000E+01 1.000000E+03 9.992202E-02
2.000000E+01 1.000000E+03 9.992202E-02
1 DISP
1
11(--, 11) 2.000000E+01 1.000000E+03 7.680080E-07
2.000000E+01 1.000000E+03 7.680080E-07
1 DISP
2
33( 5,--) 2.000000E+01 1.000000E+03 2.312926E-03
2.000000E+01 1.000000E+03 2.312926E-03
1 DISP
2
33(--, 11) 2.000000E+01 1.000000E+03 3.348117E-01
2.000000E+01 1.000000E+03 3.348117E-01
1 DISP
3
55( 5,--) 2.000000E+01 1.000000E+03 2.434351E-02
2.000000E+01 1.000000E+03 2.434351E-02
1 DISP
3
55(--, 11) 2.000000E+01 1.000000E+03 3.690138E+00
2.000000E+01 1.000000E+03 3.690138E+00

N S E )
X FOR
YMIN
3.600000E+02
3.600000E+02
1.000000E+03
1.000000E+03
6.000000E+02
6.000000E+02
9.799999E+02
9.799999E+02
1.000000E+03
1.000000E+03
1.000000E+03
1.000000E+03

YMAX-FRAME/
ALL DATA
9.992512E-02
9.992512E-02
3.828149E-04
3.828149E-04
8.446401E-01
8.446401E-01
3.599947E+02
3.599947E+02
1.624350E+00
1.624350E+00
3.599892E+02
3.599892E+02

X FOR
YMAX
2.000000E+0
2.000000E+01
3.800000E+02
3.800000E+02
3.800000E+02
3.800000E+02
2.000000E+01
2.000000E+01
3.800000E+02
3.800000E+02
2.000000E+01
2.000000E+01

WORKSHOP 8

Enforced Motion with Direct Frequency Response

14. MSC/NASTRAN Users have finished this exercise. MSC/


PATRAN Users should proceed to the next step.
15. Proceed with the Reverse Translation process, that is importing
the prob8.op2 results file into MSC/PATRAN. To do this,
return to the Analysis form and proceed as follows:
Analysis
Action:

Read Output2

Object:

Result Entities

Method:

Translate

Select Results File...


Select Results File

prob8.op2

OK
Apply
16.

Plot the results in XY plots.

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

Results
Form Type:

Advanced

Select Result Cases


(Highlight all cases.)

Get Results
Select Result

1.1-Displacements, Translational

Plot Type:

XY Plot

Plot Type Options...


Result XY Plot Types

Results Versus Global Variables

Global Var...
Global Variable:

1-Frequency

Apply
Result (Y)...
MSC/NASTRAN 102 Exercise Workbook

8-19

Results:

1.1-Displacements, Translational

Vector Component

Numerical Form for Complex


Results

Mag.

OK
Node IDs

Node 33

Apply...
New Title or Title Filter

Displacement vs Frequency
at Node 11

Rename
Apply
Figure 8.2-Displacement Response at Node 11.

8-20

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 8

Enforced Motion with Direct Frequency Response

The next step is to make the plot of Phase versus Frequency. Return to the
Results Display form. If the Curves for XY Plot form and the Result XY
Plot Options form are still open, close them by pushing the Cancel button.
Plot Type Options...
Result (Y)...
Results:

1.1-Displacements, Translational

Vector Component

Numerical Form for Complex


Results

Phase

OK
Node 11

Node IDs
Apply...
Result XY Window Name:

XYWindow2

New Title or Title Filter

Phase vs Frequency at Node 11

Rename
Apply
Figure 8.3-Phase Angle at Node 11

MSC/NASTRAN 102 Exercise Workbook

8-21

Repeat the above steps of plotting the XY plots of Grids 11 for Node 33
and 55. Once again, push Cancel to remove any miscellaneous forms until
the Results Display form.
Plot Type Options...
Result (Y)...
Numerical Form for Complex
Results

Mag.

OK
Node IDs

Node 33

Apply...
Result XY Window Name:

XYWindow3

New Title or Title Filter

Displacement vs Frequency
at Node 33

Rename
Apply

8-22

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 8

Enforced Motion with Direct Frequency Response

Figure 8.4-Displacement Response at Node 33

Plot Type Options...


Result (Y)...
Numerical Form for Complex
Results

Phase

OK
Node 33

Node IDs
Apply...
Result XY Window Name:

XYWindow4

New Title or Title Filter

Phase vs Frequency
at Node 33

Rename
Apply

MSC/NASTRAN 102 Exercise Workbook

8-23

Figure 8.5-Phase Angle at Node 33

Plot Type Options...


Result (Y)...
Numerical Form for Complex
Results

Mag.

OK
Node IDs

Node 55

Apply...
Result XY Window Name:

XYWindow5

New Title or Title Filter

Displacement vs Frequency
at Node 55

Rename
Apply

8-24

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 8

Enforced Motion with Direct Frequency Response

Figure 8.6-Displacement Response at Node 55

Plot Type Options...


Result (Y)...
Numerical Form for Complex
Results

Phase

OK
Node 55

Node IDs
Apply...
Result XY Window Name:

XYWindow6

New Title or Title Filter

Phase vs Frequency
at Node 55

Rename
Apply

MSC/NASTRAN 102 Exercise Workbook

8-25

Fig 8.7-Phase Angle at Node 55

Quit MSC/PATRAN when you have completed this exercise.

8-26

MSC/NASTRAN 102 Exercise Workbook

9a-2

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 9a

Shock Response Spectrum

Model Description:
Define the shock response of the plate due to a 2.0 in/sec2 sine pulse
applied at the clamped edge. Use modes to a frequency of 1000 Hz
with 3% critical damping. Use the SRSS option for model response
summation.
Figure 9a.1-Model Description and Loading Diagram.

x(t)

FN = 0.0

1.0

2.0

FMAX

uB(t)

Point on larger, vibrating structure.


Response computed for transient analysis.

MSC/NASTRAN 102 Exercise Workbook

9a-3

Suggested Exercise Steps:


Generate the finite element representation of the model
using (GRID) and (CMASS2) elements.
Apply loading to mass, (TLOAD2) and (DAREA).
Specify integration time steps (TSTEP).
Define frequency and damping values for the SDOF
oscillators (DTI).
Specify damping information (FREQ) and natural
frequency (FREQ1).
Define the parameter to calculate shock spectrum.
PARAM, RESPECTRA, 0
Generate an input file and submit it to the MSC/NASTRAN
solver for direct transient analysis.

9a-4

Review the results.

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 9a

Shock Response Spectrum

ID SEMINAR,PROB9A
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________

CEND
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
MSC/NASTRAN 102 Exercise Workbook

9a-5

BEGIN BULK
1

9a-6

MSC/NASTRAN 102 Exercise Workbook

10

WORKSHOP 9a

Shock Response Spectrum

10

ENDDATA
MSC/NASTRAN 102 Exercise Workbook

9a-7

Generating an input file for MSC/NASTRAN Users:


MSC/NASTRAN users can generate an input file using the data from
page 9-3 (Model Description). The result should be similar to the
output below.
1.

MSC/NASTRAN input file: prob9a.dat.

ID SEMINAR, PROB9a
SOL 109
TIME 30
CEND
TITLE= TRANSIENT RESPONSE
SUBTITLE= USING DIRECT TRANSIENT METHOD
LABEL= SHOCK SPECTRUM CALCULATION
ECHO= UNSORTED
SPC= 100
SET 111= 3000
DISPLACEMENT (SORT2)= 111 $ AT LEAST DISP AND VEL MUST APPEAR
VELOCITY (SORT2)= 111
ACCELERATION ()= 111
DLOAD= 500
TSTEP= 100
$
OUTPUT (XYPLOT)
$
$ SHOCK RESPONSE IS ONLY AVAILABLE IN PLOT OR PUNCH OUTPUT. THEREFORE,
$ THE OUTPUT(XYPLOT) SECTION OF THE CASE CONTROL MUST BE USED.
$
XGRID=YES
YGRID=YES
XYPLOT ACCE / 3000(T1)
XLOG= YES
YLOG= YES
$
$ RELATIVE SHOCK RESPONSES ARE CONTAINED IN THE IMAGINARY/PHASE
$ COMPONENTS OF THE OUTPUT
$ ABSOLUTE SHOCK RESPONSES ARE CONTAINED IN THE REAL/MAGNITUDE
$ COMPONENTS OF THE OUTPUT
$
XTITLE= FREQUENCY (CYCLES/SEC)
YTITLE= RELATIVE DISPLACEMENT
XYPLOT DISP SPECTRAL 1 / 3000 (T1IP)
YTITLE= RELATIVE VELOCITY
XYPLOT VELOCITY SPECTRAL 1 / 3000 (T1IP)
YTITLE= ABSOLUTE ACCELERATION
XYPLOT ACCELERATION SPECTRAL 1 / 3000 (T1RM)
$
9a-8

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 9a

Shock Response Spectrum

$ PUNCH SHOCK SPECTRUM FOR LATER USE


$
XYPUNCH ACCELERATION SPECTRAL 1 / 3000(T1RM)
$
BEGIN BULK
$
$ DEFINE GRID POINT
$
GRID, 3000, ,0.,0.,0., ,23456
$
$ DEFINE MASS
$
CMASS2, 100, 1.0, 3000, 1
$
$ APPLY LOADING TO MASS
$
TLOAD2, 500, 600, , 0, 0., 0.004, 250., -90.
$
DAREA, 600, 3000, 1, 1.
$
$ SPECIFY INTEGRATION TIME STEPS
$
TSTEP, 100, 100, 4.0E-4, 1
$
$ PARAMETER TO CALCULATE SHOCK SPECTRUM
$
PARAM, RSPECTRA, 0
$
$ SPECIFY FREQUENCY AND DAMPING VALUES FOR
$ THE SDOF OSCILLATORS AT GRID 3000
$
DTI, SPSEL, 0
DTI, SPSEL, 1, 111, 222, 3000
$ 1= SUBCASE... 111= DAMPING... 222= FREQUENCIES... 3000= GRID NUMBER
$
$ DAMPING INFORMATION FOR OSCILLATORS
$
FREQ, 111, 0., 0.02, 0.04
$
$ NATURAL FREQUENCIES OF OSCILLATORS
$
FREQ1, 222, 20., 20., 49
$
ENDDATA
MSC/NASTRAN 102 Exercise Workbook

9a-9

Submitting the input file:


2. Submit the input file to MSC/NASTRAN for analysis.
To submit the MSC/NASTRAN .dat file, find an available UNIX shell
window and at the command prompt enter nastran prob9a scr=yes.
Monitor the run using the UNIX ps command.
3.

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


postscript file, prob9a.ps, from the binary plot file
prob9a.plt. The nonlinear force and displacement plots
are shown on the following pages.

4.

When the run is completed, edit the prob9a.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.

Comparison of Results
5.

9a-10

Compare the plot made from the exercise with the plots on
the following pages.

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 9a

Shock Response Spectrum

Figure 9a.2

MSC/NASTRAN 102 Exercise Workbook

9a-11

Figure 9a.3

9a-12

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 9a

Shock Response Spectrum

Figure 9a.4

MSC/NASTRAN 102 Exercise Workbook

9a-13

Figure 9a.5

9a-14

MSC/NASTRAN 102 Exercise Workbook

9b-2

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 9b

Shock Response Spectrum

Model Description:
Define the shock response of the plate due to a 2.0 in/sec2 sine pulse
applied at the clamped edge. Use modes to a frequency of 1000 Hz
with 3% critical damping. Use the SRSS option for model response
summation.
Figure 9b.1-Model Description and Loading Diagram
x(t)

FN = 0.0

1.0

2.0

FMAX

uB(t)

Point on larger, vibrating structure.


Response computed for transient analysis.

MSC/NASTRAN 102 Exercise Workbook

9b-3

Suggested Exercise Steps:


Reference a previously created dynamic math model,
plate.bdf, by using the INCLUDE statement.
Modify boundary conditions for clamped modes.
Place big foundation mass (BFM) at base to simulate
clamped modes (CMASS2).
RBE mass to remaining base point (RBE2).
Identify excitation DOFs (SUPORT).
Specify damping table (TABDMP1).
Specify shock spectrum to be used (DLOAD).
Specify shock tables (DTI).
Insert punch output for shock spectrum calculation.
Specify the appropriate parameters.
PARAM, SCRSPEC, 0
PARAM, OPTION, SRSS
PARAM, LFREQ, 0.1
PARAM, HFREQ, 1000
Generate an input file and submit it to the MSC/NASTRAN
solver (SOL103).

9b-4

Review the results.

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 9b

Shock Response Spectrum

ID SEMINAR,PROB9B
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________

CEND
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
MSC/NASTRAN 102 Exercise Workbook

9b-5

BEGIN BULK
1

9b-6

MSC/NASTRAN 102 Exercise Workbook

10

WORKSHOP 9b

Shock Response Spectrum

10

ENDDATA
MSC/NASTRAN 102 Exercise Workbook

9b-7

Generating an input file for MSC/NASTRAN Users:


MSC/NASTRAN users can generate an input file using the data from
page 9-3 (Model Description). The result should be similar to the
output below.

1.

MSC/NASTRAN input file: prob9b.dat.

ID SEMINAR, PROB9B
SOL 103
TIME 30
CEND
TITLE= RESPONSE SPECTRUM ANALYSIS
SUBTITLE= USING CALCULATED SHOCK RESPONSE
LABEL= SHOCK WILL BE INPUT IN Z DIRECTION
ECHO= UNSORTED
SET 111= ALL
DISPLACEMENT= 111
SPC= 200
SUBCASE 1
METHOD= 100
SDAMP= 200
DLOAD= 500
$
BEGIN BULK
$
$ PLATE MODEL DESCRIBED IN NORMAL MODES EXAMPLE
$
INCLUDE plate.bdf
PARAM,COUPMASS,1
PARAM,WTMASS,0.00259
$
$ BOUNDARY CONDITIONS FOR `CLAMPED MODES
$
SPC1, 200, 1245, 1, 12, 23, 34, 45
$
$ PLACE BIG FOUNDATION MASS (BFM) AT BASE
$ TO STIMULATE `CLAMPED MODES
$
CMASS2, 110, 1000., 23, 3
$
$ RBE MASS TO REMAINING BASE POINTS
$
RBE2, 101, 23, 3, 1, 12, 34, 45
9b-8

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 9b

Shock Response Spectrum

$
$ SUPPORT CARD TO IDENTIFY EXCITATION DOFS
$
SUPORT, 23, 3
$
$ EIGENVALUE EXTRACTION
$ MUST BE MASS NORMALIZED (DEFAULT)
$
EIGR, 100, MGIV, 0., 1000.
$
$ TABLE TO SPECIFY DAMPING FOR USE IN THE ANALYSIS
$
TABDMP1, 200, CRIT,
, 0., 0.03, 1000., 0.03, ENDT
$
$ SPECIFICATION OF SHOCK SPECTRUM TO BE USED
$
DLOAD, 500, 1.0, 2.0, 1
$
$ DLOAD, ID, OVERALL SCALE, SCALE FOR R-SET DOF# 1, SHOCK TABLE FOR DOF# 1,
$ SCALE FOR R-SET DOF# 2, SHOCK TABLE FOR DOF# 2, ETC.
$
$ SELECT SHOCK RESPONSE CALCULATION
$
PARAM, SCRSPEC, 0
$
$ SELECT SUMMATION OPTION
$
PARAM, OPTION, SRSS
$
$ MODAL FREQUENCY RANGE CAN BE SELECTED USING
PARAM, LFREQ, 0.1
PARAM, HFREQ, 1000.
$
$ SPECIFICATION FOR SHOCK TABLES
$
DTI, SPECSEL, 0
DTI, SPECSEL, 1, , A, 2, 0., 3, 0.02,
, 4, 0.04, ENDREC
$
$ DTI, SPECSEL, SHOCK TABLE NUMBER, , [(A)CCELERATION, (V)ELOCITY, OR (D)ISP],
$ TABLED1 POINTER, DAMPING FOR TABLE, ETC.
$
MSC/NASTRAN 102 Exercise Workbook

9b-9

$ PUNCH OUTPUT FOR SHOCK SPECTRUM CALCULATION


$
$ ACCE
4
3000
3
1
$
0.000000E+00
$
TABLED1 2
20.
.038683 40.
.152539 60.
100.
.862049 120.
1.17619 140.
180.
2.11404 200.
2.36801 220.
260.
2.76275 280.
2.75073 300.
340.
2.4218 360.
2.39068 380.
420.
1.78538 440.
1.70355 460.
500.
1.22608 520.
1.20483 540.
580.
1.10048 600.
1.05582 620.
660.
.908725 680.
.859158 700.
740.
.728996 760.
.694088 780.
820.
.598496 840.
.571831 860.
900.
.528854 920.
.509281 940.
980.
.488793 1000.
.468321 ENDT
$ACCE
4
3000
3
52
$
2.0000000E-02
TABLED1 3
20.
.037708 40.
.143365 60.
100.
.80976 120.
1.10506 140.
180.
1.98167 200.
2.22217 220.
260.
2.56231 280.
2.55577 300.
340.
2.29411 360.
2.25956 380.
420.
1.68656 440.
1.61355 460.
500.
1.19796 520.
1.17707 540.
580.
1.07807 600.
1.03637 620.
660.
.900171 680.
.854434 700.
740.
.727923 760.
.691288 780.
820.
.602091 840.
.579207 860.
900.
.526973 920.
.51329 940.
980.
.474171 1000.
.459408 ENDT
$ACCE
4
3000
3
103
$
4.0000000E-02
TABLED1 4
20.
.039336 40.
.137673 60.
100.
.764891 120.
1.04406 140.
180.
1.85678 200.
2.10175 220.
260.
2.39929 280.
2.42782 300.
340.
2.17923 360.
2.14283 380.
420.
1.62279 440.
1.53417 460.
500.
1.17212 520.
1.15165 540.

9b-10

MSC/NASTRAN 102 Exercise Workbook

.33511
1.50169
2.56617
2.74632
2.24931
1.57056
1.17631
1.00818
.827667
.668602
.563072
.500534

80.
160.
240.
320.
400.
480.
560.
640.
720.
800.
880.
960.

.576059
1.82018
2.70027
2.61887
2.02296
1.40493
1.14097
.958761
.782127
.635044
.550499
.498016

.314936
1.40671
2.35249
2.58668
2.12901
1.4968
1.14947
.992124
.810016
.658039
.559362
.500403

80.
160.
240.
320.
400.
480.
560.
640.
720.
800.
880.
960.

.541342
1.69567
2.53055
2.45921
1.92605
1.35263
1.11613
.946383
.767647
.628311
.542128
.487602

.297382
1.31588
2.19165
2.44263
2.0227
1.43168
1.12513

80.
160.
240.
320.
400.
480.
560.

.511244
1.58461
2.3921
2.317
1.8407
1.30597
1.09349

WORKSHOP 9b

580.
660.
740.
820.
900.
980.

Shock Response Spectrum


1.05768
.892143
.731968
.611319
.532205
.469568

600.
680.
760.
840.
920.
1000.

1.01868
.849752
.69746
.588261
.516369
.452243

620.
700.
780.
860.
940.
ENDT

.977462
.808538
.665814
.567655
.501047

640.
720.
800.
880.
960.

.934986
.769114
.637115
.549125
.485644

$
ENDDATA

MSC/NASTRAN 102 Exercise Workbook

9b-11

Submitting the input file for analysis:


2. Submit the input file to MSC/NASTRAN for analysis.
To submit the MSC/NASTRAN .dat file, find an available UNIX shell window
and at the command prompt enter nastran prob9b scr=yes. Monitor the run
using the UNIX ps command.
3. When the run is completed, edit the prob9b.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.

9b-12

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 9b

Shock Response Spectrum

Comparison of Results
4. Compare the results obtained in the .f06 file with the results on the
following page:

MSC/NASTRAN 102 Exercise Workbook

9b-13

9b-14

*** USER INFORMATION MESSAGE 3035 FOR DATA BLOCK


SUPPORT PT.NO.
EPSILON
1
9.5369362E-15
RESPONSE SPECTRUM ANALYSIS
USING CALCULATED SHOCK RESPONSE
SHOCK WILL BE INPUT IN Z DIRECTION

KLR

STRAIN
ENERGY
2.2118911E-09

EPSILONS LARGER THAN 0.001 ARE FLAGGED WITH ASTERISKS


JULY

5, 1997

MSC/NASTRAN

1/23/97

PAGE

12

SUBCASE 1

MSC/NASTRAN 102 Exercise Workbook

*** USER

INFORMATION MESSAGE 4415,


THE FOLLOWING A-SET DEGREES OF FREEDOM HAVE EITHER NULL MASSES OR NULL MASSES AND STIFFNESS.
RESPONSE SPECTRUM ANALYSIS
JULY
5, 1997 MSC/NASTRAN
1/23/97
PAGE
13
USING CALCULATED SHOCK RESPONSE
SHOCK WILL BE INPUT IN Z DIRECTION
SUBCASE 1

R E A L
E I G E N V A L U E S
EXTRACTION
EIGENVALUE
RADIANS
CYCLES
GENERALIZED
GENERALIZED
ORDER
MASS
STIFFNESS
1
101
0.0
0.0
0.0
1.000000E+00
0.0
2
102
7.058213E+05
8.401317E+02
1.337111E+02
1.000000E+00
7.058213E+05
3
103
1.878432E+07
4.334088E+03
6.897916E+02
1.000000E+00
1.878432E+07
4
104
2.811620E+07
5.302471E+03
8.439145E+02
1.000000E+00
2.811620E+07
5
105
1.931709E+08
1.389859E+04
2.212030E+03
0.0
0.0
6
1
2.234434E+08
1.494802E+04
2.379052E+03
0.0
0.0
7
106
2.328846E+08
1.526056E+04
2.428793E+03
0.0
0.0
8
107
6.845925E+08
2.616472E+04
4.164244E+03
0.0
0.0
9
108
9.602310E+08
3.098759E+04
4.931829E+03
0.0
0.0
10
109
1.365612E+09
3.695418E+04
5.881441E+03
0.0
0.0
^^^ USER INFORMATION MESSAGE 9047 (POSTREIG) - SCALED RESPONSE SPECTRA FOR RESIDUAL STRUCTURE ONLY
RESPONSE SPECTRUM ANALYSIS
JULY
5, 1997 MSC/NASTRAN
1/23/97
PAGE
20
USING CALCULATED SHOCK RESPONSE
SHOCK WILL BE INPUT IN Z DIRECTION
SUBCASE 1
MATRIX FN
(GINO NAME 101 ) IS A DB PREC
1 COLUMN X
3 ROW RECTANG MATRIX.
COLUMN
1
ROWS
1 THRU
3
-------------------------------------------------ROW
1)
1.337111500777D+02 6.897918198043D+02 8.439147827213D+02
THE NUMBER OF NON-ZERO TERMS IN THE DENSEST COLUMN =
3
THE DENSITY OF THIS MATRIX IS 100.00 PERCENT.
MODE
NO.

JULY

-3-

1/23/97

PAGE

VALUE

POINT

VALUE

POINT

JULY

5, 1997

T A B L E
( I N T E R N A L
R
DISPLACEMENT SET
-4-5-6-

1=
23-3
SCALED SPECTRAL RESPONSE, SRSS OPTION, DLOAD =

500

CLOSE =

VALUE

POINT

MSC/NASTRAN

S E Q U E N C E ,
-7-

-8-

VALU

1/23/97

R O W
-9-

PAGE
SUBCAS
S O R
-

1.00

9b-15

RESPONSE SPECTRUM ANALYSIS


JULY
5, 1997 MSC/NASTRAN
1/23/97
PAGE
USING CALCULATED SHOCK RESPONSE
SHOCK WILL BE INPUT IN Z DIRECTION
SUBCAS
MATRIX UHVR
(GINO NAME 101 ) IS A REAL
3 COLUMN X
3 ROW SQUARE
MATRIX.
COLUMN
1
ROWS
1 THRU
3
-------------------------------------------------ROW
1)
7.6201E-08 0.0000E+00 4.8912E-10
COLUMN
2
ROWS
1 THRU
3
-------------------------------------------------ROW
1)
6.4019E-05 0.0000E+00 2.5935E-06
COLUMN
3
ROWS
1 THRU
3
-------------------------------------------------ROW
1)
5.3784E-02 0.0000E+00 1.3752E-02
THE NUMBER OF NON-ZERO TERMS IN THE DENSEST COLUMN =
2
THE DENSITY OF THIS MATRIX IS 66.67 PERCENT.
RESPONSE SPECTRUM ANALYSIS
JULY
5, 1997 MSC/NASTRAN
1/23/97
PAGE

Shock Response Spectrum

MSC/NASTRAN 102 Exercise Workbook

-2-

MSC/NASTRAN

SUBCAS

COLUMN
1
23 T3 -2.11560E-02
COLUMN
2
23 T3 -4.40398E-16
COLUMN
3
23 T3 1.18599E-02
RESPONSE SPECTRUM ANALYSIS
USING CALCULATED SHOCK RESPONSE
SHOCK WILL BE INPUT IN Z DIRECTION
U S E T
D E F I N I T I O N
-1-

5, 1997

WORKSHOP 9b

RESPONSE SPECTRUM ANALYSIS


USING CALCULATED SHOCK RESPONSE
SHOCK WILL BE INPUT IN Z DIRECTION
PSIT
POINT
VALUE
POINT

9b-16

USING CALCULATED SHOCK RESPONSE


SHOCK WILL BE INPUT IN Z DIRECTION
TIME = 0.000000E+00

SUBCASE 1
D I S P L A C E M E N T

MSC/NASTRAN 102 Exercise Workbook

POINT ID.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
.
.

TYPE
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G
G

T1
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0

T2
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0

T3
6.222643E-10
7.374101E-08
3.173752E-07
7.194814E-07
1.246589E-06
1.870825E-06
2.565942E-06
3.308983E-06
4.080410E-06
4.865109E-06
5.653596E-06
6.222643E-10
9.569825E-08
3.577567E-07
7.614155E-07
1.284521E-06
1.901587E-06
2.589113E-06
3.325147E-06
4.091029E-06
4.872190E-06
5.659176E-06
6.222643E-10
9.883527E-08
3.680177E-07
7.745986E-07
1.296574E-06
1.911807E-06
2.596946E-06
3.330809E-06

V E C T O R
R1
0.0
8.308994E-08
1.214202E-07
1.167081E-07
1.039592E-07
8.258849E-08
6.149954E-08
4.203442E-08
2.671514E-08
1.691765E-08
1.301237E-08
0.0
1.430303E-08
4.357941E-08
5.321706E-08
4.871796E-08
4.101944E-08
3.141624E-08
2.263403E-08
1.555984E-08
1.101371E-08
8.953050E-09
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0

R2
0.0
3.059527E-07
6.616859E-07
9.376687E-07
1.160147E-06
1.326412E-06
1.444516E-06
1.519611E-06
1.559985E-06
1.575028E-06
1.577945E-06
0.0
3.703413E-07
6.684602E-07
9.362424E-07
1.147647E-06
1.312728E-06
1.430590E-06
1.507938E-06
1.551409E-06
1.570698E-06
1.576591E-06
0.0
3.798185E-07
6.836093E-07
9.334775E-07
1.145700E-06
1.307858E-06
1.426281E-06
1.503975E-06

R3
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0

10-2

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 10

Random Analysis

Model Description:
For the plate model, enforce a base motion in the z-direction described
by the following power spectral density, (PSD).
Autospectra of the Base Excitation

large mass
55

33

Frequency
(Hz)

G2/Hz

20

0.1

30

100

500

0.1

1000

0.1

Use the modal method with a large mass attached to the edge with an
RBE2 entry.
Below is a finite element representation of the flat plate. It also
contains the loads and boundary constraints.
Determine:
The response displacement and acceleration PSD at the
drive location, (the large mass).
The displacement PSD at the corner and center of the free
edge, (Grids 33 and 55).
Use modal solution.
Assume a constant critical damping ratio of 3% across the
whole frequency range.

MSC/NASTRAN 102 Exercise Workbook

10-3

Suggested Exercise Steps:


Reference a previously created dynamic math model,
plate.bdf, by using the INCLUDE statement.
Attach the large mass to the edge of the plate (CONM2 and
RBE2).
Specify modal damping as a tabular function of natural
frequency (TABDMP1).
Define the frequency-varying tip load (DAREA and
RLOAD2).
Define a set of frequencies to be used in the solution (FREQ,
FREQ1, and FREQ 4).
Specify Spectral Density (RANDPS and TABRND1).
Prepare the model for a direct transient analysis (SOL 111).
Request acceleration responses at base, tip center, and
opposite corner.
Generate an input file and submit it to the MSC/NASTRAN
solver for direct transient analysis.

10-4

Review the results.

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 10

Random Analysis

ID SEMINAR,PROB10
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________

CEND
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
MSC/NASTRAN 102 Exercise Workbook

10-5

BEGIN BULK
1

10-6

MSC/NASTRAN 102 Exercise Workbook

10

WORKSHOP 10

Random Analysis

10

ENDDATA
MSC/NASTRAN 102 Exercise Workbook

10-7

Generating an input file for MSC/NASTRAN Users:


MSC/NASTRAN users can generate an input file using the data from
pages 10-3 (general model description). The result should be similar
to the output below.

1.

MSC/NASTRAN input file: prob10.dat.

ID SEMINAR, PROB10
SOL 111
TIME 30
CEND
TITLE= RANDOM ANALYSIS - BASE EXCITATION
SUBTITLE= USING THE MODAL METHOD WITH LANCZOS
ECHO= UNSORTED
SPC= 101
SET 111= 33, 55, 9999
ACCELERATION(SORT2, PHASE)= 111
METHOD= 100
FREQUENCY= 100
SDAMPING= 100
RANDOM= 100
DLOAD= 100
$
OUTPUT(XYPLOT)
XTGRID= YES
YTGRID= YES
XBGRID= YES
YBGRID= YES
YTLOG= YES
XTITLE= FREQUENCY
YTTITLE= ACCEL RESPONSE BASE, MAGNITUDE
YBTITLE= ACCEL RESPONSE AT BASE, PHASE
XYPLOT ACCEL RESPONSE / 9999 (T3RM, T3IP)
YTTITLE= ACCEL RESPONSE AT TIP CENTER, MAGNITUDE
YBTITLE= ACCEL RESPONSE AT TIP CENTER, PHASE
XYPLOT ACCEL RESPONSE / 33 (T3RM, T3IP)
YTTITLE= ACCEL RESPONSE AT OPPOSITE CORNER, MAGNITUDE
YBTITLE= ACCEL RESPONSE AT OPPOSITE CORNER, PHASE
XYPLOT ACCEL RESPONSE / 55 (T3RM, T3IP)
$
$ PLOT OUTPUT IS ONLY MEANS OF VIEWING PSD DATA
$
XGRID= YES
YGRID= YES
10-8

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 10

Random Analysis

XLOG= YES
YLOG= YES
YTITLE= ACCEL P S D AT LOADED CORNER
XYPLOT ACCEL PSDF / 9999(T3)
YTITLE= ACCEL P S D AT TIP CENTER
XYPLOT ACCEL PSDF / 33(T3)
YTITLE= ACCEL P S D AT OPPOSITE CORNER
XYPLOT ACCEL PSDF / 55(T3)
$
BEGIN BULK
PARAM,COUPMASS,1
PARAM,WTMASS,0.00259
$
INCLUDE plate.bdf
$
GRID, 9999, , 0., 0., 0.
$
RBE2, 101, 9999, 12345, 1, 12, 23, 34, 45
$
SPC1, 101, 12456, 9999
$
CONM2, 6000, 9999, , 1.0E8
$
$ EIGENVALUE EXTRACTION PARAMETERS
$
EIGRL, 100 , , 2000.
$
$ SPECIFY MODAL DAMPING
$
TABDMP1, 100, CRIT,
+, 0., .03, 10., .03, ENDT
$
$ POINT LOADING AT TIP CENTER
$
RLOAD2, 100, 600, , , 310
$
TABLED1, 310,
+, 10., 1., 1000., 1., ENDT
$
DAREA, 600, 9999, 3, 1.E8
$
$ SPECIFY FREQUENCY STEPS
$
MSC/NASTRAN 102 Exercise Workbook

10-9

FREQ,100,30.
FREQ1,100,20.,20.,50
FREQ4,100,20.,1000.,.03,5
$
$ SPECIFY SPECTRAL DENSITY
$
RANDPS, 100, 1, 1, 1., 0., 111
$
TABRND1, 111,LOG,LOG
+, 20., 0.1, 30., 1., 100., 1., 500., .1,
+, 1000., .1, ENDT
$
ENDDATA

10-10

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 10

Random Analysis

Submitting the input file:


2. Submit the input file to MSC/NASTRAN for analysis.
To submit the MSC/NASTRAN .dat file, find an available UNIX shell
window and at the command prompt enter nastran prob10 scr=yes.
Monitor the run using the UNIX ps command.
3.

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


postscript file, prob10.ps, from the binary plot file
prob10.plt. The nonlinear force and displacement plots
are shown on the following pages.

4.

When the run is completed, edit the prob10.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.

MSC/NASTRAN 102 Exercise Workbook

10-11

Comparison of Results
5.

10-12

Compare the plot made from the exercise with the plots on
the following pages.

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 10

Random Analysis

Figure 10.1

MSC/NASTRAN 102 Exercise Workbook

10-13

Figure 10.2

10-14

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 10

Random Analysis

Figure 10.3

MSC/NASTRAN 102 Exercise Workbook

10-15

Figure 10.4

10-16

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 10

Random Analysis

Figure 10.5

MSC/NASTRAN 102 Exercise Workbook

10-17

Figure 10.6

10-18

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP PROBLEM 11

Random Analysis

Objectives:
Define a frequency-varying excitation.
Produce a MSC/NASTRAN input file from a dynamic
math model created in Workshop 1.
Submit the file for analysis in MSC/NASTRAN.
Compute nodal displacements for desired frequency
domain.

MSC/NASTRAN 102 Exercise Workbook11-1

11-2

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 11

Random Analysis

Model Description:
Using the modal method, determine the displacement response
spectrum of the tip center point due to the input spectrum of the
pressure and point loads listed below. Solve using the complex matrix
representation [Sab] for the cross spectrum.
Table 11.1
Autospectra of Pressure Load

Autospectra of Corner
Load

Frequency (Hz)

psi2/Hz

Frequency (Hz)

lb2/Hz

20

0.1

20

0.5

30

30

2.5

100

500

2.5

500

0.1

1000

1000

0.1

Table 11.2
Cross-Spectrum of Pressure and Corner Loads Real/
Imaginary
Frequency (Hz)

Real Part

Imaginary Part

20

-0.099619

0.007816

100

-0.498097

0.043579

500

0.070711

-0.070711

1000

MSC/NASTRAN 102 Exercise Workbook 11-3

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


contains the loads and boundary constraints.
Figure 11.1-Loads and Boundary Conditions
Unit psi over the total surface

1.0

11-4

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 11

Random Analysis

Suggested Exercise Steps:


Reference a previously created dynamic math model,
plate.bdf, 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).
Specify spectral density (RANDPS and TABRND1).
Prepare the model for a random analysis (SOL 111).
Request displacement response at loaded corner, tip center,
and opposite corner.
Generate an input file and submit it to the MSC/NASTRAN
solver for random analysis.

Review the results, specifically the nodal displacements.

MSC/NASTRAN 102 Exercise Workbook 11-5

ID SEMINAR,PROB11
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________

CEND
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________

11-6

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 11

Random Analysis

BEGIN BULK
1

MSC/NASTRAN 102 Exercise Workbook 11-7

10

ENDDATA

11-8

MSC/NASTRAN 102 Exercise Workbook

10

WORKSHOP 11

Random Analysis

Generating an input file for MSC/NASTRAN Users:


MSC/NASTRAN users can generate an input file using the data from
pages 11-3 (general model description). The result should be similar
to the output below.
1.

MSC/NASTRAN input file: prob11.dat.

ID SEMINAR, PROB11
SOL 111
TIME 30
CEND
TITLE= FREQUENCY RESPONSE WITH PRESSURE AND POINT LOADS
SUBTITLE= USING THE MODAL METHOD WITH LANCZOS
ECHO= UNSORTED
SPC= 1
SET 111= 11, 33, 55
DISPLACEMENT(SORT2, PHASE)= 111
METHOD= 100
FREQUENCY= 100
SDAMPING= 100
RANDOM= 100
SUBCASE 1
LABEL= PRESSURE LOAD
DLOAD= 100
LOADSET= 100
SUBCASE 2
LABEL CORNER LOAD
DLOAD= 200
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
MSC/NASTRAN 102 Exercise Workbook 11-9

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
XYPLOT DISP RESPONSE / 55 (T3RM, T3IP)
$
$ PLOT OUTPUT IS ONLY MEANS OF VIEWING PSD DATA
$
XGRID= YES
YGRID= YES
XLOG= YES
YLOG= YES
YTITLE= DISP P S D AT LOADED CORNER
XYPLOT DISP PSDF / 11(T3)
YTITLE= DISP P S D AT TIP CENTER
XYPLOT DISP PSDF / 33(T3)
YTITLE= DISP P S D AT OPPOSITE CORNER
XYPLOT DISP PSDF / 55(T3)
$
BEGIN BULK
$
PARAM,COUPMASS,1
$
PARAM,WTMASS,0.00259
$
$ MODEL DESCRIBED IN NORMAL MODES EXAMPLE
$
INCLUDE plate.bdf
$
$ EIGENVALUE EXTRACTION PARAMETERS
$
EIGRL, 100, 10., 2000.
$
$ SPECIFY MODAL DAMPING
$
TABDMP1, 100, CRIT,
+, 0., .03, 10., .03, ENDT
$
$ FIRST LOADING
$
RLOAD2, 100, 300, , , 310
$
TABLED1, 310,
+, 10., 1., 1000., 1., ENDT
$

11-10

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 11

Random Analysis

$ UNIT PRESSURE LOAD TO PLATE


$
LSEQ, 100, 300, 400
$
PLOAD2, 400, 1., 1, THRU, 40
$
$ SECOND LOADING
$
RLOAD2, 200, 600, , , 310
$
$ POINT LOAD AT TIP CENTER
$
DAREA, 600, 11, 3, 1.
$
$ SPECIFY FREQUENCY STEPS
$
FREQ1, 100, 20., 20., 49
$
$ SPECIFY SPECTRAL DENSITY
$
RANDPS, 100, 1, 1, 1., 0., 100
RANDPS, 100, 2, 2, 1., 0., 200
RANDPS, 100, 1, 2, 1., 0., 300
RANDPS, 100, 1, 2, 0., 1.0, 400
$
TABRND1, 100,
+, 20., 0.1, 30., 1., 100., 1., 500., .1,
+, 1000., .1, ENDT
$
TABRND1, 200,
+, 20., 0.5, 30., 2.5, 500., 2.5, 1000., 0.,
+, ENDT
$
TABRND1, 300,
+, 20., -.099619, 100., -.498097, 500., .070711, 1000., 0.,
+, ENDT
$
TABRND1, 400,
+, 20., .0078158, 100., .0435791, 500., -.70711, 1000., 0.,
+, ENDT
$
ENDDATA

MSC/NASTRAN 102 Exercise Workbook 11-11

Submitting the input file for analysis:


2.

Submit the input file to MSC/NASTRAN for analysis.


2a. To submit the MSC/NASTRAN .dat file, find an available
UNIX shell window and at the command prompt enter
nastran prob11 scr=yes. Monitor the run using the UNIX
ps command.

3.

11-12

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


postscript file, prob11.ps, from the binary plot file
prob11.plt. Compare the results with the plots below.

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 11

Random Analysis

Figure 11.2

MSC/NASTRAN 102 Exercise Workbook 11-13

Figure 11.3

11-14

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 11

Random Analysis

Figure 11.4

MSC/NASTRAN 102 Exercise Workbook 11-15

Figure 11.5

11-16

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 11

Random Analysis

Figure 11.6

MSC/NASTRAN 102 Exercise Workbook 11-17

Figure 11.7

11-18

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 11

Random Analysis

Figure 11.8

MSC/NASTRAN 102 Exercise Workbook 11-19

Figure 11.9

11-20

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 11

Random Analysis

Figure 11.10

MSC/NASTRAN 102 Exercise Workbook 11-21

11-22

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP PROBLEM 12

Complex Modes of a Pile


Driver

EXCITER

PILE

GROUND

Objectives
Define complex eigenvalue extraction parameters.
Submit the file for analysis in MSC/NASTRAN.
Compute complex modes.

MSC/NASTRAN 102 Exercise Workbook

12-1

12-2

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 12

Complex Modes of a Pile Driver

Model Description:
The model is idealized as shown below in Figure 12.1. (Note that both a
spring element and a damper element will be created connected Grid 2 and
Grid 3.)
Figure 12.1-Model Description

m1

m2

K1

K2

C2

Table 12.1
m1

3.0 lb-sec2/in

m2

1.5 lb-sec2/in

K1

50,000 lb/in

K2

12,500 lb/in

C2

30 lb-sec/in

MSC/NASTRAN 102 Exercise Workbook

12-3

Suggested Exercise Steps


Generate an input file and submit it to the MSC/NASTRAN
solver for complex eigenvalue analysis.
Generate a finite element representation of the pile driver using
GRID, CONM2, CELAS, and CVISC elements.
Define material (MAT1), and element (PELAS) and (PVISC)
properties.
Apply x-direction boundary constraint (SPC1).
Specify complex eigenvalue extraction parameters
(CMETHOD) and (EIGC).
Prepare the model for complex eigenvalue analysis (SOL107).

12-4

Review the results, specifically the complex eigenvalues.

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 12

Complex Modes of a Pile Driver

ID SEMINAR,PROB12
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________

CEND
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
BEGIN BULK
MSC/NASTRAN 102 Exercise Workbook

12-5

12-6

MSC/NASTRAN 102 Exercise Workbook

10

WORKSHOP 12

Complex Modes of a Pile Driver

10

ENDDATA
MSC/NASTRAN 102 Exercise Workbook

12-7

Exercise Procedure:
1. Users who are not utilizing MSC/PATRAN for generating an input file
should go to Step 16, otherwise, proceed to step 2.
2. Create a new database and named prob12.db
File/New Database
New Database Name

prob12

OK
In the New Model Preference form set the following:
Tolerance

Default

Analysis code:

MSC/NASTRAN

OK
3. Create the model by the edit method in Finite Elements.
Finite Elements
Action:

Create

Object:

Node

Method:

Edit

Associate with Geometry


Node Location List

[0 0 0]

Apply
Turn on the label and increase the node size by using the Quick Pick
buttons.
Show Label

12-8

MSC/NASTRAN 102 Exercise Workbook

Node Size

WORKSHOP 12

Complex Modes of a Pile Driver

4. Similarly, create Nodes 2 and 3.


Node

Location

Node 2

[ 1 0 0]

Node 3

[2 0 0]

5. Create the Bar Element for Node 1 and Node 2.

Finite Element
Action:

Create

Object:

Element

Method:

Edit

Shape

Bar

Node 1 =

Node 1

Node 2 =

Node 2

Apply
6.

Similarly, create the 2nd bar element by:

Node 1 =

Node 2

Node 2 =

Node 3

Apply
7.

Create the 2 mass elements at Node 1 and Node 2.

Finite Element
Action:

Create

Object:

Element

Method:

Edit

Shape:

Point

Element ID List

Node 1 =

Node 1

Apply
MSC/NASTRAN 102 Exercise Workbook

12-9

Element ID List

Node 1 =

Node 2

Apply
8.

Create the damper elements connecting Node 2 and Node 3.

Finite Element
Action:

Create

Object:

Element

Method:

Edit

Shape

Bar

Node 1 =

Node 2

Node 2 =

Node 3

Apply
9.

Create Element Properties, (spring constant).

Properties
Action:

Create

Dimension:

1D

Type:

Spring

Property Set Name:

spring1

Input Properties ...


Spring Constant:

50000

DOF at Node 1:

UX

DOF at Node 2:

UX

OK
Application Region
(In the select menu, select the Beam Element
filter.)

12-10

MSC/NASTRAN 102 Exercise Workbook

Element 1

WORKSHOP 12

Complex Modes of a Pile Driver

Beam Element

Add
Apply
10.

Similarly, create the spring constant of 12,500 for the 2nd


spring element.

Property Set Name:

spring2

Input Properties ...


Spring Constant:

12500

DOF at Node 1:

UX

DOF at Node 2:

UX

OK
Element 2

Application Region
Add
Apply
11.

Create Element Properties, (damper coefficient), for the


damper element:

Properties
Action:

Create

Dimension:

1D

Type:

Damper
damper

Property Set Name:


Option(s)...

Viscous

Input Properties ...


[Ext. Viscous Coeff.]

30

OK
Application Region

Element 5
MSC/NASTRAN 102 Exercise Workbook

12-11

Add
Apply
12.

Create the mass properties of the mass elements.

Properties
Action:

Create

Dimension:

0D

Type:

Mass

Property Set Name:

mass1

Option(s):

Lumped

Input Properties ...


3

Mass:
OK

Element 3

Application Region
(In the select menu, select the Point Element
filter.)

Point Element

Add
Apply
13.

Similarly, create the mass property of the 2nd mass element:

Properties

12-12

Action:

Create

Dimension:

0D

Type:

Mass

Property Set Name:

mass2

Option:

Lumped

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 12

Complex Modes of a Pile Driver

Input Properties ...


1.5

Mass:
OK

Element 4

Application Region
Add
Apply
14.

Create the constraint at the ground, Node 3.

Load/BCs
Action:

Create

Object:

Displacement

Type:

Nodal
constraint

New Set Name:


Input Data...
Translations < T1 T2 T3 >

< 0, , >

OK
Select Application Region...

FEM
Node 3

Select Nodes:
Add
OK
Apply
15.

Create the analysis deck.

Analysis
Action:

Analyze

Object:

Entire Model

Method:

Analysis Deck

Job Name:

prob12
MSC/NASTRAN 102 Exercise Workbook

12-13

Solution Type...
Solution Type:

COMPLEX EIGENVALUES

Solution Parameters ...


Formulation

Direct

Complex Eigenvalue...
Number of Desired Roots=
OK
OK
OK
Apply

12-14

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 12

Complex Modes of a Pile Driver

Generating an input file for MSC/NASTRAN Users:


MSC/NASTRAN users can generate an input file using the data
previously stated. The result should be similar to the output below.
16. MSC/NASTRAN input file: prob12.dat
ID SEMINAR, PROB12
SOL 107
TIME 5
CEND
TITLE= TWO-DOF MODEL (IMAC 8, PG 891)
SUBTITLE= COMPLEX MODES
DISPLACEMENT= ALL $ DEFAULT= REAL, IMAGINARY
SPC= 100
CMETHOD= 99
$
BEGIN BULK
$
$ COMPLEX EIGENVALUE EXTRACTION PARAMETERS
$
EIGC, 99, HESS, , , , , 4
$
$ DEFINE GRIDS, MASSES, AND STIFFNESSES
$ GRID 1 = EXCITER (X=2, MASS=3) 50K STIFFNESS BETWEEN GRIDS 1 AND 2
$ GRID 2 = PILE (X=1, MASS=3)
12.5K STIFFNESS BETWEEN GRIDS 2 AND 3
$ GRID 3 = BASE (X=0, FIX BASE)
$
GRID, 1, , 2., 0., 0.
GRID, 2, , 1., 0., 0.
GRID, 3, , 0., 0., 0.
GRDSET, , , , , , , 23456
CELAS2, 1, 50000., 1, 1, 2, 1
CELAS2, 2, 12500., 2, 1, 3, 1
CONM2, 201, 1, , 3.0
CONM2, 202, 2, , 1.5
SPC, 100, 3, 1
$
$ DEFINE DAMPER OF 30 BETWEEN GRIDS 2 AND 3
$
CVISC, 101, 1, 2, 3
PVISC, 1, 30.
$
ENDDATA

MSC/NASTRAN 102 Exercise Workbook

12-15

Submitting the input file for analysis:


17.

Submit the input file to MSC/NASTRAN for analysis.


17a. To submit the MSC/PATRAN .bdf file, find an available UNIX
shell window. At the command prompt enter nastran
prob12.bdf scr=yes. Monitor the run using the UNIX ps
command.
17b. To submit the MSC/NASTRAN .dat file, find an available
UNIX shell window and at the command prompt enter nastran
prob12 scr=yes. Monitor the run using the UNIX ps command.

18.

When the run is completed, edit the prob12.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.

19.

While still editing prob12.f06, search for the word:


E I G E N V A L U E (spaces are necessary).

12-16

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 12

Complex Modes of a Pile Driver

Comparison of Results
20.

Compare the results obtained in the .f06 file with the results on
the following page:

MSC/NASTRAN 102 Exercise Workbook

12-17

12-18
MSC/NASTRAN 102 Exercise Workbook

1
2
3
4

EXTRACTION
ORDER
2
3
1
4

E I G E N V A L U E
EIGENVALUE
(REAL)
(IMAG)
-2.660969E+00
-4.983521E+01
-7.339031E+00
-2.360312E+02
-2.660969E+00
4.983521E+01
-7.339031E+00
2.360312E+02

C O M P L E X

T1
1.000000E+00
.0
8.514119E-01
1.591320E-02
.0
.0
.0
.0
.0
.0
.0
.0

T2
.0
.0
.0
.0
.0
.0

T3
.0
.0
.0
.0
.0
.0

T1
-4.241094E-01
-3.768431E-02
1.000000E+00
.0
.0
.0
.0
.0
.0
.0
.0
.0

T2
.0
.0
.0
.0
.0
.0

T3
.0
.0
.0
.0
.0
.0

T1
1.000000E+00
.0
8.514119E-01
-1.591320E-02
.0
.0

.0
.0
.0
.0
.0
.0

T2
.0
.0
.0
.0
.0
.0

T3
.0
.0
.0
.0
.0
.0

G
G

TYPE
G

POINT ID.
1

T1
-4.241094E-01
3.768431E-02
1.000000E+00
.0
.0
.0

.0
.0
.0
.0
.0
.0

T2
.0
.0
.0
.0
.0
.0

T3
.0
.0
.0
.0
.0
.0

COMPLEX EIGENVALUE = -7.339031E+00, 2.360312E+02


C O M P L E X
E I G E N V E C T O R
(REAL/IMAGINARY)

TYPE
G

POINT ID.
1

COMPLEX EIGENVALUE = -2.660969E+00, 4.983521E+01


C O M P L E X
E I G E N V E C T O R
(REAL/IMAGINARY)

TYPE
G

POINT ID.
1

COMPLEX EIGENVALUE = -7.339031E+00, -2.360312E+02


C O M P L E X
E I G E N V E C T O R
(REAL/IMAGINARY)

TYPE
G

POINT ID.
1

NO.

NO.

R1

NO.

R1

NO.

.0
.0
.0
.0
.0
.0

.0
.0
.0
.0
.0
.0

.0
.0
.0
.0
.0
.0

.0
.0
.0
.0
.0
.0

R2

R2

R2

R2

S U M M A R Y
FREQUENCY
(CYCLES)
7.931520E+00
3.756553E+01
7.931520E+00
3.756553E+01

R1

R1

COMPLEX EIGENVALUE = -2.660969E+00, -4.983521E+01


C O M P L E X
E I G E N V E C T O R
(REAL/IMAGINARY)

ROOT
NO.

.0
.0
.0
.0
.0
.0

.0
.0
.0
.0
.0
.0

.0
.0
.0
.0
.0
.0

.0
.0
.0
.0
.0
.0

R3

R3

R3

SUBCASE 1

SUBCASE 1

SUBCASE 1

R3

DAMPING
COEFFICIENT
1.067907E-01
6.218695E-02
1.067907E-01
6.218695E-02
SUBCASE 1

WORKSHOP 12

Complex Modes of a Pile Driver

21. MSC/NASTRAN Users have finished this exercise. MSC/


PATRAN Users should proceed to the next step.
22.

Proceed with the Reverse Translation process, that is


importing the prob12.op2 results file into MSC/PATRAN. To
do this, return to the Analysis form and proceed as follows.

Analysis
Action:

Read Output2

Object:

Result Entities

Method:

Translate

Select Results File...


Select Available Files

prob12.op2

OK
Apply
23.

View the results.

Results
Form Type:

Basic

Select Result Cases

<select one of the modes>

Select Deformation Result

2.1-Eigenvectors, Translation

MSC/NASTRAN 102 Exercise Workbook

12-19

12-20

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP PROBLEM 13

Nolins in
Linear Transient

Objectives:
Represent non-structural variables using non-structural
DOFs.
Define dynamic functions with transfer functions.
Create a nonlinear transient force.
Prepare a MSC/NASTRAN input file for a transient
analysis.
Visualize analysis results.
MSC/NASTRAN 102 Exercise Workbook

13-1

13-2

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 13

Nolins in Linear Transient

Model Description:
Figure 13.1-Car traveling over a speed bump

u
1

Velocity = 100 in/sec


3864 lbs

120

50

k:

c:

u 2.0in

197.4 lb/in

u < 2.0in

394.8 lb/in

u 0

1.88 lb/(in/sec)

u 0

1.88 lb/(in/sec) + 0.3 lb/(in/sec)2

MSC/NASTRAN 102 Exercise Workbook

13-3

Figure 13.2-Force vs Relative Velocity


F

-2.0

13-4

Urel.

MSC/NASTRAN 102 Exercise Workbook

Urel.

WORKSHOP 13

Nolins in Linear Transient

Suggested Exercise Steps:


Generate a finite element representation of the model using (GRID),
(CBAR), and (CELAS2) elements.
Define material (MAT1) and element (PBAR) properties.
Constraints to eliminate rigid-body modes (SPC1).
Define non-structural variables (CONM2).
Specify scalar damper property and connection (CDAMP2).
Define extra points (EPOINT).
Define dynamic transfer functions (TF).
Add nonlinear portion of the spring (NOLIN1).
Add nonlinear portion of the damper (NOLIN4).
Define the time-varying load (DAREA & TLOAD2).
Define time delay (DELAY).
Specify integration time step.
Prepare the model for a direct transient analysis (SOL109).
Request response in terms of nodal displacement and nonlinear load
output.
Generate an input file and submit it to the MSC/NASTRAN solver for
direct transient analysis.
Review the results, specially the xy plot of nodal displacements and
nonlinear load.

MSC/NASTRAN 102 Exercise Workbook

13-5

ID SEMINAR,PROB13
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________

CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK

13-6

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 13

Nolins in Linear Transient

MSC/NASTRAN 102 Exercise Workbook

10

13-7

ENDDATA
13-8

MSC/NASTRAN 102 Exercise Workbook

10

WORKSHOP 13

Nolins in Linear Transient

Generating an input file for MSC/NASTRAN Users:


1.

MSC/NASTRAN users can generate an input file using


the data from pages 13-3 and 13-4 (general model
description). The result should be similar to the output
below (prob13.dat):

ASSIGN OUTPUT2 = prob13.op2, UNIT=12


ID NAS102, WORKSHOP13
SOL 109
TIME 100
CEND
TITLE= SIMPLE CAR MODEL WITH NOLINEAR
SUBTITLE= SPRINGS AND DAMPERS RUNNING OVER A BUMP
LABEL= SOL 109, CONSTANT DELTA TIME
SEALL= ALL
SPC= 100
TFL= 100
NONLINEAR = 100
DLOAD = 100
TSTEP = 100
DISPLACEMENT(PLOT)= ALL
NLLOAD(PLOT)= ALL
$
OUTPUT(XYPLOT)
CSCALE=1.3
XAXIS= YES
YAXIS= YES
XGRID LINES= YES
YGRID LINES= YES
XTITLE= TIME (SEC)
YTITLE= VERTICAL DISPLACEMENT OF POINT 1
XYPLOT DISP/1(T2)
YTITLE= VERTICAL DISPLACEMENT OF POINT 2
XYPLOT DISP/2(T2)
YTITLE= VERTICAL DISPLACEMENT OF POINT 3
XYPLOT DISP/3(T2)
YTITLE= VERTICAL DISPLACEMENT OF POINT 4
XYPLOT DISP/4(T2)
YTITLE= VERTICAL DISPLACEMENT OF POINT 5
XYPLOT DISP/5(T2)
YTITLE= NONLINEAR FORCES AT POINT 1
XYPLOT NONLINEAR/1(T2)
MSC/NASTRAN 102 Exercise Workbook

13-9

YTITLE= NONLINEAR FORCES AT POINT 2


XYPLOT NONLINEAR/2(T2)
$
BEGIN BULK
PARAM,POST,-1
PARAM,PATVER,3.0
$
$ CARRIAGE POINTS
$
GRID, 1, , 0., 0., 0.
GRID, 2, , 120., 0., 0.
GRID, 5, , 60., 0., 0.
$
$ WHEEL POINTS
$
GRID, 3, , 0., -10., 0.
GRID, 4, , 120., -10., 0.
$
$ CAR CARRIAGE
$
CBAR, 5, 11, 1, 5, 0., 1., 0.
CBAR, 6, 11, 5, 2, 0., 1., 0.
PBAR, 11, 12, 10., 10., 10.
MAT1, 12, 3.0E+7, , .33
$
$ CONSTRAINTS TO ELIMINATE RIGID-BODY MODES
$
SPC1, 100, 1345, 1, 2, 5
SPC1, 100, 13456, 3, 4
$
$ SYSTEM WILL HAVE A NATURAL FREQUENCY OF 1 HZ
$ WITH CRITICAL DAMPING OF 1 PERCENT
$
CONM2, 10, 1, ,2.5
CONM2, 15, 2, ,2.5
CONM2, 20, 5, ,5.
$
CELAS2, 30, 197.4, 1, 2, 3, 2
CELAS2, 40, 197.4, 2, 2, 4, 2
$
CDAMP2, 50, 1.88, 1, 2, 3, 2
CDAMP2, 60, 1.88, 2, 2, 4, 2
$
$ DEFINE EXTRA POINTS TO HOLD DIFFERENCES
$ BETWEEN WHEELS AND CARRIAGE
13-10

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 13

Nolins in Linear Transient

$
EPOINT, 101, 102
$
$ USE TRANSFER FUNCTIONS TO TRACK DIFFERENCES
$ 101= V1 - V3
$ 102= V2 - V4
$
TF, 100, 101, 0, 1., 0., 0.,
, 1, 2, -1., 0., 0.,
, 3, 2, 1., 0., 0.
$
TF, 100, 102, 0, 1., 0., 0.,
, 2, 2, -1., 0., 0.,
, 4, 2, 1., 0., 0.
$
$ ADD NONLINEAR PORTION OF SPRINGS
$
NOLIN1, 100, 1, 2, 197.4, 101, 0, 111
NOLIN1, 100, 2, 2, 197.4, 102, 0, 111
TABLED2, 111, -2.0,
, -1., 1., 0., 0., 1., 0.,ENDT
$
$ ADD NONLINEAR PORTION OF DAMPERS
$
NOLIN4, 100, 1, 2, -0.3, 101, 10, 2.
NOLIN4, 100, 2, 2, -0.3, 102, 10, 2.
$
$ USE LAGRANGE MULTIPLIERS TO IMPOSE WHEEL DISPLACEMENT
$ 103= V3
$ 104= V4
$
EPOINT, 103, 104
$
TF, 100, 103, 0, 0., 0., 0.,
, 3, 2, 1., 0., 0.
TF, 100, 3, 2, 0., 0., 0.,
, 103, 0, 1., 0., 0.
$
TF, 100, 104, 0, 0., 0., 0.,
, 4, 2, 1., 0., 0.
TF, 100, 4, 2, 0., 0., 0.,
, 104, 0, 1., 0., 0.
$
MSC/NASTRAN 102 Exercise Workbook

13-11

$ MOVE WHEELS OVER BUMP


$
TLOAD2, 100, 222, 333, 0, 0., 0.5, 1., -90.
DAREA, 222, 103, 0, 4.
DAREA, 222, 104, 0, 4.
DELAY, 333, 104, 0, 1.2
$
$ INTEGRATION INFORMATION
TSTEP, 100, 200, .05, 1
$
ENDDATA

13-12

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 13

Nolins in Linear Transient

Submitting the input file for analysis


2.

Submit the input file to MSC/NASTRAN for analysis.


To submit the MSC/NASTRAN .dat file, find an available UNIX
shell window and at the command prompt enter nastran prob13
scr=yes. Monitor the run using the UNIX ps command.

3.

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


postscript file, prob13.ps, from the binary plot file
prob13.plt. The nonlinear force and displacement plots
are shown in figures 13.3 to 13.9.

4.

When the run is completed, edit the prob13.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.

5.

While still editing prob13.f06, search for the word:


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

MSC/NASTRAN 102 Exercise Workbook

13-13

Comparison of Results
6.

13-14

Compare the results obtained in the .f06 file with the


results on the following page:

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 13

Nolins in Linear Transient

MSC/NASTRAN 102 Exercise Workbook

13-15

13-16
MSC/NASTRAN 102 Exercise Workbook

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
ID
TYPE
NO. CURVE ID.
ALL DATA
ALL DATA
ALL DATA
YMIN
1 NONLIN
1
1( 4)
0.000000E+00 1.000002E+01 0.000000E+00 0.000000E+00
0.000000E+00 1.000002E+01 0.000000E+00 0.000000E+00
1 NONLIN
2
2( 4)
0.000000E+00 1.000002E+01 0.000000E+00 0.000000E+00
0.000000E+00 1.000002E+01 0.000000E+00 0.000000E+00
1 DISP
3
1( 4)
0.000000E+00 1.000002E+01 -2.836541E+00 2.199999E+00
0.000000E+00 1.000002E+01 -2.836541E+00 2.199999E+00
1 DISP
4
2( 4)
0.000000E+00 1.000002E+01 -2.671688E+00 1.999999E+00
0.000000E+00 1.000002E+01 -2.671688E+00 1.999999E+00
1 DISP
5
3( 4)
0.000000E+00 1.000002E+01 -4.440892E-16 5.000001E-01
0.000000E+00 1.000002E+01 -4.440892E-16 5.000001E-01
1 DISP
6
4( 4)
0.000000E+00 1.000002E+01 -1.332268E-15 1.749999E+00
0.000000E+00 1.000002E+01 -1.332268E-15 1.749999E+00
1 DISP
7
5( 4)
0.000000E+00 1.000002E+01 -2.150819E+00 2.149999E+00
0.000000E+00 1.000002E+01 -2.150819E+00 2.149999E+00

YMAX-FRAME/
X FOR
ALL DATA
YMAX
2.975151E+02 6.000001E-01
2.975151E+02 6.000001E-01
4.661461E+02 1.749999E+00
4.661461E+02 1.749999E+00
5.942877E+00 4.500000E-01
5.942877E+00 4.500000E-01
7.464218E+00 1.600000E+00
7.464218E+00 1.600000E+00
4.000000E+00 2.500000E-01
4.000000E+00 2.500000E-01
4.000000E+00 1.450000E+00
4.000000E+00 1.450000E+00
3.963917E+00 1.600000E+00
3.963917E+00 1.600000E+00

WORKSHOP 13

Nolins in Linear Transient

7.

MSC/NASTRAN Users have finished this exercise.


MSC/PATRAN Users should proceed to the next step.

8.

Proceed with the Reverse Translation process, that is


importing the prob13.op2 results file into MSC/PATRAN.
To do this, return to the Analysis form and proceed as
follows.

Analysis
Action:

Read Output2

Object:

Both

Method:

Translate

Select Results File...


Select Available Files

prob13.op2

OK
Apply
9.

Plot the results in xy plots.

The first plot is to make the Displacement versus Frequency plot at


Node 1.

Results
Form Type:
Select Result Cases

Advanced
(highlight all cases)

Get Results
Select Result

1.1-Nonlinear Applied Load

Plot Type:

XY Plot

Plot Type Options...


Result XY Plot Types

Results Versus Global Variables

Global Var...
Global Variables:

1-Time

Apply
MSC/NASTRAN 102 Exercise Workbook

13-17

Result (Y)...
Results:

1.1-Nonlinear Applied
Translational

Vector Component

Mag.

Loads,

OK
Node Ids

Node 1

Apply
New Title or Title Filter

Nonlinear Forces at Point 1

Rename
Apply
Figure 13.3-Nonlinear Forces at Point 1

13-18

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 13

Nolins in Linear Transient

Figure 13.4-Nonlinear Forces at Point 2

To get the displacement results, go all the way back to the Result
menu.
Select Result

2.1 - Displacement, Translational

Plot Type:

XY Plot

Plot Type Options...


Result XY Plot Types

Results Versus Global Variables

Global Var...
Global Variable:

1-time

Apply
Result (Y)...
Results:

2.1-Displacement, Translational

Vector Component

Component
X
Y
Z

OK
Node IDs

Node 1
MSC/NASTRAN 102 Exercise Workbook

13-19

Apply
New Title or Title Filter

Vertical Displacement of Point 1

Rename
Apply
Figure 13.5-Vertical Displacement of Point 1

Figure 13.6-Vertical Displacement of Point 2

13-20

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 13

Nolins in Linear Transient

Figure 13.7-Vertical Displacement of Point 3

Figure 13.8-Vertical Displacement of Point 4

MSC/NASTRAN 102 Exercise Workbook

13-21

Figure 13.9-Vertical Displacement of Point 5

13-22

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP PROBLEM 14a

Modal Analysis of a Beam

Objectives
Perform normal modes analysis of a cantilever beam.
Submit the file for analysis in MSC/NASTRAN.
Find the first three natural frequencies and mode shapes of
the beam.

MSC/NASTRAN 102 Exercise Workbook

14a-1

14a-2

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 14a

Modal Analysis of a Beam

Model Description:
The goal of this example is to find the first 3 modes of a beam pinned at
both ends.
Figure 14a.1 below is a finite element representation of the beam. One end
is contrained in all translation and the other is free to move in the X. Both
ends are held in the X-rotation.
Figure 14a.1-Grid Coordinates and Element Connectivities

MSC/NASTRAN 102 Exercise Workbook

14a-3

Figure 14a.2-Beam Cross Section


1.0 in
0.1 in

2.0 in

0.1 in

0.1 in
1.0 in

Table 14a.1

14a-4

Length

100 in

Height

2 in

Width

1 in

Thickness

0.100 in

Area

0.38 in2

I1

0.229 in4

I2

0.017 in4

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 14a

Modal Analysis of a Beam

Hand Calculations

K n EIg 1 2
f n = ------ ---------2 W l 4
1 10 106 ( 0.229 ) ( 386.4 )
f n = K n ------ ------------------------------------------------------ 2 ( 0.38 ) ( 0.101 ) ( 100 ) 4

12

f n = K n ( 2.417 )
* I of the strong axis is used since translational Z DOF has been constrained
by the permanent constraint.

From Theory
Mode

Kn

fn

9.87

23.85 Hz

39.5

95.46 Hz

88.8

214.59 Hz

MSC/NASTRAN 102 Exercise Workbook

14a-5

Suggested Exercise Steps


Explicitly generate a finite element representation of the beam structure.
(i.e., the grids (GRID) and element connectivities (CBAR) should be
defined manually.)
Define material (MAT1) and element (PBARL) properties.
Apply the fixed boundary constraints (SPC1).
Prepare the model for a normal modes analysis (SOL 103 and PARAMS).
PARAM, WTMASS, 0.00259
PARAM, COUPMASS, 1
EIGRL (To select Lanczos.)
Generate an input file and submit it to the MSC/NASTRAN solver for
normal modes analysis.

14a-6

Review the results, specifically the eigenvalues.

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 14a

Modal Analysis of a Beam

ID SEMINAR,PROB1
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK

MSC/NASTRAN 102 Exercise Workbook

14a-7

14a-8

MSC/NASTRAN 102 Exercise Workbook

10

WORKSHOP 14a

Modal Analysis of a Beam


4

ENDDATA

MSC/NASTRAN 102 Exercise Workbook

14a-9

10

Exercise Procedure:
1. Users who are not utilizing MSC/PATRAN for generating an input file
should go to Step 10, otherwise, proceed to step 2.
2. Create a new database named prob14a.db.
File/New Database
prob14a

New Database Name


OK

In the New Model Preference form set the following:


Tolerance

Default

Analysis Code:

MSC/NASTRAN

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

4. Create a curve.
Geometry
Action:

Create

Object:

Curve

Method:

XYZ

Vector Coordinates List

<100, 0, 0>

Origin Coordinates List

[ 0, 0, 0]

Apply

14a-10

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 14a

Modal Analysis of a Beam

5. Create the finite element model and mesh the surface.


Finite Elements
Action:

Create

Object:

Mesh

Type:

Curve

Global Edge Length

10

Curve List

Curve 1

Apply
6. Create nodal displacements.
Loads/BCs
Action:

Create

Object:

Displacement

Type:

Nodal
disp1

New Set Name


Input Data...
Translations <T1 T2 T3>

<0 0 0>

Rotations <R1 R2 R3>

<0 , , >

OK
Select Application Region...

Geometry
Select Geometry Entities

Point 1

Add
OK
Apply
disp2

New Set Name


Input Data...
Translations <T1 T2 T3>

< , 0 0>

MSC/NASTRAN 102 Exercise Workbook

14a-11

Rotations <R1 R2 R3>

<0 , , >

OK
Select Application Region...
Select Geometry Entities

Point 2

Add
OK
Apply
New Set Name

permanent_constraint

Input Data...
Translations <T1 T2 T3>

< , , 0>

Rotations <R1 R2 R3>

<0 , 0 , >

OK
Select Application Region...
Select Geometry Entities

Curve 1

Add
OK
Apply
7. Create a set of material properties for the bar.
Materials
Action:

Create

Object:

Isotropic

Method:

Manual Input

Material Name

alum

Input Properties...

14a-12

Elastic Modulus =

10.0E6

Poisson Ratio =

.3

Density =

.101

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 14a

Modal Analysis of a Beam

Apply
Cancel
8. Define the bar properties.
Properties
Action:

Create

Dimension:

1D

Type:

Beam
bar

Property Set Name


Input Properties...

m:alum

Material Name
(Select from Material Property Sets box)

Use Beam Section

<Click on Beam Library>

New Section Name

ibeam

W1

W2

0.1

t1

0.1

t2

0.1

OK
Bar Orientation

Coord 0.2

OK
Select Members

Curve 1

Add
Apply

MSC/NASTRAN 102 Exercise Workbook

14a-13

9.

Now, you will generate the input file for analysis.

Analysis
Action:

Analyze

Object:

Entire Model

Method

Analysis Deck
prob14a

Job Name
Solution Type...

NORMAL MODES

Solution Type:
Solution Parameters ...
<deselect Automatic
Constraints>

Automatic

Mass Calculation:

Coupled

Data Deck Echo:

None

Wt. -Mass Conversion =

.00259

Constraints

OK
OK
Subcase Create...
Available Subcases

Default

Subcase Parameters...
Number of Desired Roots =

OK
Apply
Cancel

<to close form>

Apply
An MSC/NASTRAN input file called prob14a.bdf will be generated. The process
of translating your model into an input file is called Forward Translation. The
Forward Translation is complete when the Heartbeat turns green. MSC/PATRAN
Users should proceed to step 11.

14a-14

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 14a

Modal Analysis of a Beam

Generating an input file for MSC/NASTRAN Users:


MSC/NASTRAN users can generate an input file using the data from
Table 14a.1. The result should be similar to the output below.
10. MSC/NASTRAN Input File: prob14a.dat
SOL 103
TIME 600
CEND
TITLE = Normal Modes Example
SUBCASE 1
METHOD = 1
SPC = 1
VECTOR=ALL
BEGIN BULK
PARAM,WTMASS,.00259
PARAM,COUPMASS,1
EIGRL
1
PBARL
1
1
I
+
A 2.
1.
1.
CBAR
1
1
1
CBAR
2
1
2
CBAR
3
1
3
CBAR
4
1
4
CBAR
5
1
5
CBAR
6
1
6
CBAR
7
1
7
CBAR
8
1
8
CBAR
9
1
9
CBAR
10
1
10
MAT1
1
1.+7
GRID
1
0.
GRID
2
10.
GRID
3
20.
GRID
4
30.
GRID
5
39.9999
GRID
6
49.9999
GRID
7
60.
GRID
8
70.
GRID
9
80.
GRID
10
90.
GRID
11
100.
SPC1
1
1234
1
SPC1
1
234
11
ENDDATA

.1
2
3
4
5
6
7
8
9
10
11
.3
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.

.1
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
.101
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.

+
.1
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.

0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
345
345
345
345
345
345
345
345
345
345
345

MSC/NASTRAN 102 Exercise Workbook

14a-15

Submit the input file for analysis


11.

Submit the input file to MSC/NASTRAN for analysis.


11a. To submit the MSC/PATRAN .bdf file for analysis, find an
available UNIX shell window. At the command prompt
enter: nastran prob14a.bdf scr=yes. Monitor the run using
the UNIX ps command.
11b. To submit the MSC/NASTRAN .dat file for analysis, find an
available UNIX shell window. At the command prompt
enter: nastran prob14a scr=yes. Monitor the run using the
UNIX ps command.

12. When the run is completed, edit the prob14a.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.
13. While still editing prob14a.f06, search for the word:
E I G E N (spaces are necessary)
What are the first three modes?
1st = __________Hz
2nd = __________Hz
3rd = __________Hz

Comparison of Results
14.

14a-16

Compare the results obtained in the .f06 file with the results on
the following page:

MSC/NASTRAN 102 Exercise Workbook

1
2
3

EXTRACTION
ORDER
1
2
3

EIGENVALUE
2.239398E+04
3.549898E+05
1.771818E+06

R E A L
E I G E N V A L U E S
RADIANS
CYCLES
1.496462E+02
5.958102E+02
1.331096E+03

2.381693E+01
9.482614E+01
2.118506E+02

GENERALIZED GENERALIZED
MASS STIFFNESS
1.000000E+00 2.239398E+04
1.000000E+00 3.549898E+05
1.000000E+00 1.771818E+06

WORKSHOP 14a

MODE
NO.

Modal Analysis of a Beam

MSC/NASTRAN 102 Exercise Workbook


14a-17

MSC/NASTRAN Users have finished this exercise. MSC/PATRAN Users


should proceed to the next step.
15. Proceed with the Reverse Translation process, that is
importing the prob14a.op2 results file into MSC/PATRAN.
To do this, return to the Analysis form and proceed as follows:
Analysis
Action:

Read Output2

Object:

Result Entities

Method

Translate

Select Results File...


Select Results File

prob14a.op2

OK
Apply
When the translation is complete bring up the Results form.
Results
Form Type:

Basic

Select Results Cases

1.1-Default, Mode 1:Freq=23.816

Select Deformation Result

1.1 Eigenvectors, Translational

Apply
To reset the graphics, click on this icon:
Reset Graphics
You can go back and select any Results Case, Fringe Results or
Deformation Results you are interested in.
Quit MSC/PATRAN when you are finished with this exercise.

14a-18

MSC/NASTRAN 102 Exercise Workbook

14b-2

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 14b

Normal Modes with Differential Stiffness

Model Description:
The goal of this example is to analyze a stiffened model. In this case, the
beam from Problem 14a. with a 500 lb force applied.
Figure 14b.1 below is a finite element representation of the beam. This is
no longer a simple normal modes analysis. Instead we will be using a
nonlinear static solution (SOL 106) with (PARAM, NMLOOP and
METHOD and EIGRL).
Below is a finite element representation of the beam. One end is pinned
in 3 translations and one rotation. The other is pinned in 2 translations and
one rotation with a 500 lb force applied.
Figure 14b.1-Grid Coordinates and Element Connectivities

MSC/NASTRAN 102 Exercise Workbook

14b-3

Figure 14b.2-Beam Cross Section

1.0 in
0.1 in

2.0 in

0.1 in

0.1 in
1.0 in

Table 14b.1

14b-4

Length

100 in

Height

2 in

Width

1 in

Thickness

0.100 in

Area

0.38 in2

I1

0.229 in4

I2

0.017 in4

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 14b

Normal Modes with Differential Stiffness

Theoretical Solution

K n EIg
1 Pl 2 1 2
f n = ------ ---------- 1 + ------ --------
2 W l 4
Kr EI
For Mode 1, Kr =9.87
6
12
1
( 500 ) ( 100 ) 2
9.87 10 10 ( 0.229 ) ( 386.4 )
f n = ---------- ------------------------------------------------------- x 1 + ---------- -----------------------------------------
9.87 ( 10 106 ) ( 0.229 )
2 ( 0.38 ) ( 0.101 ) ( 100 ) 4

f n = 26.36Hz

For Static Load

PL
= -------AE
500 ( 100 )
= --------------------------------6
0.38 ( 10 10 )
= 0.0132

MSC/NASTRAN 102 Exercise Workbook

14b-5

Suggested Exercise Steps


Open database created in Problem 1a in order to modify it, adding a load
and reanalyze.
Create 500 lb force applied at one end (FORCE).
Make sure analysis is set to nonlinear static (SOL 106).
Prepare nonlinear analysis to also analyze for normal mode (PARAM
NMLOOP, EIGRL, LGDISP, NLPARM).

14b-6

Review the results, specifically the eigenvectors.

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 14b

Normal Modes with Differential Stiffness

ID SEMINAR,PROB1
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________

CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK
MSC/NASTRAN 102 Exercise Workbook

14b-7

14b-8

MSC/NASTRAN 102 Exercise Workbook

10

WORKSHOP 14b

Normal Modes with Differential Stiffness


4

ENDDATA

MSC/NASTRAN 102 Exercise Workbook

14b-9

10

Exercise Procedure:
1. Users who are not utilizing MSC/PATRAN for generating an input file
should go to Step 6, otherwise, proceed to step 2.
2. Open database created in Problem 14a named prob14a.db.
File/Open Database
prob14a

Existing Database Name


OK

3. Activate the entity labels by selecting the Show Labels icon on the toolbar.
Show Labels
4. Create force.
Loads/BCs
Action:

Create

Object:

Force

Type:

Nodal

New Set Name

pull

Input Data...
Force <F1 F2 F3>

<500, , >

OK
Select Application Region...
Select Geometry Entities
Add
OK
Apply

14b-10

MSC/NASTRAN 102 Exercise Workbook

Point 2

Normal Modes with Differential Stiffness

WORKSHOP 14b
5.

Now, you will generate the input file for analysis.

Analysis
Action:

Analyze

Object:

Entire Model

Method

Analysis Deck
prob14b

Job Name
Solution Type...

NONLINEAR STATIC

Solution Type:
Solution Parameters ...
<deselect Automatic
Constraints>

Automatic

Mass Calculation:

Coupled

Data Deck Echo:

None

Wt. -Mass Conversion =

.00259

Contraints

OK
OK
Direct Text Input...

Case Control Section

METHOD = 10

Bulk Data Section

PARAM, NMLOOP, 5
EIGRL, 10, , , 3

OK
Subcase Create...
Available Subcases

Default

Subcase Parameters...
Number of Load Increments =

OK
Apply
Cancel
Apply
MSC/NASTRAN 102 Exercise Workbook

14b-11

An MSC/NASTRAN input file called prob14b.bdf will be generated. The


process of translating your model into an input file is called Forward
Translation. The Forward Translation is complete when the Heartbeat
turns green. MSC/PATRAN Users should proceed to step 7.

14b-12

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 14b

Normal Modes with Differential Stiffness

Generating an input file for MSC/NASTRAN Users:


MSC/NASTRAN users can generate an input file using the data from Table 14b.1.
The result should be similar to the output below.
6.

MSC/NASTRAN Input File: prob14b.dat

SOL 106
TIME 600
CEND
$
TITLE = Normal Modes with Differential Stiffness
METHOD = 10
SUBCASE 1
NLPARM = 1
SPC = 1
LOAD = 1
DISPLACEMENT=ALL
$
BEGIN BULK
PARAM
COUPMASS 1
PARAM
WTMASS .00259
PARAM
LGDISP 1
NLPARM
1
5
AUTO
5
25
+
A
.001
1.-7
PARAM,NMLOOP,5
$
EIGRL,10,,,3
PBARL
1
1
I
+
B 2.
1.
1.
.1
.1
.1
CBAR
1
1
1
2
0.
1.
CBAR
2
1
2
3
0.
1.
CBAR
3
1
3
4
0.
1.
CBAR
4
1
4
5
0.
1.
CBAR
5
1
5
6
0.
1.
CBAR
6
1
6
7
0.
1.
CBAR
7
1
7
8
0.
1.
CBAR
8
1
8
9
0.
1.
CBAR
9
1
9
10
0.
1.
CBAR
10
1
10
11
0.
1.
$
MAT1
1
1.+7
.3
.101
GRID
1
0.
0.
0.
345
GRID
2
10.
0.
0.
345
GRID
3
20.
0.
0.
GRID
4
30.
0.
0.
GRID
5
39.9999 0.
0.
GRID
6
49.9999 0.
0.
GRID
7
60.
0.
0.
345
GRID
8
70.
0.
0.
345
GRID
9
80.
0.
0.
345
GRID
10
90.
0.
0.

PW

NO

0.
0.
0.
0.
0.
0.
0.
0.
0.
0.

345
345
345
345

345

MSC/NASTRAN 102 Exercise Workbook

14b-13

GRID
LOAD
SPC1
SPC1
FORCE
ENDDATA

14b-14

11
2
1
1
1

1.
1234
234
11

100.
1.
1
11
0

0.
1

0.

500.

1.

MSC/NASTRAN 102 Exercise Workbook

345

0.

0.

WORKSHOP 14b

Normal Modes with Differential Stiffness

Submit the input file for analysis


7.

Submit the input file to MSC/NASTRAN for analysis.


7a.

To submit the MSC/PATRAN .bdf file for analysis, find an


available UNIX shell window. At the command prompt
enter: nastran prob14b.bdf scr=yes. Monitor the run
using the UNIX ps command.

7b.

To submit the MSC/NASTRAN .dat file for analysis, find


an available UNIX shell window. At the command prompt
enter: nastran prob14b scr=yes. Monitor the run using
the UNIX ps command.

8.

When the run is completed, edit the prob14b.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.

9.

While still editing prob14b.f06, search for the word:

E I G E N (spaces are necessary)


What are the first three natural frequencies?
1st = __________Hz
2nd = __________Hz
3rd = __________Hz

Comparison of Results
10.

Compare the results obtained in the .f06 file with the results
on the following page:

MSC/NASTRAN 102 Exercise Workbook

14b-15

14b-16
MODE
NO.

EXTRACTION
ORDER

MSC/NASTRAN 102 Exercise Workbook

1
2
3

1
2
3

EIGENVALUE

2.735837E+04
3.748482E+05
1.816509E+06

R E A L
E I G E N V A L U E S
RADIANS
CYCLES

1.654037E+02
6.122484E+02
1.347779E+03

2.632481E+01
9.744236E+01
2.145057E+02

GENERALIZED
MASS
1.000000E+00
1.000000E+00
1.000000E+00

GENERALIZED
STIFFNESS
2.735837E+04
3.748482E+05
1.816509E+06

WORKSHOP 14b

Normal Modes with Differential Stiffness

11. MSC/NASTRAN Usershave finished this exercise. MSC/


PATRAN Users should proceed to the next step.
12. Proceed with the Reverse Translation process, that is importing
the prob14b.op2 results file into MSC/PATRAN. To do this,
return to the Analysis form and proceed as follows:
Analysis
Action:

Read Output2

Object:

Result Entities

Method

Translate

Select Results File...


Select Results File

prob14b.op2

OK
Apply
When the translation is complete bring up the Results form.
Results
Form Type:

Basic

Select Results Cases

1.1-Default, Mode 1:Freq=23.36

Select Deformation Result

1.1 Eigenvectors, Translational

Apply
To reset the graphics, click on this icon:
Reset Graphics
You can go back and select any Results Case, Fringe Results or
Deformation Results you are interested in.
Quit MSC/PATRAN when you are finished with this exercise.

MSC/NASTRAN 102 Exercise Workbook

14b-17

14b-18

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP PROBLEM 15

Weight Minimization of a
Three Bar Truss

Objectives:
Minimize the weight of the truss.
First mode must be between 1500-1550 Hz.
Submit the file for analysis in MSC/NASTRAN.
Recover the desired objective while satisfying the
frequency requirement.

MSC/NASTRAN 102 Exercise Workbook

15-1

15-2

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 15

Weight Minimization of a Three Bar Truss

Model Description:
You must minimize the weight of the following three bar truss
problem. The first mode must be between 1500-1550 hz. The model
will have different areas for the inside and outside beams. The
structure must remain symetric.
Below is a Geometric representation of the truss. It also contains the
loads and boundary constraints.
Figure 15.1-Loads and Boundary Conditions
123456

123456

Node 1

123456

Node 3

Node 2

Area 1

Area 2

Area 1

10

3456

Node 4
10

Table 15.1
Elastic Modulus

10E6

Poissons Ratio

.33

Density

.1

Wt.-Mass Conversion

.00259

Area 1

1.0

Area 2

2.0
MSC/NASTRAN 102 Exercise Workbook

15-3

Optimization Statement
Design Variables
These are the areas of the three rod elements (A1, A2, A3).
Objective
Minimize the weight of the truss.
Subject to the follow constraints
The first mode must be between 1500-1550 Hz.
A1=A3 to impose symmetry.

15-4

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 15

Weight Minimization of a Three Bar Truss

Suggested Exercise Steps:


Generate the analysis model. The nodes (GRID) and
element connectivities (CROD) should be defined manually.
Define material (MAT1) and element (PROD) properties.
Apply fixed boundary constraints (SPC1) to the upper
nodes.
Create the appropriate design optimization model.
Define the disign variables (DESVAR).
Relate one design variable to another design variable
(DLINK).
Define design variable to analysis model parameter relations
(DVPREL).
Specify design sensitivity response quantities (DRESP1).
Define constraints (DCONSTR).
Define optimization control parameters (DOPTPRM).
Prepare the model for linear static analysis and normal
modes analysis using Lanczos Method.
PARAM, WTMASS, 0.00259
Generate an input file and submit it to MSC/NASTRAN for
structural optimization analysis.
Review the results, specifically the eigenvalues and the
design variable history.

MSC/NASTRAN 102 Exercise Workbook

15-5

ID SEMINAR,PROB5
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________

CEND
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
BEGIN BULK
15-6

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 15

Weight Minimization of a Three Bar Truss

MSC/NASTRAN 102 Exercise Workbook

10

15-7

ENDDATA

15-8

MSC/NASTRAN 102 Exercise Workbook

10

WORKSHOP 15

Weight Minimization of a Three Bar Truss

Generating an input file for MSC/NASTRAN Users:


MSC/NASTRAN users can generate an input file using the data from page 153 (Model Description). The result should be similar to the output below.
1. MSC/NASTRAN input file: prob15.dat
ASSIGN OUTPUT2=prob15.op2, UNIT=12
ID NAS102, WORKSHOP 15
TIME 10
SOL 200
$ OPTIMIZATION
CEND
TITLE= SYMMETRIC THREE BAR TRUSS DESIGN OPTIMIZATION - VARIATION OF D200X1
SUBTITLE= GOAL IS TO MIN WT WHILE KEEPING THE 1ST MODE BETWEEN 1500-1550 HZ
ECHO= SORT
SPC= 100
DISP(PLOT) ALL
DESOBJ(MIN)= 100
$ (DESIGN OBJECTIVE = DRESP ID)
DESSUB= 200
$ DEFINE CONSTRAINT SET FOR BOTH SUBCASES
SUBCASE 1
ANALYSIS= MODES
METHOD= 10
BEGIN BULK
$
$---------------------------------------------------------------------$ ANALYSIS MODEL
$---------------------------------------------------------------------$
EIGRL, 10, , , 2
PARAM, POST, -1
PARAM, PATVER, 3.0
$
$ GRID DATA
$
2
3
4
5
6
7
8
9
10
GRID,
1,
,
-10.0,
0.0, 0.0
GRID,
2,
,
0.0,
0.0, 0.0
GRID,
3,
,
10.0,
0.0, 0.0
GRID,
4,
,
0.0, -10.0, 0.0
$ SUPPORT DATA
SPC,
100,
1,
123456, ,
2,
123456
SPC,
100,
3,
123456, ,
4,
3456
$ ELEMENT DATA
CROD,
1,
11,
1,
4
CROD,
2,
12,
2,
4
MSC/NASTRAN 102 Exercise Workbook

15-9

CROD,
3,
13,
3,
4
$ PROPERTY DATA
PROD,
11,
1,
1.0
PROD,
12,
1,
2.0
PROD,
13,
1,
1.0
MAT1,
1,
1.0E+7, ,
0.33,
0.1
$
PARAM, WTMASS, .00259
$
$---------------------------------------------------------------------$ DESIGN MODEL
$---------------------------------------------------------------------$
$...DESIGN VARIABLE DEFINITION
$
$DESVAR,ID,
LABEL, XINIT, XLB,
XUB,
DELXV(OPTIONAL)
DESVAR, 1,
A1,
1.0,
0.1,
100.0
DESVAR, 2,
A2,
2.0,
0.1,
100.0
DESVAR, 3,
A3,
1.0,
0.1,
100.0
$
$...IMPOSE X3=X1 (LEADS TO A3=A1)
$
$DLINK, ID,
DDVID, CO,
CMULT, IDV1,
C1,
IDV2,
C2,
$+,
IDV3,
C3,
...
DLINK, 1,
3,
0.0,
1.0,
1
1.00
$
$...DEFINITION OF DESIGN VARIABLE TO ANALYSIS MODEL PARAMETER RELATIONS
$
$DVPREL1,ID,
TYPE,
PID,
FID,
PMIN,
PMAX,
CO,
,
$+,
DVID1, COEF1, DVID2, COEF2, ...
DVPREL1, 10,
PROD,
11,
4,
,
,
,
,
+DP1,
1,
1.0
DVPREL1, 20,
PROD,
12,
4,
,
,
,
,
+DP2,
2,
1.0
DVPREL1, 30,
PROD,
13,
4,
,
,
,
,
+DP3,
3,
1.0
$
$...STRUCTURAL RESPONSE INDENTIFICATION
$
$DRESP1 ID
LABEL
RTYPE
PTYPE
REGION ATTA
ATTB
ATT1
$+
ATT2
...
DRESP1 100
W
WEIGHT
DRESP1 210
MODE1
EIGN
1
$
$...CONSTRAINTS
15-10

MSC/NASTRAN 102 Exercise Workbook

+
+DP1
+DP2
+DP3

WORKSHOP 15

Weight Minimization of a Three Bar Truss

$
$DCONSTR,DCID, RID,
LALLOW, UALLOW
DCONSTR, 200,
210,
8.883E7, 9.485E7
$
$...OPTIMIZATION CONTROL
$
DOPTPRM, DESMAX, 30
$
$.......2.......3.......4.......5.......6.......7.......8.......9.......0
ENDDATA

MSC/NASTRAN 102 Exercise Workbook

15-11

Submitting the input file for analysis:


2. Submit the input file to MSC/NASTRAN for analysis.
To submit the MSC/NASTRAN .dat file, find an available UNIX shell
window and at the command prompt enter nastran prob15 scr=yes.
Monitor the run using the UNIX ps command.
3. When the run is completed, edit the prob15.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.
4. While still editing prob15.f06, search for the word:
DESIGN VARIABLE HISTORY

Design Variable

15-12

Initial Value

Optimization Value

Iteration Value

_________

______________

___________

_________

______________

___________

_________

______________

___________

MSC/NASTRAN 102 Exercise Workbook

WORKSHOP 15

Weight Minimization of a Three Bar Truss

Comparison of Results
5. Compare the results obtained in the .f06 file with the results on the
following page:

MSC/NASTRAN 102 Exercise Workbook

15-13

15-14
***************************************************************
SUMMARY OF DESIGN CYCLE HISTORY
***************************************************************
(HARD CONVERGENCE ACHIEVED)

MSC/NASTRAN 102 Exercise Workbook

(SOFT CONVERGENCE ACHIEVED)


NUMBER OF FINITE ELEMENT ANALYSES COMPLETED
13
NUMBER OF OPTIMIZATIONS W.R.T. APPROXIMATE MODELS 12
OBJECTIVE AND MAXIMUM CONSTRAINT HISTORY
------------------------------------------------------------------------------------------------------------------------------------------------OBJECTIVE FROM
OBJECTIVE FROM
FRACTIONAL ERROR
MAXIMUM VALUE
CYCLE
APPROXIMATE
EXACT
OF
OF
NUMBER
OPTIMIZATION
ANALYSIS
APPROXIMATION
CONSTRAINT
-----------------------------------------------------------------------------------------------------------------------------------------------INITIAL
4.828427E+00
1.922634E-01
1
4.623621E+00
4.623881E+00
-5.620305E-05
-3.995783E-03
2
3.699294E+00
3.699106E+00
5.104681E-05
-3.995698E-03
3
2.958779E+00
2.959009E+00
-7.783429E-05
-4.092440E-03
4
2.367016E+00
2.367207E+00
-8.087595E-05
-4.092440E-03
5
1.892135E+00
1.892194E+00
-3.087028E-05
-4.955952E-03
6
1.514454E+00
1.514385E+00
4.518410E-05
-4.523099E-03
7
1.210997E+00
1.210911E+00
7.078262E-05
-5.035487E-03
8
9.690463E-01
9.690970E-01
-5.227955E-05
-4.640506E-03
9
7.749816E-01
7.748774E-01
1.344586E-04
-4.126600E-03
10
6.199054E-01
6.199020E-01
5.480649E-06
-4.126516E-03
11
5.758885E-01
5.758446E-01
7.618205E-05
-2.931787E-04
12
5.758446E-01
5.758446E-01
0.000000E+00
-2.931787E-04
------------------------------------------------------------------------------------------------------------------------------------------------

MSC/NASTRAN 102 Exercise Workbook

----------------------------------------------------------------------------------------------------------------------------INTERNAL |
EXTERNAL
|
|
DV. ID. |
DV. ID.
|
LABEL |
6
:
7
:
8
:
9
:
10
:
11
----------------------------------------------------------------------------------------------------------------------------1 |
1
| A1
|
2.6187E-01 :
2.0929E-01 :
1.6756E-01 :
1.3405E-01 :
1.0724E-01 :
1.0000E-01
2 |
2
| A1
|
7.7370E-01 :
6.1896E-01 :
4.9517E-01 :
3.9573E-01 :
3.1659E-01 :
2.9300E-01
3 |
3
| A1
|
2.6187E-01 :
2.0929E-01 :
1.6756E-01 :
1.3405E-01 :
1.0724E-01 :
1.0000E-01

----------------------------------------------------------------------------------------------------------------------------INTERNAL |
EXTERNAL
|
|
DV. ID. |
DV. ID.
|
LABEL |
12
:
13
:
14
:
15
:
16
:
17
----------------------------------------------------------------------------------------------------------------------------1 |
1
| A1
|
1.0000E-01 :
2 |
2
| A1
|
2.9300E-01 :
3 |
3
| A1
|
1.0000E-01 :
-----------------------------------------------------------------------------------------------------------------------------

15-15

Weight Minimization of a Three Bar Truss

DV. ID. |
DV. ID.
|
LABEL |
INITIAL
:
1
:
2
:
3
:
4
:
5
:
----------------------------------------------------------------------------------------------------------------------------1 |
1
| A1
|
1.0000E+00 :
8.0000E-01 :
6.4000E-01 :
5.1190E-01 :
4.0952E-01 :
3.2706E-01
2 |
2
| A1
|
2.0000E+00 :
2.3611E+00 :
1.8889E+00 :
1.5111E+00 :
1.2089E+00 :
9.6712E-01
3 |
3
| A1
|
1.0000E+00 :
8.0000E-01 :
6.4000E-01 :
5.1190E-01 :
4.0952E-01 :
3.2706E-01

WORKSHOP 15

DESIGN VARIABLE HISTORY


----------------------------------------------------------------------------------------------------------------------------INTERNAL |
EXTERNAL
|
|

15-16

MSC/NASTRAN 102 Exercise Workbook

You might also like