Total NAS102
Total NAS102
WORKSHOP 1
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.
1-3
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
WORKSHOP 1
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
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
WORKSHOP 1
1-7
ID SEMINAR,PROB1
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK
1-8
WORKSHOP 1
10
1-9
ENDDATA
1-10
10
WORKSHOP 1
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
<5, 2, 0>
[ 0, 0, 0]
Apply
1-11
Surface 1.2
Surface 1.3
Y
Create
Object:
Mesh Seed
Type:
Uniform
Number of Elements
Number =
10
Curve List
Surface 1.2
Apply
5a. Change the number of mesh seeds to 4 and select the right edge.
Number =
Curve List
Surface 1.3
1-12
WORKSHOP 1
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.
Create
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
plate
Input Properties...
Material Name
m:mat_1
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
Create
WORKSHOP 1
Object:
Displacement
Type:
Nodal
fixed
<0, 0, 0>
<0, 0, >
Coord 0
OK
Select Application Region...
Select the curve or edge icon.
Curve
Surface 1.1
Add
OK
Apply
1-15
Figure 1.5
Surface 1.1
9.
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.
Analysis
1-16
Action:
Analyze
Object:
Entire Model
Method
Analysis Deck
WORKSHOP 1
Job Name
Solution Type...
NORMAL MODES
Solution Type:
Solution Parameters...
Mass Calculation:
Coupled
Unsorted
.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.
1-17
5
1
2
13
1
12
13
24
23
24
35
34
35
46
45
.3
0.
.282
0.
.5
0.
1.
0.
WORKSHOP 1
GRID
34
=,*1,=,*0.5,==
=9
GRID
45
=,*1,=,*0.5,==
=9
SPC1
1
ENDDATA
12345
1.5
0.
0.
2.
0.
0.
12
23
34
45
1-19
11a.
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
45
WORKSHOP 1
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
1-21
Comparison of Results
15.
1-22
Compare the results obtained in the .f06 file with the results on
the following page:
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
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
OK
prob1.op2
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
1.1-Default, Mode 1
Apply
The results should resemble Figure 1.6.
1-24
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
1-25
1-26
WORKSHOP PROBLEM 2
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.
2-1
2-2
WORKSHOP 2
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
2-3
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
WORKSHOP 2
2-5
ID SEMINAR,PROB2
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK
2-6
WORKSHOP 2
10
2-7
ENDDATA
2-8
10
WORKSHOP 2
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:
Object:
Model Data
Method
Translate
plate.bdf
OK
Apply
OK
4. Activate the entity labels by selecting the Show Labels icon on the toolbar.
Show Labels
2-9
Modify
Default
Assign/Prioritize Loads/BCs
Select Load/BCs to Add to
Spreadsheet
Displ_spc1.1
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
Unsorted
.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
WORKSHOP 2
OK
Subcase Create...
Available Subcases
Default
Subcase Parameters...
Extraction Method:
Automatic Givens
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.
2-11
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
WORKSHOP 2
9.
8b.
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.
Comparison of Results
11.
Compare the results obtained in the .f06 file with the results on
the following page:
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
Read Output2
Object:
Result Entities
Method
Translate
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
2-15
Basic
1.1-Default, Mode1
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
WORKSHOP PROBLEM 3
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.
3-1
3-2
WORKSHOP 3
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
3-3
3-4
WORKSHOP 3
ID SEMINAR,PROB3
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK
3-5
3-6
10
WORKSHOP 3
10
ENDDATA
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
Default
Analysis Code:
MSC/NASTRAN
OK
3. Create the model by importing an existing MSC/NASTRAN input file,
(plate.bdf).
Analysis
Action:
Object:
Model Data
Method
Translate
Select File
OK
Apply
OK
4. Activate the entity labels by selecting the Show Labels icon on the toolbar.
Show Labels
3-8
WORKSHOP 3
Create
transient_response
Time Dependent
Assign/Prioritize Loads/BCs
Select Load/BCs to Add to
Spreadsheet
Displ_spc1.1
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
WORKSHOP 3
0.04
Value
0.0
OK
Apply
8. Create the time dependent pressure.
Loads/BCs
Action:
Create
Object:
Pressure
Type:
Element Uniform
pressure
2D
Input Data...
-1
f:time_dependent_pressure
OK
Select Application Region...
FEM
Select 2D Elements or Edge
Elem 1:40
Add
OK
Apply
3-11
Create
Object:
Force
Type:
Nodal
force
<0 0 50>
Time Dependence:
f:time_dependent_force
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
WORKSHOP 3
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
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
0.06
1571
OK
OK
3-13
transient_response
Subcase Parameters...
Ending Time =
.04
100
OK
Output Requests...
Form Type:
Advanced
select DISPLACEMENT(...
Sorting:
By Freq/Time
Modify
OK
Apply
Cancel
Subcase Select...
Subcases Selected:
Default
(Click to de-select.)
OK
3-14
transient_response
WORKSHOP 3
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.
3-15
WORKSHOP 3
3-17
T3
.0024 = __________
.0052 = __________
.02
= __________
Displacement at Grid 33
Time
T3
.0024 = __________
.0052 = __________
.02
3-18
= __________
WORKSHOP 3
Displacement at Grid 55
Time
T3
.0024 = __________
.0052 = __________
.02
= __________
3-19
Comparison of Results
16.
3-20
Compare the results obtained in the .f06 file with the results
on the following page:
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
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
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
Read Output2
Object:
Result Entities
Method
Translate
prob3.op2
OK
Apply
When the translation is complete bring up the Results form.
Results
Form Type:
Advanced
Get Results
Select Result
Plot Type
XY Plot
1-Time
Apply
Result(Y)...
Results
Vector Component
1.1-Displacements, Translational
X
OK
Node IDs
3-22
Node 11
WORKSHOP 3
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
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
Rename
Apply
Figure 3.4-Displacement Response at Node 33
3-24
WORKSHOP 3
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
3-25
3-26
WORKSHOP PROBLEM 4
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.
4-1
4-2
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
25.0
Z
Y
X
4-3
4-4
WORKSHOP 4
ID SEMINAR,PROB4
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK
4-5
4-6
10
WORKSHOP 4
10
ENDDATA
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
Default
Analysis code:
MSC/NASTRAN
OK
3. Create the model by importing an existing MSC/NASTRAN input file,
(plate.bdf).
Analysis
Action:
Object:
Model Data
Method:
Translate
4. Activate the entity labels by selecting the Show Labels icon on the toolbar.
Show Labels
4-8
WORKSHOP 4
Create
transient_response
Time Dependent
Assign/Prioritize Loads/BCs
Select Load/BCs to Add to
Spreadsheet
Displ_spc1.1
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
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
32
WORKSHOP 4
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
pressure
2D
Input Data...
Top Surf Pressure
-1
Time Dependence
f:time_dependent_pressure
OK
Select Application Region ...
FEM
MSC/NASTRAN 102 Exercise Workbook
4-11
Elm 1:40
Add
OK
Apply
9. Create the time-dependent nodal force.
Loads/BCs
Action:
Create
Object:
Force
Type:
Nodal
force
<0 0 25>
Time Dependence
f:time_dependent_force
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
WORKSHOP 4
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
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
4-13
OK
OK
SDAMPING = 100
OK
Subcase Create...
Available Subcases
transient_response
Subcase Parameters...
Ending Time =
.04
100
OK
Output Requests...
Form Type:
Advanced
select DISPLACEMENT(...
Sorting:
By Freq/Time
Modify
OK
Apply
Cancel
4-14
WORKSHOP 4
default
(Click to deselect.)
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.
4-15
WORKSHOP 4
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
4-17
T3
.0064= __________
.0092= __________
.02 = __________
Displacement at Grid 33.
Time
T3
.0068= __________
.0092= __________
.02 = __________
4-18
WORKSHOP 4
T3
.0068= __________
.0092= __________
.02 = __________
4-19
Comparison of Results
16. Compare the results obtained in the .f06 file with the results on the following page:
4-20
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
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
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
Analysis
Action:
Read Output2
Object:
Result Entities
Method:
Translate
prob4.op2
OK
Apply
When the translation is complete bring up the Results form.
Results
Form Type:
Advanced
Get Results
Select Result
Plot Type
XY Plot
1. Time
Apply
Result (Y)
4-22
Results:
Vector Component
WORKSHOP 4
OK
Node 11
Node IDs
Apply
Displacement Response At
Loaded Corner
.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
Get Results
Select Result
4-23
Plot Type
XY Plot
1. Time
Apply
Result(Y)...
Results:
Result(Y)...
OK
Node 33
Node IDs
Apply
.180
.120
.0600
0.
-.0600
-.120
-.180
0.
4-24
.00700
.0140
.0210
.0280
.0350
.0420
WORKSHOP 4
20. Repeat the procedure to find the nodal displacement for Node 55.
Results
Form Type:
Advanced
Get Results
Select Result
Plot Type
XY Plot
1. Time
Apply
Result(Y)...
Results:
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.
4-25
.140
.0700
0.
-.0700
-.140
-.210
0.
.00700
.0140
.0210
.0280
.0350
.0420
4-26
5-2
WORKSHOP 5
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
5-3
5-4
WORKSHOP 5
ID SEMINAR,PROB5
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK
MSC/NASTRAN 102 Exercise Workbook
5-5
5-6
10
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
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:
Object:
Model Data
Method:
Translate
5-8
WORKSHOP 5
Load Cases
Action:
Create
frequency_response
Time Dependent
Assign/Prioritize Loads/BCs
Displ_spc1.1
OK
Apply
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
<0 0 1>
Time Dependence:
f:frequency_dependent_load
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
5-10
WORKSHOP 5
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
5-11
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
select DISPLACEMENT(...
Sorting:
By Freq/Time
Modify
OK
Apply
Cancel
Subcase Select ...
Subcases Selected:
Default
(Click to de-select.)
frequency_response
to select.)
OK
Apply
5-12
WORKSHOP 5
5-13
WORKSHOP 5
$
$ 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
5-15
5-16
WORKSHOP 5
Displacement (Y)
140
__________
380
__________
Displacement at Grid 33
Frequency (X)
Displacement (Y)
140
__________
600
__________
Displacement at Grid 55
Frequency (X)
Displacement (Y)
140
__________
1000
__________
5-17
Comparison of Results
14. Compare the results obtained in the .f06 file with the results on the
following page:
5-18
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
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
FREQUENCY
2.000000E+01
D I S P L A C E M E N T
(MAGNITUDE/PHASE)
WORKSHOP 5
POINT-ID =
Analysis
Action:
Read Output2
Object:
Result Entities
Method:
Translate
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
Get Results
Select Result
1.1-Displacements, Translational
Plot Type:
XY Plot
Global Var...
Global Variable:
Apply
Result (Y)...
5-20
1-Frequency
WORKSHOP 5
Results
1.1-Displacements, Translational
Vector Component
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
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
Phase
OK
Node IDs
Node 11
Apply...
Result XY Window Name:
XYWindow2
Rename
Apply
Figure 5.4-Phase Angle at Node 11
5-22
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
Displacement vs Frequency
at Node 33
Rename
Apply
5-23
Phase
OK
Node IDs
Node 33
Apply
Result XY Window Name:
XYWindow4
Phase vs Frequency
at Node 33
Rename
Apply
5-24
WORKSHOP 5
Mag.
OK
Node 55
Node IDs
Apply...
Result XY Window Name:
XYWindow5
Displacement vs Frequency
at Node 55
Rename
Apply
5-25
Phase
OK
Node IDs
Node 55
Apply...
Result XY Window Name:
XYWindow6
Phase vs Frequency
at Node 55
Rename
Apply
5-26
WORKSHOP 5
5-27
5-28
6-2
WORKSHOP 6
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
1.0
6-3
6-4
WORKSHOP 6
ID SEMINAR,PROB6
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK
MSC/NASTRAN 102 Exercise Workbook
6-5
6-6
10
WORKSHOP 6
10
ENDDATA
MSC/NASTRAN 102 Exercise Workbook
6-7
Exercise Procedure:
1.
2.
File/New Database
prob6
Default
Analysis Code:
MSC/NASTRAN
Analysis Type:
Structural
OK
3.
Analysis
Action:
Object:
Model Data
Method
Translate
Select File
OK
Apply
OK
4.
6-8
WORKSHOP 6
5.
Load Cases
Action:
Create
frequency_dependent
Time Dependent
Assign/Prioritize Loads/BCs
Displ_spc1.1
OK
Apply
6.
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
6-9
OK
Apply
7.
Loads/BCs
Action:
Create
Object:
Pressure
Type:
Element Uniform
pressure
2D
Input Data...
Top Surf Pressure
-0.1
Time Dependence
f:frequency_dependent_load
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
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
<0, 0 , 1>
Time Dependence:
f:frequency_dependent_load
OK
MSC/NASTRAN 102 Exercise Workbook
6-11
FEM
Select Nodes
Node 11
Add
OK
Apply
9.
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
WORKSHOP 6
SDAMPING = 100
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.)
frequency_dependent
(Click to select.)
OK
Apply
6-13
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
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
WORKSHOP 6
11.
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
WORKSHOP 6
$
DAREA, 600, 11, 3, 1.0
$
$ COMBINE LOADS
$
DLOAD, 100, 1., .1, 400, 1.0, 500
$
$ SPECIFY FREQUENCY STEPS
$
FREQ1, 100, 20., 20., 49
FREQ4, 100, 20., 1000., .03, 5
$
ENDDATA
6-17
6-18
13.
14.
WORKSHOP 6
Displacement at Grid 11
Frequency (X)
Displacement (Y)
140
__________
440
__________
Displacement at Grid 33
Frequency (X)
Displacement (Y)
140
__________
660
__________
Displacement at Grid 55
Frequency (X)
Displacement (Y)
140
__________
1000
__________
6-19
Comparison of Results
16.
6-20
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
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
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.
18.
Analysis
Action:
Read Output2
Object:
Result Entities
Method
Translate
prob6.op2
OK
Apply
19.
The first plot is the Displacement versus Frequency plot at Node 11.
Results
Form Type:
Advanced
Get Results
Select Result
1.1-Displacements, Translational
Plot Type:
XY Plot
Global Var...
Global Variable:
1-Frequency
Apply
Result (Y)...
Results:
6-22
1.1-Displacements, Translational
WORKSHOP 6
Vector Component
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
6-23
Results
Plot Type Options...
Result (Y)...
Numerical Form for Complex
Results
Mag.
OK
Node IDs
Node 33
Apply...
Result XY Window Name:
XYWindow2
Displacement Response
at Tip Center
Rename
6-24
WORKSHOP 6
XY Plot
Action:
Modify
Object:
Axis
Active Axis:
Scale...
Scale:
Logarithmic
Apply
Cancel
6-25
Results
Plot Type Options...
Result (Y)...
Numerical Form for Complex
Results
Mag.
OK
Node IDs
Node 55
Apply...
Result XY Window Name:
XYWindow3
Displacement Response
at Opposite Corner
Rename
Apply
XY Plot
Action:
Modify
Object:
Axis
Active Axis:
Scale...
Scale:
Apply
Cancel
6-26
Logarithmic
WORKSHOP 6
6-27
6-28
7-2
WORKSHOP 7
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
7-3
7-4
WORKSHOP 7
ID SEMINAR,PROB7
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK
7-5
7-6
10
WORKSHOP 7
10
ENDDATA
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
Default
Analysis code:
MSC/NASTRAN
OK
3. Create the model by importing an existing MSC/NASTRAN input file,
(plate.bdf).
Analysis
Action:
Object:
Model Data
Method:
Translate
4. Activate the entity labels by selecting the Show Labels icon on the toolbar.
Show Labels
Iso 3 View
7-8
WORKSHOP 7
Create
transient_response
Time Dependent
Assign/Prioritize loads/BCs
Displ_spc1.1
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
scalar_mass
Option(s)
Grounded
MSC/NASTRAN 102 Exercise Workbook
7-9
Input Properties...
Mass
1000
DOF at Node 1
UZ
OK
Select Members:
Elm 41
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
Auto Execute
Create Independent
Node List:
Node 23
Apply
Create Dependent
Node List:
7-10
WORKSHOP 7
DOFs
Apply
Cancel
Apply
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
Create
Object:
Force
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
spc1.1
Modify Data...
Translations < T1 T2 T3 >
< 0, 0, >
< 0, 0, >
OK
Apply
7-12
WORKSHOP 7
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
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
7-13
.06
1571
OK
OK
Subcase Create ...
Available Subcases
transient_response
Subcase Parameters...
Ending Time =
.04
200
OK
Output Request
under Select Result Type,
highlight:
Velocities
Accelerations
Default
(Click to deselect.)
transient_response
(Click to select.)
OK
Apply
7-14
WORKSHOP 7
7-15
WORKSHOP 7
7-17
7-18
WORKSHOP 7
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 = __________
7-19
T3
= __________
.02 = __________
.04 = __________
Acceleration at Grid 33
Time
.0
T3
= __________
.02 = __________
.04 = __________
7-20
WORKSHOP 7
Comparison of Results
16. Compare the results obtained in the .f06 file with the results on the following page:
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
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
Read Output2
Object:
Result Entities
Method:
Translate
Advanced
Get Results
Select Result
Plot Type
XY Plot
1. Time
Apply
Result (Y)...
Results
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
Get Results
7-24
Select Result
Plot Type
XY Plot
WORKSHOP 7
1. Time
Apply
Result (Y)...
Results
Vector Component
OK
Node 33
Node IDs:
Apply...
Result XY Window Name:
XYWindow2
Rename
Apply
The output should look similar to Figure 7.4
7-25
20. Repeat the procedure to find the nodal velocity for Node 23.
Results
Form Type:
Advanced
Get Results
Select Result
Plot Type
XY Plot
1. Time
Apply
Result (Y)...
7-26
Results
Vector Component
WORKSHOP 7
OK
Node 23
Node IDs:
Apply...
Result XY Window Name:
XYWindow3
Rename
Apply
The output should look similar to Figure 7.5
Figure 7.5-Base Velocity at Node 23
Advanced
Get Results
MSC/NASTRAN 102 Exercise Workbook
7-27
Select Result
Plot Type
XY Plot
1. Time
Apply
Result (Y)...
Results
Vector Component
OK
Node IDs:
Node 33
Apply...
Result XY Window Name:
XYWindow4
Rename
Apply
The output should look similar to Figure 7.6
7-28
WORKSHOP 7
22. Repeat the procedure to find the nodal acceleration for Node 23.
Results
Form Type:
Advanced
Get Results
Select Result
Plot Type
XY Plot
1. Time
Apply...
Result (Y)...
Results
Vector Component
7-29
OK
Node IDs:
Node 23
Apply...
Result XY Window Name:
XYWindow5
Rename
Apply
The output should be similar to Figure 7.7
Figure 7.7-Base Acceleration at Node 23
Advanced
Get Results
Select Result
7-30
WORKSHOP 7
Plot Type
XY plot
1. Time
Apply
Result (Y)...
Results:
Vector Component
OK
Node IDs:
Node 33
Apply...
Result XY Window Name:
XYWindow6
Rename
Apply
The output should look similar to Figure 7.8
7-31
7-32
8-2
WORKSHOP 8
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
8-3
8-4
WORKSHOP 8
ID SEMINAR,PROB8
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
BEGIN BULK
MSC/NASTRAN 102 Exercise Workbook
8-5
8-6
10
WORKSHOP 8
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
Default
Analysis code:
MSC/NASTRAN
OK
3. Create the model by importing an existing MSC/NASTRAN input file,
(plate.bdf).
Analysis
Action:
Object:
Model Data
Method:
Translate
4. Activate the entity labels by selecting the Show Labels icon on the toolbar.
Show Labels
8-8
WORKSHOP 8
Create
frequency response
time_dependent
Assign/Prioritize loads/BCs
Displ_spc1.1
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
scalar_mass
Option(s):
Grounded
1.0E+5
UZ
8-9
OK
Select Members:
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
Direct
Mass Calculation
Coupled
Wt.-Mass Conversion
.00259
0.06
OK
OK
Direct Text Input ...
Bulk Data Section
8-10
WORKSHOP 8
RLOAD2,500,600, , ,310
TABLED4,310,0.,1.,0.,10000.,
+,0.,0.,-39.4784,ENDT
DAREA,600,11,3,25.8799
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
frequency_response
OK
Apply
8-11
8-12
WORKSHOP 8
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
8-14
WORKSHOP 8
8-15
13.
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
__________
WORKSHOP 8
Comparison of Results
Compare the results obtained in the .f06 file with the results on the
following page:
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
Read Output2
Object:
Result Entities
Method:
Translate
prob8.op2
OK
Apply
16.
The first plot is to make the Displacement versus Frequency plot at Node
11.
Results
Form Type:
Advanced
Get Results
Select Result
1.1-Displacements, Translational
Plot Type:
XY Plot
Global Var...
Global Variable:
1-Frequency
Apply
Result (Y)...
MSC/NASTRAN 102 Exercise Workbook
8-19
Results:
1.1-Displacements, Translational
Vector Component
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
WORKSHOP 8
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
Phase
OK
Node 11
Node IDs
Apply...
Result XY Window Name:
XYWindow2
Rename
Apply
Figure 8.3-Phase Angle at Node 11
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
Displacement vs Frequency
at Node 33
Rename
Apply
8-22
WORKSHOP 8
Phase
OK
Node 33
Node IDs
Apply...
Result XY Window Name:
XYWindow4
Phase vs Frequency
at Node 33
Rename
Apply
8-23
Mag.
OK
Node IDs
Node 55
Apply...
Result XY Window Name:
XYWindow5
Displacement vs Frequency
at Node 55
Rename
Apply
8-24
WORKSHOP 8
Phase
OK
Node 55
Node IDs
Apply...
Result XY Window Name:
XYWindow6
Phase vs Frequency
at Node 55
Rename
Apply
8-25
8-26
9a-2
WORKSHOP 9a
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)
9a-3
9a-4
WORKSHOP 9a
ID SEMINAR,PROB9A
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
MSC/NASTRAN 102 Exercise Workbook
9a-5
BEGIN BULK
1
9a-6
10
WORKSHOP 9a
10
ENDDATA
MSC/NASTRAN 102 Exercise Workbook
9a-7
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
WORKSHOP 9a
9a-9
4.
Comparison of Results
5.
9a-10
Compare the plot made from the exercise with the plots on
the following pages.
WORKSHOP 9a
Figure 9a.2
9a-11
Figure 9a.3
9a-12
WORKSHOP 9a
Figure 9a.4
9a-13
Figure 9a.5
9a-14
9b-2
WORKSHOP 9b
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)
9b-3
9b-4
WORKSHOP 9b
ID SEMINAR,PROB9B
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
MSC/NASTRAN 102 Exercise Workbook
9b-5
BEGIN BULK
1
9b-6
10
WORKSHOP 9b
10
ENDDATA
MSC/NASTRAN 102 Exercise Workbook
9b-7
1.
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
WORKSHOP 9b
$
$ 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
9b-10
.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.
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
9b-11
9b-12
WORKSHOP 9b
Comparison of Results
4. Compare the results obtained in the .f06 file with the results on the
following page:
9b-13
9b-14
KLR
STRAIN
ENERGY
2.2118911E-09
5, 1997
MSC/NASTRAN
1/23/97
PAGE
12
SUBCASE 1
*** USER
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
-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
9b-16
SUBCASE 1
D I S P L A C E M E N T
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
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.
10-3
10-4
WORKSHOP 10
Random Analysis
ID SEMINAR,PROB10
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
MSC/NASTRAN 102 Exercise Workbook
10-5
BEGIN BULK
1
10-6
10
WORKSHOP 10
Random Analysis
10
ENDDATA
MSC/NASTRAN 102 Exercise Workbook
10-7
1.
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
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
WORKSHOP 10
Random Analysis
4.
10-11
Comparison of Results
5.
10-12
Compare the plot made from the exercise with the plots on
the following pages.
WORKSHOP 10
Random Analysis
Figure 10.1
10-13
Figure 10.2
10-14
WORKSHOP 10
Random Analysis
Figure 10.3
10-15
Figure 10.4
10-16
WORKSHOP 10
Random Analysis
Figure 10.5
10-17
Figure 10.6
10-18
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.
11-2
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
1.0
11-4
WORKSHOP 11
Random Analysis
ID SEMINAR,PROB11
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
11-6
WORKSHOP 11
Random Analysis
BEGIN BULK
1
10
ENDDATA
11-8
10
WORKSHOP 11
Random Analysis
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
11-10
WORKSHOP 11
Random Analysis
3.
11-12
WORKSHOP 11
Random Analysis
Figure 11.2
Figure 11.3
11-14
WORKSHOP 11
Random Analysis
Figure 11.4
Figure 11.5
11-16
WORKSHOP 11
Random Analysis
Figure 11.6
Figure 11.7
11-18
WORKSHOP 11
Random Analysis
Figure 11.8
Figure 11.9
11-20
WORKSHOP 11
Random Analysis
Figure 11.10
11-22
WORKSHOP PROBLEM 12
EXCITER
PILE
GROUND
Objectives
Define complex eigenvalue extraction parameters.
Submit the file for analysis in MSC/NASTRAN.
Compute complex modes.
12-1
12-2
WORKSHOP 12
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
12-3
12-4
WORKSHOP 12
ID SEMINAR,PROB12
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
_____________________________________________________
BEGIN BULK
MSC/NASTRAN 102 Exercise Workbook
12-5
12-6
10
WORKSHOP 12
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
[0 0 0]
Apply
Turn on the label and increase the node size by using the Quick Pick
buttons.
Show Label
12-8
Node Size
WORKSHOP 12
Location
Node 2
[ 1 0 0]
Node 3
[2 0 0]
Finite Element
Action:
Create
Object:
Element
Method:
Edit
Shape
Bar
Node 1 =
Node 1
Node 2 =
Node 2
Apply
6.
Node 1 =
Node 2
Node 2 =
Node 3
Apply
7.
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.
Finite Element
Action:
Create
Object:
Element
Method:
Edit
Shape
Bar
Node 1 =
Node 2
Node 2 =
Node 3
Apply
9.
Properties
Action:
Create
Dimension:
1D
Type:
Spring
spring1
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
Element 1
WORKSHOP 12
Beam Element
Add
Apply
10.
spring2
12500
DOF at Node 1:
UX
DOF at Node 2:
UX
OK
Element 2
Application Region
Add
Apply
11.
Properties
Action:
Create
Dimension:
1D
Type:
Damper
damper
Viscous
30
OK
Application Region
Element 5
MSC/NASTRAN 102 Exercise Workbook
12-11
Add
Apply
12.
Properties
Action:
Create
Dimension:
0D
Type:
Mass
mass1
Option(s):
Lumped
Mass:
OK
Element 3
Application Region
(In the select menu, select the Point Element
filter.)
Point Element
Add
Apply
13.
Properties
12-12
Action:
Create
Dimension:
0D
Type:
Mass
mass2
Option:
Lumped
WORKSHOP 12
Mass:
OK
Element 4
Application Region
Add
Apply
14.
Load/BCs
Action:
Create
Object:
Displacement
Type:
Nodal
constraint
< 0, , >
OK
Select Application Region...
FEM
Node 3
Select Nodes:
Add
OK
Apply
15.
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
Direct
Complex Eigenvalue...
Number of Desired Roots=
OK
OK
OK
Apply
12-14
WORKSHOP 12
12-15
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.
12-16
WORKSHOP 12
Comparison of Results
20.
Compare the results obtained in the .f06 file with the results on
the following page:
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
TYPE
G
POINT ID.
1
TYPE
G
POINT ID.
1
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
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
Analysis
Action:
Read Output2
Object:
Result Entities
Method:
Translate
prob12.op2
OK
Apply
23.
Results
Form Type:
Basic
2.1-Eigenvectors, Translation
12-19
12-20
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
WORKSHOP 13
Model Description:
Figure 13.1-Car traveling over a speed bump
u
1
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
13-3
-2.0
13-4
Urel.
Urel.
WORKSHOP 13
13-5
ID SEMINAR,PROB13
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK
13-6
WORKSHOP 13
10
13-7
ENDDATA
13-8
10
WORKSHOP 13
13-9
WORKSHOP 13
$
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
13-12
WORKSHOP 13
3.
4.
5.
13-13
Comparison of Results
6.
13-14
WORKSHOP 13
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
7.
8.
Analysis
Action:
Read Output2
Object:
Both
Method:
Translate
prob13.op2
OK
Apply
9.
Results
Form Type:
Select Result Cases
Advanced
(highlight all cases)
Get Results
Select Result
Plot Type:
XY Plot
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
Rename
Apply
Figure 13.3-Nonlinear Forces at Point 1
13-18
WORKSHOP 13
To get the displacement results, go all the way back to the Result
menu.
Select Result
Plot Type:
XY Plot
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
Rename
Apply
Figure 13.5-Vertical Displacement of Point 1
13-20
WORKSHOP 13
13-21
13-22
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.
14a-1
14a-2
WORKSHOP 14a
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
14a-3
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
WORKSHOP 14a
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
14a-5
14a-6
WORKSHOP 14a
ID SEMINAR,PROB1
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK
14a-7
14a-8
10
WORKSHOP 14a
ENDDATA
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
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
<100, 0, 0>
[ 0, 0, 0]
Apply
14a-10
WORKSHOP 14a
Create
Object:
Mesh
Type:
Curve
10
Curve List
Curve 1
Apply
6. Create nodal displacements.
Loads/BCs
Action:
Create
Object:
Displacement
Type:
Nodal
disp1
<0 0 0>
<0 , , >
OK
Select Application Region...
Geometry
Select Geometry Entities
Point 1
Add
OK
Apply
disp2
< , 0 0>
14a-11
<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>
<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
WORKSHOP 14a
Apply
Cancel
8. Define the bar properties.
Properties
Action:
Create
Dimension:
1D
Type:
Beam
bar
m:alum
Material Name
(Select from Material Property Sets box)
ibeam
W1
W2
0.1
t1
0.1
t2
0.1
OK
Bar Orientation
Coord 0.2
OK
Select Members
Curve 1
Add
Apply
14a-13
9.
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
None
.00259
Constraints
OK
OK
Subcase Create...
Available Subcases
Default
Subcase Parameters...
Number of Desired Roots =
OK
Apply
Cancel
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
WORKSHOP 14a
.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
14a-15
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:
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.
Read Output2
Object:
Result Entities
Method
Translate
prob14a.op2
OK
Apply
When the translation is complete bring up the Results form.
Results
Form Type:
Basic
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
14b-2
WORKSHOP 14b
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
14b-3
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
WORKSHOP 14b
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
PL
= -------AE
500 ( 100 )
= --------------------------------6
0.38 ( 10 10 )
= 0.0132
14b-5
14b-6
WORKSHOP 14b
ID SEMINAR,PROB1
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
__________________________________________________________
BEGIN BULK
MSC/NASTRAN 102 Exercise Workbook
14b-7
14b-8
10
WORKSHOP 14b
ENDDATA
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
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
pull
Input Data...
Force <F1 F2 F3>
<500, , >
OK
Select Application Region...
Select Geometry Entities
Add
OK
Apply
14b-10
Point 2
WORKSHOP 14b
5.
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
None
.00259
Contraints
OK
OK
Direct Text Input...
METHOD = 10
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
14b-12
WORKSHOP 14b
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
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.
345
0.
0.
WORKSHOP 14b
7b.
8.
9.
Comparison of Results
10.
Compare the results obtained in the .f06 file with the results
on the following page:
14b-15
14b-16
MODE
NO.
EXTRACTION
ORDER
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
Read Output2
Object:
Result Entities
Method
Translate
prob14b.op2
OK
Apply
When the translation is complete bring up the Results form.
Results
Form Type:
Basic
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.
14b-17
14b-18
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.
15-1
15-2
WORKSHOP 15
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
WORKSHOP 15
15-5
ID SEMINAR,PROB5
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
__________________________________________
CEND
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
___________________________________________________
BEGIN BULK
15-6
WORKSHOP 15
10
15-7
ENDDATA
15-8
10
WORKSHOP 15
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
+
+DP1
+DP2
+DP3
WORKSHOP 15
$
$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
15-11
Design Variable
15-12
Initial Value
Optimization Value
Iteration Value
_________
______________
___________
_________
______________
___________
_________
______________
___________
WORKSHOP 15
Comparison of Results
5. Compare the results obtained in the .f06 file with the results on the
following page:
15-13
15-14
***************************************************************
SUMMARY OF DESIGN CYCLE HISTORY
***************************************************************
(HARD CONVERGENCE ACHIEVED)
----------------------------------------------------------------------------------------------------------------------------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
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
15-16