1#@ *: DakotaConfig=HAVE_NPSOL
2# DAKOTA Input File: dakota_rbdo_steel_column_trsb2.in
3# Optimization under uncertainty using reliability methods within a
4# fully-analytic sequential RBDO approach.
5
6environment,
7	method_pointer = 'SBLO'
8
9method,
10	id_method = 'SBLO'
11	surrogate_based_local
12	model_pointer = 'OPTIM_M'
13	  approx_method_pointer = 'OPTIM'
14	  max_iterations = 50
15	  soft_convergence_limit = 4
16	  trust_region
17	    initial_size = 0.2
18	    contraction_factor = 0.5
19	    expansion_factor   = 1.50
20	  output verbose
21
22###########################
23# begin opt specification #
24###########################
25method,
26	id_method = 'OPTIM'
27	npsol_sqp
28#	dot_sqp
29	  convergence_tolerance = 1.e-8
30
31model,
32	id_model = 'OPTIM_M'
33	surrogate local taylor_series
34	  variables_pointer = 'OPTIM_V'
35	  responses_pointer = 'OPTIM_R_SURR'
36	  actual_model_pointer = 'OPTIM_TRUTH'
37
38variables,
39	id_variables = 'OPTIM_V'
40	continuous_design = 3
41	  initial_point   300.     20.      300.
42	  lower_bounds    200.     10.      100.
43	  upper_bounds    400.     30.      500.
44	  descriptors     'b'      'd'      'h'
45
46responses,
47# NOTE: This specifies the TOTAL RESPONSE for the optimization,
48#       as computed _from_ the surrogate model.
49	id_responses = 'OPTIM_R_SURR'
50	objective_functions = 1
51	nonlinear_inequality_constraints = 1
52	nonlinear_inequality_upper_bounds = 4000.
53	analytic_gradients
54	no_hessians
55
56##########################
57# begin TS specification #
58##########################
59
60model,
61	id_model = 'OPTIM_TRUTH'
62	nested
63	  variables_pointer  = 'OPTIM_V'
64	  sub_method_pointer = 'UQ'
65	  optional_interface_pointer  = 'OPTIONAL_I'
66	  optional_interface_responses_pointer = 'OPTIONAL_I_R'
67	  responses_pointer  = 'OPTIM_R_TRUTH'
68	  primary_variable_mapping   = 'B'    'D'    'H'
69	  secondary_variable_mapping = 'mean' 'mean' 'mean'
70# minimize p		Cases 0,1,2
71# maximize beta		Cases 3,4,5
72# maximize z		Cases 6,7,8,9,10,11
73# s.t.     Cost <= 4000.
74	  primary_response_mapping = 0. 0.  1.	#s0,#s1,#s2
75#	  primary_response_mapping = 0. 0. -1.	#s3,#s4,#s5,#s6,#s7,#s8,#s9,#s10,#s11
76
77
78responses,
79# NOTE: This specifies the TOTAL RESPONSE from the nested mapping
80#       used for constructing the surrogate.
81	id_responses = 'OPTIM_R_TRUTH'
82	objective_functions = 1
83	nonlinear_inequality_constraints = 1
84	analytic_gradients
85	quasi_hessians sr1
86
87##########################################
88# begin optional interface specification #
89##########################################
90interface,
91	id_interface = 'OPTIONAL_I'
92  direct
93	  analysis_driver = 'steel_column_cost'
94
95responses,
96# NOTE: This is the response set from the optional interface.
97	id_responses = 'OPTIONAL_I_R'
98	objective_functions = 0
99	nonlinear_inequality_constraints = 1
100	analytic_gradients
101	no_hessians
102
103##########################
104# begin UQ specification #
105##########################
106method,
107	id_method = 'UQ'
108	model_pointer = 'UQ_M'
109	local_reliability nip
110	  mpp_search x_taylor_mpp		#s0,#s3,#s6,#s9
111#	  mpp_search u_taylor_mpp		#s1,#s4,#s7,#s10
112#	  mpp_search no_approx			#s2,#s5,#s8,#s11
113	  num_response_levels = 1		#s0,#s1,#s2,#s3,#s4,#s5
114	  response_levels = 0.0			#s0,#s1,#s2,#s3,#s4,#s5
115#	  compute gen_reliabilities		#s3,#s4,#s5
116	  integration second_order
117#	  num_probability_levels = 1		#s6,#s7,#s8
118#	  probability_levels = .001373		#s6,#s7,#s8
119#	  num_gen_reliability_levels = 1	#s9,#s10,#s11
120#	  gen_reliability_levels = 2.995	#s9,#s10,#s11
121	  cumulative distribution
122
123model,
124	id_model = 'UQ_M'
125	single
126	  variables_pointer = 'UQ_V'
127	  interface_pointer = 'UQ_I'
128	  responses_pointer = 'UQ_R'
129
130variables,
131	id_variables = 'UQ_V'
132	normal_uncertain = 2
133	  means           =  30. 500000.
134	  std_deviations  =  10.  50000.
135	  descriptors     = 'F0'  'P1'
136	lognormal_uncertain = 4
137	  means          = 300.  20. 300. 400.
138	  std_deviations =   3.   2.   5.  35.
139	  descriptors    =  'B'  'D'  'H' 'Fs'
140	gumbel_uncertain = 2
141	  alphas         = 1.4250554e-5 1.4250554e-5
142	  betas          = 559496.31    559496.31
143	  descriptors    =  'P2'    'P3'
144	weibull_uncertain = 1
145	  alphas          =     5.7974
146	  betas           =	22679.4777
147	  descriptors     =    'E'
148
149interface,
150	id_interface = 'UQ_I'
151  direct
152	  analysis_driver = 'steel_column_perf'
153	  deactivate evaluation_cache restart_file
154
155responses,
156	id_responses = 'UQ_R'
157	response_functions = 1
158	analytic_gradients
159	analytic_hessians
160