NASA / T M,2000-210305

Advancement of Bi-Level Integrated

System Synthesis (BLISS)

Jaroslaw Sobieszczanski- Sobieski

Langley Research Center, Hampton, Virginia

Mark S. Emiley

George Washington University

Joint Institute for Advancement of Flight Sciences, Hampton, Virginia

Jeremy S. Agte

United States Air Force, San Antonio, Texas

Robert R. Sandusky, Jr.

George Washington University

Joint Institute for Advancement of Flight Sciences, Hampton, Virginia

December 2000

https://ntrs.nasa.gov/search.jsp?R=20010020389 2018-05-14T14:44:44+00:00Z
The NASA STI Program Office ... in Profile

Since its founding, NASA has been dedicated to

the advancement of aeronautics and space

science. The NASA Scientific and Technical

Information (STI) Program Office plays a key

part in helping NASA maintain this important

role.

The NASA STI Program Office is operated by

Langley Research Center, the lead center for

NASA's scientific and technical information. The

NASA STI Program Office provides access to the

NASA STI Database, the largest collection of

aeronautical and space science STI in the world.

The Program Office is also NASA's institutional

mechanism for disseminating the results of its

research and development activities. These

results are published by NASA in the NASA STI

Report Series, which includes the following

report types:

TECHNICAL PUBLICATION. Reports of

completed research or a major significant

phase of research that present the results of

NASA programs and include extensive

data or theoretical analysis. Includes

compilations of significant scientific and

technical data and information deemed to

be of continuing reference value. NASA

counterpart of peer-reviewed formal

professional papers, but having less

stringent limitations on manuscript length

and extent of graphic presentations.

TECHNICAL MEMORANDUM. Scientific

and technical findings that are preliminary

or of specialized interest, e.g., quick release

reports, working papers, and

bibliographies that contain minimal

annotation. Does not contain extensive

analysis.

CONTRACTOR REPORT. Scientific and

technical findings by NASA-sponsored

contractors and grantees.

CONFERENCE PUBLICATION. Collected

papers from scientific and technical

conferences, symposia, seminars, or other

meetings sponsored or co-sponsored by

NASA.

SPECIAL PUBLICATION. Scientific,

technical, or historical information from

NASA programs, projects, and missions,

often concerned with subjects having

substantial public interest.

TECHNICAL TRANSLATION. Englishlanguage translations of foreign scientific

and technical material pertinent to NASA's

mission.

Specialized services that complement the STI

Program Office's diverse offerings include

creating custom thesauri, building customized

databases, organizing and publishing research

results ... even providing videos.

For more information about the NASA STI

Program Office, see the following:

• Access the NASA STI Program Home Page

at http'//www.sti.nasa.gov

• E-mail your question via the Internet to

help@sti.nasa.gov

• Fax your question to the NASA STI Help

Desk at (301) 621-0134

• Phone the NASA STI Help Desk at

(301) 621-0390

Write to:

NASA STI Help Desk

NASA Center for AeroSpace Information

7121 Standard Drive

Hanover, MD 21076-1320
N AS A / TM-2000-210305

Advancement of Bi-Level Integrated

System Synthesis (BLISS)

Jaroslaw Sobieszczanski-Sobieski

Langley Research Center, Hampton, Virginia

Mark S. Emiley

George Washington University

Joint Institute for Advancement of Flight Sciences, Hampton, Virginia

Jeremy S. Agte

United States Air Force, San Antonio, Texas

Robert R. Sandusky, Jr.

George Washington University

Joint Institute for Advancement of Flight Sciences, Hampton, Virginia

National Aeronautics and

Space Administration

Langley Research Center

Hampton, Virginia 23681-2199

December 2000
Available from:

NASA Center for AeroSpace Information (CASI)

7121 Standard Drive

Hanover, MD 21076-1320

(301) 621-0390

National Technical Information Service (NTIS)

5285 Port Royal Road

Springfield, VA 22161-2171

(703) 605-6000
ADVANCEMENT OF BI-LEVEL INTEGRATED SYSTEM

SYNTHESIS (BLISS)

Jaroslaw Sobieszczanski-Sobieski*

NASA Langley Research Center, Hampton, Virginia

j.sobieski@larc.nasa.gov

Mark S. Emileyt

George Washington University

Joint Institute for the Advancement of Flight Sciences (JIAFS)

m.s.emiley@larc.nasa.gov

Jeremy S. Agte$

United States Air Force, San Antonio, Texas

j eremy.agte@kelly.af.mil

Robert R. Sandusky, Jr.§

George Washington University

Joint Institute for the Advancement of Flight Sciences (JIAFS)

r.r.sandusky@larc.nasa.gov

Abstract Introduction

Bi-Level Integrated System Synthesis (BLISS) is a

method for optimization of an engineering system, e.g.,

an aerospace vehicle. BLISS consists of optimizations

at the subsystem (module) and system levels to divide

the overall large optimization task into sets of smaller

ones that can be executed concurrently. In the initial

version of BLISS that was introduced and documented

in previous publications, analysis in the modules was

kept at the early conceptual design level. This paper

reports on the next step in the BLISS development in

