Design Sensitivity and Optimization in MSC/NASTRAN
Design Sensitivity and Optimization in MSC/NASTRAN
MacNeal-Schwendler
Corporation
815 Colorado Boulevard Tel: (213) 258-9111
Los Angeles, California 90041-1777 FAX: (213) 259-3838
SEMINAR NOTES
NA∗V70.5∗Z∗Z∗Z∗SM-NAS107-NT1
February 1999
DISCLAIMER
The concepts, methods, and examples presented in this text are for educational
purposes only and are not intended to be exhaustive or to apply to any particular
engineering problem or design. The MacNeal-Schwendler Corporation assumes no
liability or responsibility to any person or company for direct or indirect damages
resulting from the use of any information contained herein.
Printed in U.S.A.
1998 by The MacNeal-Schwendler Corporation
All rights reserved.
The MacNeal-Schwendler Corporation
Seminar Questionnaire
Company _____________________________________________________________________
Your feedback is very important in helping MSC meet the training needs of its clients. Your specific
written comments are especially helpful. Please use the following scale:
1 = poor 2 = marginal 3 = average 4 = good 5 = very good
Effectiveness of instructor in
explaining technical concepts.
Organization of material in a
clear and understandable
manner.
(Continued on back)
Rev 5/93
The MacNeal-Schwendler Corporation
Seminar Questionnaire (Cont.)
Impression of hotel
accommodations. At which
hotel did you stay? Any
problems?
Would you recommend this instructor for future classes? Why or why not? ______________________
__________________________________________________________________________________
Rev 5/93
TABLE OF CONTENTS
TITLE SECTION
1.0 INTRODUCTION
iii
TABLE OF CONTENTS
TITLE SECTION
4.0 EXAMPLES I
7.0 EXAMPLES II
iv
TABLE OF CONTENTS
TITLE SECTION
v
TABLE OF CONTENTS
TITLE SECTION
vi
TABLE OF CONTENTS
TITLE SECTION
DESVAR . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-19
DLINK. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-20
DOPTPRM . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-21
DRESP1. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-25
DRESP2. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-31
DSAPRT. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-34
DSCREEN . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-37
DTABLE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-39
DVBSHAP . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-40
DVGRID . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-41
DVPREL1. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-43
DVPREL2. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-46
DVSHAP . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-48
MODTRAK . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-49
vii
SECTION 1
INTRODUCTION
• Configuration evaluations
• Sensitivity analysis
1-1
BASIC FEATURES IMPLEMENTED IN MSC/NASTRAN
1-2
MSC/NASTRAN Implementation of Structural
Optimization
Initial
Design
Improved
Design
Constraint Approximate
Optimizer
Screening Model
Finite Element
Analysis
1-3
Concepts Presented in this Course
• Optimization
• Sensitivity Analysis
• Objective Function
• Constrants
• Design Model
• Approximate Problem
• Constraint Screening
• Move Limits
• Convergence
1-4
STRENGTHS OF MSC/NASTRAN STRUCTURAL
OPTIMIZATION
• Continuous enhancements
1-5
GENERAL FUNCTIONS
• Solution Sequence
• Analysis Types
Statics
Normal modes
Buckling
Direct frequency*
Modal frequency*
Modal transient*
Static aeroelastic
Aeroelastic flutter
• Design Variables
1-6
SUMMARY OF NEW CAPABILITIES
Version 69
• Mode Tracking
• Buckling
• Normal Modes
• Flutter
1-7
SUMMARY OF NEW CAPABILITIES (Cont.)
Version 70
• NVH Applications
1-8
SUMMARY OF NEW CAPABILITIES (Cont.)
Version 70.5
1-9
BASIC OPTIMIZATION PROBLEM STATEMENT
• Design variables
• Objective function:
Minimize F(X)
• Subject to:
• Inequality constraints:
Gj (X) ≤ 0 j = 1, 2, . . ., L
• Equality constraints:
Hk (X) = 0 k = 1, 2, . . ., M
• Side constraints:
x iL ≤ x i ≤ x iu i = 1 , 2 , . . ., N
1-10
EXAMPLE FORMULATIONS
Example 1:
P2
l1 = 10 cm K1 = 8N/cm
X2
P2 = 5N P1
P1 = 5N
X1
l2 = 10 cm
K2 = 1N/cm
Undeformed Deformed
1-11
EXAMPLE FORMULATIONS (Cont.)
1 1
PE = --- K 1 ( ∆l 1 ) 2 + --- K 2 ( ∆l 2 ) 2 – P 1 X 1 – P 2 X 2
2 2
1 2
PE ( X 1, X 2 ) = --- K 1 [ X 12 + ( l 1 – X 2 ) 2 – l 1 ]
2
1 2
+ --- K 2 [ X 12 + ( l 2 – X 2 ) 2 – l 2 ] – P 1 X 1 – P 2 X 2
2
1-12
EXAMPLE FORMULATIONS (Cont.)
t3
t2
t1
h
P
l
---
4
l b
σ tube ≤ σ buck
t 1 ≥ t 1min, 2t 1 ≤ b
t 2 ≥ t 2min, 2t 2 ≤ h
t 3 ≥ t 3min, t3 ≤ b
1-13
EXAMPLE FORMULATIONS (Cont.)
• Minimize:
2
ρv = ρt 3 --- + ρl [ bh – ( h – 2t 2 ) ( b – 2t 1 ) ]
l
4
T Design
where: { x } = [ b h t 1 t 2 t 3 ] Variables
σ min plate – σ plate ≤ 0
2t 1 – b ≤ 0
Design
2t 2 – h ≤ 0
Constraints
t3 – b ≤ 0
t1 – t1 ≤ 0
min
t2 – t2 ≤ 0 Side
min Constraints
t3 – t3 ≤ 0
min
1-14
EXAMPLE FORMULATIONS (Cont.)
{ Y }1 , { Y }2 , { Y }3
where { Y } i = { A 1 , A 2 , r 1 , r 2 , t, l 1 , l 2 }
i
• Proposed design:
{ Y′ } = α 1 { Y } 1 + α 2 { Y } 2 + α 3 { Y } 3
• Minimize:
W = F ( α1 , α2 , α3 )
• Subject to:
g j = g j ( α1 , α2 , α3 ) ≤ 0
hk = hk ( α1 , α2 , α3 ) = 0
1-15
ONE BAR SUBJECT TO COMPRESSIVE LOAD
• Design problem:
P = 1.0E4 N
L = 500 mm
2
Ainit = 75 mm
ρ, A, E
2
E = 7.0E4 N/mm
L
3
ρ = 2.7E–6 kg/mm
2
σmax = 200 N/mm
Optimal Design:
P
σ max = -------------
A opt
2
A opt = 1.0E4 N ⁄ 200 N ⁄ mm
2
= 50 mm
1-16
ONE BAR SUBJECT TO COMPRESSIVE LOAD (Cont.)
ID MSC, D200X9
TIME 5
SOL 200
CEND
TITLE = STRUCTURAL OPTIMIZATION TEST PROBLEM (STATIC LOADING)
SUBTITLE = SINGLE BAR, 1 DOF, SUBJECT TO COMPRESSIVE LOAD
STRESS = ALL
DISP = ALL
SPC = 1
LOAD = 1
ANALYSIS = STATICS
DESOBJ(MIN) = 20 $ OBJECTIVE FUNCTION DEFINITION
DESSUB = 100 $ CONSTRAINT SET SELECTION
BEGIN BULK
$
$ ANALYSIS MODEL
$
$ units are N, kg, mm:
$
GRID, 1, , 0.0, 0.0, 0.0
GRID, 2, , 0.0, 0.0, 500.0
SPC, 1, 1, 123456
SPC, 1, 2, 12456
CROD, 1, 1, 1, 2
PROD, 1, 1, 75.
$ mm**2
MAT1, 1, 7.0E4, , 0.3, 2.7E-6
$ N/mm**2 kg/mm**3
FORCE, 1, 2, , 1.0, 0.0, 0.0, -1.0E4
$ N
$
$ DESIGN MODEL DESCRIPTION
$
DESVAR, 1, A, 75., 10., 150.
DVPREL1,1, PROD, 1, 4, , , , , +DVP1
+DVP1, 1, 1.0
DRESP1, 10, STRESS, STRESS, PROD, , 2, , 1
DRESP1, 20, WEIGHT, WEIGHT
DCONSTR,100, 10, -200., 200.
$ N/mm**2
DSCREEN,STRESS, -0.5, 1
DOPTPRM,IPRINT, 2, DESMAX, 10, DELP, 0.5, P1, 1 +
+, P2, 15
ENDDATA
1-17
ONE BAR SUBJECT TO COMPRESSIVE LOAD (Cont.)
***** NORMAL CONVERGENCE CRITERIA SATISFIED ***** (HARD CONVERGENCE DECISION LOGIC)
**************************************************************************************
CONVERGENCE ACHIEVED BASED ON THE FOLLOWING CRITERIA
(HARD CONVERGENCE DECISION LOGIC)
***************************************************************
S U M M A R Y O F D E S I G N C Y C L E H I S T O R Y
***************************************************************
(HARD CONVERGENCE ACHIEVED)
(SOFT CONVERGENCE ACHIEVED)
NUMBER OF FINITE ELEMENT ANALYSES COMPLETED 3
NUMBER OF OPTIMIZATIONS W.R.T. APPROXIMATE MODELS 2
OBJECTIVE AND MAXIMUM CONSTRAINT HISTORY
---------------------------------------------------------------------------------------------------------------
OBJECTIVE FROM OBJECTIVE FROM FRACTIONAL ERROR MAXIMUM VALUE
CYCLE APPROXIMATE EXACT OF OF
NUMBER OPTIMIZATION ANALYSIS APPROXIMATION CONSTRAINT
---------------------------------------------------------------------------------------------------------------
INITIAL 1.012500E-01 -3.333334E-01
1 6.746810E-02 6.747905E-02 -1.623075E-04 3.104401E-04
2 6.747905E-02 6.747905E-02 0.000000E+00 3.104401E-04
---------------------------------------------------------------------------------------------------------------
1-19
1 STRUCTURAL OPTIMIZATION TEST PROBLEM (STATIC LOADING) APRIL 20, 1995 MSC/NASTRAN 8/17/94 PAGE 38
SINGLE BAR, 1 DOF, SUBJECT TO COMPRESSIVE LOAD
NUMERICAL OPTIMIZATION
GRADIENTS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-4
F = F(X)
∗ dF
extremum is defined for X at which -------- = 0
dX
′
( X 0 , F 0 ), ( X 0 , F 0 ), ( X 1 , F 1 )
2
F˜ ( X ) = a 0 + a 1 X + a 2 X
2-1
MINIMIZATION OF FUNCTIONS OF A SINGLE
VARIABLE (Cont.)
′
F 0 = a 1 + 2a 2 X 0
2
F0 = a0 + a1 X0 + a2 X0
2
F1 = a0 + a1 X1 + a2 X1
From which:
′
( F1 – F0 ) ⁄ ( ( X1 – X0 ) – F0 )
a 2 = -----------------------------------------------------------------------
( X1 – X0 )
′
a 1 = F 0 – 2a 2 X 0
2
a0 = F0 – a1 X0 – a2 X0
Therefore
˜ ( X ) = a + 2a X = 0
F′ 1 2
∗ a1
X ≅ – ----------
2a 2
F˜ ′′ ( X ) = 2a 2
∗
Minimum at X if a2 > 0
∗
Maximum at X if a2 < 0
2-2
MINIMIZATION OF FUNCTIONS OF A SINGLE
VARIABLE (Cont.)
• Unconstrained minimum:
F(X)
Fopt
∗ X
X0 X X1
• Constrained minimum:
F(X)
u
X≤X
Fopt
u X
X0 X X1
u
Side constraint X≤X
u
X∗ = X
2-3
GRADIENTS
Let:
f = f ( X 1 , X 2 , …, X n )
∂f ∂f ∂f
∇f = ---------- ⋅ e 1 + ---------- ⋅ e 2 + … + ---------- ⋅ e n
∂X 1 ∂X 2 ∂X n
th
where e 1 , e 2 , …, e n are unit vectors in the n direction.
∂ ∂ ∂
∇ = ---------- e 1 + ---------- e 2 + … + ---------- e n
∂X 1 ∂X 2 ∂X n
∂
∇ = --------e i
∂X i
2-4
GRADIENTS (Cont.)
2 2
Example: f ( X 1 , X 2 ) = X 1 + X 2
(Note that this is the equation for a family of circles in the X1, X2
plane with centers at the origin and radii of f .)
∂f ∂f
∇f ( X 1 , X 2 ) = ---------- e 1 + ---------- e 2
∂X 1 ∂X 2
= 2X 1 e 1 + 2X 2 e 2
2X 1
or ∇f ( X 1 , X 2 ) =
2X 2
Graphically: X2
∇f ( X 1, X 2 )
X1
2-5
UNCONSTRAINED MINIMIZATION METHODS
• Theoretical solution:
2
1 2 2
PE = --- K 1 X1 + ( l 1 – X2 ) – l 1
2
2
1 2 2
+ --- K 2 X1 + ( l 2 – X2 ) – l 2 – P1 X1 – P2 X2
2
minimum occurs at
∂ ( PE )
---------------
∂X
-
1
∇ ( PE ) = = 0
∂ ( )
----------------
PE
∂X
2
PEmin = −41.81 N − cm
X1∗ = 8.631 cm
X2∗ = 4.533 cm
2-6
UNCONSTRAINED MINIMIZATION METHODS (Cont.)
• Steepest descent:
2
X
8
100
80 6
60 X*
0 +
40 X
1 4
S
−40
30 −20 −30
1 20 −10
8
X 2 3 X
5
10 X 6
2 X
S 2
X 4 1
X X
−4 −2 0 2 4 6 8 10 12
−2
q q–1 q
X = X + α∗ S
2-7
UNCONSTRAINED MINIMIZATION METHODS (Cont.)
where
q–1
∂F
----------
∂X 1
∂F
q q–1 ----------
S = – ∇F ( X ) = – ∂X 2
.
.
.
∂F
----------
∂X 1
2-8
UNCONSTRAINED MINIMIZATION METHODS (Cont.)
First iteration (q = 1)
q q–1
S = – ∇F ( X )
th
q iteration (q > 1)
q q–1 q–1
S = – ∇F ( X ) + βS
2
X
100
8
3
X
80 6
60
0
40 X +
1 4 X*
S
−40
30 −20 −30
1 20 8 −10
X
S
2
10 2 5
3 3 X
X
2S
X 1
−2 X
−4 0 2 4 6 8 10 12
−2
2-9
KUHN-TUCKER CONDITIONS FOR OPTIMALITY
∗
1. X is feasible
∗ j = 1, …, m
gj ( X ) ≤ 0
∗ k = 1, …, l
hk ( X ) = 0
∗ j = 1, …, m
2. λj g ( X ) = 0
j
λj ≥ 0
m l
∗ ∗ ∗
3. ∇F ( X ) + ∑ λj ∇gj ( X )+ ∑ λk + m ∇hk ( X ) = 0
j=1 k=1
2-10
KUHN-TUCKER CONDITIONS FOR
OPTIMALITY (Cont.)
X2
F ( X ) = constant
∗
∇F ( X )
g1 ( X ) = 0
X1
∗
∇g 2 ( X )
∗
∇g 1 ( X ) g2 ( X ) = 0
∗
∇F ( X )
∗
X
∗ ∗
λ 2 ∇g 2 ( X ) λ 1 ∇g 1 ( X )
∗
–∇ F ( X )
2-11
SIMPLE CANTILEVER EXAMPLE
• Problem description
P = 2250 N
A
L = 500 cm
Section A-A
6 N
E = 1 × 10 ------------
2
cm
2-12
SIMPLE CANTILEVER EXAMPLE (Cont.)
Minimize V = B • H • L
Subject to:
Mc
σ = -------- ≤ 700 Bending Stress
l
3
PL Tip Deflection
δ = ----------- ≤ 2.54
3El
H
---- ≤ 12 Aspect Ratio
B
1 ≤ B ≤ 20
Gauge Requirements
20 ≤ H ≤ 50
2-13
SIMPLE CANTILEVER EXAMPLE (Cont.)
H/B = 12
65
60
V = 200000
175000
55 100000 150000
Height H (cm)
125000
H = 50
50
Optimum
45
δ = 2.54
40 75000
50000
σb = 700
35
2.5 3 4 5 6 7 7.5
Width B (cm)
2-14
SERIES APPROXIMATIONS
2 2
o o df ∆x
f ( x + ∆x ) = f ( x ) + ------- d
⋅ ∆x + ---------f- ⋅ ---------- + …
dx x o 2 2!
dx x o
˜f ( x o + ∆x ) = f ( x o ) + df 2
------- o ⋅ ∆x + 0 ( ∆x )
dx x
where
2 2
0 ( ∆x ) ≡ error on the order of ∆x
2-15
SERIES APPROXIMATIONS (Cont.)
• Minimize V = B • H • L
• Subject to:
MC 6PL N
σ = --------- = ------------ ≤ 700 ------------
I 2 2
BH cm
3 3
PL 4PL
σ = ----------- = ---------------- ≤ 2.54 cm
3EI 3
BH E
2-16
SERIES APPROXIMATIONS (Cont.)
• First-order approximations:
o o o o
V˜ ( B + ∆B, H + ∆H, L ) = V ( B , H , L )
∂V ∂V
- o o ⋅ ∆B + ------- o o ⋅ ∆H
+ ------
∂B B , H ∂H B , H
o o o o
σ̃ ( B + ∆B, H + ∆H, L ) = σ̃ ( B , H , L )
∂σ ∂σ
- o o ⋅ ∆B + ------- o o ⋅ ∆H
+ ------
∂B B , H ∂H B , H
o o o o
δ̃ ( B + ∆B, H + ∆H, L ) = δ̃ ( B , H , L )
∂δ ∂δ
- o o ⋅ ∆B + ------- o o ⋅ ∆H
+ ------
∂B B , H ∂H B , H
o o
• At ( B , H ) = ( 6, 45 )
o o 5 4 3
V ( B + ∆B, H + ∆H, L ) = 1.35 × 10 + 2.25 × 10 ∆B + 3.0 × 10 ∆H
o o
σ ( B + ∆B, H + ∆H, L ) = 555.56 – 92.593∆B – 24.691∆H
o o
δ ( B + ∆B, H + ∆H, L ) = 2.0576 – 0.34294∆B – 0.13717∆H
2-17
SERIES APPROXIMATIONS (Cont.)
65 H/B = 12
60
175000
55 100000 150000
125000
Height H (cm)
H = 50
50
Optimum
Approximate
Optimum
0
45 X
75000
40 δ̃
σ̃
35
2.5 3 4 5 6 7 7.5
Width B (cm)
2-18
SEQUENTIAL LINEAR PROGRAMMING
• First linearize:
o o o
F˜ ( x + ∆x ) = F ( x ) + ∇F ( x ) ⋅ ∆x
o o o
g̃ j ( x + ∆x ) = g j ( x ) + ∇g j ( x ) ⋅ ∆x j ∈ J, (J ≡ set of
active constraints)
o
∆x = x – x
2-19
SEQUENTIAL LINEAR PROGRAMMING (Cont.)
x2
o
∇F ( x )
o
x
1
x
∗
x
2
x
x1
2-20
SEQUENTIAL LINEAR PROGRAMMING (Cont.)
X2 move limits
o
F˜ ( X )
true optimum
o
g̃ ( X ) g(X) = 0
0 X1
• Move limits are required because higher order terms have been
ignored
2-21
METHOD OF FEASIBLE DIRECTIONS
X2 o
∇F ( X )
o
∇g j ( X ) F = Constant
gj ( X )
X1
∇F ⋅ S ≤ 0
∇g j ⋅ S ≤ 0
X2 Feasible
Usable/Feasible
S
Usable gj ( X )
X1
2-22
METHOD OF FEASIBLE DIRECTIONS (Cont.)
• To maintain feasibility:
∇g j ⋅ S + θ j ≤ 0 (1)
• To maintain usability:
∇F ⋅ S + β ≤ 0 (2)
2-23
METHOD OF FEASIBLE DIRECTIONS (Cont.)
H/B = 12
60
55
Height H (cm)
H = 50
50
Optimum
45 o
X
δ = 2.54
40
σ = 700
35
3 4 5 6 7
Width B (cm)
2-24
METHOD OF FEASIBLE DIRECTIONS (Cont.)
• The problem with the method is that it may tend to zigzag between
active and inactive constraints as the optimum is approached.
X2
g1 ( X )
o
X g2 ( X )
1
X
3 F ( X ) = Constant
X 4 2
X X
X1
2-25
MODIFIED METHOD OF FEASIBLE DIRECTIONS
Maximize – ∇F ⋅ s
s⋅s≤1
2-26
MODIFIED METHOD OF FEASIBLE
DIRECTIONS (Cont.)
maximize – ∇F ⋅ s
subject to:
∇g j ( X ) ⋅ s ≤ 0
for j ∈ J
s⋅s≤1
2-27
MODIFIED METHOD OF FEASIBLE
DIRECTIONS (Cont.)
q q–1 q
X = X + αS
T T –1
δX = A [ AA ] G
T
where A contains the gradients of the active constraints and G
contains the constraint values.
2-28
MODIFIED METHOD OF FEASIBLE
DIRECTIONS (Cont.)
• Features
2-29
MODIFIED METHOD OF FEASIBLE
DIRECTIONS (Cont.)
H/B = 12
60
55
Height H (cm)
H = 50
50
Optimum
o
45 X
δ = 2.54
40
σ = 700
35
3 4 5 6 7
Width B (cm)
2-30
MODIFIED METHOD OF FEASIBLE
DIRECTIONS (Cont.)
x2
F = Constant
∇F
Minimize Constraint
Violation
o
∇G x
x1
Minimize Objective
Function Compromise Direction
x2
∗
x
x1
o
x
No Feasible Design
2-31
NUMERICAL IDENTIFICATION OF ACTIVE AND
VIOLATED CONSTRAINTS
• Parameters of interest
x2
Feasible Region
gj ( x ) < 0
g j ( x ) = CT
Infeasible Region gj ( x ) = 0
gj ( x ) > 0 g j ( x ) = CTMIN
x1
2-32
NUMERICAL IDENTIFICATION OF ACTIVE AND
VIOLATED CONSTRAINTS (Cont.)
Gj(X)
Feasible Violated
Active
CTMIN
X
0
CT Numerical Constraint
Boundary
Exact Constraint
Boundary
CT = −0.03
CTMIN = 0.003
2-33
SECTION 3
DESIGN MODELING I
DESVAR . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-8
DVPREL1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-11
DRESP1. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-16
DESOBJ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-25
DCONSTR . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-28
DCONADD. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-31
DESSUB . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-33
DESGLB . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-34
DOPTPRM . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-39
WHAT IS A DESIGN MODEL? 3
It must:
3-1
HOW ARE THE DESIGN AND ANALYSIS MODELS
RELATED?
DESIGN MODEL
ANALYSIS MODEL
3-2
DESIGN AND ANALYSIS MODEL COMPARISON
h tw tf
...
...
...
...
PBAR 21 6 2.36 ...
h, b, tw, tf
3-3
DESIGN AND ANALYSIS MODEL
COMPARISON (Cont.)
A = 2∗b∗t f + ( h – 2t f )∗t w
3 (b – t )
b∗h w 3
I 1 = -------------- – --------------------∗ ( h – 2∗t f )
12 12
3-4
DESIGN MODEL DEFINITION PROCESS
3-5
DEFINING THE ANALYSIS DISCIPLINES
• Executive Section
SOL 200
STATICS – statics
MODES – normal modes
BUCK – buckling
DFREQ – direct frequency*
Analysis =
MFREQ – modal frequency*
MTRAN – modal transient*
SAERO – static aeroelasticity
FLUTTER – flutter
*Includes acoustics
3-6
DEFINING THE DESIGN VARIABLES
The DESVAR entry defines a design variable, its initial value, and
bounds. It also can provide design cycle move limits.
3-7
DESVAR Design Variable
Format:
1 2 3 4 5 6 7 8 9 10
DESVAR ID LABEL XINIT XLB XUB DELXV
Example:
Field Contents
ID Unique design variable identification number. (Integer > 0)
DELXV Fractional change allowed for the design variable during approximate
optimization. (Real > 0.0, for Default see Remark 2.)
Remarks:
1. DELXV can be used to control the change in the design variable during one optimization
cycle.
2. If DELXV is blank, the default is taken from the specification of the DELX parameter on
the DOPTPRM entry. If DELX is not specified, then the default is 1.0.
3-8
DEFINING THE DESIGN VARIABLES (Cont.)
Example:
0.01 ≤ x 10 ≤ 0.1
0.01 ≤ x 11 ≤ 0.08
T
X = [ x 1, x 2, …, x n ]
I u
xi ≤ xi ≤ xi
3-9
RELATING DESIGN VARIABLES TO PROPERTIES
pj = Co + ∑ Ci xi
i
3-10
Design Variable to Property Relation DVPREL1
DVPREL1 the relation between an analysis model property and design variables.
Defines Design Variable to Property Relation
Format:
1 2 3 4 5 6 7 8 9 10
DVPREL1 ID TYPE PID FID PMIN PMAX C0
Example:
Field Contents
ID Unique identification number. (Integer > 0)
FlD Field position of the property entry, or word position in the element property table
of the analysis model. (Integer ≠ 0)
PMIN Minimum value allowed for this property. If FID references a stress recovery
location, then the default value for PMIN is –1.0+35. PMIN must be explicitly set
to a negative number for properties that may be less than zero (for example, field
ZO on the PCOMP entry). (Real; Default = 1.0E–20)
PMAX Maximum value allowed for this property. (Real; Default = 1.0E–20)
Remarks:
1. The relationship between the analysis model property and design variables is given by:
P i = C0 + ∑ COEFi ∗ DVIDi
i
(Continued)
3-11
DVPREL1 Design Variable to Property Relation
3. PTYPE = “PBEND” is not supported, either directly through FIDs or indirectly via word po-
sitions in the element property table.
4. FID may be either a positive or a negative number. If FID > 0, it identifies the field position
on a property entry. If FID < 0, it identifies the word position of an entry in the element
property table. For example, to specify the area of a PBAR, either FID = +4 or FID = –3
can be used. However, if PTYPE = “PBEAM”, FID must be negative. See the following
element property table for the word positions for PBEAM.
1 I Property ID
2 I Material ID
3 I Number of segments
4 I Constant cross section flag (1 = yes, 0 = no)
5 I Unused
6 I Stress output request flag, SO (1 = yes, 0 = no)
7 R X/XB ratio; at end A, X/XB = 0.0
8 through 13 R A, I1, I2, I12, J, NSM
14 through 21 R C1, C2, D1, D2, E1, E2, F1, F2
22 I
23 R Repeat of words 6 through 21 for the 1st intermediate
24 through 29 R station
30 through 37 R
38 I
39 R
40 through 45 R 2nd intermediate station
46 through 53 R
54 I
55 R
56 through 61 R 3rd intermediate station
62 through 69 R
70 I
71 R
72 through 77 R 4th intermediate station
78 through 85 R
86 I
87 R
88 through 93 R 5th intermediate station
94 through 101 R
(Continued)
3-12
Design Variable to Property Relation DVPREL1
102 I
103 R
104 through 109 R 6th intermediate station
110 through 117 R
118 I
119 R
120 through 125 R 7th intermediate station
126 through 133 R
134 I
135 R
136 through 141 R 8th intermediate station
142 through 149 R
150 I
151 R
152 through 157 R 9th intermediate station
158 through 165 R
166 I
167 R
168 through 173 R End B
174 through 181 R
182 through 189 R K1, K2, S1, S2, NSI(A), NSI(B), CW(A), CW(B)
190 through 197 R M1(A), M2(A), M1(B), M2(B), N1(A), N2(A), N1(B),
N2(B)
3-13
RELATING DESIGN VARIABLES TO
PROPERTIES (Cont.)
Example:
Web
t i = 1.0 x 2
Web Cap
A i = 1.0 x 3
Plate
t i = 1.0 x 1
PSHELL, 1, 1, 0.15, 1
$
$...Define the design variables:
$
$DESVAR,ID, LABEL, XINIT, XLB, XUB, DELXV
DESVAR, 1, T-PLATE,0.15, 0.001, 10.0
$
$...Relate the design variables to analysis model properties
$ (linear relations, so use DVPREL1)
$
$DVPREL1,ID, TYPE, PID, FID, PMIN, PMAX, C0, , +
$+, DVID1, COEF1, DVID2, COEF2, ...
DVPREL1,1, PSHELL, 1, 4, 0.01, , , , +DP1
+DP1, 1, 1.0
3-14
IDENTIFYING THE DESIGN RESPONSES
3-15
DRESP1 Design Sensitivity Response Quantities
Format:
1 2 3 4 5 6 7 8 9 10
DRESP1 ID LABEL RTYPE PTYPE REGION ATTA ATTB ATT1
ATT2 -etc.-
Example:
103
Field Contents
ID Unique entry identifier. (Integer > 0)
PTYPE Element flag (PTYPE = “ELEM”) or property entry name. Used with element type
responses (stress, strain, force, etc.) to identify the property type, since property
entry IDs are not unique across property types. (Character: “ELEM”, “PBAR”,
“PSHELL”, etc.)
REGION Region identifier for constraint screening. See Remark 10 for defaults.
(Integer > 0)
ATTA, ATTB, Response attributes. See Table 1. (Integer > 0 or Real or blank)
ATTi
(Continued)
3-16
Design Sensitivity Response Quantities DRESP1
Table 1. Design Sensitivity Response Attributes.
Response Attributes
Response
Type ATTB (Integer > 0 or
(RTYPE) ATTA (Integer > 0) Real > 0.0) ATTi (Integer > 0)
FRSTRE Stress Item Code Frequency Value. (Real > Property entry (PID)
0.0) See Remark 15..
(Continued)
3-17
DRESP1 Design Sensitivity Response Quantities
Response Attributes
Response
Type ATTB (Integer > 0 or
(RTYPE) ATTA (Integer > 0) Real > 0.0) ATTi (Integer > 0)
FRFORC Force Item Code Frequency Value. (Real > Property entry (PID)
0.0) See Remark 15.
TSTRE Stress Item Code Time Value. (Real) See Property entry (PID)
Remark 16.
TFORC Force Item Code Time Value. (Real) See Property entry (PID)
Remark 16.
Remarks:
1. Stress, strain, and force item codes can be found in the MSC/NASTRAN Quick Reference
Guide, Appendix A. For stress or strain item codes that have dual meanings, such as von
Mises or maximum shear, the option specified in the Case Control Section will be used;
i.e., STRESS(VONM) or STRESS(MAXS).
2. RTYPE=“CSTRESS”, “CSTRAIN”, and “CFAILURE” are used only with the PCOMP entry.
“CSTRESS” and “CSTRAIN” item codes are described under Table 1. (Element
Stress/Strain Item Codes) in the MSC/NASTRAN Quick Reference Guide, Appendix A.
“CFAILURE” item codes are described under Table Table 2. (Element Force Item Codes)
in the MSC/NASTRAN Quick Reference Guide, Appendix A. Only force item codes that
refer to failure indices of direct stress and interlaminar shear stress are valid.
(Continued)
3-18
Design Sensitivity Response Quantities DRESP1
The CFAILURE response type requires the following specifications on the applicable
entries:
c. Stress limits in the ST, SC, and SS fields on all MATi entries.
3. ATTB is used only for responses of composite laminae, dynamics, and stability deriva-
tives. For other responses, this field must be blank.
4. All grids associated with a DRESP1 entry are considered to be in the same region for
screening purposes. Only up to NSTR displacement constraints (see DSCREEN entry)
per group per load case will be retained in the design optimization phase.
10. REGION is used for constraint screening. The NSTR field on DSCREEN entries gives the
maximum number of constraints retained for each region per load case.
If the REGION field is not blank, all the responses on this entry as well as all responses on
other DRESP1 entries that have the same RTYPE and REGION identification number will
be grouped into the same region.
(Continued)
3-19
DRESP1 Design Sensitivity Response Quantities
WEIGHT No region
VOLUME No region
LAMA No region
EIGN No region
FREQ No region
11. REGION is valid only among the same type of responses. Responses of different types
will never be grouped into the same region, even if they are assigned the same REGION
identification number by the user.
12. If RTYPE = “WEIGHT” or “VOLUME”, field ATTi = “ALL” implies total weight/volume of all
superelements except external superelements.
13. RTYPE = “STABDER” identifies a stability derivative response. ATTB is the restraint flag
for the stability derivative. ATTB = 0 means unrestrained, and ATTB = 1 means re-
strained. For example, ATTA = 4000, ATTB = 0, and ATT1 = 3 reference the unrestrained
Cz derivative for the AESTAT (or AESURF) entry ID = 4000.
(Continued)
3-20
Design Sensitivity Response Quantities DRESP1
14. RTYPE = “FLUTTER” identifies a set of damping responses. The set is specified by ATTi:
ATT3 = Identification number of an FLFACT entry that specifies a list of Mach numbers.
15. For RTYP = “FRDISP”, “FRVELO”, “FRACCL”, “FRSPCF”, “FRFORC”, and “FRSTRE”,
ATTB specifies a frequency value in cycles per unit time. If ATTB is specified, then the
responses are evaluated at the closest frequency selected by the OFREQ command. The
default for ATTB is all frequencies selected by the OFREQ command.
16. For RTYPE = “TDISP”, “TVELO”, “TACCL”, “TSPCF”, “TFORC”, and “TSTRE”, ATTB
specifies a time value. If ATTB is specified, then the responses are evaluated at the clos-
est time selected by the OTIME command. The default for ATTB is all time steps selected
by the OTIME command.
17. Intermediate station responses on CBAR elements due to PLOAD1 and/or CBARAO en-
tries may not be defined on the DRESP1 entry.
18. RTYPE = “EIGN” refers to normal modes response in terms of eigenvalue (radi-
an/time)∗∗2 while RTYPE = “FREQ” refers to normal modes response in terms of natural
frequency or units of cycles per unit time.
19. For RTYPE = LAMA, EIGN or FREQ, the response approximation used for optimization
can be individually selected. (Approximation Code = 1 = direct linearization, = 2 = Inverse
Linearization).
3-21
IDENTIFYING THE DESIGN RESPONSES (Cont.)
$
$DRESP1, ID, LABEL, RTYPE, PTYPE, REGION, ATTA, ATTB, ATT1, +
$+, ATT2, ...
$
$ TRANSLATIONAL DISPLACEMENTS AT GRIDS 501-505:
$
DRESP1, 201, TXYZ1, DISP, , , 123, , 501, +
+, 502, 503, 504, 505
$
$DRESP1, ID, LABEL, RTYPE, PTYPE, REGION, ATTA, ATTB, ATT1, +
$+, ATT2, ...
$
$ AXIAL STRESSES FOR ROD ELEMENTS IN PID GROUPS 150,160,170,180:
$
DRESP1, 250, SIG1, STRESS, PROD, , 2, , 150, +
+, 160, 170, 180
3-22
IDENTIFYING THE DESIGN RESPONSES (Cont.)
EID1
EID2
PID1 .
.
.
EIDn
PID2
DRESP1 .
.
.
PIDm
3-23
DEFINING THE OBJECTIVE FUNCTION
3-24
Design Objective DESOBJ
DESOBJ the DRESP1 or DRESP2 entry to be used as the design objective.
Selects Design Objective
Format:
MAX
DESOBJ ( ) = N
MIN
Examples:
DESOBJ = 10
DESO = 25
Describer Meaning
MIN Specifies that the objective is to be minimized.
Remarks:
1. A DESOBJ command is required for a design optimization task and is optional for a sensi-
tivity task.
2. If the DESOBJ command is specified within a SUBCASE, the identified DRESPi Bulk Data
entry uses a response only from that subcase. If DESOBJ appears above all SUBCASE
commands and there are multiple subcases, it uses a global response.
3-25
DEFINING THE OBJECTIVE FUNCTION (Cont.)
Example:
3-26
DEFINING THE DESIGN CONSTRAINTS
L U
rj ≤ rj ( x ) ≤ rj
3-27
DCONSTR Design Constraints
DCONSTR
Define design constraints. Design Constraints
Format:
1 2 3 4 5 6 7 8 9 10
DCONSTR DCID RID LALLOW UALLOW
Example:
DCONSTR 10 4 1.25
Field Contents
DCID Design constraint set identification number. (Integer > 0)
Remarks:
1. The DCONSTR entry may be selected in the Case Control Section by the DESSUB or
DESGLB command.
3. For a given DCID, the associated RID can be referenced only once.
4. The units of LALLOW and UALLOW must be consistent with the referenced response de-
fined on the DRESPi entry. If RID refers to an “EIGN” response, then the imposed bounds
2
must be expressed in units of eigenvalue, (radian/time) . If RID refers to a “FREQ” re-
sponse, then the imposed bounds must be expressed in cycles/time.
5. LALLOW and UALLOW are unrelated to the stress limits specified on the MATi entry.
LALLOW – r
g = -------------------------------- for lower bound constraints
GNORM
r – UALLOW
g = --------------------------------- for upper bound constraints
GNORM
where r is the response defined on the DRESPi entry and if the RID refers to a “FREQ”
response, then the imposed bounds must be expressed in cycles/time.
(Continued)
3-28
Design Constraints DCONSTR
7. As Remark 6 indicates, small values of UALLOW and LALLOW require special processing
and should be avoided. Bounds of exactly zero are particularly troublesome. This can be
avoided by using a DRESP2 entry that offsets the constrained response from zero.
3-29
DEFINING THE DESIGN CONSTRAINTS (Cont.)
3-30
Design Constraint Set Combination DCONADD
DCONADDthe design constraints for a subcase as a union of DCONSTRDesign
Defines entries.
Constraint Set Combination
Format:
1 2 3 4 5 6 7 8 9 10
DCONADD DCID DC1 DC2 DC3 DC4 DC5 DC6 DC7
DC8 -etc.-
Example:
DCONADD 10 4 12
Field Contents
DCID Design constraint set identification number. (Integer > 0)
Remarks:
1. The DCONADD entry is selected by a DESSUB or DESGLB Case Control command.
3-31
DEFINING THE DESIGN CONSTRAINTS (Cont.)
3-32
Design Constraints Request at the Subcase Level DESSUB
DESSUBthe design constraints to be used in a design optimization
Select task for the
Design Constraints current
Request subcase.
at the Subcase Level
Format:
DESSUB = n
Examples:
DESSUB = 10
DESS = 25
Describer Meaning
n Set identification of a DCONSTR or DCONADD Bulk Data entry identification
number. (Integer ≥ 0)
Remark:
1. A DESSUB command is required for every subcase for which constraints are to be
applied.
3-33
DESGLB Request Design Constraints at the Global Level
Format:
DESGLB = n
Examples:
DESGLB = 10
DESG = 25
Describer Meaning
n Set identification of a DCONSTR or DCONADD Bulk Data entry identification
number. (Integer > 0)
Remarks:
1. A DESGLB command is optional and invokes constraints that are to be applied indepen-
dent of a particular subcase. These constraints could be based on responses that are
independent of subcases (e.g., WEIGHT or VOLUME).
2. The DESGLB command can be used to invoke constraints that are not a function of
DRESP1 entries; e.g., DRESP2 responses that are not functions of DRESP1 responses
are subcase independent.
3-34
DEFINING THE DESIGN CONSTRAINTS (Cont.)
$
$ IN CASE CONTROL:
$
SUBCASE 20
ANALYSIS = STATICS
DESSUB = 100
(ETC.)
$
$ IN BULK DATA:
$
$DRESP1, ID, LABEL, RTYPE, PTYPE, REGION, ATTA, ATTB, ATT1, +
$+, ATT2, ...
$
$ MAXIMUM STRESS AT END A:
DRESP1, 110, SAMAX, STRESS, PBAR, , 7, , 2010, +
+ 2020, 2030
$
$ MAXIMUM STRESS AT END B:
DRESP1, 120, SBMAX, STRESS, PBAR, , 14, , 2010, +
+ 2020, 2030
$
$ BOUNDS APPLIED TO BOTH ENDS A & B:
$DCONSTR,DCID, RID, LALLOW, UALLOW
DCONSTR,100, 110, -85., +85. $ N/mm**2 untreated struct. steel,
DCONSTR,100, 120, -85., +85. $ FS=2 on fatigue limit
$
3-35
DEFINING THE DESIGN CONSTRAINTS (Cont.)
EID1
EID2
PID1 .
.
.
PID2
EIDn
DCONSTR DRESP1 .
.
.
PIDm
3-36
DEFINING THE DESIGN CONSTRAINTS (Cont.)
L
rj – rj ( x )
g 2j – 1 ( x ) = ------------------------- ≤ 0
L
rj
U
rj ( x ) – rj
g 2j ( x ) = -------------------------- ≤ 0
U
rj
Normalization provides:
3-37
BASIC DESIGN PROCESS CONTROL
3-38
Design Optimization Parameters DOPTPRM
DOPTPRM default values of parameters used in design optimization.
Overrides Design Optimization Parameters
Format:
1 2 3 4 5 6 7 8 9 10
DOPTPRM PARAM1 VAL1 PARAM2 VAL2 PARAM3 VAL3 PARAM4 VAL4
Example:
Field Contents
PARAMi Name of the design optimization parameter. Allowable names are given in
Table 1. (Character)
VALi Value of the parameter. (Real or Integer, see Table 1.)
Remarks:
1. Only one DOPTPRM entry is allowed in the Bulk Data Section.
CONVDV Relative convergence criterion on design variables. (Real > 0.0; Default = 0.001)
CONVPR Relative convergence criterion on properties. (Real > 0.0; Default = 0.001)
(Continued)
3-39
DOPTPRM Design Optimization Parameters
DELB Relative finite difference move parameter. (Real > 0.0; Default = 0.0001)
DELP Fractional change allowed in each property during any optimization design cycle.
This provides constraints on property moves. (Real > 0.0; Default = 0.2)
DELX Fractional change allowed in each design variable during any optimization cycle.
(Real > 0.0; Default = 1.0)
DOBJ1 Relative change in objective attempted on the first optimization iteration. Used to
estimate initial move in the one-dimensional search. Updated as the optimization
progresses. (Real > 0.0; Default = 0.1)
DOBJ2 Absolute change in objective attempted on the first optimization iteration. (Real
> 0.0; Default = 0.2 ∗ (F0))
DPMIN Minimum move limit imposed. (Real > 0.0; Default = 0.01)
DX1 Maximum relative change in a design variable attempted on the first optimization
iteration. Used to estimate the initial move in the one dimensional search. Updated
as the optimization progresses. (Real > 0.0; Default = 0.01)
DX2 Absolute change in a design variable attempted on the first optimization iteration.
(Real > 0.0; Default = 0.2 ∗ MAX[X(I)])
DXMIN Minimum design variable move limit (Real > 0.0; Default = 0.05).
GMAX Maximum constraint violation allowed at the converged optimum. (Real > 0.0;
Default = 0.005)
GSCAL Constraint normalization factor. See Remarks under the DSCREEN and
DCONSTR entries. (Real > 0.0; Default = 0.001)
(Continued)
3-40
Design Optimization Parameters DOPTPRM
IGMAX If IGMAX = 0, only gradients of active and violated constraints are calculated. If
IGMAX > 0, up to NCOLA gradients are calculated including active, violated, and
near active constraints. (Integer > 0; Default = 0)
IPRINT Print control during approximate optimization phase. Increasing values represent
increasing levels of optimizer information. (0 ≤ Integer ≤ 7; Default = 0)
0 no output (Default)
1 internal optimization parameters, initial information, and results
2 same, plus objective function and design variables at each iterations
3 same, plus constraint values and identification of critical constraints
4 same, plus gradients
5 same, plus search direction
6 same, plus scaling factors and miscellaneous search information
7 same, plus one dimensional search information
IPRNT1 If IPRNT1 = 1, print scaling factors for design variable vector. (Integer 0 or 1;
Default = 0)
ISCAL Design variables are rescaled every ISCAL iterations. Set ISCAL= –1 to turn off
scaling. (Integer; Default = NDV (number of design variables))
ITMAX Maximum number of iterations allowed at optimizer level during each design cycle.
(Integer; Default = 40)
ITRMOP Number of consecutive iterations for which convergence criteria must be satisfied
to indicate convergence at the optimizer level. (Integer; Default = 2)
ITRMST Number of consecutive iterations for which convergence criteria must be met at
the optimizer level to indicate convergence in the Sequential Linear Programming
Method. (Integer > 0; Default = 2)
JTMAX Maximum number of iterations allowed at the optimizer level for the Sequential
Linear Programming Method. This is the number of linearized subproblems
solved. (Integer ≥ 0; Default = 20)
JPRINT Sequential Linear Programming subproblem print. If JPRINT > 0, IPRINT is turned
on during the approximate linear subproblem. (Default = 0)
JWRITE If JWRITE > 0, file number on which iteration history will be written. (Integer > 0;
Default = 0)
(Continued)
3-41
DOPTPRM Design Optimization Parameters
P1 Print control items specified for P2. (Integer ≥ 0; Default = 0) Initial results are
always printed prior to the first approximate optimization. If an optimization task
is performed, final results are always printed for the final analysis unless
PARAM,SOFTEXIT,YES is specified. These two sets of print are not controllable.
n: Print at every n-th design cycle.
STPSCL Scaling factor for shape finite difference step sizes, to be applied to all shape
design variables. (Real > 0.0; Default = 1.0)
3-42
BASIC DESIGN PROCESS CONTROL (Cont.)
PARAM,NASPRT,n
3-43
BASIC DESIGN PROCESS CONTROL (Cont.)
PARAM,DESPCH,n
n<0 never
3-44
SECTION 4
EXAMPLES I
• Constraints:
4-1
3-BAR TRUSS (Cont.)
10" 10"
2 3
1
x
A1 A2
3
1 10"
2
A1
Subcase 2 Subcase 1
20,000 lbs 20,000 lbs
x: –16,000 lbs x: 16,000 lbs
y: –12,000 lbs y: –12,000 lbs
4-2
3-BAR TRUSS (Cont.)
ID MSC, D200X1
TIME 10 $
SOL 200 $ OPTIMIZATION
CEND
TITLE = SYMMETRIC THREE BAR TRUSS DESIGN OPTIMIZATION - D200X1
SUBTITLE = BASELINE - 2 CROSS SECTIONAL AREAS AS DESIGN VARIABLES
ECHO = SORT
SPC = 100
DISP = ALL
STRESS = ALL
DESOBJ(MIN) = 20 $ (DESIGN OBJECTIVE = DRESP ID)
DESSUB = 21 $ DEFINE CONSTRAINT SET FOR BOTH SUBCASES
ANALYSIS = STATICS
SUBCASE 1
LABEL = LOAD CONDITION 1
LOAD = 300
SUBCASE 2
LABEL = LOAD CONDITION 2
LOAD = 310
BEGIN BULK
$
$------------------------------------------------------------------------
$ ANALYSIS MODEL
$------------------------------------------------------------------------
$
$ GRID DATA
$ 2 3 4 5 6 7 8 910
GRID, 1, , -10.0 , 0.0, 0.0
GRID, 2, , 0.0 , 0.0, 0.0
GRID, 3, , 10.0 , 0.0, 0.0
GRID, 4, , 0.0 , -10.0, 0.0
$ SUPPORT DATA
SPC, 100, 1, 123456, , 2, 123456
SPC, 100, 3, 123456, , 4, 3456
$ ELEMENT DATA
CROD, 1, 11, 1, 4
CROD, 2, 12, 2, 4
CROD, 3, 11, 3, 4
$ PROPERTY DATA
PROD, 11, 1, 1.0
PROD, 12, 1, 2.0
MAT1, 1, 1.0E+7, , 0.33, 0.1
$ EXTERNAL LOADS DATA
FORCE, 300, 4, , 20000., 0.8, -0.6
FORCE, 310, 4, , 20000., -0.8, -0.6
$
$------------------------------------------------------------------------
$ DESIGN MODEL
$------------------------------------------------------------------------
$
4-3
3-BAR TRUSS (Cont.)
4-4
3-BAR TRUSS (Cont.)
***************************************************************
S U M M A R Y O F D E S I G N C Y C L E H I S T O R Y
***************************************************************
(HARD CONVERGENCE ACHIEVED)
NUMBER OF FINITE ELEMENT ANALYSES COMPLETED 6
NUMBER OF OPTIMIZATIONS W.R.T. APPROXIMATE MODELS 5
OBJECTIVE AND MAXIMUM CONSTRAINT HISTORY
-------------------------------------------------------------------------------------------------
OBJECTIVE FROM OBJECTIVE FROM FRACTIONAL ERROR MAXIMUM VALUE
CYCLE APPROXIMATE EXACT OF OF
NUMBER OPTIMIZATION ANALYSIS APPROXIMATION CONSTRAINT
--------------------------------------------------------------------------------------------------
INITIAL 4.828427E+00 -3.234952E-01
1 3.007897E+00 3.008492E+00 -1.977251E-04 -3.737402E-03
2 2.821953E+00 2.821638E+00 1.118734E-04 -1.967246E-02
3 2.734469E+00 2.734299E+00 6.217039E-05 -7.241016E-03
4-5
--------------------------------------------------------------------------------------------------
• Constraints:
360" 360"
5 3 1
1 2
7 10
5 6
360"
8 9
3 4
6 4 2
4-6
10-BAR TRUSS (Cont.)
ID MSC, D200X2
TIME 10
SOL 200 $ OPTIMIZATION
CEND
$
TITLE = OPTIMIZATION TEST CASE 2
SUBTITLE = TEN BAR TRUSS D200X2
DESOBJ = 11 $ OBJECTIVE FUNCTION DEFINITION
DESSUB = 1 $ CONSTRAINT SET SELECTION
ANALYSIS = STATICS
DISP = ALL
STRESS = ALL
SPC = 100
LOAD = 300
BEGIN BULK
$
$------------------------------------------------------------------------
$ ANALYSIS MODEL:
$------------------------------------------------------------------------
$
GRDSET, , , , , , , 3456
GRID, 1, , 720., 360., 0.
GRID, 2, , 720., 0., 0.
GRID, 3, , 360., 360., 0.
GRID, 4, , 360., 0., 0.
GRID, 5, , 0., 360., 0.
GRID, 6, , 0., 0., 0.
CROD, 1, 1, 5, 3
CROD, 2, 2, 3, 1
CROD, 3, 3, 6, 4
CROD, 4, 4, 4, 2
CROD, 5, 5, 3, 4
CROD, 6, 6, 1, 2
CROD, 7, 7, 5, 4
CROD, 8, 8, 6, 3
CROD, 9, 9, 3, 2
CROD, 10, 10, 4, 1
PROD, 1, 1, 5.0
PROD, 2, 1, 5.0
PROD, 3, 1, 5.0
PROD, 4, 1, 5.0
PROD, 5, 1, 5.0
PROD, 6, 1, 5.0
PROD, 7, 1, 5.0
PROD, 8, 1, 5.0
PROD, 9, 1, 5.0
PROD, 10, 1, 5.0
MAT1, 1, 1.0E7, , , 0.1, , , , +
+, 25000., 25000.
FORCE, 300, 2, , 100000.,0.0, -1.0, 0.0
FORCE, 300, 4, , 100000.,0.0, -1.0, 0.0
4-7
10-BAR TRUSS (Cont.)
4-8
10-BAR TRUSS (Cont.)
DCONSTR,1, 5, -25000.,25000.
DCONSTR,1, 6, -25000.,25000.
DCONSTR,1, 7, -25000.,25000.
DCONSTR,1, 8, -25000.,25000.
DCONSTR,1, 9, -25000.,25000.
DCONSTR,1, 10, -25000.,25000.
DOPTPRM,IPRINT, 2, DESMAX, 20, DELP, 0.8
$DSCREEN STRESS -100. 1
$DSCREEN DISP 1
$.......2.......3.......4.......5.......6.......7.......8.......9.......0
ENDDATA
4-9
10-BAR TRUSS RESULTS
***** NORMAL CONVERGENCE CRITERIA SATISFIED ***** (HARD CONVERGENCE DECISION LOGIC)
**************************************************************************************
CONVERGENCE ACHIEVED BASED ON THE FOLLOWING CRITERIA
(HARD CONVERGENCE DECISION LOGIC)
***************************************************************
S U M M A R Y O F D E S I G N C Y C L E H I S T O R Y
***************************************************************
----------------------------------------------------------------------------------------------------------------------------------
1 | 1 | A1 | 7.6879E+00 : 7.7337E+00 : 7.8566E+00 : 7.8665E+00 : 7.9267E+00 :
2 | 2 | A2 | 4.6169E-01 : 3.6299E-01 : 2.4340E-01 : 1.8105E-01 : 1.1949E-01 :
3 | 3 | A3 | 8.3136E+00 : 8.2719E+00 : 8.1506E+00 : 8.1431E+00 : 8.0690E+00 :
4 | 4 | A4 | 3.6911E+00 : 3.7331E+00 : 3.8545E+00 : 3.8673E+00 : 3.9281E+00 :
5 | 5 | A5 | 1.4644E-01 : 8.7864E-02 : 5.2718E-02 : 3.1658E-02 : 1.8995E-02 :
6 | 6 | A6 | 4.6169E-01 : 3.6299E-01 : 2.4340E-01 : 1.8105E-01 : 1.1949E-01 :
7 | 7 | A7 | 6.1055E+00 : 6.0343E+00 : 5.8732E+00 : 5.8515E+00 : 5.7556E+00 :
8 | 8 | A8 | 5.2173E+00 : 5.2757E+00 : 5.4508E+00 : 5.4628E+00 : 5.5563E+00 :
9 | 9 | A9 | 5.2201E+00 : 5.2795E+00 : 5.4454E+00 : 5.4669E+00 : 5.5591E+00 :
10 | 10 | A10 | 6.4001E-01 : 5.1250E-01 : 3.4406E-01 : 2.5765E-01 : 1.6734E-01 :
*** USER INFORMATION MESSAGE 6464 (DOM12E)
RUN TERMINATED DUE TO HARD CONVERGENCE TO AN OPTIMUM AT CYCLE NUMBER = 10.
10-BAR TRUSS RESULTS (Cont.)
Objective Function
4-13
10-BAR TRUSS RESULTS (Cont.)
4-14
BEAM VIBRATION FREQUENCY AS A DESIGN
CONSTRAINT (M. J. Turner 1970)
A1 A2 A3
6" t1 t2 t3 X
• Objective
• Constraint
• Design variables
4-15
BEAM VIBRATION FREQUENCY AS A DESIGN
CONSTRAINT (Cont.)
nastran oldq4k
ID MSC, D200X6 $
TIME 10
SOL 200 $ OPTIMIZATION
CEND
TITLE = VIBRATION OF A BEAM. D200X6
ECHO = UNSORT
DISP = ALL
STRESS = ALL
METHOD = 1
ANALYSIS = MODES
DESOBJ(MIN) = 1 $ OBJECTIVE FUNCTION DEFINITION
DESSUB = 10 $ CONSTRAINT SET SELECTION
$
BEGIN BULK
$
$------------------------------------------------------------------------
$ ANALYSIS MODEL:
$------------------------------------------------------------------------
$
EIGRL, 1, , ,
GRID, 1, , 0.0, 0.0, -3.0, , 123456
GRID, 2, , 20.0, 0.0, -3.0, , 2456
GRID, 3, , 40.0, 0.0, -3.0, , 2456
GRID, 4, , 60.0, 0.0, -3.0, , 2456
GRID, 5, , 0.0, 0.0, 3.0, , 123456
GRID, 6, , 20.0, 0.0, 3.0, , 2456
GRID, 7, , 40.0, 0.0, 3.0, , 2456
GRID, 8, , 60.0, 0.0, 3.0, , 2456
$
CROD, 1, 201, 5, 6
CROD, 2, 202, 6, 7
CROD, 3, 203, 7, 8
CROD, 7, 201, 1, 2
CROD, 8, 202, 2, 3
CROD, 9, 203, 3, 4
PROD, 201, 1, 1.0, 0.0
PROD, 202, 1, 1.0, 0.0
PROD, 203, 1, 1.0, 0.0
$
CQUAD4, 4, 204, 1, 2, 6, 5
CQUAD4, 5, 205, 2, 3, 7, 6
CQUAD4, 6, 206, 3, 4, 8, 7
PSHELL, 204, 1, 0.2
PSHELL, 205, 1, 0.2
PSHELL, 206, 1, 0.2
$
4-16
BEAM VIBRATION FREQUENCY AS A DESIGN
CONSTRAINT (Cont.)
4-17
BEAM VIBRATION FREQUENCY AS A DESIGN
CONSTRAINT (Cont.)
4-18
BEAM VIBRATION FREQUENCY AS A DESIGN CONSTRAINT (Cont.)
***** NORMAL CONVERGENCE CRITERIA SATISFIED ***** (HARD CONVERGENCE DECISION LOGIC)
**************************************************************************************
CONVERGENCE ACHIEVED BASED ON THE FOLLOWING CRITERIA
(HARD CONVERGENCE DECISION LOGIC)
***************************************************************
S U M M A R Y O F D E S I G N C Y C L E H I S T O R Y
***************************************************************
----------------------------------------------------------------------------------------------------------------------------------
INTERNAL | EXTERNAL | |
DV. ID. | DV. ID. | LABEL | 6 : 7 : 8 : 9 : 10 : 11 :
4-21
----------------------------------------------------------------------------------------------------------------------------------
1 | 1 | A1 | 8.2998E-01 :
2 | 2 | A2 | 4.3973E-01 :
3 | 3 | A3 | 1.1371E-01 :
4 | 4 | T1 | 4.9550E-02 :
5 | 5 | T2 | 4.2818E-02 :
6 | 6 | T3 | 2.6718E-02 :
*** USER INFORMATION MESSAGE 6464 (DOM12E)
RUN TERMINATED DUE TO HARD CONVERGENCE TO AN OPTIMUM AT CYCLE NUMBER = 6.
BEAM VIBRATION FREQUENCY AS A DESIGN
CONSTRAINT RESULTS
4-22
BEAM VIBRATION FREQUENCY AS A DESIGN
CONSTRAINT RESULTS (Cont.)
4-23
BEAM VIBRATION FREQUENCY AS A DESIGN
CONSTRAINT RESULTS (Cont.)
SECTION 5
DLINK . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5-2
5-1
DLINK Multiple Design Variable Linking
Relates
DLINK one design variable to one or more other design variables. Multiple Design Variable Linking
Format:
1 2 3 4 5 6 7 8 9 10
DLINK ID DDVID C0 CMULT IDV1 C1 IDV2 C2
IDV3 C3 -etc.-
Example:
8 7.0
Field Contents
ID Unique identification number. (Integer > 0)
Remarks:
1. DLINK defines the relationship
2. This capability provides a means of linking physical design variables such as element
thicknesses to nonphysical design variables such as the coefficients of interpolating
functions.
3. CMULT provides a simple means of scaling the Ci. For example if Ci = 1/7, 2/7, 4/7, etc.
is desired, then CMULT = 1/7 and Ci = 1, 2, 4, etc., may be input.
4. An independent IDVi must not occur on the same DLINK entry more than once.
5-2
DESIGN VARIABLE LINKING (Cont.)
t1 t2 t3 t4
x
0 10. 20. 30. 40.
t4 – t1
t i = ---------------- x i + t 1
30
2 1
so, t 2 = --- t 1 + --- t 4
3 3
1 2
t 3 = --- t 1 + --- t 4
3 3
5-3
DESIGN VARIABLE LINKING (Cont.)
The only independent variables are t1 and t4. t2 and t3 are the
dependent quantities. The design space is now two-dimensional,
rather than four.
5-4
REDUCED BASIS FORMULATIONS
x=0
t1
t2
t10
x=L
5-5
REDUCED BASIS FORMULATIONS (Cont.)
Basis Functions
F(x/L)
F1
1.0
F2
F3
x/L
0.5 1.0
F 1 --- = C 1
x
L
F 2 --- = C 1 1 – ---
x x
L L
x 2
F 3 --- = C 1 1 – ---
x
L L
C 1 = 1.0
5-6
REDUCED BASIS FORMULATIONS (Cont.)
t1
1.0 1.0 1.0
t2
1.0 0.9 0.81
t
3 1.0 0.8 0.64
t
4 1.0 0.7 0.49
t x1
5 1.0 0.6 0.36
= x2
t6 1.0 0.5 0.25
x3
t7 1.0 0.4 0.16
t8 1.0 0.3 0.09
t 1.0 0.2 0.04
9
t 1.0 0.1 0.01
10
{ p } M = [ T ] MXN { x } N
where M » N
5-7
REDUCED BASIS FORMULATIONS (Cont.)
5-8
DESIGN RESPONSES AND CASE CONTROL
Exceptions:
• Dynamic Responses
Example:
SUBCASE 10
ANALYSIS = DFREQ
SET 200 = 1000,1001,1003
DISPLACEMENT(PHASE) = 200
5-9
DESIGN RESPONSES AND CASE CONTROL (Cont.)
• Stress/Strain
Example:
SUBCASE 20
ANALYSIS = STATICS
SET 200 = 1000,1001,1003
STRESS(MAXS) = 200
5-10
MODE TRACKING FEATURES
t
Φi Mi Φi – 1 = ti
$
$ Mode Tracking has been performed successfully. Updated
DRESP1 entries are:
$
DRESP1 301LFREQ FREQ 2
DRESP1 303HFREQ FREQ 4
5-11
MODTRAK Mode Tracking Request
Format:
MODTRAK = n
Example:
MODTRAK = 100
Describer Meaning
n Set identification of a MODTRAK Bulk Data entry. (Integer > 0)
Remark:
1. Selection of a MODTRAK Bulk Data entry with the MODTRAK Case Control command
activates mode tracking for the current subcase. This request is limited to normal modes
subcases (ANALYSIS = MODES) in design optimization (SOL 200).
5-12
Mode Tracking Parameters MODTRAK
MODTRAK parameters for mode tracking in design optimization (SOL 200).
Specifies Mode Tracking Parameters
Format:
1 2 3 4 5 6 7 8 9 10
MODTRAK SID LOWRNG HIGHRNG MTFILTER
Example:
Field Contents
SID Sets identification number that is selected in the Case Control Section with the
MODTRAK command. (Integer; No Default) See Remark 1.
HIGHRNG Highest mode number in range to search. See Remark 2. (Integer > 0,
Default = number of eigenvalues extracted. If nonzero, LOWRNG < HIGHRNG.)
Remarks:
1. Only the designed modes for the subcase will be tracked. A designed mode is one that is
used in the design model (in connection with either objective or constraints) and, there-
fore, identified on a DRESP1 entry.
2. The range of modes LOWRNG through HIGHRNG, inclusive, will be used to track the de-
signed modes. If LOWRNG and HIGHRNG are both blank, then all computed modes will
be used to search for the designed modes. Since large numbers of computed modes will
result in higher computational costs, limiting the search range with LOWRNG and HIGH-
RNG is recommended.
3. Modes are considered to correlate if their mass normalized cross orthogonalities are
greater than MTFILTER.
5-13
EQUALITY CONSTRAINTS
Equality constraints are not provided for directly, although they can
be formulated in MSC/NASTRAN.
The simplest way is to provide lower and upper bounds that are
equivalent:
r ≤ r(x) ≤ r
$
$ EIGENVALUE FOR MODE 2:
DRESP1, 15, LAMA1, EIGN, , , 2
$
$ BOUNDS:
DCONSTR,100, 15, 3.5531E4,3.5531E4 $ 30 Hz
$
5-14
FREQUENCY MATCHING
• Design objective:
• Approach 1
Λi – ε ≤ λi ≤ Λi + ε
• Approach 2
5-15
5-16
SECTION 6
DESIGN MODELING II
DVPREL2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6-3
DEQATN . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6-5
DTABLE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6-9
DRESP2. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6-13
6-1
USER-DEFINED DESIGN VARIABLE-TO-PROPERTY
RELATIONS
PL
Analysis = Po + C2 X DVPREL1
Model =
Parameters PNL
= f(X,C) DVPREL2
6-2
Design Variable to Property Relation DVPREL2
Defines the relation between an analysis model property and design Design
DVPREL2 variables withtoaProperty
Variable user-supplied
Relation
equation.
Format:
1 2 3 4 5 6 7 8 9 10
DVPREL2 ID TYPE PID FID PMIN PMAX EQID
Example:
DESVAR 4 11 13 5
DTABLE PI YM
Field Contents
ID Unique identification number. (Integer > 0)
FID Field position of the property in the analysis model entry. (Integer ≠ 0)
PMIN Minimum value allowed for this property. If FID references a stress recovery
location field, then the default value for PMIN is –1.0+35. PMIN must be
explicitly set to a negative number for properties that may be less than zero (for
example, field ZO on the PCOMP entry). (Real; Default = 1.E–20)
PMAX Maximum value allowed for this property. (Real; Default = 1.0E20)
“DESVAR” DESVAR flag. Indicates that the IDs of DESVAR entries follow. (Character)
“DTABLE” DTABLE flag. Indicates that the IDs for the constants in a DTABLE entry follow.
This field may be omitted if there are no constants involved in this relation.
(Character)
(Continued)
6-3
DVPREL2 Design Variable to Property Relation
Remarks:
1. The variables identified by DVIDi and LABLi correspond to variable names (x1, x2, etc.)
listed in the left-hand side of the first equation on the DEQATN entry identified by EQID.
The variable names x1 through xN (where N = m + n) are assigned in the order DVID1,
DVID2, ..., DVIDn, LABL1, LABL2, ..., LABLm.
2. If both “DESVAR” and “DTABLE” are specified in field 2, “DESVAR” must appear first.
3. FID may be either a positive or a negative number. If FID > 0, it identifies the field position
on a property entry. If FID < 0, it identifies the word position of an entry in EPT. For
example, to specify the area of a PBAR, either FID = +4 or FID = –3 may be used. How-
ever, if PTYPE = “PBEAM”, FID must be negative. See Remark 4 on the DVPREL1 entry
description for specification of the PBEAM element property tables.
4. PTYPE = “PBEND” is not supported, either directly through FIDs or indirectly via word po-
sitions in the element property table.
6-4
Design Equation Definition DEQATN
DEQATN one or more equations for use in design sensitivity or p-element analysis.
Defines Design Equation Definition
Format:
1 2 3 4 5 6 7 8 9 10
DEQATN EQID EQUATION
EQUATION (Cont.)
Example:
+ A ∗ ∗ 2/(B – C); F = A + B – F1 ∗ D
Field Contents
EQID Unique equation identification number. (Integer > 0)
Remarks:
1. EQUATION is a single equation or a set of nested equations and is specified in fields 3
through 9 on the first entry and may be continued on fields 2 through 9 on the continuation
entries. On the continuation entries, no commas can appear in columns 1 through 8. All
data in fields 2 through 9 must be specified in columns 9 through 72. The large-field for-
mat is not allowed.
(Continued)
6-5
DEQATN Design Equation Definition
2
3 A
1 = A + B ∗ C – ( D + 10 ) + sin ( PI ( 1 ) ∗ R ) + --------------
B–C
F = A + B – F1 ∗ D
2. EQUATION may contain embedded blanks. EQUATION must contain less than 12,500
nonblank characters. This is equivalent to approximately 195 continuation entries.
3. The syntax of the expressions follows FORTRAN language standards. The allowable
arithmetic operations are shown in Table 1 in the order of execution precedence. Paren-
thesis are used to change the order of precedence. Operations within parentheses are
performed first with the usual order of precedence being maintained within the
parentheses.
4. The expressions may contain intrinsic functions. Table 2 contains the format and descrip-
tions of functions that may appear in the expressions. The use of functions that may be
discontinuous must be used with caution because they can cause discontinuous deriva-
tives. These are ABS, DIM, MAX, MIN, and MOD. For examples and further details see
the MSC/NASTRAN DMAP Module Dictionary.
(Continued)
6-6
Design Equation Definition DEQATN
Table 2. DEQATN Entry Functions.
ACOS(x) arccosine -1
cos x
ASIN(x) arcsine -1
sin x
ATAN(x) arctangent -1
tan x
EXP(x) exponential x
e
6-7
DEQATN Design Equation Definition
b. DRESP2 entry, then xi represents the DVIDj, LABLk, NRm, and Gp fields in that
order.
7. The DMAP logical operators NOT, AND, OR < XOR, and XQV cannot be used as Xi
names.
8. Input errors on the DEQATN entry often result in poor messages. Substituting a “[” for a
parenthesis or violating the restriction against large field format are examples. Known
messages are UFM 215, SFM 233 and UFM 5199. If any of these messages are encoun-
tered then review the DEQATN entry input.
6-8
Table Constants DTABLE
DTABLE a table of real constants that are used in equations (see DEQATN entry).
Defines Table Constants
Format:
1 2 3 4 5 6 7 8 9 10
DTABLE LABL1 VALU1 LABL2 VALU2 LABL3 VALU3 LABL4 VALU4
Example:
G 5.5E5 B 100.
Field Contents
LABLi Label for the constant. (Character)
Remarks:
1. Only one DTABLE entry may be specified in the Bulk Data Section.
6-9
EXAMPLE – DEFINING A RECTANGULAR
SECTION BAR
b = 0.3 1
2 C 1 = – --- b
c h = 0.4 2
A = b⋅h 1
C 2 = --- h
2
h 1 3
bh
I 1 = ----------
12
3
hb
b I 2 = ----------
12
6-10
USER-DEFINED RESPONSES
DRESP1
• Design Variables
• Table Constants
• Type 1 Response(s)
• Grid Coordinates
• DVPREL1 Entries
6-11
USER-DEFINED RESPONSES (Cont.)
• Examples:
2 2 2
U = ux + uy + uz
6-12
Design Sensitivity Equation Response Quantities DRESP2
DRESP2 equation responses that are used in the design, either
Defines asSensitivity
Design constraints or asResponse
Equation an objective.
Quantities
Format:
1 2 3 4 5 6 7 8 9 10
DRESP2 ID LABEL EQID REGION
DVID8 -etc.-
LABL8 -etc.-
NR8 -etc.-
“DNODE” G1 C1 G2 C2 G3 C3
G4 C4 -etc.-
Example:
DRESP2 1 LBUCK 5 3
201
DTABLE PI YM L
DRESP1 14 1 4 22 6 33 2
DNODE 14 1 4 1 22 3
2 1 43 1
Field Contents
ID Unique identification number. (Integer > 0)
(Continued)
6-13
DRESP2 Design Sensitivity Equation Response Quantities
REGION Region identifier for constraint screening. See Remark 5. (Integer > 0)
“DTABLE” Flag indicating that the labels for the constants in a DTABLE entry follow.
(Character)
“DNODE” Flag signifying that the following fields are “designed grid points.” See Remark 7.
(Character)
Remarks:
1. DRESP2 entries may only reference DESVAR, DTABLE, DRESP1, and DVPREL1 en-
tries. They may not reference other DRESP2 entries.
3. DRESP2 entries must have unique identification numbers with respect to DRESP1
entries.
4. The “DESVAR”, “DTABLE”, “DRESP1”, “DNODE”, and “DVPREL1” flags in field 2 must
appear in the order given above. Any of these words, along with the identification num-
bers associated with them, may be omitted if they are not involved in this DRESP2
relationship. However, at least one of these four types of arguments must exist.
5. The REGION field follows the same rules as for the DRESP1 entries. DRESP1 and
DRESP2 responses will never be contained in the same region, even if they are assigned
the same REGION identification number. The default is to put all responses referenced
by one DRESP2 entry in the same region.
(Continued)
6-14
Design Sensitivity Equation Response Quantities DRESP2
6. The variables identified by DVIDi, LABLj, NRk, the Gm, Cm pairs and DPIPi are assigned
(in that order) to the variable names (x1, x2, x3, etc.) specified in the left-hand side of the
first equation on the DEQATN entry referenced by EQID. The variable names x1 through
xN (N = m + n + p + q) are assigned in the order DVID1, ..., DVIDm, LABL1, ..., LABLn,
NR1, ..., NRp, G1, ..., Gq, DPIP1, ..., DPIPr. In the example below,
DRESP2 1 LBUCK 5 3
DESVAR 101 3
DTABLE PI YM
DNODE 14 1
7. (Gm, Cm) refer to a designed grid component. Depending on the scheme used in gener-
ating basis vectors, a designed grid component can be one of the following:
b. A grid component that is free to move in a full auxiliary model when the external dis-
placement fields are used.
c. When the auxiliary model boundary shape method is used, a designed grid compo-
nent can be either a loaded grid component in the auxiliary model or a grid
component with motion that is obtained from interpolation of boundary shapes.
d. When the interface using geometric boundary shapes is used, a designed grid com-
ponent is either a grid component defined on a DVGRID entry or one with motion that
is obtained from interpolation of boundary shapes.
6-15
EXAMPLE – DEFORMATION CONTROL WITH
TYPE-2 RESPONSES
desired shape
deformed shape
u1 u2 u3 u4
1 2 3 4
original shape
Design objective:
2u 1 + u 4
ũ 2 = -----------------------
3
and
u 1 + 2u 4
ũ 3 = -----------------------
3
6-16
EXAMPLE – DEFORMATION CONTROL WITH
TYPE-2 RESPONSES (Cont.)
2 2
DEF = ( u 2 – ũ 2 ) + ( u 3 – ũ 3 )
$ IN CASE CONTROL:
$
DESOBJ(MIN) = 21
$
$ IN BULK DATA:
$
$...IDENTIFY THE FIRST-LEVEL DISPLACEMENT RESPONSES:
$
$DRESP1,ID, LABEL, RTYPE, PTYPE, REGION, ATTA, ATTB, ATT1, +
$+, ATT2, ...
DRESP1, 11, UY1, DISP, , , 2, , 1
DRESP1, 12, UY2, DISP, , , 2, , 2
DRESP1, 13, UY3, DISP, , , 2, , 3
DRESP1, 14, UY4, DISP, , , 2, , 4
$
$DRESP2,ID, LABEL, EQID, REGION, , , , , +
$+, DRESP1, NR1, NR2, ...
DRESP2, 21, DEF2, 100, , , , , , +
+, DRESP1, 11, 12, 13, 14
$
DEQATN 100 F1(A,B,C,D) = (2.*A + D)/3.0 ; +
+ F2 = (A + 2.*D)/3.0 ; +
+ F = SQRT( (B-F1)**2 + (C-F2)**2 )
$
$.......2.......3.......4.......5.......6.......7.......8.......9.......0
6-17
DRESP2 DEPENDENCY ON DRESP1 RESPONSES
EID1
EID2
PID1 .
.
.
PID2
EIDn
DRESP2 DRESP1 .
.
.
PIDm
6-18
EXAMPLE – BAR AXIAL PLUS BENDING STRESSES
Stress recovery locations for a BAR element yield only bending stress
components. DRESP2 entries can be used to include axial stress
component, yielding total stress.
$ ANALYSIS:
$
CBAR, 10, 100, 1, 2, 1., 1., 0.
CBAR, 11, 100, 2, 3, 1., 1., 0.
CBAR, 12, 100, 3, 4, 1., 1., 0.
CBAR, 13, 100, 4, 5, 1., 1., 0.
$
PBAR, 100, 101, 1.25E3, 2.6042E5,6.5104E4,, , , +
+, 12.5, 25., 12.5, -25., -12.5, -25., -12.5, 25.
$
$ DESIGN:
$
DRESP1, 110, SA1, STRESS, PBAR, , 2, , 100
DRESP1, 111, SA2, STRESS, PBAR, , 3, , 100
DRESP1, 112, SA3, STRESS, PBAR, , 4, , 100
DRESP1, 113, SA4, STRESS, PBAR, , 5, , 100
DRESP1, 114, AXIAL, STRESS, PBAR, , 6, , 100
$
DRESP2, 210, SR1, 303, , , , , , +
+, DRESP1, 110, 114
DRESP2, 211, SR2, 303, , , , , , +
+, DRESP1, 111, 114
DRESP2, 212, SR3, 303, , , , , , +
+, DRESP1, 112, 114
DRESP2, 213, SR4, 303, , , , , , +
+, DRESP1, 113, 114
$
DEQATN 303 SRI(POINT,AXIAL) = POINT+AXIAL
6-19
RESTRICTIONS IN FORMING SYNTHETIC
RESPONSES
6-20
SECTION 7
EXAMPLES II
Fz
10404 20204
Web Cap
Fx
10204
X
10203
Y
20200
10004
10400
10200
Web
Plate
10000
Load cases:
7-1
STIFFENED PANEL WEIGHT MINIMIZATION (Cont.)
• Objective
• Weight minimization
• Constraints
• Stress:
• Displacement:
7-2
STIFFENED PANEL WEIGHT MINIMIZATION (Cont.)
• Design variables
• Thickness of web
• Thickness of plate
A 2 A 2
I 1 = I 1 ------- , I 2 = I 2 -------
0 A 0 0 A 0
1
---
2 2
I 12 = I 12 ------- , C j = C j -------
A A
0 A 0 0 A 0
7-3
STIFFENED PANEL WEIGHT MINIMIZATION (Cont.)
ID MSC, D200X7 $
TIME 10
SOL 200 $ OPTIMIZATION
CEND
$
TITLE = STATIC ANALYSIS OF A STIFFENED PLATE D200X7
ECHO = UNSORT
DISP = ALL
STRESS = ALL
SPC = 1
ANALYSIS = STATICS
DESOBJ(MIN) = 15 $ OBJECTIVE FUNCTION DEFINITION
$ (MIN IS THE DEFAULT)
SUBCASE 1
LABEL = LOAD CONDITION 1
LOAD = 1
DESSUB = 100 $ CONSTRAINT DEFININITION
SUBCASE 2
LABEL = LOAD CONDITION 2
LOAD = 2
DESSUB = 200 $ CONSTRAINT DEFININITION
BEGIN BULK
$
$------------------------------------------------------------------------
$ ANALYSIS MODEL:
$------------------------------------------------------------------------
$
GRID, 10000, , 0.0, 0.0, 0.0
GRID, 10001, , 2.5, 0.0, 0.0
GRID, 10002, , 5.0, 0.0, 0.0
GRID, 10003, , 7.5, 0.0, 0.0
GRID, 10004, , 10.0, 0.0, 0.0
GRID, 10100, , 0.0, 2.5, 0.0
GRID, 10101, , 2.5, 2.5, 0.0
GRID, 10102, , 5.0, 2.5, 0.0
GRID, 10103, , 7.5, 2.5, 0.0
GRID, 10104, , 10.0, 2.5, 0.0
GRID, 10200, , 0.0, 5.0, 0.0
GRID, 10201, , 2.5, 5.0, 0.0
GRID, 10202, , 5.0, 5.0, 0.0
GRID, 10203, , 7.5, 5.0, 0.0
GRID, 10204, , 10.0, 5.0, 0.0
GRID, 10300, , 0.0, 7.5, 0.0
GRID, 10301, , 2.5, 7.5, 0.0
GRID, 10302, , 5.0, 7.5, 0.0
GRID, 10303, , 7.5, 7.5, 0.0
GRID, 10304, , 10.0, 7.5, 0.0
GRID, 10400, , 0.0, 10.0, 0.0
GRID, 10401, , 2.5, 10.0, 0.0
7-4
STIFFENED PANEL WEIGHT MINIMIZATION (Cont.)
7-5
STIFFENED PANEL WEIGHT MINIMIZATION (Cont.)
$
MAT1, 1, 1.0E+7, , 0.33, 0.283
$
FORCE, 1, 10004, , 2000.0, 1.0, 0.0, 0.0
FORCE, 1, 10104, , 2000.0, 1.0, 0.0, 0.0
FORCE, 1, 10204, , 2000.0, 1.0, 0.0, 0.0
FORCE, 1, 10304, , 2000.0, 1.0, 0.0, 0.0
FORCE, 1, 10404, , 2000.0, 1.0, 0.0, 0.0
FORCE, 2, 10203, , 10000.0,0.0, 0.0, 1.0
PLOAD2, 1, 50., 1, THRU, 16
$
SPC1, 1, 1236, 10000
SPC1, 1, 136, 10100, 10300, 10400
SPC1, 1, 36, 10001, 10002, 10003, 10004, 10104
SPC1, 1, 36, 10401, 10402, 10403, 10404, 10304
SPC1, 1, 3, 10204
SPC1, 1, 13, 10200
SPC1, 1, 5, 20100, 20101, 20102, 20103, 20104
SPC1, 1, 6, 10101, 10102, 10103, 10104
SPC1, 1, 6, 10301, 10302, 10303, 10304
$
PARAM, GRDPNT, 1
PARAM, WTMASS, 0.00259
PARAM, AUTOSPC, YES
$
$------------------------------------------------------------------------
$ DESIGN MODEL:
$------------------------------------------------------------------------
$
$...Define the design variables:
$
$DESVAR,ID, LABEL, XINIT, XLB, XUB, DELXV
DESVAR, 1, T-PLATE,0.15, 0.001, 10.0
DESVAR, 2, T-WEB, 0.20, 0.001, 10.0
DESVAR, 3, A-BAR, 0.144, 0.001, 10.0
$
$...Relate the design variables to analysis model properties
$ (linear relations, so use DVPREL1)
$
$...Express shell thicknesses as functions of x1, x2:
$DVPREL1,ID, TYPE, PID, FID, PMIN, PMAX, C0, , +
$+, DVIDD1, COEF1, DVID2, COEF2, ...
DVPREL1,1, PSHELL, 1, 4, 0.01, , , , +DP1
+DP1, 1, 1.0
DVPREL1,2, PSHELL, 2, 4, 0.01, , , , +DP2
+DP2, 2, 1.0
$
$...Express bar cross sectional area as a function of x3:
DVPREL1,3, PBAR, 3, 4, 0.01, , , , +DP3
7-6
STIFFENED PANEL WEIGHT MINIMIZATION (Cont.)
+DP3, 3, 1.0
$
$...Proportionally relate bar’s I1, I2, I12 to changes in area (x3):
$ (nonlinear relations require use of DVPREL2+DEQATN)
$DVPREL2,ID, TYPE, PID, FID, PMIN, PMAX, EQID, , +
$+, DESVAR,DVID1, DVID2, ..., , , , , +
$+, DTABLE,CID1, CID2, ...
DVPREL2,11, PBAR, 3, 5, 1.0-6, , 101, , +
+, DESVAR, 3, , , , , , , +
+, DTABLE, X3INIT, I1INIT
DEQATN 101 I1(X3,X3INIT,I1INIT) = I1INIT ;
DELTA = SQRT(X3/X3INIT) ;
I1NEW = I1*DELTA**4
DVPREL2,12, PBAR, 3, 6, 1.0-6, , 102, , +
+, DESVAR, 3, , , , , , , +
+, DTABLE, X3INIT, I2INIT
DEQATN 102 I2(X3,X3INIT,I2INIT) = I2INIT ;
DELTA = SQRT(X3/X3INIT) ;
I2NEW = I2*DELTA**4
DVPREL2,13, PBAR, 3, 7, 1.0-6, , 103, , +
+, DESVAR, 3, , , , , , , +
+, DTABLE, X3INIT, I12INIT
DEQATN 103 I12(X3,X3INIT,I12INIT) = I12INIT ;
DELTA = SQRT(X3/X3INIT) ;
I12NEW = I12*DELTA**4
$
$...Modify stress recovery points accordingly:
DVPREL2,14, PBAR, 3, 12, , , 104, , +
+, DESVAR, 3, , , , , , , +
+, DTABLE, X3INIT, CYINIT
DVPREL2,15, PBAR, 3, 13, , , 104, , +
+, DESVAR, 3, , , , , , , +
+, DTABLE, X3INIT, CZINIT
DVPREL2,16, PBAR, 3, 14, , , 104, , +
+, DESVAR, 3, , , , , , , +
+, DTABLE, X3INIT, DYINIT
DVPREL2,17, PBAR, 3, 15, -1.0, , 104, , +
+, DESVAR, 3, , , , , , , +
+, DTABLE, X3INIT, DZINIT
DVPREL2,18, PBAR, 3, 16, -1.0, , 104, , +
+, DESVAR, 3, , , , , , , +
+, DTABLE, X3INIT, EYINIT
DVPREL2,19, PBAR, 3, 17, -1.0, , 104, , +
+, DESVAR, 3, , , , , , , +
+, DTABLE, X3INIT, EZINIT
DVPREL2,20, PBAR, 3, 18, -1.0, , 104, , +
+, DESVAR, 3, , , , , , , +
+, DTABLE, X3INIT, FYINIT
DVPREL2,21, PBAR, 3, 19, , , 104, , +
7-7
STIFFENED PANEL WEIGHT MINIMIZATION (Cont.)
+, DESVAR, 3, , , , , , , +
+, DTABLE, X3INIT, FZINIT
$
$...Equation for stress recovery points:
DEQATN 104 NEWPOINT(X3,X3INIT,POINT) = POINT*SQRT(X3/X3INIT)
$
$...Table constants for all equations:
DTABLE, X3INIT, 0.144, I1INIT, 1.728-4,I2INIT, 1.728-2,I12INIT,1.745-2,+
+, CYINIT, 0.06, CZINIT, 0.6, DYINIT, 0.06, DZINIT, -0.6, +
+, EYINIT, -0.06, EZINIT, -0.6, FYINIT, -0.06, FZINIT, 0.6
$
$...Identify the design responses:
$
$DRESP1,ID, LABEL, RTYPE, PTYPE, REGION, ATTA, ATTB, ATT1, +
$+, ATT2, ...
DRESP1, 1, SBARA, STRESS, PBAR, , 7, , 3
DRESP1, 2, SBARB, STRESS, PBAR, , 14, , 3
DRESP1, 3, S13, STRESS, PSHELL, , 9, , 1
DRESP1, 6, S16, STRESS, PSHELL, , 17, , 1
DRESP1, 9, S23, STRESS, PSHELL, , 9, , 2
DRESP1, 12, S26, STRESS, PSHELL, , 17, , 2
DRESP1, 13, D1, DISP, , , 3, , 10302
DRESP1, 14, D2, DISP, , , 3, , 10203
DRESP1, 15, W, WEIGHT
$
$...Place bounds on the responses:
$
$DCONSTR,DCID, RID, LALLOW, UALLOW
DCONSTR, 10, 1, -25000.,25000.
DCONSTR, 10, 2, -25000.,25000.
DCONSTR, 10, 3, -25000.,25000.
DCONSTR, 10, 6, -25000.,25000.
DCONSTR, 10, 9, -25000.,25000.
DCONSTR, 10, 12, -25000.,25000.
DCONSTR, 20, 13, -0.1, 0.1
DCONSTR, 30, 14, -0.03, 0.03
$
$DCONADD,DCID, DC1, DC2, ...
DCONADD, 100, 10, 20 $ summed constraint set for subcase 1
DCONADD, 200, 10, 30 $ summed constraint set for subcase 2
$
$...Optional override of optimization control parameters:
$
DOPTPRM,IPRINT, 1, DESMAX, 20, DELP, 0.5, P1, 1, +
+, P2, 15
$ (DELP=0.5 allows larger moves, thus overcoming constraint
$ violations quicker)
$
ENDDATA
7-8
STIFFENED PANEL WEIGHT MINIMIZATION (Cont.)
***** NORMAL CONVERGENCE CRITERIA SATISFIED ***** (HARD CONVERGENCE DECISION LOGIC)
**************************************************************************************
CONVERGENCE ACHIEVED BASED ON THE FOLLOWING CRITERIA
(HARD CONVERGENCE DECISION LOGIC)
RELATIVE CHANGE IN OBJECTIVE 0.0000E+00 MUST BE LESS THAN 1.0000E-03
OR ABSOLUTE CHANGE IN OBJECTIVE 0.0000E+00 MUST BE LESS THAN 1.0000E-02
--- AND ---
MAXIMUM CONSTRAINT VALUE 4.1453E-04 MUST BE LESS THAN 5.0000E-03
(CONVERGENCE TO A FEASIBLE DESIGN)
--- OR ---
MAXIMUM OF RELATIVE PROP. CHANGES 0.0000E+00 MUST BE LESS THAN 1.0000E-03
AND MAXIMUM OF RELATIVE D.V. CHANGES 0.0000E+00 MUST BE LESS THAN 1.0000E-03
(CONVERGENCE TO A BEST COMPROMISE INFEASIBLE DESIGN)
**************************************************************************************
7-9
STIFFENED PANEL WEIGHT MINIMIZATION (Cont.)
***************************************************************
S U M M A R Y O F D E S I G N C Y C L E H I S T O R Y
***************************************************************
(HARD CONVERGENCE ACHIEVED)
---------------------------------------------------------------------------------------------------------------
OBJECTIVE FROM OBJECTIVE FROM FRACTIONAL ERROR MAXIMUM VALUE
7-10
STATIC ANALYSIS OF A STIFFENED PLATE D200X7 JULY 26, 1994 MSC/NASTRAN 7/25/94 PAGE 92 DESIGN VARIABLE HISTORY
------------------------------------------------------------------------------------------------------------------------
INTERNAL | EXTERNAL | |
DV. ID. | DV. ID. | LABEL | INITIAL : 1 : 2 : 3 : 4 : 5 :
-------------------------------------------------------------------------------------------------------------------------
1 | 1 | T-PLATE | 1.5000E-01 : 3.0000E-01 : 4.4321E-01 : 2.8106E-01 : 1.4038E-01 : 1.0879E-01 :
2 | 2 | T-WEB | 2.0000E-01 : 4.0000E-01 : 6.1265E-01 : 8.4962E-01 : 8.2809E-01 : 7.5169E-01 :
3 | 3 | A-BAR | 1.4400E-01 : 1.6592E-01 : 1.2682E-01 : 6.4814E-02 : 3.2407E-02 : 1.6165E-02 :
-------------------------------------------------------------------------------------------------------------------------
INTERNAL | EXTERNAL | |
DV. ID. | DV. ID. | LABEL | 6 : 7 : 8 : 9 : 10 : 11 :
-------------------------------------------------------------------------------------------------------------------------
1 | 1 | T-PLATE | 1.1169E-01 : 1.1169E-01 :
2 | 2 | T-WEB | 8.3887E-01 : 8.3887E-01 :
3 | 3 | A-BAR | 1.6165E-02 : 1.6165E-02 :
-------------------------------------------------------------------------------------------------------------------------
*** USER INFORMATION MESSAGE 6464 (DOM12E)
7-11
6 3
4
5
10
7 Y
8
9 X
7-12
25-BAR TRUSS (Cont.)
• Three-dimensional
• Symmetric with respect to the x-z plane and the y-z plane
• Material
• E = 1.0E7 psi
3
• Weight density = 0.1 lbs/in
• Design variables
7-13
25-BAR TRUSS (Cont.)
• Constraints
σ cr ≤ F s ∗σ
Rearranging,
F s ∗σ
-------------- ≤ 1
σ cr
P cr 2
π EI
σ cr = --------- ; P cr = – ------------
A 2
L
7-14
25-BAR TRUSS (Cont.)
or,
2
π EI
σ cr = – ------------
2
AL
With:
A = 2πR∗t
π 2
= ------D
10
3
I = πR ∗t
π 4
= ------D
80
(2)
Our synthetic responses r are:
2
(2) σ 8L F s σ
r = --------- = – --------------------- ≤ 1
σ cr ( πD ) E
2
7-15
25-BAR TRUSS (Cont.)
ID MSC, D200X3 $
TIME 10
SOL 200 $ OPTIMIZATION
CEND
TITLE = OPTIMAL SIZING OF A 25-BAR TRUSS - D200X3
SUBTITLE = EIGHT INDEPENDENT CROSS SECTIONAL AREAS,
SUBTITLE = PLUS EULER BUCKLING
ECHO = UNSORT
OLOAD = ALL
DISP = ALL
SPCFORCE = ALL
ELFORCE = ALL
STRESS = ALL
SPC = 100
ANALYSIS = STATICS $
DESOBJ(MIN) = 11 $ OBJECTIVE FUNCTION DEFINITION
DESSUB = 12 $ CONSTRAINT DEFININITION
SUBCASE 1
LABEL = LOAD CONDITION 1
LOAD = 300
SUBCASE 2
LABEL = LOAD CONDITION 2
LOAD = 310
BEGIN BULK
$
$-----------------------------------------------------------------------------
$ ANALYSIS MODEL
$-----------------------------------------------------------------------------
$
GRDSET, , , , , , , 456
MAT1, 1, 10.0E6, , , 0.1
SPC1, 100, 123, 7, THRU, 10
GRID, 1 , , -37.5, 0.0, 200.0
GRID, 2 , , 37.5, 0.0, 200.0
GRID, 3 , , -37.5, 37.5, 100.0
GRID, 4 , , 37.5, 37.5, 100.0
GRID, 5 , , 37.5, -37.5, 100.0
GRID, 6 , , -37.5, -37.5, 100.0
GRID, 7 , , -100.0, 100.0, 0.0
GRID, 8 , , 100.0, 100.0, 0.0
GRID, 9 , , 100.0, -100.0, 0.0
GRID, 10, , -100.0, -100.0, 0.0
CROD, 1 , 1, 1, 2
CROD, 2 , 2, 1, 4
CROD, 3 , 2, 2, 3
CROD, 4 , 2, 1, 5
CROD, 5 , 2, 2, 6
CROD, 6 , 3, 2, 4
7-16
25-BAR TRUSS (Cont.)
CROD, 7 , 3, 2, 5
CROD, 8 , 3, 1, 3
CROD, 9 , 3, 1, 6
CROD, 10, 4, 3, 6
CROD, 11, 4, 4, 5
CROD, 12, 5, 3, 4
CROD, 13, 5, 5, 6
CROD, 14, 6, 3, 10
CROD, 15, 6, 6, 7
CROD, 16, 6, 4, 9
CROD, 17, 6, 5, 8
CROD, 18, 7, 4, 7
CROD, 19, 7, 3, 8
CROD, 20, 7, 5, 10
CROD, 21, 7, 6, 9
CROD, 22, 8, 6, 10
CROD, 23, 8, 3, 7
CROD, 24, 8, 4, 8
CROD, 25, 8, 5, 9
$
PROD, 1, 1, 2.0, 0.0
PROD, 2, 1, 2.0, 0.0
PROD, 3, 1, 2.0, 0.0
PROD, 4, 1, 2.0, 0.0
PROD, 5, 1, 2.0, 0.0
PROD, 6, 1, 2.0, 0.0
PROD, 7, 1, 2.0, 0.0
PROD, 8, 1, 2.0, 0.0
$
FORCE, 300, 1, , 1.0, 1000., 10000.,-5000.
FORCE, 300, 2, , 1.0, 0., 10000.,-5000.
FORCE, 300, 3, , 1.0, 500., 0., 0.
FORCE, 300, 6, , 1.0, 500., 0., 0.
FORCE, 310, 1, , 1.0, 0., 20000.,-5000.
FORCE, 310, 2, , 1.0, 0., -20000.,-5000.
$
$-----------------------------------------------------------------------------
$ DESIGN MODEL
$-----------------------------------------------------------------------------
$
$...Define the design variables
$ average (mean) diameters of tubes:
$DESVAR,ID, LABEL, XINIT, XLB, XUB, DELXV
$
DESVAR, 1, D1, 2.5, 0.01, 100.0
DESVAR, 2, D2, 2.5, 0.01, 100.0
DESVAR, 3, D3, 2.5, 0.01, 100.0
DESVAR, 4, D4, 2.5, 0.01, 100.0
7-17
25-BAR TRUSS (Cont.)
7-18
25-BAR TRUSS (Cont.)
7-19
25-BAR TRUSS (Cont.)
7-20
25-BAR TRUSS (Cont.)
***** NORMAL CONVERGENCE CRITERIA SATISFIED ***** (HARD CONVERGENCE DECISION LOGIC)
**************************************************************************************
CONVERGENCE ACHIEVED BASED ON THE FOLLOWING CRITERIA
(HARD CONVERGENCE DECISION LOGIC)
7-24
25-BAR TRUSS (Cont.)
7-25
CANTILEVER PLATE (REDUCED BASIS
FORMULATION)
p1
Y
p2
t1 t2 t3 t4 t5 t6 t7 t8 10"
40"
• Design objective
• Weight minimization
7-26
CANTILEVER PLATE (REDUCED BASIS
FORMULATION) (Cont.)
• Design constraints
Load case 2: 60.44 psi uniform pressure acting on the upper surface
7-27
CANTILEVER PLATE (REDUCED BASIS
FORMULATION) (Cont.)
ID MSC, D200X5 $
TIME 10
SOL 200 $ OPTIMIZATION
CEND
TITLE = CANTILEVERED PLATE - D200X5 D200X5
SUBTITLE = REDUCED BASIS FORMULATION
SPC = 100
OLOAD = ALL
DISP = ALL
STRESS = ALL
DESOBJ(MIN) = 35 $ OBJECTIVE FUNCTION DEFINITION
SUBCASE 1
ANALYSIS = STATICS
LABEL = LOAD CONDITION 1
LOAD = 300
DESSUB = 10
SUBCASE 2
ANALYSIS = STATICS
LABEL = LOAD CONDITION 2
LOAD = 310
DESSUB = 10
BEGIN BULK
$
$------------------------------------------------------------------------
$ ANALYSIS MODEL:
$------------------------------------------------------------------------
$
MAT1, 51, 1.0E+7, , 0.33, 0.1, , , , +M2
+M2, 50000., 50000., 29000.
SPC1, 100, 123456, 1, 11, 21
GRID, 1 , , 0., -5., 0.
GRID, 2 , , 5., -5., 0.
GRID, 3 , , 10., -5., 0.
GRID, 4 , , 15., -5., 0.
GRID, 5 , , 20., -5., 0.
GRID, 6 , , 25., -5., 0.
GRID, 7 , , 30., -5., 0.
GRID, 8 , , 35., -5., 0.
GRID, 9 , , 40., -5., 0.
GRID, 11, , 0., 0., 0.
GRID, 12, , 5., 0., 0.
GRID, 13, , 10., 0., 0.
GRID, 14, , 15., 0., 0.
GRID, 15, , 20., 0., 0.
GRID, 16, , 25., 0., 0.
GRID, 17, , 30., 0., 0.
GRID, 18, , 35., 0., 0.
GRID, 19, , 40., 0., 0.
7-28
CANTILEVER PLATE (REDUCED BASIS
FORMULATION) (Cont.)
7-29
CANTILEVER PLATE (REDUCED BASIS
FORMULATION) (Cont.)
7-30
CANTILEVER PLATE (REDUCED BASIS
FORMULATION) (Cont.)
7-31
CANTILEVER PLATE (REDUCED BASIS
FORMULATION) (Cont.)
7-32
CANTILEVER PLATE (REDUCED BASIS FORMULATION) (Cont.)
***** NORMAL CONVERGENCE CRITERIA SATISFIED ***** (HARD CONVERGENCE DECISION LOGIC)
**************************************************************************************
CONVERGENCE ACHIEVED BASED ON THE FOLLOWING CRITERIA
(HARD CONVERGENCE DECISION LOGIC)
***************************************************************
S U M M A R Y O F D E S I G N C Y C L E H I S T O R Y
***************************************************************
---------------------------------------------------------------------------------------------------------------
OBJECTIVE FROM OBJECTIVE FROM FRACTIONAL ERROR MAXIMUM VALUE
CYCLE APPROXIMATE EXACT OF OF
7-34
• May specify cross section dimensions (e.g., height and width for a
rectangular section) rather than properties (e.g., I, A, etc.)
8-1
BEAM CROSS-SECTION LIBRAY OPTIMIZATION
• For BEAM element, you still have to write the DVPREL2 and
DEQATN entries to optimize section dimension.
8-2
COMPARISON BETWEEN DESIGNING PBAR AND PBARL
1
--- B, --- H
1
2 2
C D
1 H
F E
B
*************************************
* Input Data for Designing PBAR *
*************************************
$
PBAR 9 1 5.0 10.4167 0.4167
-0.5 2.5 0.5 2.5 0.5 -2.5 -0.5 -2.5
C
$ ID LABEL INIT MIN MAX
DESVAR 1 HEIGHT 5.0 0.01 10.
DESVAR 2 WIDTH 1.0 0.01 10.
$
$ Explict expression for cross-section area
$ ID TYPE PID FID PMIN PMAX EQID
DVPREL2 1 PBAR 9 4 0.01 100
DESVAR 1 2
DEQATN 100 A(H,B)=H*B
$
$ Explict expression for first moment intertia
DVPREL2 2 PBAR 9 5 0.01 200
DESVAR 1 2
DEQATN 200 A(H,B)=B*H**3/12.
$
$ Explict expression for stress data recovery point D (0.5B, 0.5H)
$ ID TYPE PID FID PMIN PMAX C0
DVPREL1 11 PBAR 9 14 0.01
1 0.5
DVPREL1 12 PBAR 9 15 0.01
2 0.5
8-3
COMPARISON BETWEEN DESIGNING PBAR AND
PBARL (Cont.)
*************************************
* Input Data for Designing PBARL *
*************************************
$
PBARL, 9, 1, , BAR
1.0 5.0
$ ID LABEL INIT MIN MAX
DESVAR 1 HEIGHT5.0 0.01 10.
DESVAR 2 WIDTH 1.0 0.01 10.
$
$ Dimension Height references to design variable 1
$ ID TYPE PID FID PMIN PMAX C0
DVPREL1 11 BARL 9 13 0.01
1 1
$
$ Dimension Width references to design variable 2
DVPREL1 12 PBARL 9 12 0.01
1 2
8-4
DVPREL1 INCLUDED ON DRESP2
Format:
1 2 3 4 5 6 7 8 9 10
DRESP2 ID LABEL EQID REGION
“DESVAR” DVID1 DVID2 DVID3 DVID4 DVID5 DVID6 DVID7
DVID8 -etc.-
“DTABLE” LABL1 LABL2 LABL3 LABL4 LABL5 LABL6 LABL7
LABL8 -etc.-
“DRESP1” NR1 NR2 NR3 NR4 NR5 NR6 NR7
NR8 -etc.-
“DNODE” G1 C1 G2 C2 G3 C3
G4 C4 -etc.-
“DVPREL1” DPIP1 DPIP2 DPIP3 DPIP4 DPIP5 DPIP6 DPIP7
NEW!
8-5
DESIGN VARIABLE LINKING
The DLINK entry is the only way to implement design variable linking,
when DVPREL2’s are involved.
2
2 A = b ⋅ h = 5.0E-5 m
3
bh 4
1 h I 1 = ---------- = 1.042E-10 m
12
3
hb 4
b I 2 = ---------- = 4.167E-10 m
12
8-6
DESIGN VARIABLE LINKING (Cont.)
8-7
DESIGN VARIABLE LINKING (Cont.)
8-8
CONSTRAINT EQUATION SCALING
• Examples:
2 2
π E 2 2 8L
σ c + ---------- ( D + T ) ≥ 0 → ------------------------------------ σ c ≥ – 1.0
2 2 2 2
8L π E(D + T )
Synthetic Response
0.4ET D
σ c + ---------------- ≥ 0 → ---------------- σ c ≥ – 1.0
D 0.4ET
Synthetic Response
8-9
CONSTRAINT EQUATION SCALING (Cont.)
σ1 + σ2
------------------- ≤ σ max
2
σ1 + σ2
------------------- – σ max ≤ 0
2
σ1 + σ2
------------------- – σ max – ε
2
g u = ------------------------------------------------------ ≤ 0
ε
∂g u 1
---------- = ------
∂σ 1 2ε
• Then
g
1
------
2ε
CTMIN
σ1
σ1
CT
Active
8-10
CONSTRAINT EQUATION SCALING (Cont.)
σ1 + σ2
------------------- ≤ σ max
2
• Then
σ1 + σ2
------------------- –σ
2 max
g u = ---------------------------------------------- ≤ 0
σ max
∂g u 1
---------- = -----------------
∂σ 1 σ max
1
---------------------
2 σ max
CTMIN
σ1
σ1
CT
Active
8-11
CONSTRAINT EQUATION SCALING (Cont.)
σ1 + σ2
------------------- ≤ 1
2σ max
8-12
BEAM DESIGN – GENERAL ISSUES
End A End B
Words Words
6-21 166-181
8-13
BEAM DESIGN – GENERAL ISSUES (Cont.)
8-14
BEAM DESIGN – GENERAL ISSUES (Cont.)
8-15
BEAM DESIGN – SIMPLE EXAMPLE
b1
h1
h2
B
b2
b
1 1.0
h
1 2.0
{x} = =
b2 0.5
h2 1.0
8-16
BEAM DESIGN – SIMPLE EXAMPLE (Cont.)
$>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
$ BEAM MODELING TEST, VERSION 68
$
$ TAPERED BEAM WITH RECTANGULAR SECTION, DESIGN SENSITIVITY ANALYSIS. SINCE
$ BOTH ENDS A AND B ARE SPECIFIED ON THE PBEAM ENTRY, DESIGN CHANGES MUST
$ BE SPECIFIED NOT ONLY FOR THE ENDS, BUT ALSO FOR THE FIRST INTERMEDIATE
$ STATION, WHICH CONTAINS A COPY OF END B DATA.
$
$
$ PLANE 2
$ |
$ ----*----
$ | |
$ H | | --- PLANE 1
$ | |
$ ----*----
$ B * = STRESS RECOVERY LOCATIONS
$
$
$>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
$
TIME 10
SOL 200
CEND
$
$ CC FOR ANALYSIS:
DISP = ALL
STRESS = ALL
SPC = 100
LOAD = 300
$ CC FOR OPTIMIZATION:
ANALYSIS = STATICS
DESOBJ(MIN) = 8 $ OBJECTIVE FUNCTION DEFINITION
DESSUB = 10 $ CONSTRAINT SET SELECTION
$
BEGIN BULK
$PARAM, OPTIM, NO
PARAM, OPTEXIT, 4
$-------------------------------------------------------------------------------
$ ANALYSIS MODEL
$-------------------------------------------------------------------------------
MAT1 110 10.0E6 0.33 0.1 +M1
+M1 50000. 50000. 29000.
GRDSET 4
GRID 1 0.0 0.0 0.0
GRID 2 20. 0.0 0.0
CBEAM 1 100 1 2 0.0 1.0 0.0
$
$ PBEAM ENTRY INPUT WITH SLIGHT ‘ERROR’ IN TERMS. THIS HELPS VALIDATE THE DE-
$ SIGN MODEL BECAUSE USER WARNING MESSAGE WILL BE ISSUED, CONFIRMING OVERRIDE.
8-17
BEAM DESIGN – SIMPLE EXAMPLE (Cont.)
$
PBEAM 100 110 2.01 .167 .667 +P11
+P11 0.0 1.01 0.0 -1.01 +P12
+P12 YES 1.0 0.51 .0104 .042 +P13
+P13 0.0 0.51 0.0 -0.51
$
SPC1 100 123456 1
FORCE 300 2 20000.0 0.0 0.0 -1.0
$-------------------------------------------------------------------------------
$ DESIGN MODEL
$-------------------------------------------------------------------------------
$
$...END A DATA: (A,I1,I2,C2,D2)
$
$DESVAR,ID, LABEL, XINIT, XLB, XUB, DELXV
DESVAR, 1, B1, 1.0, 0.1, 10.0
DESVAR, 2, H1, 2.0, 0.2, 20.0
$
$DVPREL2,ID, TYPE, PID, FID, PMIN, PMAX, EQID, , +
$+, DESVAR, DVID1, DVID2, ..., , , , , +
$+, DTABLE, CID1, CID2, ...
DVPREL2,1, PBEAM, 100, -8, , , 101, , +
+, DESVAR, 1, 2
DVPREL2,2, PBEAM, 100, -9, , , 102, , +
+, DESVAR, 1, 2
DVPREL2,3, PBEAM, 100, -10, , , 103, , +
+, DESVAR, 1, 2
DVPREL2,4, PBEAM, 100, -15, , , 104, , +
+, DESVAR, 2
DVPREL2,5, PBEAM, 100, -17, , , 105, , +
+, DESVAR, 2
$
$...END B DATA: (A,I1,I2,C2,D2)
$
$DESVAR,ID, LABEL, XINIT, XLB, XUB, DELXV
DESVAR, 3, B2, 0.5, 0.05, 10.0
DESVAR, 4, H2, 1.0, 0.1, 20.0
$
$DVPREL2,ID, TYPE, PID, FID, PMIN, PMAX, EQID, , +
$+, DESVAR, DVID1, DVID2, ..., , , , , +
$+, DTABLE, CID1, CID2, ...
DVPREL2,6, PBEAM, 100, -168, , , 101, , +
+, DESVAR, 3, 4
DVPREL2,7, PBEAM, 100, -169, , , 102, , +
+, DESVAR, 3, 4
DVPREL2,8, PBEAM, 100, -170, , , 103, , +
+, DESVAR, 3, 4
DVPREL2,9, PBEAM, 100, -175, , , 104, , +
+, DESVAR, 4
DVPREL2,10, PBEAM, 100, -177, , , 105, , +
+, DESVAR, 4
$
$...FIRST INTERMEDIATE STATION (COPY OF END B DATA): (A,I1,I2,C2,D2)
$
8-18
BEAM DESIGN – SIMPLE EXAMPLE (Cont.)
8-19
BEAM DESIGN – SIMPLE EXAMPLE (Cont.)
As proof that the override took place, one can refer to the comparison table between analysis and
design models:
----- COMPARISON BETWEEN INPUT PROPERTY VALUES FROM ANALYSIS AND DESIGN MODELS -----
---------------------------------------------------------------------------------------------------------------
PROPERTY PROPERTY FIELD ANALYSIS DESIGN LOWER UPPER DIFFERENCE
TYPE ID ID VALUE VALUE BOUND BOUND FLAG
---------------------------------------------------------------------------------------------------------------
PBEAM 100 -177 -5.100000E-01 -5.000000E-01 -1.000000E+35 1.000000E+20 WARNING
PBEAM 100 -175 5.100000E-01 5.000000E-01 -1.000000E+35 1.000000E+20 WARNING
PBEAM 100 -170 4.200000E-02 4.166667E-02 1.000000E-03 1.000000E+20 WARNING
PBEAM 100 -169 1.040000E-02 1.041667E-02 1.000000E-03 1.000000E+20 WARNING
PBEAM 100 -168 5.100000E-01 5.000000E-01 1.000000E-03 1.000000E+20 WARNING
8-20
1. IF FIELD ID IS LESS THAN ZERO, IT IDENTIFIES THE WORD POSITION OF AN ENTRY IN EPT.
2. IF FIELD ID IS GREATER THAN ZERO, IT IDENTIFIES THE FIELD POSITION ON A PROPERTY BULK DATA ENTRY.
3. THE DIFFERENCE FLAG IS USED TO CHARACTERIZE DIFFERENCES BETWEEN ANALYSIS AND DESIGN MODEL PROPERTIES:
IF THE FLAG IS NONE, THEN THERE IS NO SIGNIFICANT DIFFERENCE BETWEEN THE TWO VALUES.
IF THE FLAG IS WARNING, THEN THE USER IS ADVISED THAT DIFFERENCES EXIST.
IF THE FLAG IS FATAL, THEN THE DIFFERENCES ARE GREATER THAN 1.00000E+35 AND THE RUN WILL BE TERMINATED.
SECTION 9
DSCREEN . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9-18
DSAPRT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9-49
FUNDAMENTALS OF STRUCTURAL
OPTIMIZATION (Cont.)
9-1
STRUCTURAL OPTIMIZATION –
THE DIFFICULTIES (Cont.)
9-2
STRUCTURAL OPTIMIZATION – THE SOLUTION
Experienced Approximation
Engineers Concepts
9-3
STRUCTURAL OPTIMIZATION –
THE SOLUTION (Cont.)
• Repeat if necessary
9-4
STRUCTURAL OPTIMIZATION –
THE SOLUTION (Cont.)
Improved Designs
Finite
Numerical
Element
Optimization
Analysis
Program
Program
Approximate
Design
Model
9-5
STRUCTURAL OPTIMIZATION –
THE SOLUTION (Cont.)
Initial
Design
Improved
Design
Constraint Approximate
Optimizer
Screening Model
Finite Element
Analysis
9-6
DESIGN VARIABLE LINKING
360" 360"
5 3 1
1 2
7 10
5 6
360"
8 9
3 4
6 4 2
9-7
DESIGN VARIABLE LINKING (Cont.)
A1
1.0 0.0 0.0
A2
1.0 0.0 0.0
A
3 1.0 0.0 0.0
A
4 1.0 0.0 0.0
A x1
5 0.0 1.0 0.0
= x2
A6 0.0 1.0 0.0
x3
A7 0.0 1.0 1.0
A8 0.0 0.0 1.0
A 0.0 0.0 1.0
9
A 0.0 0.0 1.0
10
9-8
DESIGN VARIABLE LINKING (Cont.)
360" 360"
5 3 1
1 2
7 10
5 6
360"
8 9
3 4
6 4 2
1 1 1 1
• Proposal 1: A 1, A 2, A 3, …, A 10
2 2 2 2
• Proposal 2: A 1, A 2, A 3, …, A 10
3 3 3 3
• Proposal 3: A 1, A 2, A 3, …, A 10
9-9
DESIGN VARIABLE LINKING (Cont.)
A1 A1
1
A1
2
A1
3
A 1 2 3
2 A2 A2 A2
A3 1 2 3
A3 A3 A3
A4 A4
1
A4
2
A4
3
x
A 1 2 3 1
5 A5 A5 A5
= x2
A A
1
A6
2
A6
4
6 6
1 2 3 x3
A7 A7 A7 A7
1 2 4
A8 A8 A8 A8
A A9
1
A9
2
A9
5
9
1 2 6
A A 10 A 10 A 10
10
The reduced basis method used here will yield a design consisting of
ten unique cross-sectional areas. Contrast this to the previous
example where linking of the properties (member size linking) yielded
3 characteristic cross-sectional areas. Though neither approach is
generally capable of yielding the global optimum, the techniques
often provide a convenient way of expressing manufacturability, or
other process constraints.
9-10
DESIGN VARIABLE LINKING IN SHAPE REDESIGN
x1 x2
1 2
3 4
h1
h2
5 6
x1 = –x2
x5 = –x6
1
x 4 = – ------ [ x 6 ( h 1 – h 2 ) + x 2 h 2 ]
h1
x3 = –x4
T
{ x }I = [ x2 x6 ]
T
{ x }D = [ x1 x3 x4 x5 ]
9-11
DESIGN VARIABLE LINKING IN SHAPE
REDESIGN (Cont.)
9-12
DESIGN VARIABLE LINKING IN SHAPE
REDESIGN (Cont.)
Cubic Polynomial
3 Grid Points
1 Slope
9-13
DESIGN VARIABLE LINKING IN SHAPE
REDESIGN (Cont.)
F ( ξ ) = x 1 F 1 ( ξ ) + x 2 F 2 ( ξ ) + x 3 F 3 ( ξ )…
F1 ( ξ )
F1 ( ξ )
F2 ( ξ )
F3 ( ξ )
.
.
.
Basis Functions
9-14
GENERATION OF CONSTRAINTS IN MSC/NASTRAN
rj ( x ) (DRESP1, 2)
r jL ≤ r j ( x ) ≤ r jU
r jL – r j ( x )
g 2j – 1 ( x ) = ------------------------- ≤ 0
r jL
(automatic)
r j ( x ) – r jU
g 2j ( x ) = -------------------------- ≤ 0
r jU
9-15
CONSTRAINT SCREENING
• Purpose:
Step I - Deletion:
Gj(X)
x x x x x x x x x
TRS
–1.0 j
List of Constraints
9-16
CONSTRAINT SCREENING (Cont.)
Step II - Regionalization:
Assume NSTR = 2
G(X)
x x x x x x x x x
TRS
–1.0
9-17
DSCREEN Design Constraint Screening Data
Format:
1 2 3 4 5 6 7 8 9 10
DSCREEN RTYPE TRS NSTR
Example:
Field Contents
RTYPE Response type for which the screening criteria apply. See Table 1. (Character)
NSTR Maximum number of constraints to be retained per region per load case. See
Remark 3. (Integer > 0; Default = 20)
Remarks:
1. Displacement and displacement derivative constraints associated with one particular load
case are grouped by the specification of DRESP1 entries. From each group, a maximum
of NSTR constraints are retained per load case.
(Continued)
9-18
Design Constraint Screening Data DSCREEN
2. Stress-strain constraints are grouped by the property; i.e., all elements belonging to the
set of PIDs specified under ATTi on a DRESPi entry are regarded as belonging to the
same region. In superelement sensitivity analysis, if the property (PID) is defined in more
than one superelement, then separate regions are defined. A particular stress constraint
specification may be applied to many elements in a region generating many stress con-
straints, but only up to NSTR constraints per load case will be retained.
4. If a certain type of constraint exists but no corresponding DSCREEN entry is specified, all
the screening criteria used for this type of constraint will be furnished by the default
values.
5. Constraints can be retained only if they are greater than TRS. See the Remarks under
the DCONSTR entry for a definition of constraint value.
9-19
CONSTRAINT SCREENING (Cont.)
• Recommendations
• Constraints that switch from the active and inactive sets from one
design cycle to the next
9-20
FORMAL APPROXIMATIONS
9-21
FORMAL APPROXIMATIONS (Cont.)
2
∆x
f ( x + ∆x ) = f ( x ) + f′ ( x )∆x + f″ ( x ) ---------- + …
2!
n
∂f
∑ ------- (x – x )
o o
f(x) ≅ f(x ) +
∂x i x o i i
i=1
9-22
FORMAL APPROXIMATIONS (Cont.)
Yi = Yi ( xi )
n
∂f
∑ -------
- ( Y ( x ) – Y i ( x i ) )
o o
f(x) ≅ f(x ) +
∂Y i x o i i
i=1
1
now, let Y i = ----
xi
∂Y i 1
then -------- = – ------
∂x i xi
2
n
∂f o 2 1 1
∑ -------
o
and fR ( x ) = f ( x ) – ( x ) ---- – ------
∂x x o i x o
i i x
i=1 i
n o
∂f x i o
∑ -------
o
= f(x ) – ------ ( x – x i )
∂x x o x i
i i
i=1
9-23
FORMAL APPROXIMATIONS (Cont.)
A
P
Stress = P/A
DISPL = PL/AE
σ̃
A˚ A
9-24
FORMAL APPROXIMATIONS (Cont.)
1/A˚ 1/A
9-25
FORMAL APPROXIMATIONS (Cont.)
Minimize F(x)
subject to:
gj ( x ) ≤ 0 j = 1,m
l u
xi ≤ xi ≤ xi i = 1,n
n
∂F
∑ -------
- ( Y i ( x i ) – Y i ( x i ) )
o o
Minimize: F˜ ( x ) = F ( x ) +
∂Y i x o
i=1
subject to:
n
∂g j
∑
o o
g˜ j ( x ) = g j ( x ) + -------- Y i ( x i ) – Y i ( x i ) j = 1,m
∂Y i x o
i=1
Move limits:
l u
xi ≤ xi ≤ xi i = 1,n
Yi = 1
----- for reciprocal approximation
Xi
9-26
FORMAL APPROXIMATIONS (Cont.)
• Weight
• Volume
• Internal force
• Displacement
• Stress
• Strain
• Eigenvalue/Natural Frequency
9-27
APRCOD – SELECTION OF APPROXIMATE FORM
N
∂g
g̃ ( X ) = g ( X o ) + ∑ -------
∂x i
( x i – x io )
i=1 x = xo
9-28
APRCOD – SELECTION OF APPROXIMATE
FORM (Cont.)
if g D ( X ) – g R ( X ) ≥ 0, use g D ( X )
if g D ( X ) – g R ( X ) < 0, use g R ( X )
Thus:
N
g(X ) + ∂g ∂g
∑
------- ( x i – x io ) if x i ⋅ ------- ≥ 0
o ∂x i ∂x i
i=1 x = xo
x = xo
g̃ ( X ) =
N
∂g ∂g
∑ x io ---- – --------
2 1 1
g(X ) – ------- if x i ⋅ ------- < 0
o ∂x x x ∂x i
i=1 i x=x i io
x = xo
o
9-29
SENSITIVITY ANALYSIS
KU = P
∂K ∂U ∂P
------- U + K ------- = -------
∂x i ∂x i ∂x i
∂U
• Solve for -------
∂x i
∂U ∂P ∂K
K ------- = ------- – ------- ⋅ U
∂x i ∂x i ∂x i
dr j ∆r j r j ( X + ∆x j , U + ∆u ) – r j ( X, U )
-------- ≅ -------- = ---------------------------------------------------------------------------
-
dx i ∆x i ∆x i
∂U
where ∆U = ------- ⋅ ∆x i
∂x i
9-30
SEMIANALYTIC SENSITIVITY ANALYSIS
∂P ∂K
• The difficulty in practice arises from the fact that ------- and ------- in the
∂x i ∂x i
pseudo-load vector are not generally available analytically. Thus,
we must introduce approximations for these derivatives. The
result is a semianalytical approximation for the displacement
sensitivity components
e o
e K ( x + ∆x ) – K ( x ) e o
∂K i
---------- ≅ ------------------------------------------------------------
∂x i ∆x i
• Central differences:
e o
e K ( x + ∆x ) – K ( x – ∆x ) e o
∂K i i
---------- ≅ ---------------------------------------------------------------------------
∂x i 2∆x i
2
the error is reduced to the order of ∆x . Central Differences are
the default for shape optimization and forward differences are the
default for property optimization. PARAM,CDIF,YES or NO can be
used to override the defaults.
9-31
SEMIANALYTIC SENSITIVITY ANALYSIS (Cont.)
• Eigenvalue sensitivity:
[ K – λ j M ]φ j = 0
T ∂K ∂M
φ j ------- – λ j ------- φ j
∂λ ∂x i ∂x i
------ = ---------------------------------------------
∂x i T
φ j Mφ j
• Buckling sensitivity
Kφ k + λ k K d φ k = 0
T ∂K ∂K d
– φ k ------- + λ k ---------- φ k
∂λ k ∂x i ∂x i
--------- = ------------------------------------------------------
∂x i T
( φk Kd φk )
9-32
DELB – FINITE DIFFERENCE STEP SIZE
∂k ∂p j
e e
∂k
--------- = --------- ⋅ --------
∂x i ∂p j ∂x i
e o o e o
∂k
e k ( p j + DELB ⋅ P j ) – k ( p j ) ∂p j
--------- ≅ ------------------------------------------------------------------------------ ⋅ --------
∂x i DELB ⋅ p j
o ∂x i x o
for DVPREL1,
pj = co + c1 x1 + … + ci xi + …
∂p j
-------- = c i
∂x i
9-33
DELB – FINITE DIFFERENCE STEP SIZE (Cont.)
for DVPREL2,
p j = f ( x 1, x 2, …, x n )
∂p j p j ( x 1, x 2, …, x i + ∆x i, …, x n ) – p j ( x 1, x 2, x i – ∆x i, …, x n )
-------- = -------------------------------------------------------------------------------------------------------------------------------------------------------
∂x i 2∆x i
• Displacement perturbation:
∂U j o
∆U j = -------- ⋅ DELB ⋅ p i
∂p i
• Stress gradient:
∂σ ∂σ ∂p j
------- = - ⋅ --------
-------
∂x i ∂p j ∂x i
o o o
σ ( u + ∆u, p j + DELB ⋅ p j ) – σ ( u, p j ) ∂p j
≅ --------------------------------------------------------------------------------------------------- ⋅ --------
DELB ⋅ p j
o ∂x i x o
o o
• If DELB ⋅ p j ≤ DELBM , then DELB ⋅ p j in the equations given
above will be replaced by DELBM = 1.0E–5.
9-34
SENSITIVITY ANALYSIS – IMPLEMENTATION
∂K ∂M
• In MSC/NASTRAN sensitivity analysis ------- and ------- are never formed
∂x i ∂x i
as G-size matrices
∂K ∂M
• Instead, vectors ------- U and ------- U are formed as the sums of the
∂x i ∂x i
products of the element level matrices and displacement vectors
M e
∂K ∂k m e
------- U =
∂x i ∑ ----------- ⋅ U
∂x i
m=1
M e
∂M ∂m m e
------- U =
∂x i ∑ ------------- ⋅ U
∂x i
m=1
9-35
SENSITIVITY ANALYSIS – IMPLEMENTATION (Cont.)
9-36
ADJOINT SENSITIVITY
Motivation
• Four subcases
• 100,000 DOF
9-37
ADJOINT SENSITIVITY THEORY
2
[ – ω M + iωB + K ] { u } = { P }
Grid response (r) is extracted from the global displacement vector (u)
r = f(u)
T
dr ⁄ dx = ∂f ⁄ ∂ u ∂u ⁄ ∂x
2 2
[ – ω M + iωB + K ] { du ⁄ dx } = – [ – ω dM ⁄ dx + iωdB ⁄ dx + dK ⁄ dx ] { u }
2 T
[ – ω M + iωB + K ] { λ } = { ∂ f ⁄ ∂ u }
Defining
2
[ FAC ] = [ – ω M + iωB + K ]
9-38
ADJOINT SENSITIVITY THEORY (Cont.)
Then
T 2
dr ⁄ dx = – { λ } [ – ω dM ⁄ dx + iωdB ⁄ dx + dK ⁄ dx ] { u }
–1
[ FAC ] [ FAC ] = [ I ]
9-39
ADJOINT SENSITIVITY
9-40
ADJOINT SENSITIVITY (Cont.)
Model
Design Task
nfreq
∑
2
φ = u si
i=1
Constraints
Weight ≤ 8.5
f 1 and f 2 ≤ 100 Hz
9-41
ADJOINT SENSITIVITY (Cont.)
Model Statistics
Number of grids: 2335
Number of elements: 2129
Number of subcases: 1
Number of excitation frequencies: 101
Number of design variables: 13
Number of responses: 91
Performance Results
9-42
ADJOINT SENSITIVITY (Cont.)
9-44
ADJOINT SENSITIVITY (Cont.)
Model Statistics
Number of grids: 2615
Number of elements: 2647
Number of subcases: 2
Number of design variables: 1075
Number of responses: 2
Performance Results
9-45
ADJOINT SENSITIVITY (Cont.)
Model Statistics
Number of grids: 14990
Number of elements: 16174
Number of subcases: 1
Number of frequencies: 1
Number of modes: 67
Number of design variables: 2514
Number of responses: 1
Performance Results
9-46
DESIGN SENSITIVITY COEFFICIENT MATRICES
∂r j Response
λ ij = -------
∂x i Design Variable
∆r j ∆r j ∆p k
--------- =
∆x i ∑ ---------- ⋅ ----------
∆p k ∆x i
k
Sensitivity DVPREL1,
Analysis DVPREL2
where pk = properties
9-47
DESIGN SENSITIVITY COEFFICIENT MATRICES (Cont.)
∆r j
---------
∆x i
r1 rj r1 rj
.
.
.
.
.
.
.
.
.
.
.
.
Subcase 1 Subcase 1 Subcase 2 Subcase 2
x1 ∆r 1j ∆r j
----------- ----------
∆x 1 1 ∆x 1 1
x2
. . . .
. . . .
.
.
.
. . . .
xn ∆r j
----------
∆x n k
9-48
Design Sensitivity Output Parameters DSAPRT
Specifies
DSAPRT design sensitivity output parameters. Design Sensitivity Output Parameters
Format:
FORMATTED ALL
NOEXPORT
DSAPRT ( UNFORMATTED , , [ START = i ], [ BY = j ], [ END = k ]) = n
EXPORT
NOPRINT NONE
Examples:
DSAPRT(FORMATTED,EXPORT)
DSAPRT(FORMATTED,START = FIRST,BY = 3,END = LAST) = 101
DSAPRT(UNFORMATTED,START = FIRST)
DSAPRT(UNFORMATTED,EXPORT)
DSAPRT(FORMATTED,END = 4) = ALL
DSAPRT(UNFORMATTED,END = SENS) = ALL
DSAPRT(NOPRINT,EXPORT)
Describers Meaning
FORMATTED Output will be presented with headings and labels.
UNFORMATTED Output will be printed as a matrix print (see description of the MATPRN module
in the MSC/NASTRAN DMAP Module Dictionary).
START = i Specifies the first design cycle for output. (Integer > 0 or Character: “FIRST” or
“LAST”; Default = 1 or “FIRST”)
BY = j Specifies the design cycle interval for output. (Integer ≤ 1, or > 0; Default = 0)
See Remark 2.
END = k Specifies the last design cycle for output. (Integer > 0 or Character: “FIRST”,
“LAST”, or “SENS”; Default = “LAST”)
ALL All design responses (defined in DRESP1 and DRESP2 entries) will be output.
(Continued)
9-49
DSAPRT Design Sensitivity Output Parameters
Remarks:
1. Only one DSAPRT may appear in the Case Control Section and should appear above all
SUBCASE commands.
2. Sensitivity data will be output at design cycles i, i + j, i + 2j, ..., k. Note that the BY = 0
default implies no sensitivity analysis at the intermediate design cycles.
3. END = SENS requests design sensitivity analysis, and no optimization will be performed.
4. Prior to Version 69, PARAM,OPTEXIT,4, –4, or 7 requested design sensitivity output and
may still be used in Version 69. However, if both DSAPRT and PARAM,OPTEXIT,4, –4,
or 7 are specified, then DSAPRT overrides PARAM,OPTEXIT,4, –4, or 7.
PARAM,OPTEXIT values and the equivalent DSAPRT commands are as follows:
–4 DSAPRT(NOPRINT,EXPORT,END = SENS)
7 DSAPRT(UNFORMATTED,START = LAST)
9-50
FORMATTED SENSITIVITY PRINT EXAMPLE
**********************************************************************
* *
* D E S I G N S E N S I T I V I T Y M A T R I X O U T P U T *
* *
* *
* R E S P O N S E S E N S I T I V I T Y C O E F F I C I E N T S *
* *
**********************************************************************
-----------------------------------------------------------------------------------------------------------------------
DRESP1 ID= 20 RESPONSE TYPE= WEIGHT SEID= 0
RESP VALUE DESIGN VARIABLE COEFFICIENT DESIGN VARIABLE COEFFICIENT
-----------------------------------------------------------------------------------------------------------------------
4.8284e+00 1 A1 2.8286e+00 2 A2 1.0006e+00
SYMMETRIC THREE BAR TRUSS DESIGN OPTIMIZATION - D200X1 JANUARY 18, 1996 MSC/NASTRAN 1/17/96 PAGE 24
9-51
-----------------------------------------------------------------------------------------------------------------------
DRESP1 ID= 23 RESPONSE TYPE= STRESS ELEM ID= 1 COMP NO= 2 SEID= 0
SUBCASE RESP VALUE DESIGN VARIABLE COEFFICIENT DESIGN VARIABLE COEFFICIENT
-----------------------------------------------------------------------------------------------------------------------
1 1.3530E+04 1 A1 -1.1895E+04 2 A2 -8.1543E+02
2 -9.0973E+03 1 A1 1.0742E+04 2 A2 -8.2031E+02
-----------------------------------------------------------------------------------------------------------------------
DRESP1 ID= 25 RESPONSE TYPE= STRESS ELEM ID= 3 COMP NO= 2 SEID= 0
SUBCASE RESP VALUE DESIGN VARIABLE COEFFICIENT DESIGN VARIABLE COEFFICIENT
-----------------------------------------------------------------------------------------------------------------------
1 -9.0973E+03 1 A1 1.0742E+04 2 A2 -8.2031E+02
2 1.3530E+04 1 A1 -1.1895E+04 2 A2 -8.1543E+02
MOVE LIMITS
true function
F(X)
linear approximation
o
X = current design
o X
X
allowable
move limits
• For sizing problems, use relatively large move limits (such as 0.5–
1.0 initially), reducing them if convergence becomes a problem.
9-52
MOVE LIMITS IMPOSED ON ANALYSIS MODEL
PROPERTIES
o o
• For the initial design x , analysis model properties p j , j = 1, 2, , ...,
J are determined. DELP specifies the bounds within which pj may
be modified as:
l u
pj ≤ pj ( x ) ≤ pj
where l = o o
pj p j – p j ⋅ DELP
u = o o
pj p j + p j ⋅ DELP
l o u
0 pj pj pj
9-53
MOVE LIMITS IMPOSED ON ANALYSIS MODEL
PROPERTIES (Cont.)
l l o
p j = max min [ p j ( DELP ), ( p j – DPMIN ) ], PMIN
u u o
p j = min max [ p j ( DELP ), ( p j + DPMIN ) ], PMAX
l u
pj pj
l pjo u
pj pj
0 1
• Note that PMIN and PMAX from the DVPREL1 and 2 entries have
been accounted for as well.
9-54
MOVE LIMITS IMPOSED ON DESIGN VARIABLES
o o o o
xi – xi ⋅ DELTA ≤ x i ≤ x i + xi ⋅ DELTA
or
or
1.0 by default
o
if x i ⋅ DELTA < DXMIN then
o o
xi – DXMIN ≤ x i ≤ xi + DXMIN
l u
However, x i and x i , provided on the DESVAR entry, are never
exceeded.
9-55
AUTOMATIC UPDATES OF MOVE LIMITS
• Every time the move limits are changed, a message is printed and
the user is advised to include a revised DOPTPRM entry reflecting
the changes in the restart Bulk Data
9-56
MOVE LIMITS UPDATE
^^^ USER INFORMATION MESSAGE 9052 (FEA)
^^^ STATIC ANALYSIS COMPLETED. DESIGN CYCLE NUMBER= 4
********************************************************
* *
* USER WARNING MESSAGE *
* *
* IF YOU WANT TO CONTINUE THE DESIGN AFTER THIS JOB *
* IS COMPLETED, YOU MUST INCLUDE A REVISED DOPTPRM *
* BULK DATA ENTRY IN THE BULK DATA SECTION WITH THE *
* FOLLOWING ITEMS MODIFIED AS SHOWN: *
* *
* DELP = 1.0000E-01 *
* DPMIN = 5.0000E-03 *
* DELX = 5.0000E-01 *
* DXMIN = 2.5000E-02 *
* *
********************************************************
********************************************************
9-57
* *
* IF A DELXV IS SPECIFIED ON A DESVAR BULK DATA ENTRY, *
* UPDATED DELXV VALUES ARE PRESENT ON THE DESVAR *
* ENTRIES CONTAINED IN THE PUNCH FILE. *
* *
* (NOTE: THERE MAY BE MORE THAN ONE MESSAGE LIKE *
* THIS. THE LAST ONE IN THIS RUN SHOULD *
* BE LOCATED AND USED.) *
* *
********************************************************
OPTEXIT
Value Description
0 Do not exit. Proceed with optimization.
1 Exit after the initialization of the analysis and design model
but before finite element analysis begins.
2 Exit after finite element analysis and initial design response
and shape basis vector processing.
3 Exit after design constraint evaluation and screening.
4 Exit after design sensitivity analysis and print the matrix of
design sensitivity coefficients (DSCM2). This is equivalent to
the DSAPRT (UNFORM,END+SENS) Case Control
command.
–4 Exit after design sensitivity analysis and write the data blocks
related to sensitivity coefficients (DSCM2 and DSCMCOL) to
an external file using the OUTPUT2 and OUTPUT4 modules.
This is equivalent to the DSAPRT (NOPRINT,EXPORT
END+SENS) Case Control command. See related
parameters ITAPE, IUNIT, and OMAXR.
5 Exit after the first approximate optimization of the design
model.
6 Exit after the first update of the analysis model based on the
first approximate optimization of the design model.
7 Compute and output design sensitivity coefficients at the end
of normal program termination: hard convergence, soft
convergence, or maximum design cycles. This is equivalent
to the DSAPRT (UNFORM,START+LAST) Case Control
command.
9-58
FLOW CHART FOR SOLUTION 200
NO Does a Design
Model Exist?
YES
1 EXITOPT
9-59
FLOW CHART FOR SOLUTION 200 (Cont.)
NO
If AMBS AMBS ≡ Auxiliary Model
Boundary Shapes
YES
Static Analysis for All Auxiliary Models,
Loop over Auxiliary Models, Boundary
Conditions (FEAOPT = Analysis)
FEA
NO If AMBS
GMBS ≡ Geometry Model
or GMBS
Boundary Shapes
YES
Collect Analysis Results, Generate
Enforced Boundary Displacements
for Auxiliary Structure Analysis
AXMDRV, DSAJ BNDSHP
Loop on Design Cycles
NO If AMBS
or GMBS
YES
Partition Out Shape Basis Vector
Solutions, Generate Basis Vectors
YES
Begin Design Sensitivity
and Optimization
9-60
FLOW CHART FOR SOLUTION 200 (Cont.)
Maximum
YES
Design Cycles EXIT EXITOPT 7
Loop on Design Cycles
Completed?
NO
3 EXITOPT
If Hard YES
Convergence
Data Recovery and
NO Exit (FEAOPT =
DATARECOVERY”) FEA,
EXITOPT
*
Generate Pseudo-load Vectors for 7
Sensitivity Analysis Across Analysis
Disciplines, Subcases, and Superelements DESAERDR
DSABO, DSAE, DSAM, DSAN, ,PSLGDV,
DSAP, DSAR, DSVG1, DSVG2 SEDRDR
*OPTEXIT=7 includes computation of sensitivities for the current design cycle (all operations up to DOM9).
9-61
FLOW CHART FOR SOLUTION 200 (Cont.)
NO 7
End Design
Cycle Loop
9-62
CONVERGENCE AT THE DESIGN CYCLE LEVEL
IFP
Initialization
PREDOM
FEA
9-63
CONVERGENCE AT THE DESIGN CYCLE
LEVEL (Cont.)
P (P) – P (P – 1)
max i i
CHGPRP
1 ≤ i ≤ NPROP
----------------------------------------
(P – 1) 0.001 CONVPR
Pi
x (P) – x (P – 1)
max i i
CHGDV 1 ≤ i ≤ NDV ---------------------------------------
(P – 1) 0.001 CONVDV
xi
9-64
CONVERGENCE CRITERIA
CHGOBJ ≤ CONV1
No
or
ACHOBJ ≤ CONV2
YES
No
CHGPRP ≤ CONVPR
YES
CONMAX ≤ GMAX No
or
CHGDV ≤ CONVDV
YES
9-65
CONVERGENCE CRITERIA (Cont.)
Yes
Yes No
CONMAX < GMAX
CHGPRP<CONVPR
Yes and No
CHGDV<CONVDV
9-66
SUMMARY OF PARAMETERS FOR DESIGN
OPTIMIZATION
9-67
SUMMARY OF PARAMETERS FOR DESIGN
OPTIMIZATION (Cont.)
9-68
SECTION 10
SHAPE OPTIMIZATION
DVGRID . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10-10
DVSHAP . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10-19
AUXCASE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10-30
AUXMODEL . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10-31
BNDGRID . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10-32
DVBSHAP . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10-33
where {∆G} = { G }
i + 1 – { G }i
= {x}
i + 1 – { x }i
{∆x}
i = current design
∆h = 1.0
1’ 2’
y
2
1
x 4’ h’
3’
1 ⁄ 2 ∆h
h
3 4
5 6
{ ∆G } = { G′ } – { G } = ∆x { T }
T
{ T } = 0. 1. 0. 0. 1. 0. 0. 0.5 0. 0. 0.5 0. 0. 0. 0. 0. 0. 0.
GRID 1 GRID 2 GRID 3 GRID 4 GRID 5 GRID 6
10-2
BASIS VECTORS – EXAMPLES (Cont.)
6
4
θ
5 3
R 5’
3’
∆R = -1.0
1’ 1 2
{ ∆G } = { G′ } – { G } = ∆x { T }
T
{ T } = – 1. 0. 0. 0. 0. 0. – 1. 0. 0. 0. 0. 0. – 1. 0. 0. 0. 0. 0.
GRID 1 GRID 2 GRID 3
10-3
BASIS VECTORS – EXAMPLES (Cont.)
Example 3
Consider a plate with a hole in the middle with unequal edge traction
forces. A quarter model is pictured below:
Ty = 1.0E4 N/m
0.15m
Material:
aluminum, 7075–T6 sheet
2
E = 7.2E10 N/m
µ = 0.33
3
ρ = 2.8E3 kg/m
10-4
BASIS VECTORS – AN EXAMPLE (Cont.)
x1 = a
x2 = b
so that,
∆x 1 ∆a
{ ∆G } = [ { T 1 } { T 2 } ]
∆x 2 ∆b
x-direction y-direction
variation variation
10-5
SHAPE BASIS VECTORS IN THE DESIGN MODEL
10-6
DESIGN MODELING INPUT FOR SHAPE
OPTIMIZATION
10-7
MANUAL GRID VARIATION
{∆G}i
z z
{G}i
y
y x
∆G ix N1
{ ∆G } i = ∆G iy = COEFF N2 ⋅ ∆x DVID
∆G iz N3
1 2 3 4 5 6 7 8 9 10
DVGRID DVID GRID CID COEFF N1 N2 N3
10-8
MANUAL GRID VARIATION (Cont.)
10-9
DVGRID Design Variable to Grid Point Relation
DVGRID the relationship between design variables and grid point locations.
Defines Design Variable to Grid Point Relation
Format:
1 2 3 4 5 6 7 8 9 10
DVGRID DVID GID CID COEFF N1 N2 N3
Example:
Field Contents
DVID DESVAR entry identification number. (Integer > 0)
GID Grid point (GRID) or geometric point (POINT) identification number. (Integer > 0)
Remarks:
1. A CID of zero or blank (the default) references the basic coordinate system.
2. Multiple references to the same grid ID and design variable result in vectorial addition of
the participation vectors defined by CID, COEFF, and Ni. There is no restriction on the
number of DVGRID entries that may reference a given grid (GID) or design variable
(DVID).
T
where { g } i is the location of the i-th grid, [ g x g y g z ] .
T
The vector { N } = [ N x N y N z ] is determined from CID and Ni. Note that it is a change
0
in a design variable from its initial value X , and not the absolute value of the design
0
variable itself, that represents a change in a grid point location, { g } i – { g } i .
(Continued)
10-10
Design Variable to Grid Point Relation DVGRID
4. The DVGRID entry defines the participation coefficients (basis vectors) of each design
variable for each of the coordinates affected by the design process in the relationship
{ ∆g } i = ∑ { T }ij ∗ ∆X j
j
5. DVGRID entries that reference grid points on MPCs or RSSCON entries produce incorrect
sensitivities. Often the sensitivities are 0.0 which may result in a warning message indi-
cating zero gradients which may be followed by UFM 6499. Other rigid elements produce
correct results.
10-11
AUXILIARY MODELS IN SHAPE OPTIMIZATION
10-12
AUXILIARY MODELS IN SHAPE OPTIMIZATION (Cont.)
{ ∆G } = [ T ] { ∆x }
10-13
EXAMPLE – DIRECT INPUT OF SHAPES
The figure below is the initial design and von Mises stress resultants
for a two-dimensional, symmetric road support structure (culvert)
From the above stress distribution, we can observe that the culvert
interior profile is not optimal.
10-14
EXAMPLE – DIRECT INPUT OF SHAPES (Cont.)
10-15
EXAMPLE – DIRECT INPUT OF SHAPES (Cont.)
10-16
EXAMPLE – DIRECT INPUT OF SHAPES (Cont.)
CQUAD4, 1,101, 1, 2, 6, 5
CQUAD4, 2,101, 2, 3, 7, 6
CQUAD4, 3,101, 3, 4, 8, 7
. . . . . . .
. . . . . . .
(see optimization input file)
. . . . . . .
. . . . . . .
CQUAD4, 25,101, 30, 31, 38, 37
CQUAD4, 26,101, 31, 32, 39, 38
CQUAD4, 27,101, 32, 33, 40, 39
PSHELL,101,102,.44
MAT1,102,2.+7,,.3
$
$ Additional CBAR elements maintain smoothness of the circular boundary
$
CBAR,31,1,13,14,,1.0
CBAR,32,1,14,15,,1.0
CBAR,33,1,15,16,,1.0
CBAR,34,1,16, 9,,1.0
CBAR,35,1, 9, 5,,1.0
CBAR,36,1,5 , 1,,1.0
PBAR 1 102 20.0 1.0 1.0
$
$ Seven load cases
$
FORCE,100,13,0,1.e5,0.,1.,0.
FORCE,101,14,0,1.e5,0.259,.9659
FORCE,102,15,0,1.e5,0.5,0.866,0.0
FORCE,103,16,0,1.e5,1.,1.,0.
FORCE,104,9,0,1.e5,0.866,0.5,0.0
FORCE,105,5,0,1.e5,0.9659,0.259
FORCE,106,1,0,1.e5,1.,0.,0.
$
$ Boundary conditions satisfy functional and manufacturing requirements
$
SPC1,25,345,1,THRU,40
SPC1,25,6,2,THRU,4
SPC1,25,6,6,THRU,8
SPC1,25,6,10,THRU,12
SPC1,25,6,17,THRU,19
SPC1,25,6,20,THRU,26
SPC1,25,6,27,THRU,33
SPC1,25,6,34,THRU,40
SPC1,25,12,33,THRU,40
SPC1,25,12,4,8,12,19,26
SPC1,25,1,13,20,27
SPC1,25,2,1,2,3
ENDDATA
10-17
EXAMPLE – DIRECT INPUT OF SHAPES (Cont.)
Optimization input
10-18
Design Variable to Basis Vector(s) DVSHAP
Defines a shape basis vector by relating a design variable identificationDesign
DVSHAP number (DVID)
Variable to columns
to Basis Vector(s)
of a displacement matrix.
Format:
1 2 3 4 5 6 7 8 9 10
DVSHAP DVID COL1 SF1 COL2 SF2 COL3 SF3
Example:
Field Contents
DVID Design variable identification number on the DESVAR entry. (Integer > 0)
SFi Scaling factor applied to the COLi-th column of the displacement matrix. (Real;
Default = 1.0)
Remarks:
1. DVID must be defined on a DESVAR entry.
3. Multiple references to the same DVID and/or COLi will result in a linear combination of
displacement vectors. In the example above, the shape basis vector is a linear combina-
tion of the fourth column and twice the second column.
4. The displacement matrix must have been created by MSC/NASTRAN and be available on
a database, which is attached via the DBLOCATE FMS statement shown below:
10-19
EXAMPLE – DIRECT INPUT OF SHAPES (Cont.)
$
$ FMS section for retrieving the auxiliary displacement matrix
$
assign f1_aux=’culvert1.MASTER’
dblocate datablk=(ug/ugd,geom1/geom1d,geom2/geom2d) ,
logical=f1_aux
SOL 200 $
TIME 100
CEND
TITLE=CULVERT EXAMPLE USING EXTERNAL AUXILIARY STRUCTURE
SUBTITLE=THE PRIMARY STRUCTURE
ANALYSIS = STATICS
SPC=25
LOAD=1
DISP=ALL
STRESS=all
DESSUB = 10
desobj = 5
BEGIN BULK
PARAM,POST,0
$ PARAM,optexit,4
PARAM,NEWSEQ,-1
GRID, 1,, 3.00000, 0.00000,.00
GRID, 2,, 4.00000, 0.00000,.00
GRID, 3,, 5.00000, 0.00000,.00
GRID, 4,, 6.00000, 0.00000,.00
GRID, 5,, 2.89464, 0.78478,.00
GRID, 6,, 3.79369, 0.75885,.00
GRID, 7,, 4.69274, 0.73293,.00
GRID, 8,, 5.59178, 0.70700,.00
GRID, 9,, 2.60164, 1.49178,.00
GRID, 10,, 3.46229, 1.46585,.00
GRID, 11,, 4.32293, 1.43993,.00
GRID, 12,, 5.18357, 1.41400,.00
GRID, 13,, 0.00000, 3.00000,.00
GRID, 14,, 0.78478, 2.89464,.00
GRID, 15,, 1.49178, 2.60164,.00
GRID, 16,, 2.12100, 2.12100,.00
GRID, 17,, 3.00578, 2.12100,.00
GRID, 18,, 3.89057, 2.12100,.00
GRID, 19,, 4.77535, 2.12100,.00
10-20
EXAMPLE – DIRECT INPUT OF SHAPES (Cont.)
10-21
EXAMPLE – DIRECT INPUT OF SHAPES (Cont.)
10-22
EXAMPLE – DIRECT INPUT OF SHAPES (Cont.)
Results:
10-23
AUXILIARY BOUNDARY MODELS
Bar Elements
10-24
AUXILIARY BOUNDARY MODELS
δ δ
Constant
Linear
Quadratic
Cubic
10-25
EXAMPLE – ANALYTIC BOUNDARY SHAPES
165
154
143
132
121
110
Z
Y
10-26
EXAMPLE – ANALYTIC BOUNDARY SHAPES (Cont.)
Z
Y
To change the shape, we will modify the upper and lower planes of
the cantilever.
10-27
EXAMPLE – ANALYTIC BOUNDARY SHAPES (Cont.)
Z
Y
x1
[ T ] = [ 1.0U SUBCASE200 1.0U SUBCASE300 ]
x2
10-28
EXAMPLE – ANALYTIC BOUNDARY SHAPES (Cont.)
AUXCASE
AUXMODEL
BNDGRID
DVBSHAP
10-29
AUXCASE Auxiliary Model Case Control Delimiter
Format:
AUXCASE
Examples:
AUXCAS
AUXC
Remarks:
1. AUXCASE indicates the beginning of Case Control commands for an auxiliary model.
AUXCASE must follow the primary model Case Control commands.
2. All Case Control commands following this entry are applicable until the next AUXCASE or
BEGIN BULK command. Commands from preceding Case Control Sections are ignored.
3. Each auxiliary model Case Control must be delimited with the AUXCASE command.
4. The AUXMODEL command is used to associate the auxiliary model Case Control with a
particular auxiliary model.
10-30
Auxiliary Model Identification Number AUXMODEL
AUXMODEL an auxiliary model for generation of boundary shapes in Auxiliary
References shape optimization.
Model Identification Number
Format:
AUXMODEL = n
Examples:
AUXMODEL = 4
AUXM = 4
Describer Meaning
n Auxiliary model identification number. (Integer > 0)
Remarks:
1. AUXMODEL references a particular auxiliary model for analysis and may only be specified
in the auxiliary model Case Control Section.
2. See the BEGIN BULK command for the Bulk Data definition of an auxiliary model.
10-31
BNDGRID Boundary Grid Points
Format:
1 2 3 4 5 6 7 8 9 10
BNDGRID C GP1 GP2 GP3 GP4 GP5 GP6 GP7
GP8 -etc.-
Example:
BNDGRID 123 41 42 43 44 45 46 47
49
Field Contents
C Component number (any unique combination of integers 1 through 6 with no
embedded blanks). See Remark 1.
GPi Shape boundary grid point identification number. (0 < Integer < 1000000; For
THRU option, GP1< GP2)
Remarks:
1. C specifies the components for the listed grid points for which boundary motion is
prescribed.
2. Multiple BNDGRID entries may be used to specify the shape boundary grid point identifi-
cation numbers.
3. Both fixed and free shape boundary grid point identification numbers are listed on this
entry.
4. The degrees of freedom specified on BNDGRID entries must be sufficient to statically con-
strain the model.
5. Degrees of freedom specified on this entry form members of the mutually exclusive s-set.
They may not be specified on other entries that define mutually exclusive sets. See the
MSC/NASTRAN Quick Reference Guide, Appendix B for a list of these entries.
10-32
Design Variable to Boundary Shapes DVBSHAP
Associates a design variable identification number to a linear combination
DVBSHAP of to
Design Variable boundary shape
Boundary Shapes
vectors from a particular auxiliary model.
Format:
1 2 3 4 5 6 7 8 9 10
DVBSHAP DVID AUXMOD COL1 SF1 COL2 SF2 COL3 SF3
Example:
DVBSHAP 4 1 1 1.6
Field Contents
DVID Design variable identification number of a DESVAR entry. (Integer > 0)
COLi Load sequence identification number from AUXMODEL Case Control command.
(Integer > 0)
SFi Scaling factor for load sequence identification number. (Real; Default = 1.0)
Remarks:
1. Design variable DVID must be defined on a DESVAR entry.
2. Multiple references to the same DVID and/or COLi will result in the vector addition of the
referenced boundary shape vectors.
10-33
EXAMPLE – ANALYTIC BOUNDARY SHAPES (Cont.)
Modeling input:
TIME 600 $
SOL 200 $
CEND
TITLE = CANTILEVERED BEAM - HEXA **** D200AM3 **** D200AM3
DESOBJ = 15
DESSUB = 100
SUBCASE 100
ANALYSIS = STATICS
SPC = 1
LOAD = 1
DISPLACEMENT = ALL
OUTPUT(PLOT)
SET 1 = ALL
VIEW 90.0,0.0,90.0
PLOT SET 1
$SET 2 = ALL
$VIEW 34.0, 24.0, 0.0
$PLOT SET 2
$PLOT STATIC DEFORMATION SET 2
AUXCASE
TITLE = AUXILIARY MODEL 1
AUXMODEL = 1
SUBCASE 200
SPC = 200
LOAD = 220
LABEL = UPPER
SUBCASE 300
SPC = 300
LOAD = 330
LABEL = LOWER
BEGIN BULK
$
$------------------------------------------------------------------------
$ ANALYSIS MODEL:
$------------------------------------------------------------------------
$
PARAM AUTOSPC YES
PARAM POST 0
PARAM GRDPNT 0
PARAM MAXRATIO1.0E+8
$
CORD2S* 2 0 0.0 0.0+1A 2
*1A 2 0.0 0.0 0.0 1000.000+1B 2
*1B 2 1000.000 0.0 0.0 +1C 2
*1C 2
CORD2C* 1 0 0.0 0.0+1A 1
*1A 1 0.0 0.0 0.0 1000.000+1B 1
*1B 1 1000.000 0.0 0.0 +1C 1
*1C 1
10-34
EXAMPLE – ANALYTIC BOUNDARY SHAPES (Cont.)
10-35
EXAMPLE – ANALYTIC BOUNDARY SHAPES (Cont.)
$
$DLINK, ID, DDVID, CO, CMULT, IDV1, C1, IDV2, C2, +
$+, IDV3, C3, ...
DLINK 1 2 1.0 1 1.0
$
$ BOUNDARY CONDITIONS FOR SHAPE INTERPOLATIONS:
$
$ ---TOP SURFACE:
$
BNDGRID,C, GP1, GP2, GP3, GP4, GP5, GP6, GP7, +
$+, GP8, ...
BNDGRID 123 133 134 135 136 137 138 139
140 141 142 143 144 145 146 147
148 149 150 151 152 153 154 155
156 157 158 159 160 161 162 163
164 165
$
$ ---BOTTOM SURFACE:
BNDGRID 123 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 31
32 33
$
$ ---EXTERIOR SURFACES - INTERPOLATION IN X&Z DIRECTION ONLY:
BNDGRID 2 34 35 36 37 38 39 40
41 42 43 44
BNDGRID 2 56 57 58 59 60 61 62
63 64 65 66
BNDGRID 2 67 68 69 70 71 72 73
74 75 76 77
BNDGRID 2 89 90 91 92 93 94 95
96 97 98 99
BNDGRID 2 100 101 102 103 104 105 106
107 108 109 110
BNDGRID 2 122 123 124 125 126 127 128
129 130 131 132
$
$ ---TIP END:
BNDGRID 1 11 22 33 44 55 66
BNDGRID 1 77 88 99 110 121 132
BNDGRID 1 143 154 165
$
$ ---FIXED END:
BNDGRID 123 1 12 23 34 56 67 89
100 122 133 144 155
BNDGRID 1 45 78 111
$
$ FORMULATE WEIGHT-BASED SYNTHETIC RESPONSE: F = 1.E5*W
DRESP1 1 WEIGHT WEIGHT
DRESP2 15 WE1000 1
+ DRESP1 1
DEQATN 1 F(A)=100000.*A
10-36
EXAMPLE – ANALYTIC BOUNDARY SHAPES (Cont.)
$
$ CONSTRAINTS ON VON MISES STRESSES:
DRESP1 2 STRESS STRESS PSOLID 13 1
DSCREEN STRESS -1.0 10
DCONSTR 100 2 200.
$
$ OVERRIDE OF OPTIMIZATION PARAMETERS
DOPTPRM DESMAX 9 P1 1 P2 15
$
$------------------------------------------------------------------------
$ AUXILIARY BOUNDARY MODEL(S):
$------------------------------------------------------------------------
$
BEGIN BULK AUXMODEL=1
PARAM,PRGPST,NO
PARAM MAXRATIO1.0E+8
PARAM,AUTOSPC,YES
$
$ LOWER SURFACE:
CQUAD4 1000 2 1 2 13 12 0.0
CQUAD4 1001 2 2 3 14 13 0.0
CQUAD4 1002 2 3 4 15 14 0.0
CQUAD4 1003 2 4 5 16 15 0.0
CQUAD4 1004 2 5 6 17 16 0.0
CQUAD4 1005 2 6 7 18 17 0.0
CQUAD4 1006 2 7 8 19 18 0.0
CQUAD4 1007 2 8 9 20 19 0.0
CQUAD4 1008 2 9 10 21 20 0.0
CQUAD4 1009 2 10 11 22 21 0.0
CQUAD4 1010 2 12 13 24 23 0.0
CQUAD4 1011 2 13 14 25 24 0.0
CQUAD4 1012 2 14 15 26 25 0.0
CQUAD4 1013 2 15 16 27 26 0.0
CQUAD4 1014 2 16 17 28 27 0.0
CQUAD4 1015 2 17 18 29 28 0.0
CQUAD4 1016 2 18 19 30 29 0.0
CQUAD4 1017 2 19 20 31 30 0.0
CQUAD4 1018 2 20 21 32 31 0.0
CQUAD4 1019 2 21 22 33 32 0.0
$
$ UPPER SURFACE:
CQUAD4 950 2 133 144 145 134 0.0
CQUAD4 951 2 134 145 146 135 0.0
CQUAD4 952 2 135 146 147 136 0.0
CQUAD4 953 2 136 147 148 137 0.0
CQUAD4 954 2 137 148 149 138 0.0
CQUAD4 955 2 138 149 150 139 0.0
CQUAD4 956 2 139 150 151 140 0.0
CQUAD4 957 2 140 151 152 141 0.0
CQUAD4 958 2 141 152 153 142 0.0
CQUAD4 959 2 142 153 154 143 0.0
CQUAD4 960 2 144 155 156 145 0.0
CQUAD4 961 2 145 156 157 146 0.0
CQUAD4 962 2 146 157 158 147 0.0
10-37
EXAMPLE – ANALYTIC BOUNDARY SHAPES (Cont.)
10-38
EXAMPLE – ANALYTIC BOUNDARY SHAPES (Cont.)
Optimization results:
***************************************************************
S U M M A R Y O F D E S I G N C Y C L E H I S T O R Y
***************************************************************
(HARD CONVERGENCE ACHIEVED)
(SOFT CONVERGENCE ACHIEVED)
NUMBER OF FINITE ELEMENT ANALYSES COMPLETED 6
NUMBER OF OPTIMIZATIONS W.R.T. APPROXIMATE MODELS 5
OBJECTIVE AND MAXIMUM CONSTRAINT HISTORY
---------------------------------------------------------------------------------------------------------------
OBJECTIVE FROM OBJECTIVE FROM FRACTIONAL ERROR MAXIMUM VALUE
CYCLE APPROXIMATE EXACT OF OF
10-39
Final Shape:
Z
Y
Z
Y
10-40
SUMMARY OF SHAPE BASIS VECTOR
GENERATION METHODS
10-41
SUMMARY OF SHAPE BASIS VECTOR GENERATION
METHODS (Cont.)
10-42
GUIDELINES, RECOMMENDATIONS AND
LIMITATIONS
• Use smaller move limits for a shape design variable (delxv = 0.5 or
less)
10-43
GUIDELINES, RECOMMENDATIONS AND
LIMITATIONS (Cont.)
10-44
SECTION 11
• Direct Frequency
• Modal Frequency
• Modal Transient
• Displacement
• Velocity
• Acceleration
• SPC Force
• Stress
• Element force
• Equations (DRESP2+DEQATN)
• Weight
• Volume
11-1
GENERAL CONSIDERATIONS
Design responses are computed just for the output frequencies and
output time steps:
⇒ OFREQ
• Property changes
• Shape changes
• Superelements
In short, all of the design modeling tools available for other analysis
disciplines are available for dynamic optimization as well.
11-2
BASIC EQUATIONS
• Equation of motion
2
[ – ω M + iωB + K ]u = F
∂
∆ = -------
∂x i
2 2
[ – ω M + iωB + K ]∆u = ∆F – [ – ω ∆M + iω∆B + ∆K ]u
• Modal transformation,
u = φξ
Differentiate
∂
-------
∂x i
and apply
[ ∆u = ∆φξ + φ∆ξ = ∆η ]
11-3
BASIC EQUATIONS (Cont.)
2 T
[ – ω φ Mφ + iωφ T Bφ + φ T Kφ ]∆η
T T 2
= φ ∆F – φ [ – ω ∆M + iω∆B + ∆K ]U
• Equation of motion
Mu̇˙ + Bu̇ + Ku = F
T T T
[ φ Mφ ]∆η̇˙ + [ φ Bφ ]∆η̇ + [ φ Kφ ]∆η
T T
= φ ∆F – φ ( ∆Mu̇˙ + ∆Bu̇ + ∆Ku )
11-4
BASIC EQUATIONS (Cont.)
• Thermal loading
K = K 1 + igK 1 + i ∑ ge ke + K2
e
• Limitations:
11-5
EXAMPLE – MODAL FREQUENCY RESPONSE
OPTIMIZATION
A simple pressure loaded flat plate is shown in the figure below. The
pressure loading is a constant 1.0 lbf/in**2 applied over a frequency
range of 20.0 to 200.0 Hz. We would like to minimize the transverse
mean square displacement of the free edge midside grid over the
same frequency range.
Grid 1110
X
(Uniform Pressure)
11-6
EXAMPLE – MODAL FREQUENCY RESPONSE
OPTIMIZATION (Cont.)
Symmetric
y Boundary Conditions
Grid 1110
t1 = x1 Free
Edge
100 200
2 2
∑ ∑
i 2i
minφ = ( u z, 1110 ) + 2 ( u z, 1110 )
i = 20 i = 51
11-7
EXAMPLE – MODAL FREQUENCY RESPONSE
OPTIMIZATION (Cont.)
Input file:
11-8
EXAMPLE – MODAL FREQUENCY RESPONSE
OPTIMIZATION (Cont.)
11-9
EXAMPLE – MODAL FREQUENCY RESPONSE
OPTIMIZATION (Cont.)
11-10
EXAMPLE – MODAL FREQUENCY RESPONSE
OPTIMIZATION (Cont.)
11-11
EXAMPLE– MODAL FREQUENCY RESPONSE
OPTIMIZATION (Cont.)
$
DRESP1 201 VOLUME VOLUME
DRESP1 20 g1110L FRDISP 3 20.0 1110
= *(1) = = = = = *(1.0) =
=79
DRESP1 102 G1110H FRDISP 3 102.0 1110
= *(2) = = = = = *(2.0) =
=48
$
DRESP2 1 UZ2 1
DRESP1 20 21 22 23 24 25 26
27 28 29 30 31 32 33
34 35 36 37 38 39 40
41 42 43 44 45 46 47
48 49 50 51 52 53 54
55 56 57 58 59 60 61
62 63 64 65 66 67 68
69 70 71 72 73 74 75
76 77 78 79 80 81 82
83 84 85 86 87 88 89
90 91 92 93 94 95 96
97 98 99 100 102 104 106
108 110 112 114 116 118 120
122 124 126 128 130 132 134
136 138 140 142 144 146 148
150 152 154 156 158 160 162
164 166 168 170 172 174 176
178 180 182 184 186 188 190
192 194 196 198 200
$
DEQATN 1 UZ2(U20,u21,u22,u23,u24,u25,U26,U27,U28,U29,U30,
U31,U32,U33,U34,U35,U36,U37,U38,U39,U40,
U41,U42,U43,U44,U45,U46,U47,U48,U49,U50,
U51,U52,U53,U54,U55,U56,U57,U58,U59,U60,
U61,U62,U63,U64,U65,U66,U67,U68,U69,U70,
U71,U72,U73,U74,U75,U76,U77,U78,U79,U80,
U81,U82,U83,U84,U85,U86,U87,U88,U89,U90,
U91,U92,U93,U94,U95,U96,U97,U98,U99,U100,
U102,U104,U106,U108,U110,U112,U114,U116,U118,U120,
U122,U124,U126,U128,U130,U132,U134,U136,U138,U140,
U142,U144,U146,U148,U150,U152,U154,U156,U158,U160,
U162,U164,U166,U168,U170,U172,U174,U176,U178,U180,
U182,U184,U186,U188,U190,U192,U194,U196,U198,U200)
= u20**2 + u21**2 + U22**2 + U23**2+ U24**2 +
U25**2 + U26**2 + U27**2 + U28**2 + U29**2 +U30**2 +
U31**2 + U32**2 + U33**2 + U34**2 + U35**2 +
U36**2 + U37**2 + U38**2 + U39**2 + U40**2 +
U41**2 + U42**2 + U43**2 + U44**2 + U45**2 +
U46**2 + U47**2 + U48**2 + U49**2 + U50**2 +
U51**2 + U52**2 + U53**2 + U54**2 + U55**2 +
u56**2 + u57**2 + u58**2 + u59**2 + u60**2 +
U61**2 + U62**2 + U63**2 + U64**2 + U65**2 +
U66**2 + U67**2 + U68**2 + U69**2 + U70**2 +
11-12
EXAMPLE – MODAL FREQUENCY RESPONSE
OPTIMIZATION (Cont.)
11-13
EXAMPLE – MODAL FREQUENCY RESPONSE OPTIMIZATION (Cont.)
Results:
***************************************************************
S U M M A R Y O F D E S I G N C Y C L E H I S T O R Y
***************************************************************
----------------------------------------------------------------------------------------------------------------------------------
INTERNAL | EXTERNAL | |
DV. ID. | DV. ID. | LABEL | 6 : 7 : 8 : 9 : 10 : 11 :
----------------------------------------------------------------------------------------------------------------------------------
1 | 1 | T1 | 9.2234E-02 : 1.0245E-01 :
2 | 2 | T2 | 8.9215E-02 : 7.8416E-02 :
3 | 3 | T3 | 6.6226E-02 : 6.4560E-02 :
4 | 4 | T4 | 3.6059E-02 : 3.2658E-02 :
5 | 5 | T5 | 5.5572E-02 : 5.8001E-02 :
6 | 6 | T6 | 7.6552E-02 : 7.8469E-02 :
7 | 7 | T7 | 9.6333E-02 : 9.7656E-02 :
8 | 8 | T8 | 1.3458E-01 : 1.3406E-01 :
9 | 9 | T9 | 1.5716E-01 : 1.5932E-01 :
10 | 10 | T10 | 1.4968E-01 : 1.5128E-01 :
*** USER INFORMATION MESSAGE 6464 (DOM12E)
RUN TERMINATED DUE TO HARD CONVERGENCE TO AN OPTIMUM AT CYCLE NUMBER = 7.
EXAMPLE – MODAL FREQUENCY RESPONSE
OPTIMIZATION (Cont.)
3.622 @ 58 Hz
2.929 @ 52 Hz
11-16
EXAMPLE – MODAL FREQUENCY RESPONSE
OPTIMIZATION (Cont.)
Thickness
1 2 3 4 5 6 7 8 9 10
11-17
EXAMPLE – ACOUSTIC OPTIMIZATION
The figure below shows an acoustic box with an acoustic source and
a transducer located on opposite walls of the box (prop 3 ends). We
would like to modify the box plate thicknesses to minimize the
response peaks within the fluid (air). This is to be done without
changing the total weight of the box.
11-18
EXAMPLE – ACOUSTIC OPTIMIZATION (Cont.)
P(f)
c1 c2 c3
βx
f
f1 f2 f3
minimize βx
subject to P ( f1 ) – βx ≤ 0
P ( f2 ) – βx ≤ 0
P ( f3 ) – βx ≤ 0
11-19
EXAMPLE – ACOUSTIC OPTIMIZATION (Cont.)
Input File:
ID MSC-XL, MSC-NASTRAN
TIME 100
SOL 200 $ modal frequency response
CEND
SUBTITLE = acoustic and structural elements
LABEL = boxae1.dat
$
$ analysis case control:
set 20 = 11280
echo = sort(param,eigc,eigrl,freq,desvar,dconstr,dresp1,dresp2,deqatn,
dvprel1)
spc = 1
DISP(phase) = 20
method(structure) = 20
method(fluid) = 30
cmethod = 10
frequency = 200
dload = 100
partn = 20
$ optimization case control:
ANALYSIS = MFREQ
DESGLB = 5
DESSUB = 10
DESOBJ = 100
$
BEGIN BULK
$
EIGRL 20 200. max
eigrl ,30, 15., 155., 9, 0, , 105., max
eigc,10,clan,max,,,,
,15.,300.,,,,,4
$
$ sound pressure level
param,rms,yes
$ reference pressure for dB and dBA
param,prefdb,2.-5
$
PARAM AUTOSPC no
$
$ fluid/structure interface
acmodl,diff , , , ,0.01
$
$ structural damping:
param,g,0.02
$
11-20
EXAMPLE – ACOUSTIC OPTIMIZATION (Cont.)
$-------2-------3-------4-------5-------6-------7-------8-------9-------10----
rload1,100,101,,,102
darea,101,1288,3,100.
tabled1,102
,0.,1.,1000.,1.,endt
$freq2, 200, 1.0, 200.0, 50
$-------2-------3-------4-------5-------6-------7-------8-------9-------10----
freq 200 40. 50. 60. 70. 80. 90. 95.
100. 105. 110. 120. 130. 140. 150. 160.
170. 180. 190. 200. 97.5 102.5
$-------2-------3-------4-------5-------6-------7-------8-------9-------10----
$
$
$ Define the Design Variables:
$
$ESVAR, ID, LABEL, XINIT, XLB, XUB, DELXV
DESVAR, 1, P1, 0.02493,0.0001, 1.
DESVAR, 2, P2, 0.01953,0.0001, 1.
DESVAR, 4, P4, 0.02047,0.0001, 1.
DESVAR, 5, P5, 0.02596,0.0001, 1.
DESVAR, 6, P6, 0.02175,0.0001, 1.
DESVAR, 7, P7, 0.02426,0.0001, 1.
DESVAR, 8, BETA 1.0, 0.001
$
$ Relate the Design Variables to changes in plate thicknesses:
$
$VPREL1,ID, TYPE, PID, FID, PMIN, PMAX, C0, , +
$+, DVID1, COEF1, DVID2, COEF2, ...
$
DVPREL1,1, PSHELL, 1, 4, 0.0001, , , , +
+, 1, 1.
DVPREL1,2, PSHELL, 2, 4, 0.0001, , , , +
+, 2, 1.
DVPREL1,4, PSHELL, 4, 4, 0.0001, , , , +
+, 4, 1.
DVPREL1,5, PSHELL, 5, 4, 0.0001, , , , +
+, 5, 1.
DVPREL1,6, PSHELL, 6, 4, 0.0001, , , , +
+, 6, 1.
DVPREL1,7, PSHELL, 7, 4, 0.0001, , , , +
+, 7, 1.
$
$ Define the synthetic Objective as a function of X8:
$
$RESP2, ID, LABEL, EQID, REGION, , , , , +
$+, DESVAR, DVID1, DVID2, ..., , , , , +
DRESP2 100 BETA 100
DESVAR 8
DEQATN 100 OBJ(BETA) = 10000.0 * BETA
$
11-21
EXAMPLE – ACOUSTIC OPTIMIZATION (Cont.)
11-22
EXAMPLE – ACOUSTIC OPTIMIZATION (Cont.)
11-23
EXAMPLE – ACOUSTIC OPTIMIZATION (Cont.)
11-24
EXAMPLE – ACOUSTIC OPTIMIZATION (Cont.)
Results:
***************************************************************
S U M M A R Y O F D E S I G N C Y C L E H I S T O R Y
***************************************************************
----------------------------------------------------------------------------------------------------------------------------------
INTERNAL | EXTERNAL | |
DV. ID. | DV. ID. | LABEL | INITIAL : 1 : 2 : 3 : 4 : 5 :
----------------------------------------------------------------------------------------------------------------------------------
1 | 1 | P1 | 2.4930E-02 : 3.4858E-02 : 3.2397E-02 : 2.8322E-02 : 3.4406E-02 : 3.5523E-02 :
2 | 2 | P2 | 1.9530E-02 : 2.6421E-02 : 2.6187E-02 : 3.3966E-02 : 3.3395E-02 : 3.2716E-02 :
3 | 4 | P4 | 2.0470E-02 : 1.0487E-02 : 1.2030E-02 : 1.3776E-02 : 7.8943E-03 : 6.4859E-03 :
4 | 5 | P5 | 2.5960E-02 : 2.0480E-02 : 2.1393E-02 : 1.3477E-02 : 2.3186E-02 : 2.2668E-02 :
5 | 6 | P6 | 2.1750E-02 : 1.8127E-02 : 2.1821E-02 : 2.3308E-02 : 2.1010E-02 : 1.9540E-02 :
6 | 7 | P7 | 2.4260E-02 : 1.9210E-02 : 1.5891E-02 : 1.4294E-02 : 9.1583E-03 : 1.0095E-02 :
7 | 8 | BETA | 1.0000E+00 : 5.3345E-01 : 1.2462E-01 : 1.8099E-01 : 2.0766E-01 : 2.0756E-01 :
----------------------------------------------------------------------------------------------------------------------------------
INTERNAL | EXTERNAL | |
DV. ID. | DV. ID. | LABEL | 6 : 7 : 8 : 9 : 10 : 11 :
11-26
----------------------------------------------------------------------------------------------------------------------------------
1 | 1 | P1 | 3.4322E-02 : 2.9154E-02 : 2.4363E-02 : 2.7226E-02 : 2.4670E-02 : 2.4619E-02 :
2 | 2 | P2 | 3.1484E-02 : 2.6989E-02 : 3.1301E-02 : 2.5696E-02 : 2.3196E-02 : 2.3115E-02 :
3 | 4 | P4 | 5.6109E-03 : 6.5327E-03 : 5.8686E-03 : 5.4315E-03 : 4.3781E-03 : 3.8628E-03 :
4 | 5 | P5 | 2.4430E-02 : 3.1121E-02 : 3.3386E-02 : 3.8369E-02 : 4.0770E-02 : 4.0748E-02 :
5 | 6 | P6 | 1.8545E-02 : 1.7878E-02 : 1.7679E-02 : 1.4659E-02 : 1.6339E-02 : 1.6324E-02 :
6 | 7 | P7 | 1.3480E-02 : 1.9997E-02 : 2.1885E-02 : 2.5364E-02 : 2.7816E-02 : 2.7732E-02 :
7 | 8 | BETA | 2.0764E-01 : 2.1034E-01 : 2.1049E-01 : 2.1332E-01 : 2.4781E-01 : 1.8586E-01 :
----------------------------------------------------------------------------------------------------------------------------------
INTERNAL | EXTERNAL | |
DV. ID. | DV. ID. | LABEL | 12 : 13 : 14 : 15 : 16 : 17 :
----------------------------------------------------------------------------------------------------------------------------------
1 | 1 | P1 | 2.4258E-02 : 2.6005E-02 : 2.6021E-02 : 2.6021E-02 :
2 | 2 | P2 | 2.3058E-02 : 2.1841E-02 : 2.1834E-02 : 2.1834E-02 :
3 | 4 | P4 | 3.3257E-03 : 2.7535E-03 : 2.4141E-03 : 2.4141E-03 :
4 | 5 | P5 | 4.3248E-02 : 4.3566E-02 : 4.3587E-02 : 4.3587E-02 :
5 | 6 | P6 | 1.4020E-02 : 1.1520E-02 : 1.1632E-02 : 1.1632E-02 :
6 | 7 | P7 | 2.7829E-02 : 2.9923E-02 : 2.9964E-02 : 2.9964E-02 :
7 | 8 | BETA | 1.3939E-01 : 1.4246E-01 : 1.0685E-01 : 1.0685E-01 :
140.5 dB @ 100 Hz
115.9 dB @ 100 Hz
11-27
EXAMPLE: ACOUSTIC OPTIMIZATION (Cont.)
Structural
Eigenfrequencies
Mode Fluid
No. Initial Final Eigenfrequencies
1 75.99 7.861680E+01 50.01
2 95.29 8.777968E+01 100.41
3 104.16 9.964240E+01 100.41
4 130.78 1.187981E+02 100.41
5 133.35 1.262226E+02 112.19
6 143.74 1.309330E+02 112.19
7 153.16 1.448834E+02 142.00
8 173.69 1.457052E+02 142.00
9 – – 142.00
11-28
SECTION 12
SUPERELEMENT OPTIMIZATION
Design variables:
12-1
GENERAL CONSIDERATIONS (Cont.)
Design responses:
Supported Superelements:
12-2
CASE CONTROL SPECIFICATION
Single Discipline:
• Condensed (SUPER=ALL)
• Expanded
Multidisciplinary
ID MSC, D200SE3 $
SOL 200 $ OPTIMIZATION
CEND
TITLE=STATIC ANALYSIS OF A 25-BAR TRUSS D200SE3
DISP=ALL
SPC=100
ANALYSIS=STATICS
SUPER=ALL
DESOBJ=15
SUBCASE 1
DESSUB=1
LABEL=LOAD CONDITION 1
LOAD=300
SUBCASE 2
DESSUB=2
LABEL=LOAD CONDITION 2
LOAD=310
BEGIN BULK
12-3
CASE CONTROL SPECIFICATION (Cont.)
• Design constraints called out for each superelement and load case
ANALYSIS=STATICS
DESOBJ=15
SUBCASE 10
DESSUB=1
SUPER=1,1 $ SUPERELEMENT 1, LOAD SEQ 1
LOAD=300
SUBCASE 20
DESSUB=2
SUPER=1,2 $ SUPERELEMENT 1, LOAD SEQ 2
LOAD=310
SUBCASE 30
DESSUB=3
LOAD=300
SUBCASE 40
DESSUB=4
LOAD=310
BEGIN BULK
12-4
EXAMPLE – TWENTY FIVE BAR TRUSS,
SUPERELEMENT OPTIMIZATION
The figure below shows a twenty-five bar truss structure that has
been partitioned into an upstream superelement 1, and a residual
structure, superelement 0. (Grids 1 and 2 have been defined as
interior to superelement 1.)
2
Superelement 1
6 3
5 4
10
7 Y
Superelement 0
X 8
9
Three-dimensional truss
Symmetric with respect to x-y plane and y-z plane
Weight density = 0.1 lbs/in3
Materials: E = 1.0E7 psi
Two distinct loading conditions
12-5
EXAMPLE: TWENTY FIVE BAR TRUSS,
SUPERELEMENT OPTIMIZATION (Cont.)
1 π Ei
Pb 2
σ b = ------- = ---- – ------------
A
L
A 2
D
A ≅ πDt , where ---- ≅ 100
t
2 2
πDt ( D + t )
I ≅ ----------------------------------
8
12-6
EXAMPLE – TWENTY FIVE BAR TRUSS,
SUPERELEMENT OPTIMIZATION (Cont.)
2
AE π ( 100 + 1 ) D
σ b = – --------∗ ------------------------------- where ---- = 100
2 8∗100 t
L
σ ≥ σb
σ
------ ≤ 1
σb
12-7
EXAMPLE – TWENTY FIVE BAR TRUSS,
SUPERELEMENT OPTIMIZATION (Cont.)
Input data:
ID MSC, D200X3
TIME 10
SOL 200 $ OPTIMIZATION
CEND
super = all
ECHO = UNSORT
OLOAD = ALL
DISP = ALL
SPCFORCE = ALL
ELFORCE = ALL
STRESS = ALL
SPC = 100
ANALYSIS = STATICS $
DESOBJ(MIN) = 15 $ OBJECTIVE FUNCTION DEFINITION
DESSUB = 12 $ CONSTRAINT DEFININITION
SUBCASE 1
LABEL = LOAD CONDITION 1
LOAD = 300
SUBCASE 2
LABEL = LOAD CONDITION 2
LOAD = 310
BEGIN BULK
$
$-------------------------------------------------------------------------
$ ANALYSIS MODEL
$-------------------------------------------------------------------------
$
GRDSET, , , , , , , 456
SESET, 1, 1, 2
MAT1, 1, 10.0E6, , , 0.1, , , , +M1
+M1, 25000., 25000.
SPC1, 100, 123, 7, THRU, 10
GRID, 1 , , -37.5, 0.0, 200.0
GRID, 2 , , 37.5, 0.0, 200.0
GRID, 3 , , -37.5, 37.5, 100.0
GRID, 4 , , 37.5, 37.5, 100.0
GRID, 5 , , 37.5, -37.5, 100.0
GRID, 6 , , -37.5, -37.5, 100.0
GRID, 7 , , -100.0, 100.0, 0.0
GRID, 8 , , 100.0, 100.0, 0.0
GRID, 9 , , 100.0, -100.0, 0.0
GRID, 10, , -100.0, -100.0, 0.0
CROD, 1 , 1, 1, 2
CROD, 2 , 2, 1, 4
12-8
EXAMPLE – TWENTY FIVE BAR TRUSS,
SUPERELEMENT OPTIMIZATION (Cont.)
CROD, 3 , 2, 2, 3
CROD, 4 , 2, 1, 5
CROD, 5 , 2, 2, 6
CROD, 6 , 3, 2, 4
CROD, 7 , 3, 2, 5
CROD, 8 , 3, 1, 3
CROD, 9 , 3, 1, 6
CROD, 10, 4, 3, 6
CROD, 11, 4, 4, 5
CROD, 12, 5, 3, 4
CROD, 13, 5, 5, 6
CROD, 14, 6, 3, 10
CROD, 15, 6, 6, 7
CROD, 16, 6, 4, 9
CROD, 17, 6, 5, 8
CROD, 18, 7, 4, 7
CROD, 19, 7, 3, 8
CROD, 20, 7, 5, 10
CROD, 21, 7, 6, 9
CROD, 22, 8, 6, 10
CROD, 23, 8, 3, 7
CROD, 24, 8, 5, 9
CROD, 25, 8, 4, 8
$
PROD, 1, 1, 2.0, 0.0
PROD, 2, 1, 2.0, 0.0
PROD, 3, 1, 2.0, 0.0
PROD, 4, 1, 2.0, 0.0
PROD, 5, 1, 2.0, 0.0
PROD, 6, 1, 2.0, 0.0
PROD, 7, 1, 2.0, 0.0
PROD, 8, 1, 2.0, 0.0
$
FORCE, 300, 1, , 1.0, 1000., 10000.,-5000.
FORCE, 300, 2, , 1.0, 0., 10000.,-5000.
FORCE, 300, 3, , 1.0, 500., 0., 0.
FORCE, 300, 6, , 1.0, 500., 0., 0.
FORCE, 310, 1, , 1.0, 0., 20000.,-5000.
FORCE, 310, 2, , 1.0, 0., -20000.,-5000.
$
$-----------------------------------------------------------------------------
$ DESIGN MODEL
$-----------------------------------------------------------------------------
$
$...Define the design variables
$DESVAR,ID, LABEL, XINIT, XLB, XUB, DELXV
$
12-9
EXAMPLE – TWENTY FIVE BAR TRUSS,
SUPERELEMENT OPTIMIZATION (Cont.)
12-10
EXAMPLE – TWENTY FIVE BAR TRUSS,
SUPERELEMENT OPTIMIZATION (Cont.)
12-11
EXAMPLE – TWENTY FIVE BAR TRUSS,
SUPERELEMENT OPTIMIZATION (Cont.)
$...Table constants
$DTABLE,LABEL1, VALUE1, LABEL2, VALUE2, LABEL3, VALUE3, LABEL4, VALUE4 +
$+, LABEL5, VALUE5, ...
$
DTABLE, L1, 75.00, L2, 130.50, L3, 106.80,L4, 75.00, +
+, L5, 75.00, L6, 181.14,L7, 181.14, L8, 133.46
$
$...Define the design constraints
$CONSTR,DCID, RID, LALLOW, UALLOW
$
DCONSTR,10, 1 , -40000.,40000.
DCONSTR,10, 2 , -40000.,40000.
DCONSTR,10, 3 , -40000.,40000.
DCONSTR,10, 4 , -40000.,40000.
DCONSTR,10, 5 , -40000.,40000.
DCONSTR,10, 6 , -40000.,40000.
DCONSTR,10, 7 , -40000.,40000.
DCONSTR,10, 8 , -40000.,40000.
DCONSTR,10, 9 , -0.35 ,0.35
DCONSTR,10, 10, -0.35 ,0.35
DCONSTR,10, 11, -0.35 ,0.35
DCONSTR,10, 12, -0.35 ,0.35
DCONSTR,10, 13, -0.35 ,0.35
DCONSTR,10, 14, -0.35 ,0.35
$
DCONSTR,11, 16, -1.0E10, 1.0
DCONSTR,11, 17, -1.0E10, 1.0
DCONSTR,11, 18, -1.0E10, 1.0
DCONSTR,11, 19, -1.0E10, 1.0
DCONSTR,11, 20, -1.0E10, 1.0
DCONSTR,11, 21, -1.0E10, 1.0
DCONSTR,11, 22, -1.0E10, 1.0
DCONSTR,11, 23, -1.0E10, 1.0
$
$...Combine the two constraint sets
$ (equivalent to just putting all into the same set to begin with)
DCONADD,12, 10, 11
$
$...Override optimization parameter defaults:
$
DOPTPRM,IPRINT, 3, DESMAX, 15, DELP, 0.5, p1, 1, +
+, p2, 15
ENDDATA
$.......2.......3.......4.......5.......6.......7.......8.......9.......0
12-12
EXAMPLE – TWENTY FIVE BAR TRUSS, SUPERELEMENT
OPTIMIZATION (Cont.)
Results:
***** NORMAL CONVERGENCE CRITERIA SATISFIED ***** (HARD CONVERGENCE DECISION LOGIC)
**************************************************************************************
CONVERGENCE ACHIEVED BASED ON THE FOLLOWING CRITERIA
(HARD CONVERGENCE DECISION LOGIC)
***************************************************************
S U M M A R Y O F D E S I G N C Y C L E H I S T O R Y
***************************************************************
---------------------------------------------------------------------------------------------------------------------------
INTERNAL | EXTERNAL | |
DV. ID. | DV. ID. | LABEL | 6 : 7 : 8 : 9 : 10 : 11 :
---------------------------------------------------------------------------------------------------------------------------
1 | 1 | X1 | 5.8125E-02 : 2.9063E-02 : 2.8529E-02 :
2 | 2 | X2 | 2.0443E+00 : 2.0427E+00 : 2.0427E+00 :
3 | 3 | X3 | 3.0237E+00 : 3.0059E+00 : 3.0061E+00 :
4 | 4 | X4 | 5.3750E-02 : 2.6875E-02 : 2.5837E-02 :
5 | 5 | X5 | 5.3789E-02 : 5.4376E-02 : 5.3374E-02 :
6 | 6 | X6 | 6.8365E-01 : 6.8306E-01 : 6.8303E-01 :
7 | 7 | X7 | 1.6118E+00 : 1.6216E+00 : 1.6215E+00 :
8 | 8 | X8 | 2.6649E+00 : 2.6690E+00 : 2.6691E+00 :
---------------------------------------------------------------------------------------------------------------------------
*** USER INFORMATION MESSAGE 6464 (DOM12E)
RUN TERMINATED DUE TO HARD CONVERGENCE TO AN OPTIMUM AT CYCLE NUMBER = 8.
SECTION 13
AEROELASTIC OPTIMIZATION
Aeroelastic
Optimization
13-1
STATIC AEROELASTICITY DESIGN CONDITIONS
• Stress/strain
• Displacement
• Internal force
• α ≤ 10.0°
• – 20.0° ≤ δ e ≤ 20.0°
13-2
STATIC AEROELASTICITY SENSITIVITY ANALYSIS
13-3
FLUTTER DESIGN
• Given p:
ω = Imag ( p )
γ = Real ( p ) ⁄ ω
• Velocities
• Densities
• Mach numbers
• Mode numbers
13-4
FLUTTER RESPONSE PROPERTIES
13-5
FLUTTER RESPONSE SENSITIVITY
∂γ 1 ∂p R ∂p l
------ = ---- ---------- – γ --------
∂x ω ∂x ∂x
13-6
AEROELASTIC DESIGN EXAMPLE (HA200A and B)
• Aerodynamic panels
2000 1000
Canard
2100
2007 1007
2103 1131
1100
2115 1119
Aileron
3100
2131
3115
Wing
3103
Rudder
13-7
AEROELASTIC DESIGN EXAMPLE (Cont.)
(HA200A and B)
• Structural model
97 97
(a) Plan View (b) Side View
98 Designed Bars 98
90 90
221 99 121
99
220 120
222 122
211 111
210 110
310
312
13-8
ANALYSIS CONDITIONS FOR EXAMPLE
13-9
DESIGN CONDITIONS FOR EXAMPLE
3
W = ∑ ρi Ai l i
i=1
A 1.5000
I
1 0.1736
= X i
2
I 2.0000
0.4630
J i
and ρ1 = 3
38.490 lb/ft to give a weight of 333 lbs when X1 = 1.0
ρ2 = 3
38.490 lb/ft to give a weight of 667 lbs when X2 = 1.0
ρ3 = 3
5.7735 lb/ft to give a weight of 50 lbs when X3 = 1.0
13-10
DESIGN CONDITIONS FOR EXAMPLE (Cont.)
Subject to:
6
• σt = 50,000 psi = 7.20 × 10 psf limit allowable tensile stress
6
• σc = 40,000 psi = 5.76 × 10 psf limit allowable compressive stress
( u5 ) – ( u5 )
120 100
1.0°, M = 0.9 , Level flight, q = 1200 psf
• ≤
( u5 ) – ( u5 )
100
0.5°, M = 1.2 , Level flight, q = 863 psf
220
• g ≤ 0.00 at V = 1000, 1300 and 1500 ft/s for M = 0.9 and 1.2
13-11
HA200B OBJECTIVE RESULTS
4000
3500
Wing Structural Weight (lbs)
3000
2500
2000
1500
1000
500
0
0 1 2 3 4 5 6 7
Iteration Number
13-12
HA200B DESIGN VARIABLE RESULTS
2.0
= DV10
Design Variable
= DV20
= DV30
1.0
0.0
0 1 2 3 4 5 6 7
Iteration Number
13-13
HA200B FLUTTER RESULTS
1.00
Optimized Design
0.75
Initial Design
0.50
Damping (g)
0.25
0.00
–0.25 VREQ
–0.50
–0.75
–1.00
6.0
5.0
Frequency (Hz)
4.0
3.0
1.0
0.0 1000 1200 1400 1600 1800
Velocity (ft/sec)
13-14
HA200A INPUT FILE
ID MSC, HA200A
$$$$$$$$ HANDBOOK FOR AEROELASTIC ANALYSIS EXAMPLE HA200A $$$$$$$$
$$
$ MODEL DESCRIPTION FULL SPAN 30 DEG FWD SWEPT WING $
$ WITH AILERON, CANARD AND AFT SWEPT $
$ VERTICAL FIN AND RUDDER. $
$ BAR MODEL WITH DUMBBELL MASSES.$
$$
$ SOLUTION QUASI-STEADY AEROELASTIC ANALYSIS$
$ AND UNSTEADY FLUTTER ANALYSIS USING$
$ DOUBLET-LATTICE METHOD$
$ AERODYNAMICS AT MACH NO. 0.9.$
$$
$ OUTPUT STANDARD AEROELASTIC OUTPUT PLUS$
$ A TABLE IDENTIFYING RESPONSES$
$ FOR WHICH SENSITIVITY RESULTS ARE$
$ AVAILABLE FOLLOWED BY A MATRIX OF$
$ SENSITIVITY VALUES.$
$ $
$$
$$$$$$$$ $$$$$$$$
TIME 30 $ CPU TIME IN MINUTES
SOL 200 $ OPTIMIZATION WITH AEROELASTICITY
CEND
TITLE = EXAMPLE HA200A: 30 DEG FWD SWEPT WING WITH CANARD AND FIN
SUBTI = DEMONSTRATION OF AEROELASTIC SENSITIVITY ANALYSIS
ECHO = BOTH
SPC = 1
DESOBJ = 10
DISP = ALL $
STRESS = ALL $
FORCE = ALL $
AEROF = ALL $
APRES = ALL $
SUBCASE 1
LABEL = SUBSONIC SYMMETRIC PULLOUT
ANALYSIS = SAERO
DESSUB = 1
TRIM = 1 $
SUBCASE 2
LABEL = SUPERSONIC SYMMETRIC PULLOUT
ANALYSIS = SAERO
DESSUB = 2
TRIM = 2 $
SUBCASE 3
LABEL = HIGH SPEED ROLLING PULLOUT
ANALYSIS = SAERO
DESSUB = 200
TRIM = 3 $
13-15
HA200A INPUT FILE (Cont.)
SUBCASE 4
LABEL = HIGH SPEED PULLUP WITH ABRUPT ROLL
ANALYSIS = SAERO
DESSUB = 200
TRIM = 4 $
SUBCASE 5
LABEL = SUBSONIC ENTRY INTO SNAP ROLL
ANALYSIS = SAERO
DESSUB = 200
TRIM = 5 $
SUBCASE 6
LABEL = SUBSONIC FLUTTER ANALYSIS
ANALYSIS = FLUTTER
SET 10 = 1,THRU,100000
PARAM OPPHIPA,1
DISP = 10
STRESS = NONE $
FORCE = NONE $
AEROF = NONE $
APRES = NONE $
DESSUB = 6
METHOD = 20
FMETHOD = 30
SUBCASE 7
LABEL = SUPERSONIC FLUTTER ANALYSIS
ANALYSIS = FLUTTER
DISP = NONE $
STRESS = NONE $
FORCE = NONE $
AEROF = NONE $
APRES = NONE $
DESSUB = 7
METHOD = 20
FMETHOD = 40
BEGIN BULK
13-16
DESIGN MODEL AND BULK DATA INPUT
$
$ * * * *
$
$ * * * OPTIMIZATION * * *
$
$
$ * * *
$
$ * THE DESIGN MODEL *
$
$ DEFINITION OF THE DESIGN VARIABLES
$
$ THE DESVAR ENTRY DEFINES A DESIGN VARIABLE FOR DESIGN
$ OPTIMIZATION. LISTED ARE A UNIQUE DESIGN VARIABLE ID
$ NUMBER, A USER SUPPLIED NAME FOR PRINTING PURPOSES, AN
$ INITIAL VALUE, A LOWER BOUND, AND AN UPPER BOUND.
$
$DESVAR ID LABEL XINIT XLB XUB
DESVAR 10 PBAR101 1.0 0.001 100.0
DESVAR 20 PBAR102 1.0 0.001 100.0
DESVAR 30 PBAR103 0.1 0.001 100.0
$
$ RELATIONSHIP OF DESIGN VARIABLES TO ANALYSIS MODEL PROPERTIES
$
$ THE DVPREL1 ENTRY EXPRESSES AN AN ANALYSIS MODEL PROPERTY AS
$ A LINEAR FUNCTION OF DESIGN VARIABLES. IT LISTS A UNIQUE ID,
$ AN ANALYSIS MODEL ENTRY TYPE ID STRING, A PROPERTY ENTRY ID,
$ THE FIELD POSITION OF THE PROPERTY ENTRY OR WORD POSITION IN
$ THE ELEMENT PROPERTY TABLE OF THE ANALYSIS MODEL, THE MINIMUM
$ AND MAXIMUM VALUES ALLOWED FOR THIS PROPERTY DURING OPTIMIZATION,
$ A CONSTANT TERM OF RELATION, A DESIGN VARIABLE ENTRY (DESVAR)
$ ID, AND A COEFFICIENT OF LINEAR RELATION. THE EQUATION IS PI
$ = CO + CJXJ AND THE CONTINUATION ENTRY CAN BE USED TO LIST
$ MORE THAN ONE CJ. IN THIS CASE, THERE IS ONLY ONE CJ FOR EACH
$ DVPREL1 ENTRY.
$
$DVPREL1ID TYPE PID FID PMIN PMAX C0 XXXX DVPREL1
DVPREL1 1014 PBAR 101 4
$ DVID1 COEF1 DVID2 COEF2 DVID3 COEF3 ....
10 1.5
DVPREL1 1015 PBAR 101 5
10 0.173611
DVPREL1 1016 PBAR 101 6
10 2.0
DVPREL1 1017 PBAR 101 7
10 0.462963
DVPREL1 1024 PBAR 102 4
20 1.5
DVPREL1 1025 PBAR 102 5
20 0.173611
13-17
DESIGN MODEL AND BULK DATA INPUT (Cont.)
13-18
DESIGN MODEL AND BULK DATA INPUT (Cont.)
$ STRESS CONSTRAINTS
$
$ THE DCONSTR ENTRY DEFINES DESIGN CONSTRAINTS. LISTED ARE
$ A CONSTRAINT SET ID, THE DRESPI ENTRY ID AND THE LOWER AND
$ UPPER BOUND IMPOSED ON THIS RESPONSE QUANTITY. THE FOLLOWING
$ TWO ENTRIES IMPOSE LIMITS ON THE ALLOWABLE STRESSES IN THE
$ BARS IN UNITS OF POUNDS/FOOT**2.
$
$DCONSTR DCID RID LALLOW UALLOW
DCONSTR 200 1001 -5.76+6 7.20+6
DCONSTR 200 1002 -5.76+6 7.20+6
$
$
$ TIP TWIST RESPONSES
$
DRESP1 101 RTIPROT DISP 5 120
DRESP1 201 LTIPROT DISP 5 220
DRESP1 100 RTROT DISP 5 100
$
$ SECOND LEVEL RESPONSES FOR TIP TWIST
$
$ THE DRESP2 ENTRY DEFINES THE INPUT ARGUMENTS TO USER-SUPPLIED
$ EQUATIONS. THESE SECOND-LEVEL RESPONSES CAN THEN EITHER BE
$ USED AS CONSTRAINTS OR AS AN OBJECTIVE FUNCTION. INPUT MAY
$ CONSIST OF DESIGN VARIABLES (DESVAR), FIRST-LEVEL (DRESP1)
$ RESPONSES, TABLE CONSTANTS (DTABLE), AND GRID COORDINATES
$ (DVGRID). IT LISTS AN ID, A USER DEFINED LABEL, THE DEQATN
$ ENTRY ID, A REGION IDENTIFIER FOR CONSTRAINT SCREENING, A
$ STRING INDICATING DESVAR ID NUMBERS, A DESVAR ID, A STRING
$ INDICATING THAT THE LABELS FOR THE CONSTANTS IN A DTABLE ENTRY
$ FOLLOW, THE LABELS OF CONSTANTS IN THE DTABLE INPUT,
$ A STRING INDICATING DRESP1 ID NUMBERS, DRESP1 IDS, A STRING
$ SIGNIFYING THAT THE IDS AND DIRECTIONS FOLLOWING ARE NODE
$ NUMBERS AND CARTESIAN DIRECTION COMPONENTS, NODE NUMBERS,
$ AND CARTESIAN DIRECTIONS.
$
$DRESP2 ID LABEL EQID REGION XXXX XXXX XXXX XXXX
DRESP2 5 RHSTWIST5 +DR25
$ DRESP1 NR1 NR2 NR3 ETC
+DR25 DRESP1 101 100
$
DRESP2 6 LHSTWIST5 +DR26
+DR26 DRESP1 201 100
$
$ EQUATION DEFINING SECOND LEVEL RESPONSE FOR TIP TWIST
$
$ THE DEQATN ENTRY DEFINES THE EQUATION(S) USED IN THE DESIGN
$ PROCESS. IT LISTS A UNIQUE EQUATION ID, AND THE EQUATION IS
$ WRITTEN IN FORTRAN LIKE SYNTAX FOLLOWING THE RULES IN DEFINING
$ DMAP ASSIGNMENTS AND FUNCTIONS.
13-19
DESIGN MODEL AND BULK DATA INPUT (Cont.)
13-20
DESIGN MODEL AND BULK DATA INPUT (Cont.)
$
$ SELECTION OF FLUTTER MODES FOR OPTIMIZATION
$
SET1 88 4 THRU 7
$
$ SELECTION OF VELOCITIES FOR IMPOSING FLUTTER CONSTRAINTS
$
FLFACT 4 1000.0 1300.0 1500.0
$
$ EQUATION FOR SECOND LEVEL FLUTTER RESPONSE TO CONSTRAIN SYSTEM DAMPING
$
DRESP2 4 GDAMP 4 +DR24
+DR24 DRESP1 1
DEQATN 4 F(A) = (A - 0.03)/0.1
$
$ CONSTRAINT ON AEROELASTIC SYSTEM DAMPING
$
DCONSTR 6 4 -1.0+20 -0.3
$
$
$ RESPONSE FOR SUPERSONIC FLUTTER
$
$DRESP1 ID LABEL FLUTTER XXXX REGION SID XXXX ID_MODE +DR
DRESP1 11 FLUTTER FLUTTER 89 +DR11
$ ID_DENS ID_MACH ID_VEL
+DR11 11 12 14
$
SET1 89 4 THRU 7
$
FLFACT 14 1000.0 1300.0 1500.0
$
DRESP2 21 GDAMP 4 +DR21
+DR21 DRESP1 11
$
DCONSTR 7 21 -1.0+20 -0.3
$
$ * * *
$
$ * OPTIMIZATION CONTROL PARAMETERS *
$
$ THE DOPTPRM ENTRY IS USED TO OVERRIDE DEFAULT OPTIMIZATION
$ PARAMETERS. THE DOPTPRM ENTRY IS USED
$ HERE TO SET THE MAXIMUM NUMBER OF DESIGN CYCLES TO 25 AND TO
$ REQUEST DETAILED RESULTS FOLLOWING EACH OPTIMIZATION CYCLE.
$
$ PARAM1 VAL1 PARAM2 VAL2 ETC
DOPTPRM DESMAX 25 P1 2 P2 15 DELB 0.01
$
13-21
DESIGN MODEL AND BULK DATA INPUT (Cont.)
13-22
HA200A SELECTED RESULTS
EXAMPLE HA200A: 30 DEG FWD SWEPT WING WITH CANARD AND FIN PAGE 31
DEMONSTRATION OF AEROELASTIC SENSITIVITY ANALYSIS
----- COMPARISON BETWEEN INPUT PROPERTY VALUES FROM ANALYSIS AND DESIGN MODELS -----
---------------------------------------------------------------------------------------------------------------
PROPERTY PROPERTY FIELD ANALYSIS DESIGN LOWER UPPER DIFFERENCE
TYPE ID ID VALUE VALUE BOUND BOUND FLAG
---------------------------------------------------------------------------------------------------------------
PBAR 101 4 1.500000E+00 1.500000E+00 1.000000E-03 1.000000E+20 NONE
PBAR 101 5 1.736110E-01 1.736110E-01 1.000000E-03 1.000000E+20 NONE
PBAR 101 6 2.000000E+00 2.000000E+00 1.000000E-03 1.000000E+20 NONE
PBAR 101 7 4.629630E-01 4.629630E-01 1.000000E-03 1.000000E+20 NONE
PBAR 102 4 1.500000E+00 1.500000E+00 1.000000E-03 1.000000E+20 NONE
PBAR 102 5 1.736110E-01 1.736110E-01 1.000000E-03 1.000000E+20 NONE
PBAR 102 6 2.000000E+00 2.000000E+00 1.000000E-03 1.000000E+20 NONE
13-23
----------------------------------------------------------
COLUMN DRESP2 SUB FREQ/
NO. ENTRY ID CASE TIME
----------------------------------------------------------
86 5 1
87 6 1
88 2401 1
89 5 2
90 6 2
91 2401 2
92 4 6
93 4 6
94 4 6
95 4 6
96 4 6
13-28
97 4 6
98 4 6
99 4 6
100 4 6
101 4 6
102 4 6
103 4 6
104 21 7
105 21 7
106 21 7
107 21 7
108 21 7
109 21 7
110 21 7
111 21 7
112 21 7
113 21 7
114 21 7
115 21 7
^^^ DMAP INFORMATION MESSAGE 9029 (DESOPT) - DESIGN SENSITIVITY COEFFICIENT MATRIX FOR DIRECT
AND SYNTHETIC RESPONSES - GRADIENTS OF RESPONSES WITH RESPECT TO INDEPENDENT DESIGN VARIABLES
HA200A SELECTED RESULTS (Cont.)
MATRIX DSCM2 (GINO NAME 101 ) IS A REAL 115 COLUMN X 3 ROW RECTANG MATRIX.
COLUMN 1 ROWS 1 THRU 3 --------------------------------------------------
ROW
1) 6.6667E+02 1.3333E+03 5.0000E+01
COLUMN 2 ROWS 1 THRU 3 --------------------------------------------------
ROW
1) 2.3115E-04 5.5879E-06 5.3085E-05
COLUMN 3 ROWS 1 THRU 3 --------------------------------------------------
ROW
1) -5.4311E-03 -4.8257E-03 8.5682E-05
COLUMN 4 ROWS 1 THRU 3 --------------------------------------------------
ROW
1) -5.4313E-03 -4.8257E-03 8.5682E-05
COLUMN 5 ROWS 1 THRU 3 --------------------------------------------------
ROW
13-29
---------------------------------------------------------------------------------------------------------------
OBJECTIVE FROM OBJECTIVE FROM FRACTIONAL ERROR MAXIMUM VALUE
CYCLE APPROXIMATE EXACT OF OF
NUMBER OPTIMIZATION ANALYSIS APPROXIMATION CONSTRAINT
---------------------------------------------------------------------------------------------------------------
INITIAL 2.005000E+03 1.460522E+01
---------------------------------------------------------------------------------------------------------------
DESIGN VARIABLE HISTORY
----------------------------------------------------------------------------------------------------------------------------------
INTERNAL | EXTERNAL | |
DV. ID. | DV. ID. | LABEL | INITIAL : 1 : 2 : 3 : 4 : 5 :
----------------------------------------------------------------------------------------------------------------------------------
1 | 10 | PBAR101 | 1.0000E+00 :
2 | 20 | PBAR102 | 1.0000E+00 :
3 | 30 | PBAR103 | 1.0000E-01 :
----------------------------------------------------------------------------------------------------------------------------------
*** USER INFORMATION MESSAGE 6464 (DOM12E)
RUN TERMINATED DUE TO PARAMETER OPTEXIT = 4.
HA200B SELECTED RESULTS
EXAMPLE HA200B: 30 DEG FWD SWEPT WING WITH CANARD AND FIN PAGE 159
DEMONSTRATION OF AEROELASTIC OPTIMIZATION
SUPERSONIC FLUTTER ANALYSIS SUBCASE 7
THIS IS THE FIRST ANALYSIS - NO CONVERGENCE CHECK
-------------------------------------------------------
MAXIMUM VALUE OF CONSTRAINTS : 1.4605E+01
-------------------------------------------------------
********************************************************************
* *
* *
* D E S I G N O P T I M I Z A T I O N *
* *
* *
********************************************************************
13-31
*******************************************
* *
* D E S I G N C Y C L E 1 *
* *
*******************************************
---------------------------------------------------------------------------------------------------------------------------
INTERNAL DESVAR LOWER INPUT OUTPUT UPPER
ID ID LABEL BOUND VALUE VALUE BOUND
---------------------------------------------------------------------------------------------------------------------------
1 10 PBAR101 1.0000E-03 1.0000E+00 2.0000E+00 1.0000E+02
2 20 PBAR102 1.0000E-03 1.0000E+00 1.6390E+00 1.0000E+02
3 30 PBAR103 1.0000E-03 1.0000E-01 2.0000E-01 1.0000E+02
---------------------------------------------------------------------------------------------------------------------------
PROPERTY PROPERTY FIELD TYPE OF LOWER INPUT OUTPUT UPPER
TYPE ID ID PROPERTY BOUND VALUE VALUE BOUND
13-32
---------------------------------------------------------------------------------------------------------------------------
PBAR 101 4 DVPREL1 1.0000E-03 1.5000E+00 3.0000E+00 1.0000E+20
PBAR 101 5 DVPREL1 1.0000E-03 1.7361E-01 3.4722E-01 1.0000E+20
PBAR 101 6 DVPREL1 1.0000E-03 2.0000E+00 4.0000E+00 1.0000E+20
PBAR 101 7 DVPREL1 1.0000E-03 4.6296E-01 9.2593E-01 1.0000E+20
PBAR 102 4 DVPREL1 1.0000E-03 1.5000E+00 2.4586E+00 1.0000E+20
PBAR 102 5 DVPREL1 1.0000E-03 1.7361E-01 2.8455E-01 1.0000E+20
PBAR 102 6 DVPREL1 1.0000E-03 2.0000E+00 3.2781E+00 1.0000E+20
PBAR 102 7 DVPREL1 1.0000E-03 4.6296E-01 7.5881E-01 1.0000E+20
PBAR 103 4 DVPREL1 1.0000E-03 1.5000E-01 3.0000E-01 1.0000E+20
PBAR 103 5 DVPREL1 1.0000E-03 1.7361E-02 3.4722E-02 1.0000E+20
PBAR 103 6 DVPREL1 1.0000E-03 2.0000E-01 4.0000E-01 1.0000E+20
PBAR 103 7 DVPREL1 1.0000E-03 4.6296E-02 9.2593E-02 1.0000E+20
HA200B SELECTED RESULTS (Cont.)
---------------------------------------------------------------------------------------------------------------------------
INTERNAL INTERNAL
INTERNAL DCONSTR RESPONSE RESPONSE L/U REGION SUBCASE INPUT OUTPUT
ID ID ID TYPE FLAG ID ID VALUE VALUE
---------------------------------------------------------------------------------------------------------------------------
1 50 2 EQUA UPPER 6 1 -3.9543E-01 -6.6547E-01
2 50 1 EQUA UPPER 5 1 -3.9543E-01 -6.6440E-01
3 50 3 EQUA LOWER 2401 1 1.7848E-02 -3.4503E-02
4 60 4 EQUA LOWER 2401 2 1.1241E-02 -3.9465E-02
5 6 5 EQUA UPPER 4 6 -4.3507E-01 -3.6965E-01
6 6 6 EQUA UPPER 4 6 1.4605E+01** 1.3839E+00**
13-33
---------------------------------------------------------------------------------------------------------------------------
INTERNAL PROPERTY FIELD L/U CYCLE INPUT OUTPUT
ID ID FLAG FLAG LIMIT VALUE VALUE
---------------------------------------------------------------------------------------------------------------------------
9 101 4 LOWER 1.2000E+00 -2.5000E-01 -1.5000E+00
10 101 5 LOWER 1.3889E-01 -2.5000E-01 -1.5000E+00
11 101 6 LOWER 1.6000E+00 -2.5000E-01 -1.5000E+00
12 101 7 LOWER 3.7037E-01 -2.5000E-01 -1.5000E+00
13 102 4 LOWER 1.2000E+00 -2.5000E-01 -1.0488E+00
HA200B SELECTED RESULTS (Cont.)
---------------------------------------------------------------------------------------------------------------------------
INTERNAL PROPERTY FIELD L/U CYCLE INPUT OUTPUT
ID ID FLAG FLAG LIMIT VALUE VALUE
---------------------------------------------------------------------------------------------------------------------------
14 102 5 LOWER 1.3889E-01 -2.5000E-01 -1.0488E+00
15 102 6 LOWER 1.6000E+00 -2.5000E-01 -1.0488E+00
16 102 7 LOWER 3.7037E-01 -2.5000E-01 -1.0488E+00
17 103 4 LOWER 1.2000E-01 -2.5000E-01 -1.5000E+00
18 103 5 LOWER 7.3611E-03 -1.3585E+00 -3.7170E+00
19 103 6 LOWER 1.6000E-01 -2.5000E-01 -1.5000E+00
20 103 7 LOWER 3.6296E-02 -2.7551E-01 -1.5510E+00
21 101 4 UPPER 1.8000E+00 -1.6667E-01 6.6667E-01
22 101 5 UPPER 2.0833E-01 -1.6667E-01 6.6667E-01
13-34
----------------------------------------------------------------------
| R E S P O N S E S IN D E S I G N M O D E L |
----------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------
INTERNAL DRESP1 RESPONSE LOWER INPUT OUTPUT UPPER
ID ID LABEL BOUND VALUE VALUE BOUND
---------------------------------------------------------------------------------------------------------------------------
1 10 WEIGHT N/A 2.0050E+03 3.5277E+03 N/A
D E S I G N C Y C L E = 1 S U B C A S E = 1
13-35
---------------------------------------------------------------------------------------------------------------------------
INTERNAL DRESP1 RESPONSE GRID COMPONENT LOWER INPUT OUTPUT UPPER
ID ID LABEL ID NO. BOUND VALUE VALUE BOUND
---------------------------------------------------------------------------------------------------------------------------
2 100 RTROT 100 5 N/A 1.1291E-02 1.1403E-02 N/A
3 101 RTIPROT 120 5 N/A 2.1841E-02 1.7259E-02 N/A
4 201 LTIPROT 220 5 N/A 2.1841E-02 1.7240E-02 N/A
---------------------------------------------------------------------------------------------------------------------------
RESTRAINED / AESTAT /
INTERNAL DRESP1 RESPONSE UNRESTRAINED COMPONENT AESURF LOWER INPUT OUTPUT UPPER
ID ID LABEL FLAG NO. ID BOUND VALUE VALUE BOUND
---------------------------------------------------------------------------------------------------------------------------
5 1401 CLDELTA 0 4 517 N/A 2.6340E-01 2.6922E-01 N/A
6 1402 CLP 0 4 513 N/A -4.4697E-01 -4.3374E-01 N/A
HA200B SELECTED RESULTS (Cont.)
D E S I G N C Y C L E = 1 S U B C A S E = 2
---------------------------------------------------------------------------------------------------------------------------
RESTRAINED / AESTAT /
INTERNAL DRESP1 RESPONSE UNRESTRAINED COMPONENT AESURF LOWER INPUT OUTPUT UPPER
ID ID LABEL FLAG NO. ID BOUND VALUE VALUE BOUND
---------------------------------------------------------------------------------------------------------------------------
7 1401 CLDELTA 0 4 517 N/A 2.1655E-01 2.2531E-01 N/A
8 1402 CLP 0 4 513 N/A -5.0933E-01 -5.0408E-01 N/A
D E S I G N C Y C L E = 1 S U B C A S E = 6
------------------------------------------------------------------------------------------------------------------------
INTERNAL DRESP1 RESPONSE MODE MACH LOWER INPUT OUTPUT UPPER
ID ID LABEL NO. VELOCITY NO. BOUND VALUE VALUE BOUND
---------------------------------------------------------------------------------------------------------------------------
9 1 FLUTTER 6 1.0000E+03 9.0000E-01 N/A -1.3052E-02 -1.1089E-02 N/A
10 1 FLUTTER 6 1.5000E+03 9.0000E-01 N/A 4.3816E-01 4.1516E-02 N/A
D E S I G N C Y C L E = 1 S U B C A S E = 7
---------------------------------------------------------------------------------------------------------------------------
INTERNAL DRESP1 RESPONSE MODE MACH LOWER INPUT OUTPUT UPPER
ID ID LABEL NO. VELOCITY NO. BOUND VALUE VALUE BOUND
---------------------------------------------------------------------------------------------------------------------------
11 11 FLUTTER 5 1.5000E+03 1.2000E+00 N/A 2.4369E-01 -2.6562E-01 N/A
12 11 FLUTTER 6 1.0000E+03 1.2000E+00 N/A -1.4145E-02 -1.2566E-02 N/A
HA200B SELECTED RESULTS (Cont.)
---------------------------------------------------------------------------------------------------------------------------
INTERNAL DRESP2 RESPONSE EQUATION LOWER INPUT OUTPUT UPPER
ID ID LABEL ID BOUND VALUE VALUE BOUND
---------------------------------------------------------------------------------------------------------------------------
1 5 RHSTWIST 5 N/A 1.0550E-02 5.8562E-03 1.7450E-02
2 6 LHSTWIST 5 N/A 1.0550E-02 5.8376E-03 1.7450E-02
3 2401 ROLLEFF 103 6.0000E-01 5.8929E-01 6.2070E-01 N/A
4 2401 ROLLEFF 103 4.3000E-01 4.2517E-01 4.4697E-01 N/A
5 4 GDAMP 4 N/A -4.3052E-01 -4.1089E-01 -3.0000E-01
6 4 GDAMP 4 N/A 4.0816E+00 1.1516E-01 -3.0000E-01
7 21 GDAMP 4 N/A 2.1369E+00 -2.9562E+00 -3.0000E-01
8 21 GDAMP 4 N/A -4.4145E-01 -4.2566E-01 -3.0000E-01
****************************************************************************************
13-37
INSPECTION OF CONVERGENCE DATA FOR THE OPTIMAL DESIGN WITH RESPECT TO APPROXIMATE MODELS
(SOFT CONVERGENCE DECISION LOGIC)
****************************************************************************************
***** NORMAL CONVERGENCE CRITERIA SATISFIED ***** (HARD CONVERGENCE DECISION LOGIC)
**************************************************************************************
CONVERGENCE ACHIEVED BASED ON THE FOLLOWING CRITERIA
(HARD CONVERGENCE DECISION LOGIC)
**************************************************************************************
HA200B SELECTED RESULTS (Cont.)
***************************************************************
S U M M A R Y O F D E S I G N C Y C L E H I S T O R Y
***************************************************************
AUXCASE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-4
AUXMODEL . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-5
DESGLB . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-6
DESOBJ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-7
DESSUB . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-8
MODTRAK. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-9
BNDGRID . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-11
DCONADD. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-12
DCONSTR . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-13
DEQATN . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-15
DESVAR . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-19
DLINK . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-20
DOPTPRM . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-21
DRESP1. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-25
DRESP2. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-31
DSAPRT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-34
APPENDIX A
DSCREEN . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-37
DTABLE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-39
DVBSHAP . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-40
DVGRID . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-41
DVPREL1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-43
DVPREL2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-46
DVSHAP . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-48
MODTRAK. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-49
DESIGN MODELING INPUT DATA
• Design variables
A-1
DESIGN MODELING INPUT DATA (Cont.)
• DTABLE Constants
A-2
CASE CONTROL COMMANDS
A-3
AUXCASE Auxiliary Model Case Control Delimiter
Format:
AUXCASE
Examples:
AUXCAS
AUXC
Remarks:
1. AUXCASE indicates the beginning of Case Control commands for an auxiliary model.
AUXCASE must follow the primary model Case Control commands.
2. All Case Control commands following this entry are applicable until the next AUXCASE or
BEGIN BULK command. Commands from preceding Case Control Sections are ignored.
3. Each auxiliary model Case Control must be delimited with the AUXCASE command.
4. The AUXMODEL command is used to associate the auxiliary model Case Control with a
particular auxiliary model.
5.
6.
A-4
Auxiliary Model Identification Number AUXMODEL
AUXMODEL an auxiliary model for generation of boundary shapes in Auxiliary
References shape optimization.
Model Identification Number
Format:
AUXMODEL = n
Examples:
AUXMODEL = 4
AUXM = 4
Describer Meaning
n Auxiliary model identification number. (Integer > 0)
Remarks:
1. AUXMODEL references a particular auxiliary model for analysis and may only be specified
in the auxiliary model Case Control Section.
2. See the BEGIN BULK command for the Bulk Data definition of an auxiliary model.
A-5
DESGLB Request Design Constraints at the Global Level
Format:
DESGLB = n
Examples:
DESGLB = 10
DESG = 25
Describer Meaning
n Set identification of a DCONSTR or DCONADD Bulk Data entry identification
number. (Integer > 0)
Remarks:
1. A DESGLB command is optional and invokes constraints that are to be applied indepen-
dent of a particular subcase. These constraints could be based on responses that are
independent of subcases (e.g., WEIGHT or VOLUME).
2. The DESGLB command can be used to invoke constraints that are not a function of
DRESP1 entries; e.g., DRESP2 responses that are not functions of DRESP1 responses
are subcase independent.
A-6
Design Objective DESOBJ
DESOBJ the DRESP1 or DRESP2 entry to be used as the design objective.
Selects Design Objective
Format:
MAX
DESOBJ ( ) = N
MIN
Examples:
DESOBJ = 10
DESO = 25
Describer Meaning
MIN Specifies that the objective is to be minimized.
Remarks:
1. A DESOBJ command is required for a design optimization task and is optional for a sensi-
tivity task.
2. If the DESOBJ command is specified within a SUBCASE, the identified DRESPi Bulk Data
entry use a response only from that subcase. If DESOBJ appears above all SUBCASE
commands and there are multiple subcases, it uses a global response.
A-7
DESSUB Design Constraints Request at the Subcase Level
Format:
DESSUB = n
Examples:
DESSUB = 10
DESS = 25
Describer Meaning
n Set identification of a DCONSTR or DCONADD Bulk Data entry identification
number. (Integer ≥ 0)
Remark:
1. A DESSUB command is required for every subcase for which constraints are to be
applied.
A-8
MODTRAK Mode Tracking Request
Format:
MODTRAK = n
Example:
MODTRAK = 100
Describer Meaning
n Set identification of a MODTRAK Bulk Data entry. (Integer > 0)
Remark:
1. Selection of a MODTRAK Bulk Data entry with the MODTRAK Case Control command
activates mode tracking for the current subcase. This request is limited to normal modes
subcases (ANALYSIS = MODES) in design optimization (SOL 200).
A-9
BULK DATA ENTRIES
A-10
Boundary Grid Points BNDGRID
BNDGRIDa list of grid point identification numbers on design boundaries or surfaces
Specifies for shape
Boundary opti-
Grid Points
mization (SOL 200).
Format:
1 2 3 4 5 6 7 8 9 10
BNDGRID C GP1 GP2 GP3 GP4 GP5 GP6 GP7
GP8 -etc.-
Example:
BNDGRID 123 41 42 43 44 45 46 47
49
Field Contents
C Component number (any unique combination of integers 1 through 6 with no
embedded blanks). See Remark 1.
GPi Shape boundary grid point identification number. (0 < Integer < 1000000; For
THRU option, GP1< GP2)
Remarks:
1. C specifies the components for the listed grid points for which boundary motion is
prescribed.
2. Multiple BNDGRID entries may be used to specify the shape boundary grid point identifi-
cation numbers.
3. Both fixed and free shape boundary grid point identification numbers are listed on this
entry.
4. The degrees of freedom specified on BNDGRID entries must be sufficient to statically con-
strain the model.
5. Degrees of freedom specified on this entry form members of the mutually exclusive s-set.
They may not be specified on other entries that define mutually exclusive sets. See the
MSC/NASTRAN Quick Reference Guide, Appendix B for a list of these entries.
A-11
DCONADD Design Constraint Set Combination
Format:
1 2 3 4 5 6 7 8 9 10
DCONADD DCID DC1 DC2 DC3 DC4 DC5 DC6 DC7
DC8 -etc.-
Example:
DCONADD 10 4 12
Field Contents
DCID Design constraint set identification number. (Integer > 0)
Remarks:
1. The DCONADD entry is selected by a DESSUB or DESGLB Case Control command.
A-12
Design Constraints DCONSTR
DCONSTR
Define design constraints. Design Constraints
Format:
1 2 3 4 5 6 7 8 9 10
DCONSTR DCID RID LALLOW UALLOW
Example:
DCONSTR 10 4 1.25
Field Contents
DCID Design constraint set identification number. (Integer > 0)
Remarks:
1. The DCONSTR entry may be selected in the Case Control Section by the DESSUB or
DESGLB command.
3. For a given DCID, the associated RID can be referenced only once.
4. The units of LALLOW and UALLOW must be consistent with the referenced response de-
fined on the DRESPi entry. If RID refers to an “EIGN” response, then the imposed bounds
must be expressed in units of eigenvalue, (radian/time)2. If RID refers to a “FREQ” re-
sponse, then the imposed bounds must be expressed in cycles/time.
5. LALLOW and UALLOW are unrelated to the stress limits specified on the MATi entry.
LALLOW – r
g = -------------------------------- for lower bound constraints
GNORM
r – UALLOW
g = --------------------------------- for upper bound constraints
GNORM
where r is the response defined on the DRESPi entry and if the RID refers to a “FREQ”
response, then the imposed bounds must be expressed in cycles/time.
(Continued)
A-13
DCONSTR Design Constraints
7. As Remark 6 indicates, small values of UALLOW and LALLOW require special processing
and should be avoided. Bounds of exactly zero are particularly troublesome. This can be
avoided by using a DRESP2 entry that offsets the constrained response from zero.
A-14
Design Equation Definition DEQATN
DEQATN one or more equations for use in design sensitivity or p-element analysis.
Defines Design Equation Definition
Format:
1 2 3 4 5 6 7 8 9 10
DEQATN EQID EQUATION
EQUATION (Cont.)
Example:
+ A ∗ ∗ 2/(B – C); F = A + B – F1 ∗ D
Field Contents
EQID Unique equation identification number. (Integer > 0)
Remarks:
1. EQUATION is a single equation or a set of nested equations and is specified in fields 3
through 9 on the first entry and may be continued on fields 2 through 9 on the continuation
entries. On the continuation entries, no commas can appear in columns 1 through 8. All
data in fields 2 through 9 must be specified in columns 9 through 72. The large-field for-
mat is not allowed.
(Continued)
A-15
DEQATN Design Equation Definition
2
3 A
1 = A + B ∗ C – ( D + 10 ) + sin ( PI ( 1 ) ∗ R ) + --------------
B–C
F = A + B – F1 ∗ D
2. EQUATION may contain embedded blanks. EQUATION must contain less than 12,500
nonblank characters. This is equivalent to approximately 195 continuation entries.
3. The syntax of the expressions follows FORTRAN language standards. The allowable
arithmetic operations are shown in Table 1 in the order of execution precedence. Paren-
thesis are used to change the order of precedence. Operations within parentheses are
performed first with the usual order of precedence being maintained within the
parentheses.
4. The expressions may contain intrinsic functions. Table 2 contains the format and descrip-
tions of functions that may appear in the expressions. The use of functions that may be
discontinuous must be used with caution because they can cause discontinuous deriva-
tives. These are ABS, DIM, MAX, MIN, and MOD. For examples and further details see
the MSC/NASTRAN DMAP Module Dictionary.
(Continued)
A-16
Design Equation Definition DEQATN
Table 2. DEQATN Entry Functions.
EXP(x) exponential ex
LOG(x) natural logarithm loge x
(Continued)
A-17
DEQATN Design Equation Definition
b. DRESP2 entry, then xi represents the DVIDj, LABLk, NRm, and Gp fields in that
order.
7. The DMAP logical operators NOT, AND, OR < XOR, and XQV cannot be used as Xi
names.
8. Input errors on the DEQATN entry often result in poor messages. Substituting a “[” for a
parenthesis or violating the restriction against large field format are examples. Known
messages are UFM 215, SFM 233 and UFM 5199. If any of these messages are encoun-
tered then review the DEQATN entry input.
A-18
Design Variable DESVAR
Defines
DESVAR a design variable for design optimization. Design Variable
Format:
1 2 3 4 5 6 7 8 9 10
DESVAR ID LABEL XINIT XLB XUB DELXV
Example:
Field Contents
ID Unique design variable identification number. (Integer > 0)
DELXV Fractional change allowed for the design variable during approximate
optimization. (Real > 0.0, for Default see Remark 2.)
Remarks:
1. DELXV can be used to control the change in the design variable during one optimization
cycle.
2. If DELXV is blank, the default is taken from the specification of the DELX parameter on
the DOPTPRM entry. If DELX is not specified, then the default is 1.0.
A-19
DLINK Multiple Design Variable Linking
Relates
DLINK one design variable to one or more other design variables. Multiple Design Variable Linking
Format:
1 2 3 4 5 6 7 8 9 10
DLINK ID DDVID C0 CMULT IDV1 C1 IDV2 C2
IDV3 C3 -etc.-
Example:
8 7.0
Field Contents
ID Unique identification number. (Integer > 0)
Remarks:
1. DLINK defines the relationship
3. CMULT provides a simple means of scaling the Ci. For example if Ci = 1/7, 2/7, 4/7, etc.
is desired, then CMULT = 1/7 and Ci = 1, 2, 4, etc., may be input.
4. An independent IDVi must not occur on the same DLINK entry more than once.
A-20
Design Optimization Parameters DOPTPRM
DOPTPRM default values of parameters used in design optimization.
Overrides Design Optimization Parameters
Format:
1 2 3 4 5 6 7 8 9 10
DOPTPRM PARAM1 VAL1 PARAM2 VAL2 PARAM3 VAL3 PARAM4 VAL4
Example:
Field Contents
PARAMi Name of the design optimization parameter. Allowable names are given in
Table 1. (Character)
VALi Value of the parameter. (Real or Integer, see Table 1.)
Remarks:
1. Only one DOPTPRM entry is allowed in the Bulk Data Section.
CONVDV Relative convergence criterion on design variables. (Real > 0.0; Default = 0.001)
CONVPR Relative convergence criterion on properties. (Real > 0.0; Default = 0.001)
(Continued)
A-21
DOPTPRM Design Optimization Parameters
DELB Relative finite difference move parameter. (Real > 0.0; Default = 0.0001)
DELP Fractional change allowed in each property during any optimization design cycle.
This provides constraints on property moves. (Real > 0.0; Default = 0.2)
DELX Fractional change allowed in each design variable during any optimization cycle.
(Real > 0.0; Default = 1.0)
DOBJ1 Relative change in objective attempted on the first optimization iteration. Used to
estimate initial move in the one-dimensional search. Updated as the optimization
progresses. (Real > 0.0; Default = 0.1)
DOBJ2 Absolute change in objective attempted on the first optimization iteration. (Real
> 0.0; Default = 0.2 ∗ (F0))
DPMIN Minimum move limit imposed. (Real > 0.0; Default = 0.01)
DX1 Maximum relative change in a design variable attempted on the first optimization
iteration. Used to estimate the initial move in the one dimensional search. Updated
as the optimization progresses. (Real > 0.0; Default = 0.01)
DX2 Absolute change in a design variable attempted on the first optimization iteration.
(Real > 0.0; Default = 0.2 ∗ MAX[X(I)])
DXMIN Minimum design variable move limit (Real > 0.0; Default = 0.05).
GMAX Maximum constraint violation allowed at the converged optimum. (Real > 0.0;
Default = 0.005)
GSCAL Constraint normalization factor. See Remarks under the DSCREEN and
DCONSTR entries. (Real > 0.0; Default = 0.001)
(Continued)
A-22
Design Optimization Parameters DOPTPRM
IGMAX If IGMAX = 0, only gradients of active and violated constraints are calculated. If
IGMAX > 0, up to NCOLA gradients are calculated including active, violated, and
near active constraints. (Integer > 0; Default = 0)
IPRINT Print control during approximate optimization phase. Increasing values represent
increasing levels of optimizer information. (0 ≤ Integer ≤ 7; Default = 0)
0 no output (Default)
1 internal optimization parameters, initial information, and results
2 same, plus objective function and design variables at each iterations
3 same, plus constraint values and identification of critical constraints
4 same, plus gradients
5 same, plus search direction
6 same, plus scaling factors and miscellaneous search information
7 same, plus one dimensional search information
IPRNT1 If IPRNT1 = 1, print scaling factors for design variable vector. (Integer 0 or 1;
Default = 0)
ISCAL Design variables are rescaled every ISCAL iterations. Set ISCAL= –1 to turn off
scaling. (Integer; Default = NDV (number of design variables))
ITMAX Maximum number of iterations allowed at optimizer level during each design cycle.
(Integer; Default = 40)
ITRMOP Number of consecutive iterations for which convergence criteria must be satisfied
to indicate convergence at the optimizer level. (Integer; Default = 2)
ITRMST Number of consecutive iterations for which convergence criteria must be met at
the optimizer level to indicate convergence in the Sequential Linear Programming
Method. (Integer > 0; Default = 2)
JTMAX Maximum number of iterations allowed at the optimizer level for the Sequential
Linear Programming Method. This is the number of linearized subproblems
solved. (Integer ≥ 0; Default = 20)
JPRINT Sequential Linear Programming subproblem print. If JPRINT > 0, IPRINT is turned
on during the approximate linear subproblem. (Default = 0)
JWRITE If JWRITE > 0, file number on which iteration history will be written. (Integer > 0;
Default = 0)
(Continued)
A-23
DOPTPRM Design Optimization Parameters
P1 Print control items specified for P2. (Integer ≥ 0; Default = 0) Initial results are
always printed prior to the first approximate optimization. If an optimization task
is performed, final results are always printed for the final analysis unless
PARAM,SOFTEXIT,YES is specified. These two sets of print are not controllable.
n: Print at every n-th design cycle.
STPSCL Scaling factor for shape finite difference step sizes, to be applied to all shape
design variables. (Real > 0.0; Default = 1.0)
A-24
Design Sensitivity Response Quantities DRESP1
Defines a set of structural responses that is used in the designDesign
DRESP1 either as constraints
Sensitivity Responseor as an
Quantities
objective.
Format:
1 2 3 4 5 6 7 8 9 10
DRESP1 ID LABEL RTYPE PTYPE REGION ATTA ATTB ATT1
ATT2 -etc.-
Example:
103
Field Contents
ID Unique entry identifier. (Integer > 0)
PTYPE Element flag (PTYPE = “ELEM”) or property entry name. Used with element type
responses (stress, strain, force, etc.) to identify the property type, since property
entry IDs are not unique across property types. (Character: “ELEM”, “PBAR”,
“PSHELL”, etc.)
REGION Region identifier for constraint screening. See Remark 10 for defaults.
(Integer > 0)
ATTA, ATTB, Response attributes. See Table 1. (Integer > 0 or Real or blank)
ATTi
(Continued)
A-25
DRESP1 Design Sensitivity Response Quantities
Response Attributes
Response
Type ATTB (Integer > 0 or
(RTYPE) ATTA (Integer > 0) Real > 0.0) ATTi (Integer > 0)
FRSTRE Stress Item Code Frequency Value. (Real > Property entry (PID)
0.0) See Remark 15..
(Continued)
A-26
Design Sensitivity Response Quantities DRESP1
Response Attributes
Response
Type ATTB (Integer > 0 or
(RTYPE) ATTA (Integer > 0) Real > 0.0) ATTi (Integer > 0)
FRFORC Force Item Code Frequency Value. (Real > Property entry (PID)
0.0) See Remark 15.
TSTRE Stress Item Code Time Value. (Real) See Property entry (PID)
Remark 16.
TFORC Force Item Code Time Value. (Real) See Property entry (PID)
Remark 16.
Remarks:
1. Stress, strain, and force item codes can be found in the MSC/NASTRAN Quick Reference
Guide, Appendix A. For stress or strain item codes that have dual meanings, such as von
Mises or maximum shear, the option specified in the Case Control Section will be used;
i.e., STRESS(VONM) or STRESS(MAXS).
2. RTYPE=“CSTRESS”, “CSTRAIN”, and “CFAILURE” are used only with the PCOMP entry.
“CSTRESS” and “CSTRAIN” item codes are described under Table 1. (Element
Stress/Strain Item Codes) in the MSC/NASTRAN Quick Reference Guide, Appendix A.
“CFAILURE” item codes are described under Table Table 2. (Element Force Item Codes)
in the MSC/NASTRAN Quick Reference Guide, Appendix A. Only force item codes that
refer to failure indices of direct stress and interlaminar shear stress are valid.
(Continued)
A-27
DRESP1 Design Sensitivity Response Quantities
The CFAILURE response type requires the following specifications on the applicable
entries:
c. Stress limits in the ST, SC, and SS fields on all MATi entries.
3. ATTB is used only for responses of composite laminae, dynamics, and stability deriva-
tives. For other responses, this field must be blank.
4. All grids associated with a DRESP1 entry are considered to be in the same region for
screening purposes. Only up to NSTR displacement constraints (see DSCREEN entry)
per group per load case will be retained in the design optimization phase.
10. REGION is used for constraint screening. The NSTR field on DSCREEN entries gives the
maximum number of constraints retained for each region per load case.
If the REGION field is not blank, all the responses on this entry as well as all responses on
other DRESP1 entries that have the same RTYPE and REGION identification number will
be grouped into the same region.
(Continued)
A-28
Design Sensitivity Response Quantities DRESP1
Table 2. Default Regions for Design Sensitivity Response Types.
WEIGHT No region
VOLUME No region
LAMA No region
EIGN No region
FREQ No region
11. REGION is valid only among the same type of responses. Responses of different types
will never be grouped into the same region, even if they are assigned the same REGION
identification number by the user.
12. If RTYPE = “WEIGHT” or “VOLUME”, field ATTi = “ALL” implies total weight/volume of all
superelements except external superelements.
13. RTYPE = “STABDER” identifies a stability derivative response. ATTB is the restraint flag
for the stability derivative. ATTB = 0 means unrestrained, and ATTB = 1 means re-
strained. For example, ATTA = 4000, ATTB = 0, and ATT1 = 3 reference the unrestrained
Cz derivative for the AESTAT (or AESURF) entry ID = 4000.
(Continued)
A-29
DRESP1 Design Sensitivity Response Quantities
14. RTYPE = “FLUTTER” identifies a set of damping responses. The set is specified by ATTi:
ATT3 = Identification number of an FLFACT entry that specifies a list of Mach numbers.
15. For RTYP = “FRDISP”, “FRVELO”, “FRACCL”, “FRSPCF”, “FRFORC”, and “FRSTRE”,
ATTB specifies a frequency value in cycles per unit time. If ATTB is specified, then the
responses are evaluated at the closest frequency selected by the OFREQ command. The
default for ATTB is all frequencies selected by the OFREQ command.
16. For RTYPE = “TDISP”, “TVELO”, “TACCL”, “TSPCF”, “TFORC”, and “TSTRE”, ATTB
specifies a time value. If ATTB is specified, then the responses are evaluated at the clos-
est time selected by the OTIME command. The default for ATTB is all time steps selected
by the OTIME command.
17. Intermediate station responses on CBAR elements due to PLOAD1 and/or CBARAO en-
tries may not be defined on the DRESP1 entry.
18. RTYPE = “EIGN” refers to normal modes response in terms of eigenvalue (radi-
an/time)∗∗2 while RTYPE = “FREQ” refers to normal modes response in terms of natural
frequency or units of cycles per unit time.
19. For RTYPE = LAMA, EIGN or FREQ, the response approximation used for optimization
can be individually selected. (Approximation Code = 1 = direct linearization, = 2 = Inverse
Linearization).
A-30
Design Sensitivity Equation Response Quantities DRESP2
DRESP2 equation responses that are used in the design, either
Defines asSensitivity
Design constraints or asResponse
Equation an objective.
Quantities
Format:
1 2 3 4 5 6 7 8 9 10
DRESP2 ID LABEL EQID REGION
DVID8 -etc.-
LABL8 -etc.-
NR8 -etc.-
“DNODE” G1 C1 G2 C2 G3 C3
G4 C4 -etc.-
Example:
DRESP2 1 LBUCK 5 3
201
DTABLE PI YM L
DRESP1 14 1 4 22 6 33 2
DNODE 14 1 4 1 22 3
2 1 43 1
Field Contents
ID Unique identification number. (Integer > 0)
(Continued)
A-31
DRESP2 Design Sensitivity Equation Response Quantities
REGION Region identifier for constraint screening. See Remark 5. (Integer > 0)
“DTABLE” Flag indicating that the labels for the constants in a DTABLE entry follow.
(Character)
“DNODE” Flag signifying that the following fields are “designed grid points.” See Remark 7.
(Character)
Remarks:
1. DRESP2 entries may only reference DESVAR, DTABLE, DRESP1, and DVPREL1 en-
tries. They may not reference other DRESP2 entries.
3. DRESP2 entries must have unique identification numbers with respect to DRESP1
entries.
4. The “DESVAR”, “DTABLE”, “DRESP1”, “DNODE”, and “DVPREL1” flags in field 2 must
appear in the order given above. Any of these words, along with the identification num-
bers associated with them, may be omitted if they are not involved in this DRESP2
relationship. However, at least one of these four types of arguments must exist.
5. The REGION field follows the same rules as for the DRESP1 entries. DRESP1 and
DRESP2 responses will never be contained in the same region, even if they are assigned
the same REGION identification number. The default is to put all responses referenced
by one DRESP2 entry in the same region.
(Continued)
A-32
Design Sensitivity Equation Response Quantities DRESP2
6. The variables identified by DVIDi, LABLj, NRk, the Gm, Cm pairs and DPIPi are assigned
(in that order) to the variable names (x1, x2, x3, etc.) specified in the left-hand side of the
first equation on the DEQATN entry referenced by EQID. The variable names x1 through
xN (N = m + n + p + q) are assigned in the order DVID1, ..., DVIDm, LABL1, ..., LABLn,
NR1, ..., NRp, G1, ..., Gq, DPIP1, ..., DPIPr. In the example below,
DRESP2 1 LBUCK 5 3
DESVAR 101 3
DTABLE PI YM
DNODE 14 1
7. (Gm, Cm) refer to a designed grid component. Depending on the scheme used in gener-
ating basis vectors, a designed grid component can be one of the following:
b. A grid component that is free to move in a full auxiliary model when the external dis-
placement fields are used.
c. When the auxiliary model boundary shape method is used, a designed grid compo-
nent can be either a loaded grid component in the auxiliary model or a grid
component with motion that is obtained from interpolation of boundary shapes.
d. When the interface using geometric boundary shapes is used, a designed grid com-
ponent is either a grid component defined on a DVGRID entry or one with motion that
is obtained from interpolation of boundary shapes.
A-33
DSAPRT Design Sensitivity Output Parameters
Specifies
DSAPRT design sensitivity output parameters. Design Sensitivity Output Parameters
Format:
FORMATTED ALL
NOEXPORT
DSAPRT ( UNFORMATTED , , [ START = i ], [ BY = j ], [ END = k ] ) = n
EXPORT
NOPRINT NONE
Examples:
DSAPRT(FORMATTED,EXPORT)
DSAPRT(FORMATTED,START = FIRST,BY = 3,END = LAST) = 101
DSAPRT(UNFORMATTED,START = FIRST)
DSAPRT(UNFORMATTED,EXPORT)
DSAPRT(FORMATTED,END = 4) = ALL
DSAPRT(UNFORMATTED,END = SENS) = ALL
DSAPRT(NOPRINT,EXPORT)
Describers Meaning
FORMATTED Output will be presented with headings and labels.
UNFORMATTED Output will be printed as a matrix print (see description of the MATPRN module
in the MSC/NASTRAN DMAP Module Dictionary).
START = i Specifies the first design cycle for output. (Integer > 0 or Character: “FIRST” or
“LAST”; Default = 1 or “FIRST”)
BY = j Specifies the design cycle interval for output. (Integer ≤ 1, or > 0; Default = 0)
See Remark 2.
END = k Specifies the last design cycle for output. (Integer > 0 or Character: “FIRST”,
“LAST”, or “SENS”; Default = “LAST”)
ALL All design responses (defined in DRESP1 and DRESP2 entries) will be output.
(Continued)
A-34
Design Sensitivity Output Parameters DSAPRT
Remarks:
1. Only one DSAPRT may appear in the Case Control Section and should appear above all
SUBCASE commands.
2. Sensitivity data will be output at design cycles i, i + j, i + 2j, ..., k. Note that the BY = 0
default implies no sensitivity analysis at the intermediate design cycles.
3. END = SENS requests design sensitivity analysis, and no optimization will be performed.
4. Prior to Version 69, PARAM,OPTEXIT,4, –4, or 7 requested design sensitivity output and
may still be used in Version 69. However, if both DSAPRT and PARAM,OPTEXIT,4, –4,
or 7 are specified, then DSAPRT overrides PARAM,OPTEXIT,4, –4, or 7.
PARAM,OPTEXIT values and the equivalent DSAPRT commands are as follows:
–4 DSAPRT(NOPRINT,EXPORT,END = SENS)
7 DSAPRT(UNFORMATTED,START = LAST)
A-35
A-36
DSCREEN Design Constraint Screening Data
Format:
1 2 3 4 5 6 7 8 9 10
DSCREEN RTYPE TRS NSTR
Example:
Field Contents
RTYPE Response type for which the screening criteria apply. See Table 1. (Character)
NSTR Maximum number of constraints to be retained per region per load case. See
Remark 3. (Integer > 0; Default = 20)
Remarks:
1. Displacement and displacement derivative constraints associated with one particular load
case are grouped by the specification of DRESP1 entries. From each group, a maximum
of NSTR constraints are retained per load case.
(Continued)
A-37
Design Constraint Screening Data DSCREEN
2. Stress-strain constraints are grouped by the property; i.e., all elements belonging to the
set of PIDs specified under ATTi on a DRESPi entry are regarded as belonging to the
same region. In superelement sensitivity analysis, if the property (PID) is defined in more
than one superelement, then separate regions are defined. A particular stress constraint
specification may be applied to many elements in a region generating many stress con-
straints, but only up to NSTR constraints per load case will be retained.
4. If a certain type of constraint exists but no corresponding DSCREEN entry is specified, all
the screening criteria used for this type of constraint will be furnished by the default
values.
5. Constraints can be retained only if they are greater than TRS. See the Remarks under
the DCONSTR entry for a definition of constraint value.
A-38
DTABLE Table Constants
DTABLE a table of real constants that are used in equations (see DEQATN entry).
Defines Table Constants
Format:
1 2 3 4 5 6 7 8 9 10
DTABLE LABL1 VALU1 LABL2 VALU2 LABL3 VALU3 LABL4 VALU4
Example:
G 5.5E5 B 100.
Field Contents
LABLi Label for the constant. (Character)
Remarks:
1. Only one DTABLE entry may be specified in the Bulk Data Section.
A-39
Design Variable to Boundary Shapes DVBSHAP
Associates a design variable identification number to a linear combination
DVBSHAP of to
Design Variable boundary shape
Boundary Shapes
vectors from a particular auxiliary model.
Format:
1 2 3 4 5 6 7 8 9 10
DVBSHAP DVID AUXMOD COL1 SF1 COL2 SF2 COL3 SF3
Example:
DVBSHAP 4 1 1 1.6
Field Contents
DVID Design variable identification number of a DESVAR entry. (Integer > 0)
COLi Load sequence identification number from AUXMODEL Case Control command.
(Integer > 0)
SFi Scaling factor for load sequence identification number. (Real; Default = 1.0)
Remarks:
1. Design variable DVID must be defined on a DESVAR entry.
2. Multiple references to the same DVID and/or COLi will result in the vector addition of the
referenced boundary shape vectors.
A-40
DVGRID Design Variable to Grid Point Relation
DVGRID the relationship between design variables and grid point locations.
Defines Design Variable to Grid Point Relation
Format:
1 2 3 4 5 6 7 8 9 10
DVGRID DVID GID CID COEFF N1 N2 N3
Example:
Field Contents
DVID DESVAR entry identification number. (Integer > 0)
GID Grid point (GRID) or geometric point (POINT) identification number. (Integer > 0)
Remarks:
1. A CID of zero or blank (the default) references the basic coordinate system.
2. Multiple references to the same grid ID and design variable result in vectorial addition of
the participation vectors defined by CID, COEFF, and Ni. There is no restriction on the
number of DVGRID entries that may reference a given grid (GID) or design variable
(DVID).
T
where { g } i is the location of the i-th grid, [ g x g y g z ] .
T
The vector { N } = [ N x N y N z ] is determined from CID and Ni. Note that it is a change
0
in a design variable from its initial value X , and not the absolute value of the design
0
variable itself, that represents a change in a grid point location, { g } i – { g } i .
(Continued)
A-41
Design Variable to Grid Point Relation DVGRID
4. The DVGRID entry defines the participation coefficients (basis vectors) of each design
variable for each of the coordinates affected by the design process in the relationship
{ ∆g } i = ∑j { T }ij ∗ ∆X j
5. DVGRID entries that reference grid points on MPCs or RSSCON entries produce incorrect
sensitivities. Often the sensitivities are 0.0 which may result in a warning message indi-
cating zero gradients which may be followed by UFM 6499. Other rigid elements produce
correct results.
A-42
DVPREL1 Design Variable to Property Relation
DVPREL1 the relation between an analysis model property and design variables.
Defines Design Variable to Property Relation
Format:
1 2 3 4 5 6 7 8 9 10
DVPREL1 ID TYPE PID FID PMIN PMAX C0
Example:
Field Contents
ID Unique identification number. (Integer > 0)
FlD Field position of the property entry, or word position in the element property table
of the analysis model. (Integer ≠ 0)
PMIN Minimum value allowed for this property. If FID references a stress recovery
location, then the default value for PMIN is –1.0+35. PMIN must be explicitly set
to a negative number for properties that may be less than zero (for example, field
ZO on the PCOMP entry). (Real; Default = 1.0E–20)
PMAX Maximum value allowed for this property. (Real; Default = 1.0E–20)
Remarks:
1. The relationship between the analysis model property and design variables is given by:
P i = C0 + ∑i COEFi ∗ DVIDi
2. The continuation entry is required.
(Continued)
A-43
Design Variable to Property Relation DVPREL1
3. PTYPE = “PBEND” is not supported, either directly through FIDs or indirectly via word po-
sitions in the element property table.
4. FID may be either a positive or a negative number. If FID > 0, it identifies the field position
on a property entry. If FID < 0, it identifies the word position of an entry in the element
property table. For example, to specify the area of a PBAR, either FID = +4 or FID = –3
can be used. However, if PTYPE = “PBEAM”, FID must be negative. See the following
element property table for the word positions for PBEAM.
(Continued)
A-44
DVPREL1 Design Variable to Property Relation
A-45
Design Variable to Property Relation DVPREL2
Defines the relation between an analysis model property and design Design
DVPREL2 variables withtoaProperty
Variable user-supplied
Relation
equation.
Format:
1 2 3 4 5 6 7 8 9 10
DVPREL2 ID TYPE PID FID PMIN PMAX EQID
Example:
DESVAR 4 11 13 5
DTABLE PI YM
Field Contents
ID Unique identification number. (Integer > 0)
FID Field position of the property in the analysis model entry. (Integer ≠ 0)
PMIN Minimum value allowed for this property. If FID references a stress recovery
location field, then the default value for PMIN is –1.0+35. PMIN must be
explicitly set to a negative number for properties that may be less than zero (for
example, field ZO on the PCOMP entry). (Real; Default = 1.E–20)
PMAX Maximum value allowed for this property. (Real; Default = 1.0E20)
“DESVAR” DESVAR flag. Indicates that the IDs of DESVAR entries follow. (Character)
“DTABLE” DTABLE flag. Indicates that the IDs for the constants in a DTABLE entry follow.
This field may be omitted if there are no constants involved in this relation.
(Character)
(Continued)
A-46
DVPREL2 Design Variable to Property Relation
Remarks:
1. The variables identified by DVIDi and LABLi correspond to variable names (x1, x2, etc.)
listed in the left-hand side of the first equation on the DEQATN entry identified by EQID.
The variable names x1 through xN (where N = m + n) are assigned in the order DVID1,
DVID2, ..., DVIDn, LABL1, LABL2, ..., LABLm.
2. If both “DESVAR” and “DTABLE” are specified in field 2, “DESVAR” must appear first.
3. FID may be either a positive or a negative number. If FID > 0, it identifies the field position
on a property entry. If FID < 0, it identifies the word position of an entry in EPT. For
example, to specify the area of a PBAR, either FID = +4 or FID = –3 may be used. How-
ever, if PTYPE = “PBEAM”, FID must be negative. See Remark 4 on the DVPREL1 entry
description for specification of the PBEAM element property tables.
4. PTYPE = “PBEND” is not supported, either directly through FIDs or indirectly via word po-
sitions in the element property table.
A-47
DVSHAP Design Variable to Basis Vector(s)
Format:
1 2 3 4 5 6 7 8 9 10
DVSHAP DVID COL1 SF1 COL2 SF2 COL3 SF3
Example:
Field Contents
DVID Design variable identification number on the DESVAR entry. (Integer > 0)
SFi Scaling factor applied to the COLi-th column of the displacement matrix. (Real;
Default = 1.0)
Remarks:
1. DVID must be defined on a DESVAR entry.
3. Multiple references to the same DVID and/or COLi will result in a linear combination of
displacement vectors. In the example above, the shape basis vector is a linear combina-
tion of the fourth column and twice the second column.
4. The displacement matrix must have been created by MSC/NASTRAN and be available on
a database, which is attached via the DBLOCATE FMS statement shown below:
A-48
Mode Tracking Parameters MODTRAK
MODTRAK parameters for mode tracking in design optimization (SOL 200).
Specifies Mode Tracking Parameters
Format:
1 2 3 4 5 6 7 8 9 10
MODTRAK SID LOWRNG HIGHRNG MTFILTER
Example:
Field Contents
SID Sets identification number that is selected in the Case Control Section with the
MODTRAK command. (Integer; No Default) See Remark 1.
HIGHRNG Highest mode number in range to search. See Remark 2. (Integer > 0,
Default = number of eigenvalues extracted. If nonzero, LOWRNG < HIGHRNG.)
Remarks:
1. Only the designed modes for the subcase will be tracked. A designed mode is one that is
used in the design model (in connection with either objective or constraints) and, there-
fore, identified on a DRESP1 entry.
2. The range of modes LOWRNG through HIGHRNG, inclusive, will be used to track the de-
signed modes. If LOWRNG and HIGHRNG are both blank, then all computed modes will
be used to search for the designed modes. Since large numbers of computed modes will
result in higher computational costs, limiting the search range with LOWRNG and HIGH-
RNG is recommended.
3. Modes are considered to correlate if their mass normalized cross orthogonalities are
greater than MTFILTER.
A-49