HyperStudy

ARSM-based Sequential Optimization and Reliability Assessment (SORA_ARSM)

ARSM-based Sequential Optimization and Reliability Assessment (SORA_ARSM)

Previous topic Next topic Expand/collapse all hidden text  

ARSM-based Sequential Optimization and Reliability Assessment (SORA_ARSM)

Previous topic Next topic JavaScript is required for expanding text JavaScript is required for the print function  

Reliability and robustness based optimization methods require many design evaluations, therefore improving their efficiency is one of the issues. SORA-ARSM attempts to address this issue by using Adaptive Response Surface Method (ARSM). In this process, response surfaces are created and an optimization is carried out on the surfaces. During deterministic optimization and reliability analysis, response surfaces are adaptively updated for increased accuracy.

 

Usability Characteristics

SORA-ARSM is the most efficient of the three RBDO methods available in HyperStudy, but it is also the least accurate.
An extension of SORA_ARSM is implemented in HyperStudy to allow for robust design optimization. Robust design optimization attempts to minimize the objective variance in order to reduce its sensitivity to design variations and consequently increase the design's robustness. The implementation in HyperStudy is based on the use of percentiles for the objective function and is turned on via the Robust Optimization setting in the Specification step.
SORA-ARSM terminates if one of the conditions below are met:
-One of the three convergence criteria are met.
oThe absolute objective change is less than a convergence tolerance value (SOREPS) and there is no constraint violation (GMAX) in the last design.
oThe relative objective change is less than a convergence tolerance value (SOREPS) and there is no constraint violation (GMAX) in the last design.
oThe absolute change and relative change of each input variable (SOREPS) is less than a convergence tolerance value. Also, there is no constraint violation (GMAX) in the last design.
-The maximum number of allowable iterations (MAXDES) is reached.

An exception is when the absolute change and relative change of each input variable is less than this value in the last two consecutive designs or when we have found feasible designs and the best feasible design has not been improved during the last two consecutive iterations. When this occurs, SORA-ARSM will be terminated.

The reliability analysis is carried out by searching for the most probable point (MPP). Issues such as non-uniqueness of the MPP and highly non-linear output response functions can reduce the accuracy of the reliability calculation.

 

The flowchart below illustrates the different phases of the SORA_ARSM process.

Note:The flowchart for SORA_ARSM is the same as SORA flowchart, except that the deterministic optimization and reliability analysis are done based on response surfaces, which is updated with newly evaluated designs.

SORA_ARSM_flowchart

Settings

In the Specifications step, you can change the settings of SORA_ARSM from the following tabs:

hmtoggle_arrow1Settings

In the Settings tab, you can access the settings listed below. Please note that for most applications the default settings work optimally, and you may only need to change the Maximum Iterations and Robust Optimization.

Setting

Default

Range

Description

Maximum Iterations

(MAXDES)

25

> 0

Maximum number of iterations allowed.

Angle convergence tol.

(MPPTOL)

0.25

>0.0

Angle convergence tolerance for inverse MPP search, in unit of degrees. If the angle between the vector of hst2_u (design point in standard normal distribution space) and the negative gradient falls within the tolerance, then inverse MPP search is regarded as converged.

A smaller value favors a higher precision of reliability analysis, but more computational effort is needed.

Robust Optimization

(ROBUST)

Off

Off, On

To decide whether this is robust optimization or not; 0 means that you do not use robust optimization; 1 means that robust optimization is used.

Constraint violation tol.

(GMAX)

0.1

> 0.0

Global maximum allowable percentage constraint violation.  Constraints must not be violated by more than this value in the converged design.

Robust Min %

(ROBMINPV)

95.0

>50,

<100

Define percentile value of robust optimization for minimization objective.

Robust Max %

(ROBMAXPV)

5.0

>0,

<50

Define percentile value of robust optimization for maximization objective.

 

hmtoggle_arrow1More

In the More tab, you can access the setting listed below. Please note that for most applications, the default settings work optimally.

Setting

Default

Range

Description

Termination criteria

(SOREPS)

1.0E-4

>0.0

If the absolute or relative change of the objective value is less than this value, or the absolute or relative change of the input variables is less than this value, and the constraint violation is not larger than this value, then SORA_ARSM will be terminated. Also, there must not be any constraints with an allowable violation that has been exceeded in the last design.

An exception is when the absolute change and relative change of each input variable is less than this value in the last two consecutive designs or when we have found feasible designs and the best feasible design has not been improved during the last two consecutive iterations, SORA_ARSM will be terminated.

sora_arsm2

Where, x is the input variable; f is the objective; n is the number of input variables; k is the current iteration number; arsm_cmax is the maximum constraint violation; arsm_gmazx is the allowable constraint violation; mslm_epsilon is the value of the termination criteria.

sora_arsm1

Move limit friction

(MOVE)

0.15

0.0 < MOVE < 1.0

Move limit fraction.

Smaller values allow more steady convergence (smaller fluctuation of the output response values), but more computational effort could be consumed.  The value will be adaptively updated during the optimization process.

Initial linear move

(MPERT)

1

1 or 3

1

By DV initial.

Initial move = PERT * MOVE * abs(INI)

(Default when initial value of input variable is non-zero)

An exception is that initial move will be set to minimum move if it is less than minimum move.

Minimum move = MINMOVE * (UB-LB) if (UB-LB) is less than 1.

Minimum move = MINMOVE if (UB-LB) is not less than 1 and absolute value of INI is less than 1.

Minimum move = MINMOVE * min((UB-LB),abs(INI)) if (UB-LB) is not less than 1 and absolute value of INI is not less than 1.

3

By DV bounds.

Initial move = PERT * MOVE * (UB-LB)

(Default when initial value of input variable is zero)

INI: Initial input variable value.

LB, UB: Lower and upper bounds on input variable.

Minimal Move Factor

(MINMOVE)

0.1

0.0 < MINMOVE < MOVE

Minimal move factor.

See the usage of MINMOVE in MPERT.

Initial DV perturbation

(PERT)

1.1

≠ 0.0

Initial input variable perturbation value.

Larger values result in a wider spread of the initial bign designs (bign is the number of input variables; the bign designs together with the start design can determine a linear response surface).  ARSM will search the design space more widely.

Use SVD

(SVDRSM)

False

False or True

If SVDRSM is 0, ARSM is terminated in case of soft convergence (when the current design is the same or little change (1.0e-15) to one of the history designs).

If SVDRSM is 1, Singular Value Decomposition is activated to build the response surfaces and the optimization process will be continued.

The above can be formulated as:

sora_arsm3

Where, k is the current iteration number; n is the number of input variables; xj is one of the history designs.

Revision

A

A or B

This parameter is used to help when there is a convergence difficulty. By default, "A" is selected meaning the legacy algorithm.

Constraint threshold

(EPSCON)

1.0e-4

> 0.0

This parameter is used for constraint value calculation.  In general, constraint value is normalized to its bound value.  One exception is that, constraint value is not normalized if its absolute bound value is less than this parameter.  Recommended range is 1.0e-6 ~ 1.0.