which the fidelity of the aerodynamic drag and

structural stress and displacement analyses were

upgraded while the method's satisfactory convergence

rate was retained.

* Manager, Computational AeroSciences, and Multidisciplinary

Research Coordinator, NASA Langley Research Center, MS 139,

Hampton, Virginia, AIAA Fellow

t Graduate Student, NASA Langley Research Center,

Multidisciplinary Optimization Branch, MS 139,Member AIAA

1st Lt United States Air Force, San Antonio Air Logistics Center,

TX, Member AIAA

§ Professor, George Washington University, Joint Institute for

Advancement of Flight Sciences, Fellow AIAA

BLISS, for Bi-Level Integrated System Synthesis, is an

optimization method for engineering a modular system,

e.g., an aerospace vehicle, in which it is beneficial to

separate the design variables and constraints local to a

module from those that govern the entire system. Such

separation fosters development of a broad workfront of

people and computers, hence it fits well the current

trends for massively parallel processing in computer

technology and the concurrent engineering style of the

workforce organization.

The focus on dividing the optimization into the

suboptimizations within each module (subsystem, also

called the black box) and a coordinating optimization at

the system level places BLISS in the Multidisciplinary

Design Optimization (MDO)toolbox, in the company

of a few other methods that have the same focus as

BLISS but differ in approach. Representative examples

of these methods are the Collaborative Optimization

(CO) (Braun et al, 1965), the Concurrent SubSpace

Optimization (CSSO) (Sobieszczanski-Sobieski, 1988,

Bloebaum et al, 1992), and the Concurrent Design

Optimization (Wujek et al, 1995).

The distinguishing features of BLISS are the use of the

system objective (e.g., the aircraft range) as the

optimization objective in each of the subsystems and at

the system level, and coupling between the

optimizations at the system and subsystem levels by

the optimum sensitivity derivatives with respect to

parameters.
TheoverallarchitectureofBLISSasamethoddoesnot

dependonthefidelityoftheanalysesperformedin each

module.Consequently, in principleat least,BLISS

maybeusedin anydesignphasefromconceptual,

throughpreliminaryto detailed,providedthat

appropriatelvelof analysisis implementedin the

modules.

TheBLISSmethodwasintroducedin (SobieszczanskiSobieskiet al, 1998a)anddocumentedin detailin

(Sobieszczanski-Sobiesk! etal, 1998b*).In thepaper,

thatprototypeis referredto astheoriginalBLISS.In

itsoriginalformBLISSmoduleswerekeptverysimple

correspondingto the earlyconceptual designphase.

SatisfactoryresultsfromtheinitialtrialsofBLISSona

testcaseof abusinessjet encouraged nextstepin the

BLISSdevelopment- upgradingits structural analysis

andaerodynamic draganalysismodules- andvalidating

onthesametestcase.

ThispapereportsontheaboveBLISSupgradeand

resultsof thetestingthatadvancethemethodstoward

becomingatoolsuitableforpracticalapplications. The

reportprovidesa synopsisof the BLISSmethod,

describesthe salientfeaturesof the two upgraded

modules,presentssatisfactoryconvergence results,and

summarizestheBLISSdevelopment statusandthe

futuredevelopment direction.

Notation

ARm - tail aspect ratio

ARw - wing aspect ratio

BB_ - black box

CD- coefficient of drag

Cf- skin friction coefficient

D - drag

ESF - engine scale factor

h - altitude

k - safety factor

L - lift

L/D - lift to drag ratio

LHT- horizontal tail location, % mean aerodynamic

chord (% MAC)

Lw - wing location, % MAC

M - Mach number

Nz - maximum load factor

R - range

SFC - specific fuel consumption

Sm- horizontal tail surface area

SREr-- wing surface area

T - throttle

t/c - thickness to chord ratio

*The 1998a and 1998b references are also available at

http://techreports.larc.nasa.gov/ltrs/

t_- wingbox sandwich face sheets thicknesses

ts,_- wingbox sandwich caliper thicknesses

X_ - design variables local to BB_

XL, XU - lower and upper bounds on X, sideconstraints

W_ - engine weight

WF -- fuel weight

WT -- total weight

Y_.j - behavior variables output from BB_ and sent as

inputs BBj

Z - system-level design variables

9_- taper ratio

AHT-- horizontal tail sweep

Aw - wing sweep

