STK_TESTCASE_TRUSS3 provides information about the 'truss3' test case CALL: TC = stk_testcase_truss3 () returns a structure TC that describes the 'truss3' test case, borrowed from [1, chapter 9]. This structure contains two fields: * .constants: all the numerical constants for this problem, * .search_domain: an stk_hrect object that specifies the search domain of the optimization problem. TEST CASE OVERVIEW The system considered in this test case is the following 3-bar truss: <--------- D -----------> <--- w ---> ------A==========B==============C------ \_ | __/ ^ \_ | (2) __/ | \_ | __/ L (1) \_ | __/ (3) | \_P_/ v Nodes A, B and C are fixed (pin joints). Node P is submitted to both an horizontal load F1 (e.g., wind) and a vertical load F2 (suspended load). The design variables are the cross-sections a1, a2 and a3 of the three bars, and the horizontal position w of the vertical bar. The quantities of interest are the total volume of the structure, the mechanical (tensile) stress in the bars, and the displacement of P. Various formulations of optimization problems can be considered, depending on which quantities are selected as contraints and objectives, and whether or not uncertainties are taken into account (robust formulations). NUMERICAL CONSTANTS The numerical values borrowed from [1] have been converted to SI units. The fields of TC.constants are: * .D: truss width [m], * .L: length of the vertical bar [m], * .E: Young's modulus [Pa], * .a_min: minimal cross-section [m^2], * .a_max: maximal cross-section [m^2], * .w_min: minimal value of the position of the vertical bar [m], * .w_max: maximal value of the position of the vertical bar [m], * .F1_mean: mean (nominal) value of the horizontal load [N], * .F1_std: standard deviation of the horizontal load [N], * .F2_mean: mean (nominal) value of the vertical load [N] * .F2_std: standard deviation of the vertical load [N]. The standard deviations .F1_std and .F2_std are used in the formulation of robust optimization problems related to this test case [see 1, chap 11]. NUMERICAL FUNCTIONS Two numerical functions are provided to compute the quantities of interest of this test case: * stk_testfun_truss3_vol: computes the total volume of the structure, * stk_testfun_truss3_bb: computes the tensile stress in the bars and the displacement of P. Both functions have the same syntax: V = stk_testfun_truss3_vol (X, CONST) Z = stk_testfun_truss3_bb (X, CONST) where CONST is a structure containing the necessary numerical constants. To use the constants from [1], pass TC.constants as second input argument. Both function accept as first input argument an N x D matrix (or data frame) where D is either 4 or 6: * columns 1--3: cross-section a1, a2 and a3, * column 4: position w of the horizontal bar, * column 5-6 (optional): horizontal and vertical loads F1, F2. The second function is named 'bb' for 'black box', as it plays the role of a (supposedly expensive to evaluate) black box computer model for this test case. The output Z has five columns, corresponding to: * columns 1--2: horizontal and vertical displacement y1, y2 of P, * columns 3--5: tensile stress sigma_j in bars j = 1, 2 and 3. EXAMPLE tc = stk_testcase_truss3 (); n = 5; % Draw 5 points uniformly in the 4D input domain ("design space") xd = stk_sampling_randunif (n, [], tc.search_domain) % Compute the volumes v = stk_testfun_truss3_vol (xd, tc.constants) % Compute displacements and stresses for nominal loads z = stk_testfun_truss3_bb (xd, tc.constants) % Draw loads from normal distributions F = stk_dataframe (zeros (n, 2), {'F1' 'F2'}); F(:, 1) = tc.constants.F1_mean + tc.constants.F1_std * randn (n, 1); F(:, 2) = tc.constants.F2_mean + tc.constants.F2_std * randn (n, 1); % Compute displacements and stresses for the random loads x = [xd F] z = stk_testfun_truss3_bb (x, tc.constants) REFERENCE [1] Indraneel Das, Nonlinear Multicriteria Optimization and Robust Optimality. PhD thesis, Rice University, 1997. [2] Juhani Koski, Defectiveness of weighting method in multicriterion optimization of structures. Int. J. for Numerical Methods in Biomedical Engineering, 1(6):333-337, 1985. See also: stk_testfun_truss3_vol, stk_testfun_truss3_bb
Package: stk