1#@ *: DakotaConfig=HAVE_DOT
2# DAKOTA INPUT FILE: dakota_sbouu3_tbch.in
3#
4# OUU with surrogates: nested containing surrogate
5#
6
7environment,
8	  method_pointer = 'OPTIM'
9
10###########################
11# begin opt specification #
12###########################
13method,
14	id_method = 'OPTIM'
15	model_pointer = 'OPTIM_M'
16        dot_sqp
17	  convergence_tolerance = 1.e-6
18
19model,
20	id_model = 'OPTIM_M'
21	nested
22	  variables_pointer  = 'OPTIM_V'
23	  sub_method_pointer = 'UQ'
24	  optional_interface_pointer  = 'OPTIONAL_I'
25	  optional_interface_responses_pointer = 'OPTIONAL_I_R'
26	  responses_pointer  = 'OPTIM_R'
27	  primary_response_mapping   = 0. 0. 1. 0. 0. 0. 0. 0. 1.
28	  secondary_response_mapping = 0. 0. 0. 1. 3. 0. 0. 0. 0.
29
30variables,
31	id_variables = 'OPTIM_V'
32	continuous_design = 2
33	  initial_point    1.8    1.0
34	  upper_bounds     2.164  4.0
35	  lower_bounds     1.5    0.0
36	  descriptors      'd1'   'd2'
37
38responses,
39# minimize f + p_fail_r1 + p_fail_r3
40# s.t.     g_1/2/3 <= 0
41#          mean_r2 + 3*sigma_r2 <= 1.6e+05
42# NOTE: This specifies the TOTAL RESPONSE for the optimization,
43#       which is a combination of nested & interface responses.
44	id_responses = 'OPTIM_R'
45	objective_functions = 1
46	nonlinear_inequality_constraints = 4
47	nonlinear_inequality_upper_bounds = 0.0 0.0 0.0 1.6e+5
48	numerical_gradients
49	  method_source dakota
50	  interval_type central
51	  fd_gradient_step_size = 1.e-2
52	no_hessians
53
54##########################################
55# begin optional interface specification #
56##########################################
57interface,
58	id_interface = 'OPTIONAL_I'
59  direct
60	  analysis_driver = 'cyl_head'
61
62responses,
63# NOTE: This is the response set from the optional interface.
64	id_responses = 'OPTIONAL_I_R'
65	objective_functions = 1
66	nonlinear_inequality_constraints = 3
67	no_gradients
68	no_hessians
69
70##########################
71# begin UQ specification #
72##########################
73method,
74	id_method = 'UQ'
75	model_pointer = 'UQ_M'
76        sampling,
77	  output quiet
78	  samples = 5000 sample_type lhs
79	  seed = 12347
80#	  fixed_seed					#s1
81	  response_levels = 4.2e+11 1.6e+05 7.5e+05
82	  complementary distribution
83
84model,
85	id_model = 'UQ_M'
86	surrogate global
87	  variables_pointer = 'UQ_V'
88	  responses_pointer = 'UQ_R'
89	  dace_method_pointer = 'DACE'
90	  polynomial quadratic
91
92variables,
93	id_variables = 'UQ_V'
94	continuous_design = 2
95	normal_uncertain = 2
96	  means          =  248.89, 593.33
97	  std_deviations =   12.4,  29.7
98# 3 sigma bounds for DDACE:
99#	  lower_bounds =  211.7,  504.2
100#	  upper_bounds =  286.1,  682.4
101	  descriptors  =  'nuv1'  'nuv2'
102	uniform_uncertain = 2
103	  lower_bounds =  199.3,  474.63
104	  upper_bounds =  298.5,  712.
105	  descriptors  =  'uuv1'  'uuv2'
106	weibull_uncertain = 2
107	  alphas       =    12., 30.
108	  betas        =   250., 590.
109#	  lower_bounds =   225., 540.
110#	  upper_bounds =   275., 640.
111	  descriptors  =  'wuv1'  'wuv2'
112
113responses,
114	id_responses = 'UQ_R'
115	response_functions = 3
116	no_gradients
117	no_hessians
118
119############################
120# begin DACE specification #
121############################
122method,
123	id_method = 'DACE'
124	model_pointer = 'DACE_M'
125        sampling
126	  seed = 12347
127	  sample_type lhs #samples = 68
128#	dace lhs seed = 5
129#	  samples = 50
130
131model,
132	id_model = 'DACE_M'
133	single
134	  variables_pointer = 'DACE_V'
135	  interface_pointer = 'UQ_ACTUAL'
136	  responses_pointer = 'UQ_R'
137
138variables,
139	id_variables = 'DACE_V'
140	active all
141	continuous_design = 2
142	normal_uncertain = 2
143	  means          =  248.89, 593.33
144	  std_deviations =   12.4,  29.7
145# 3 sigma bounds for DDACE:
146#	  lower_bounds =  211.7,  504.2
147#	  upper_bounds =  286.1,  682.4
148	  descriptors  =  'nuv1'  'nuv2'
149	uniform_uncertain = 2
150	  lower_bounds =  199.3,  474.63
151	  upper_bounds =  298.5,  712.
152	  descriptors  =  'uuv1'  'uuv2'
153	weibull_uncertain = 2
154	  alphas       =    12., 30.
155	  betas        =   250., 590.
156#	  lower_bounds =   225., 540.
157#	  upper_bounds =   275., 640.
158	  descriptors  =  'wuv1'  'wuv2'
159
160interface,
161	id_interface = 'UQ_ACTUAL'
162  direct
163#	system asynch evaluation_concurrency = 5
164	  analysis_driver = 'text_book_ouu'
165