(9 - wing twist

Synopsis of BLISS

A synopsis of BLISS that also appeared in Agte et al,

1999 is as follows.

BLISS is a method for optimization of engineering

systems that separates the system-level optimization

from potentially numerous autonomous subsystem

optimizations. As shown in Figure 1, it utilizes a

system architecture in which design and behavior

variables are split into three categories. X-variables are

those design variables optimized at the local level and

are unique to each particular subsystem. Behavior

variables that are output from one subsystem and input

to another are designated Y, and the system-level

design variables are specified as Z. System-level

variables are those shared by at least two subsystems.

Z - system level

(1 _--------[ 12

Subsystem 1 i_g, S _l ubsystem 2

Subsystem 3 _-------@ubsystem 4

I I x4

Figure 1: BLISS system structure

After a best guess initialization, the first step in the

BLISS begins with the system analysis and sensitivity

analysis in which Y and the derivatives of Y with

respect to Z and X are computed. A linear

approximation to the system objective (an element of

Y) as a function of Z and X is established using the

above derivatives. That approximation is adopted as

the objective function in subdomain optimizations that
follownext. In eachsubdomain(module,or black

box), the Z andY variablesarefrozenand an

improvement itheobjectivefunctionis soughtbythe

localoptimizationsthatuselocalX separatelyin each

module.ThefrozenZandYareconstant parameters in

eachmoduleoptimizationandthemoduleoptimization

is followedby computationof thederivativesof the

optimumwithrespectto theseparameters. Thesecond

stepachievesimprovement throughthe system-level

variablesZ andis linkedto the first stepby the

derivativesof optimumwithrespectto parametersZ

andY. Thederivativesareusedto extrapolateeach

subdomainoptimumasa functionof Z andY. The

functionalrelationY=Y(Z) is approximatedby

extrapolationbasedonthesystemsensitivityanalysis.

Thesestepsalternateuntilconvergence. A flowchartof

themethodisshowninFigure2.

Sle413_ Slea2_

Opportmlty for Concun'ent Processing

initialize X & Z

System Analysis

and Updale

[Sensitivity Analysis I Variables

X = X o + AXoI, r

Z =Z0 + AZovr X = X0 + AXovr

4 BLISS CYCLE Z =Z° + AZ°pT l

Figure 2: BLISS Cycle

Note that the output of step 1 is an optimum change in

the local design variables, AXoPT, in the presence of

constant Z, and the output of step 2 is an optimum

change in system design variables, AZoPT.

In the original version of BLISS the modules shown

generically in Fig. 2 are Propulsion, Aerodynamics,

Structures, and Performance whose detailed

input/output variables are identified later. The common

denominator of these modules was the extreme

simplicity of analyses that employed closed-form

expressions for input-to-output mapping. This was so

because of the need to test the overall procedure

organization and the two-level algorithm convergence at

the initial development stage of a new method without

being encumbered by long turn-around times in the

modules. The next logical step in the BLISS

development is to upgrade the fidelity of the modules

while holding the overall procedure organization

unchanged.

Upgrades in the BLISS Structures and

Aerodynamics Modules

The modularity of BLISS permits replacing or adding

black boxes to refine or alter the optimization and

analysis tools in each modules allowing the engineer

the flexibility to exercise his judgment. Having tools

of different level of fidelity in the modules enables

applications of BLISS in different design phases. The

advanced BLISS method incorporates two new modules

that can be used in lieu of previous black boxes. The

structures module now can use the Equivalent

Laminated Plate Solution (ELAPS, Giles, 1986) and

the aerodynamics module can use a code called

AWAVE (Harris, 1964) to perform wave drag analysis.

Integration of ELAPS

In the previous application example, BLISS employed

a skin-stringer representation of the internal wing box

bays. This model broke the wing down into a three

bay wing box whose geometry varied with the taper

ratio, wing sweep, thickness to chord ratio, wingspan,

and aspect ratio, all manipulated as design variables in

the system-level optimization. The displacements,

e.g., the wing twist, and stresses, were computed using

simple, thin-walled box-beam formulas (e.g., Bruhn,

1965)

In the BLISS application shown herein, the level of

accuracy in this module is raised by substituting the

previous model with the Equivalent Laminated Plate

Solution (ELAPS) computer code. This code designed

with preliminary design stage calculations in mind is

capable of modeling aircraft wing structures with

multiple trapezoidal segments. The wing structure is

represented as a plate whose stiffness is set equivalent

to that of the original, built-up, structural box of the

wing. ELAPS employs a set of displacement fimctions

defined over each trapezoidal segment and made

compatible in regard to translations and rotations at the

segment junctions. Minimization of the strain energy

based on the Ritz method leads to equations from

which to calculate static deflections and internal forces.

The latter are then converted to stresses taking into

account the details of the wing box built-up crosssection.

The accuracy of the results of ELAPS has been found to

be somewhat below that of finite element codes (Giles,

1986) but the ELAPS input is much simpler and faster

to develop. The computation time for an ELAPS

model is more than an order of magnitude faster than

that of an equivalent finite element model - an

important feature for a tool to be integrated into an

optimization procedure.
Integratedin BLISS,ELAPSreceivesits inputfroma

we-wocessor routine that generates an input file with

the skin thickness, aspect ratio, taper ratio, thickness to

chord ratio, sweep, reference area, and aircraft weight.

The model used by ELAPS analyzes stress along the

same three bay wingbox configuration used as an

example in the original application of BLISS. Each

wingbox consists of the top and bottom sandwich

panels of different thicknesses and sandwich webs

identical in the front and rear of the wingbox. The

front spar of the wing box is located at 10% of the

chord length and the rear spar lies at 70% of the chord

length. Figure 3 depicts the configuration of the

ELAPS model used by BLISS.

y

Wingbox

/

d

at 0.5 chord

Ill X

Figure 3: Wing Model

The top and bottom panels as well as the webs have the

thickness of the sandwich face sheets (t) and the

sandwich caliper thickness (ts) as design variables, as

depicted in Fig.4. ELAPS models such a built-up

structure by representing each face and the core as

separate elements linked in a common coordinate grid.

[, L ;[

O + O

_1--___-- t_2 t_z-'_ ]['_'-

o ¢ _ Ct j i _ o

_ T

Figure 4: Wingbox Model

As it was done in the original BLISS implementation,

the aerodynamic loads are being generated within the

structures module in the pre-processor to structural

analysis. To calculate the lift loads on the wing, the

we-processor routine averages spanwise between an

elliptical lift distribution and a linear distribution that

reflects the wing chord taper. The elliptical and taper

ratio based lift distributions for the wing are each

normalized to contain an area of unity as illustrated in

Figure 5. The averaged, spanwise load distribution is

multiplied by the lift required from the wing and

distributed chordwise. The chordwise distribution is a

typical supersonic one with the center of pressure

located at 50% of the chord. The aerodynamic load

distribution would be expected to be calculated by an

aerodynamics module using a higher fidelity analysis,

e.g., a computational fluid dynamics code. Thus, the

present aerodynamic loads generation is merely a

placeholder for a real aerodynamic loads analysis in a

future BLISS upgrade.

In summary, the structural module employs ELAPS to

calculate the stresses in the wing box for the given

configuration, lift distribution, and corresponding

constraints. It also outputs the wing twist and weight

and the objective function for the local optimization.

The aerodynamics module accepts the output and

models its influence on the aerodynamic response.

Spanwise Location

I • _aper onape . • • ¢,=pt,ca¢ _ Moo,n_ I

Figure 5: ELAPS Lift Distribution

Integration of AWAVE

The cruise segment of the test case mission is

supersonic. The original model (SobieszczanskiSobieski et al, 1998) used an approximation relying on

the span efficiency factor. That approximation was

replaced herein with a code, AWAVE, that is a

streamlined version of the far-field wave drag program

(Harris, 1964). There are two versions of the Harris

wave drag program in common use at LaRC. The

original version, described in the reference, treats lifting

surfaces as a series of 3-dimensional solid elements. A

much faster but slightly less accurate version treats

lifting surfaces as 2-dimensional panels with finite

thicknesses. Due to compensating errors at positive

and negative roll angles of the Mach cutting plane, the

4
panelversiongivesexcellentresultsfor wavedrag

coefficients.Theobjectiveofthelast(AWAVE)effort

wastodevelopaversionofthewavedragprogramwith

theaccuracyof thesolidelementprogramthatis faster

than the panel version. The AWAVE code

implementedcomputesthewavedragonthebasisof

the aircraftcross-sectiondistributionalong the

centerline,henceit requiresdataaboutthe entire

configurationgeometryto enablethearearuling of

supersonicbodydesign.

Similarlyto theintegrationof ELAPS,integrationof

AWAVEwasaccomplishedbycreatingapre-processor

togeneratehenecessaryinput.Theinputprovidesthe

currentdesign'saspectratio,taperatio,thicknessto

chordratio,sweepangle,wingreferencearea,horizontal

tail sweepangle,horizontaltail aspectratio, and

horizontaltail referencearea.Thepre-processor also

createsandplacesthewingandtailairfoilsaccordingto

the designconfigurationvariables. TheAWAVE

outputis thewavedragcoefficient to beaddedto the

otherdragcomponentswhosecalculationremainsthe

sameasintheoriginalBLISS.

Numerical Implementation

Compared to the original application of BLISS to the

supersonic business jet case, incorporation of ELAPS

and AWAVE in BLISS required some changes to

constraints and allocation of the design variables to the

system and subsystem levels.

In the original BLISS, the taper ratio was a local

variable of the structures module. With the integration

of ELAPS and AWAVE, the taper ratio affects both the

aerodynamics and structures module. While the

aerodynamics module optimization may tend toward a

taper ratio to reduce induced drag, the structures

Z - Variables

hREF,S HT ,AR HT, _

z_-----_ Wvo, W o, N z, WBL. CD[V_ N,M<I,]]H

Constants

Range

T -throttle

Am,- tail sweep

Lw-see Figure 1

L m.-see Figure 1

[t]-thickness array,

size lx9

[ts]-thickness array,

size lx9

)_-tap er ratio

D-drag

ESF-eng. scale fact.

L-lift

Nz-max. load fact.

R-range

SFC-spec. fuel cons.

O-wing twist

WE-engine weight

W F-fuel weight

WT-total weight

AR w- wing aspect ratio

ARnT- tail aspect ratio

h-altitude

M- Mach #

SREF-wing surf area

SnT-tail surf area

t/c-thickness/daord

Aw-wing sweep

I I I I I I

X Y Z

Figure 6: Data Dependencies for Business Jet

Model

module may generate a different taper value to reduce

stresses. To resolve this trade-off, the taper ratio was

raised to a system variable, capable of influencing both

modules. Figure 6 shows the current black box and

variable interactions.

In this model there are nine system-level Z-variables,

each influencing a minimum of two of the subsystems.

The local variables of each subsystem are manipulated

only in the optimization local to that subsystem. The

propulsion module has the throttle as its sole local

variable. In the present state of BLISS, the range

module is an exception as it performs no optimization.

It only evaluates the Breguet range formula. The

aerodynamics module optimizes the local variables of

the horizontal tail sweep as well as the variables that

place the wing and tail along the fuselage axis. The

structural subsystem optimization operates on the

sandwich face sheet and caliper thicknesses for the wing

cover panels and the webs of the three wingbox bays.

The ten Y-variables noted in the off-diagonal boxes in

Figure 6, represent couplings of the black boxes and are

computed in the system analysis.
BLISSwasoriginallyimplementedin MATLAB5.3.0

at both the systemand subsystemlevels. The

MATLAB Optimization Toolbox was used as an

optimizer in the subsystem and system optimizations.

In the version of BLISS reported herein, the use of

MATLAB continued as above with the exception of the

structures and aerodynamics modules that incorporated

ELAPS and AWAVE, both written in FORTRAN 77.

MATLAB provides a facility to invoke FORTRAN

from a MATLAB code. To exploit that facility, the

preprocessors to both ELAPS and AWAVE were

written in FORTRAN and converted into MEX-Files

using the MATLAB mex-function (Appendix). Both

AWAVE and ELAPS were then directly called from

within the BLISS modules. On the output side,

simple post-processing generated outputs in a format

acceptable to the parts of BLISS that remained being

coded in MATLAB for further analysis. Because of the

MATLAB ability to invoke FORTRAN codes, the

BLISS upgrading process may continue by adding

FORTRAN-coded modules wherever required while

retaining the MATLAB core that executes the method

logic illustrated by the flowchart in Fig. 1.

Results

BLISS iterations terminate when the change in the

aircraft range objective varies less than ten nautical

miles. This took seven passes through the flowchart in

Figure 2. The system-level design variables converged

within the first few passes. Further optimizations

focused primarily on the local variables. Most of the

changes occurred within the structures module where

the new ELAPS-based optimization kept refining the

variables searching for the best solution. The majority

of the computational time was spent in this module.

Table 1 shows the variable progression through the

optimization process.

The table reflects the major trade-offs that occur

between the wing sweep angle, airfoil thickness ratio,

and the wing aspect ratio, all of which govern the

structural weight and drag that, in turn, influence the

range. Ultimately, influences of these variables on the

range differ in sign, therefore, the procedure seeks a

compromise. For example, the wing sweep initially

increases to approximately 70 degrees and then falls to

40 while the taper ratio decreases to 0.1. The wing

reference area rapidly reduces to 200 square feet as the

wing aspect ratio is brought down first to 2.5 and then

increased to 2.607. The wing position is briefly

changed in the fourth cycle but quickly returns to its

initial value. The wing configuration progression is

var_cycle 1 2 3 4 5 6 7 6

Range(NM) 1051 449; 393E 354_ 339 c 2432 249: 249:

tl(inner) 2 0.617 1.32(_ 0.94z 0.643 0.192 0.19; 0.19;

tt(middle) 2 0.301 0.62,1 0.63; 1.028 0.663 0.66" 0.66"

tt(outer) 2 0.383 0.53 0.32; 0.45 0.044 0.04h 0.04_

t2(inner) 2 1; 7.44_ 5.57 3.98c_ 0.012 0.01; 0.01;

t2(middle) 2 1; 2.803 0.90 C. 8.86, 0.012 0.01" 0.01;

tz(outer) 2 12 0.E 0.01; 12 8.115 8.11, = 8.11,'

t3(inner) 2 0.4E 0.72_ c 0.82L 0.173 1.075 1.07,' 1.07,'

t3(middle) 2 1.34E 0.56; 0.511 0.548 0.517 0.517 0.51;

t3(outer ) 2 0.115 0.267 0.24E 0.265 0.012 0.01; 0.01;

t,_(inner) 4 0.617 1.32_ 0.94z 1.985 2.205 2.20,' 2.20_

tsl(middle) 4 0.868 0.624 0.637 1.028 0.6631 0.66," 0.66_

ts_(outer) 4 0.383 0.53 0.32_ 0.617 0.233 0.23_ 0.23,"

t,2(inner) 4 24 24 23.17 12A6 3.46E 3.46_ 3.46_

ts2(middle) 4 24 0.464 0.66_ 23.02 14.31 14.31 14.31

t,z(outer ) 4 24 2.583 1.80; 24 2.2_ 2.22_ 2.22£

t,3(inner) 4 0.48 0.37 0.677 0.177 1.078 1.07E 1.07_

t_3(middle) 4 1.348 0.562 0.51 0.058 0.391 0.391 0.391

t,3(outer ) 4 0.115 0.267 0.24E 0.38 0.251 0.251 0.251

ANT (°) 60 70 7C 7£ 40 40 4C 4£

Lw (%MAC) 10 1 1 2 1 1 1 1

LHT (%MAC) 250 35C 35C 10C 350 10C 10C 10C

1 0.35 0.319 0.236 0.241 0.255 0.281 0.31 0.31

t/c 0.05 0.058 0.058 0.058 0.058 0.058 0.058 0.058

h (ft) 55000 60000 6000C 6000C 6000( 60000 6000C 6000C

M 1.8 2 2 2 2 2 2

ARw 4 2.5 2.607 2.607 2.60; 2.607 2.607 2.607

Aw(° ) 55 40 40.63 40.63 40.6" 40.63 40.63 40.63

SRE F (ftz) 400 200 200 20C 20( 200 200 200

SHT (ft') 150 150 150 150 15( 150 150 150

ARHT 6.5 8.5 8.5 8.5 8._ 8.5 8.! 8.5

'_w 0.2 0.1 0.1 0.1 0._ 0.1 0." 0.1i

Table 1: Supersonic Business Jet Results

depicted in Figure 7. The aircraft finds its optimal

cruise conditions after the first cycle of Mach 2.0 at

60,000 feet.

2O

18

12

_,o

e

_5

Wing Ranform Progression

Cycle 2 Cycles 3_ ycI__

2 4 6 _ 10 12 14 1_ _ 20 22 _4 2_ _a a0 _ _

Distance from Wing Leading Edge (ft)

Figure 7: Wing Planform Progression

The horizontal tail position and geometry stabilize after

the main wing variables reach their settling points. The

tail position varies significantly but settles at a value of

100 percent of the mean aerodynamic chord. The tail

sweep ends up almost matching the wing sweep but has

a significantly larger aspect ratio. Further analysis of

the tail may involve incorporating an ELAPS model of
the tail to increasefidelity of analysisin that

component.

Theskinthicknesseschangethroughoutthe process

seekingthemaximumof thestructurecontributionto

the rangeunderthestressconstraintsforthe given

configurationgeometry,thelattergovernedbytheZvariables.TheresultinghistogramisseeninFigure8.

4_

35

3

'_25

_2

u

1

o5

o

ts3

__ _t3(hner) I

- - t3(rdddle) _

= = 't3(outer)

_ts3 (inner) I

m . ts3 (mddle) ---

- - ,ts3 (outer}

• . -. .........

2 3 4 5 6 7 S

BLISS Cycle Number

Figure 8: Plot of Skin Thickness Variation

Though the first cycle was able to converge to

reasonable thicknesses, the second through fifth cycles

were unable to satisfy all constraints given the systemlevel configuration. Then, by the sixth cycle the

optimizer had found a solution that allowed all

constraints to be met and in the seventh cycle it found

the optimal configuration. Figure 9 shows the

progression of the aircraft Take-Off Gross Weight and

its components of empty weight and fuel.

Aircraft Total and Empty Weigh

_Z

2 _ 4 _ 6 7 S

BLISS Cycle Number

Figure 9: Aircraft Weight

Figure 10 depicts a histogram of the aircraft range. It

starts off at a feasible design point. The cycles two

through five did not lie within the design space, but

BLISS returned to the design space and settled on a

feasible design with optimized range.

The last implementation of BLISS to the supersonic

business jet test case (Agte, 1999) yielded a range of

2,189 nautical miles. With the addition of AWAVE

and ELAPS, the more refined analysis increased the

range to 2,493 nautical miles.

Distribution of elapsed computing time over the BLISS

modules is displayed in Table 2. It is evident that

most of the elapsed time is spent in ELAPS but that

would change drastically if a CFD-level analysis were

used in the aerodynamics module. If BLISS were

grown to the point where all the major modules would

consume about equal amount of the elapsed time, then

distributed ,execution on concurrently operating

machines (or processors within a multiprocessor

machine) would radically compress the elapsed time of

the entire BLISS execution.

I BLISS I ELAPS AWAVE I

I Percent of Time 8.36% 18978% 1.87%

Table 2: Processor Time Use

The next step in the BLISS development is to

incorporate additional modules to increase the analysis

fidelity. The largest refinement would be expected

from adding a computational fluid dynamics code to

perform the aerodynamic analysis, including the loads.

The propulsion data quality would benefit from

replacing the current response surface fitted to a look-up

table with a comprehensive engine analysis. Also, the

simple Breguet formula for the aircraft range would

need to be replaced by a complete performance analysis

Conclusions and Remarks

Integration of ELAPS and AWAVE into BLISS

demonstrated the modular nature of the method and its

ability to accommodate refinements. Used in a limited

test case of a supersonic business jet design, the twolevel optimization in BLISS was effective in satisfying

the system-level and local constraints while attaining a

system-level objective within a reasonable number of

iterations. Separation of the system-level design

variables from the local ones enabled optimization for a

system-level objective while providing autonomy of the

design decision and tool choice within disciplines

represented in the modules.

The method is open to further upgrades in terms of the

fidelity of analysis and optimization techniques

employed in the modules. In this regard, it is up to the

user to decide on the variety of tools to be integrated in

BLISS as needed by the multidisciplinary optimization

task at hand.

Further advancement of BLISS from its present status

of a method concept demonstrator to a tool useful in
actualapplicationscallsforinsertingaCFDcodein the

aerodynamics module,addinga comprehensive engine

analysisto thepropulsionmodule,andextendingthe

performanceanalysismoduleto includemorethanjust

thecruisephaseofamission.

Finally,asthe increasedfidelity of analysesin the

moduleswill exactits pricein terms of the computing

elapsed time, concurrent execution of the modular

analyses and optimizations will become an attractive

option.

Acknowledgement

Michael Parks of the George Washington University

JIAFS program provided critical help with MATLAB

refinements and FORTRAN integration. His

contributions were gratefully appreciated.

References

Agte, J.; Sobieszczanski-Sobieski, J.; and Sandusky,

R.: "Supersonic Business Jet Design Through Bi-Level

Integrated System Synthesis," 1999 World Aviation

Conference, SAE 1999-01-5622.

AIAA/UTC/Pratt & Whitney Undergraduate Individual

Aircraft Design Competition, "Supersonic Cruise

Business Jet RFP." 1995/1996.

Bloebaum, C. L.; Hajela, P.; and SobieszczanskiSobieski, J.: "Non-Hierarchic System decomposition in

Structural Optimization; Engineering Optimization,"

19, pp 171-186, 1992.

Braun, R. D., and Kroo, I.M.: "Development and

Applications of the Collaborative Optimization

Architecture in a Multidisciplinary Design

Environment," SIAM Journal on Optimization,

Alexandrov, N, and Hussaini, M. Y., (eds), 1996.

Bruhn, E. F.: "Analysis and Design of Flight Vehicle

Structures," Tri-State Offset Co., 1965.

Giles, G. L.: "Equivalent Plate Analysis of Aircraft

Wing Box Structures with General Planform

Geometry," NASA TM-87697, March 1986.

Harris, R. V.: "An Analysis and Correlation of Aircraft

Wave Drag," NASA TMX-947, March 1964.

iSIGHT Designers and Developers Manual, Version

3.1, Engineous Software Inc., Morrisville, North

Carolina, 1998.

MATLAB Manual, the MathWorks, Inc., July 1993,

Version 5.0, 1997, and MATLAB Optimization

Toolbox, User's Guide, 1996.

Raymer, D. P.: Aircraft Design: A Conceptual

Approach, AIAA Education Series, 1992.

Sobieszczanski-Sobieski, J.: "Optimization by

Decomposition: A Step from Hierarchic to Nonhierarchic Systems," 2nu NASA/AF Symposium on

Recent Advances in Multidisciplinary Analysis and

Optimization; Hampton; Virginia, Sept. 1988; NASA

CP-3031, pp.51-78; also in NASA TM 101494.

Sobieszczanski-Sobieski, J.; Agte, J.; and Sandusky.

R.: "Bi-Level Integrated System Synthesis," Proc. 7tfi

AIAA/USA/NASA/ISSMO Symposium on

Multidisciplinary Analysis and Optimization, AIAA

Paper No. 98-4916. Also NASA TM-1998-208715,

Sept. 1998a; (to appear in AIAA J. Jan. 2000)

Sobieszczanski-Sobieski, J.; Agte, J.; and Sandusky,

R.: "Bi-Level Integrated System Synthesis (BLISS)",

NASA/TM-1998-208715, August 1998b.

Wujek, B. A.; Renaud, J.E.; Batill, S. M.; and

Brockman, J.B.: "Concurrent Subspace Optimization

using Design Variable Sharing in a Distributed Design

Environment," Azarm, S. (ed), Proc. Design

Engineering Technical Conf., Advances in Design

Automation, ASME DE Vol. 82, pp.181-188, 1995.

Appendix

The Appendix explains the details of the technique used

to integrate a FORTRAN code into MATLAB. It

describes the steps required for replacing sections of

BLISS code with new FORTRAN analysis.

MATLAB-FORTRAN Integration

There are a few steps required for integrating a

FORTRAN code into the MATLAB environment that

BLISS is currently programmed in. The user must

locate the place in BLISS where a call is made to

analysis that is to be replaced. Then he must examine

the input and output of the new and old analysis to

ensure that the remainder of BLISS is capable of

supplying input the new analysis requires and that the

new analysis produces all the output expected. Then a

pre-processor routine must be created to present the

variable information to the FORTRAN code in an

appropriate format. Finally the data must be harvested

and returned to the MATLAB module from which the

FORTRAN code is being called in a format compatible

with MATLAB.

8
Pre-processing Data Collection

In pre-processing, BLISS must pass the needed

variables to a routine which will manipulate them into

a form that the FORTRAN code will accept. The

programmer must first ensure that the module in which

he is pre-processing has access to the required variables

(i.e., the subsystem must not be using X variables

assigned to other subsystems).

The two FORTRAN codes integrated required

formatted text input files. The easiest way to prepare

these files was to use a FORTRAN subroutine to take

the design variables, configure them in the way needed

to represent the geometry or conditions needed by the

analysis, and write the formatted input file for analysis.

In order to pass variables from the MATLAB

environment to the FORTRAN pre-processing

subroutine, MATLAB's mex-function was invoked. In

this function, a standard gateway subroutine is added to

the FORTRAN pre-processing subroutine. This new

routine collects an array of variables passed in from

MATLAB and assigns them to a FORTRAN array.

The gateway routine sends these variables into the preprocessing subroutine.

In MATLAB, the user compiles the FORTRAN code

including the gateway routine and the pre-processing

routine using the mex command. This creates a mexfile which is treated as a MATLAB function requiring

an input array and an output array. The user then places

his variables for the pre-processing function into an

array and puts this array into the new mex-file. This

sends the variables to the gateway routine which

assigns them to the variables used in the FORTRAN

pre-processing subroutine. The input file is generated

and control returns to MATLAB.

Program Insertion

Having prepared the data for analysis by the program,

the programmer must locate the section of BLISS that

he wishes to upgrade. The previous analysis must be

removed and the code must be placed such that BLISS

will have performed the new analysis and have data

ready for later analysis that the user is not replacing.

Having located the desired calling spot and removed the

replaced analysis, the user simply calls the program

from within BLISS. By previously compiling the

FORTRAN code in question, the user calls the program

by typing .tprogram_name in the BLISS code where

programname is the command that runs the program

from the operating system. The program then processes

the prepared input file and returns to BLISS.

The final step that the programmer must perform in

order for BLISS to carry on its optimization is

harvesting the data produced by the new program.

ELAPS and AWAVE both created output files with

data required for BLISS. There are two basic ways to

collect data produced by FORTRAN codes.

The first way is to use a post-processing technique

similar to that of pre-processing. The user would create

a search algorithm to locate and collect the data from

the output file. This in turn would be harvested by

using the mex-function to create a gateway between the

FORTRAN data collection routine and the BLISS

variables. The programmer would compile the mexfunction gateway routine combined with his data

collection routine, run the new mex-file with an array

prepared to collect the output of the routine, and extract

his data to the array. The programmer would then need

to assign the array variables to the proper variables in

the BLISS code. This technique would be best for

cases where the user did not have access to the source

code of the FORTRAN program.

In cases where the programmer does have the

FORTRAN program's source code and a fair knowledge

of how the program works, he can edit the code to

output the needed results in a format compatible with

MATLAB. If the user can locate the sections of code

that write out the results to the output file, they can

change the code to output to a file with a .m extension.

Files in this format are recognized as MATLAB

programs that can be called without FORTRAN

interaction. By creating .m files with MATLAB

variable assignments corresponding to the data that the

user wishes to collect, the programmer simply runs

these MATLAB files after completion of the main

program call and the output is already in MATLAB

format. This avoids unnecessary data file searching and

reduces MATLAB-FORTRAN interaction. After the

data are collected and assigned to the proper variables in

BLISS, the procedure would then continue.
REPORT DOCUMENTATION PAGE FormApproved

OMB No. 0704-0188

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data

sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other

aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and

Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188),

Washin_lton, DC 20503.

1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED

December 2000 Technical Memorandum

4. TITLE AND SUBTITLE

Advancement of Bi-Level Integrated System Synthesis (BLISS)

6. AUTHOR(S)

Jaroslaw Sobieszczanski-Sobieski, Mark S. Emiley, Jeremy S. Agte, and

Robert R. Sandusky, Jr.

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES)

NASA Langley Research Center

Hampton, VA 23681-2199

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES)

National Aeronautics and Space Administration

Washington, DC 20546-0001

5. FUNDING NUMBERS

WU 509-10-31-03

8. PERFORMING ORGANIZATION

REPORT NUMBER

L-18000

10. SPONSORING/MONITORING

AGENCY REPORTNUMBER

NASA/TM-2000-210305

11. SUPPLEMENTARY NOTES

Sobieski: Langley Research Center; Emiley and Sandusky: George Washington University, JIAFS; Agte: U.S.

Air Force. Paper presented at the 38th AIAA Aerospace Sciences Meeting and Exhibit, January 10-13, 2000,

Reno, NV, Paper No. AIAA 2000-0421.

12a. DISTRIBUTION/AVAILABILITY STATEMENT

Unclassified-Unlimited

Subject Category 05 Distribution: Nonstandard

Availability: NASA CASI (301) 621-0390

12b. DISTRIBUTION CODE

13. ABSTRACT (Maximum 200 words)

Bi-Level Integrated System Synthesis (BLISS) is a method for optimization of an engineering system, e.g., an

aerospace vehicle. BLISS consists of optimizations at the subsystem (module) and system levels to divide the

overall large optimization task into sets of smaller ones that can be executed concurrently. In the initial version

of BLISS that was introduced and documented in previous publications, analysis in the modules was kept at the

early conceptual design level. This paper reports on the next step in the BLISS development in which the

fidelity of the aerodynamic drag and structural stress and displacement analyses were upgraded while the

method's satisfactory convergence rate was retained.

14. SUBJECT TERMS

Parallel processing; Optimization by decomposition; Engineering system design

17. SECURITY CLASSIFICATION

OF REPORT

Unclassified

NSN 7540-01-280-5500

18. SECUHITY CLASSIFICATION

OF THIS PAGE

Unclassified

19. SECURITY CLASSIFICATION

OF ABSTRACT

Unclassified

15. NUMBER OF PAGES

14

16. PRICE CODE

A03

20. LIMITATION

OF ABSTRACT

UL

StanOara Form 298 (Hey. 2-59)

Prescribed by ANSI Std. Z-39-18

298-102