1#@ s*: Label=FastTest
2#@ *: DakotaConfig=HAVE_DOT
3# DAKOTA Input File: dakota_ouu1_cantilever.in
4# Optimization under uncertainty
5# using the cantilever test function.
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#	npsol_sqp
18	  convergence_tolerance = 1.e-6
19	scaling
20
21model,
22	id_model = 'OPTIM_M'
23	nested
24	  variables_pointer  = 'OPTIM_V'
25	  sub_method_pointer = 'UQ'
26	  responses_pointer  = 'OPTIM_R'
27# Approach 1: constrain the probabilities < .00135
28#	  primary_response_mapping   = 1. 0. 0. 0. 0. 0. 0. 0. 0.
29#	  secondary_response_mapping = 0. 0. 0. 0. 0. 1. 0. 0. 0.
30#                                      0. 0. 0. 0. 0. 0. 0. 0. 1.
31# Approach 2 (better for sampling): constrain mean + 3 sigma
32	  primary_response_mapping   = 1. 0. 0. 0. 0. 0.
33	  secondary_response_mapping = 0. 0. 1. 3. 0. 0.
34                                       0. 0. 0. 0. 1. 3.
35
36variables,
37	id_variables = 'OPTIM_V'
38	continuous_design = 2
39	  initial_point    2.5    2.5
40	  upper_bounds    10.0   10.0
41	  lower_bounds     1.0    1.0
42	  descriptors 'w' 't'
43	scale_type = 'auto'
44
45responses,
46# minimize mean Weight
47# s.t.     p_fail_S <= .00135 (beta = 3.0) [Approach 1]
48#          p_fail_D <= .00135 (beta = 3.0) [Approach 1]
49# or
50#          mean_S + 3 sigma_S <= 0         [Approach 2]
51#          mean_D + 3 sigma_D <= 0         [Approach 2]
52#
53# NOTE: This specifies the TOTAL RESPONSE for the optimization,
54#       which is a combination of nested & interface responses.
55	id_responses = 'OPTIM_R'
56	objective_functions = 1
57	nonlinear_inequality_constraints = 2
58# Approach 1 allowables
59#	nonlinear_inequality_upper_bounds = .00135 .00135
60# Approach 2 uses default upper bounds (mean + 3 sigma < 0.0)
61	numerical_gradients
62	  method_source dakota
63	  interval_type central
64	  fd_gradient_step_size = 1.e-2
65	no_hessians
66
67##########################
68# begin UQ specification #
69##########################
70method,
71	id_method = 'UQ'
72	model_pointer = 'UQ_M'
73        sampling
74	  samples = 50 sample_type lhs
75	  seed = 12347
76#	  fixed_seed					#s1
77	  output quiet
78
79model,
80	id_model = 'UQ_M'
81	single
82	  variables_pointer = 'UQ_V'
83	  interface_pointer = 'UQ_I'
84	  responses_pointer = 'UQ_R'
85
86variables,
87	id_variables = 'UQ_V'
88	continuous_design = 2
89	normal_uncertain = 4
90	  means             =  40000. 29.E+6 500. 1000.
91	  std_deviations    =  2000. 1.45E+6 100. 100.
92	  descriptors       =  'R' 'E' 'X' 'Y'
93
94interface,
95	id_interface = 'UQ_I'
96	direct
97	  analysis_driver = 'cantilever'
98	  deactivate evaluation_cache restart_file
99
100responses,
101	id_responses = 'UQ_R'
102	response_functions = 3
103	no_gradients
104	no_hessians
105