Chapter 5 Robust Performance Tailoring with Tuning - SSL - MIT
Chapter 5 Robust Performance Tailoring with Tuning - SSL - MIT Chapter 5 Robust Performance Tailoring with Tuning - SSL - MIT
defined by the user and indicate the number of designs that must be accepted or rejected, respectively, at a given temperature. When either A = kA or R = kR the temperature is reduced according to a user-defined cooling constant, kC: Tk = Tk−1 (1 − kC) (2.45) and the A and R counters are reset to zero. The search continues in this manner until the termination conditions are met. The termination conditions are chosen by the user based on the specific problem statement. In the PT SCI search the algorithm terminates when kR designs are rejected at three temperatures in a row. Equation 2.44 is called the Boltzman probability, and is the critical feature of the SA algorithm. At high temperatures the probability of accepting a design that increases the cost is high. By randomly accepting less favorable designs the algorithm allows itself to climb out of areas of local minimum and more thoroughly search the space. As the temperature decreases the acceptance probability, pA, becomes lower, and the search eventually converges in the neighborhood of a favorable design. The minimum cost and corresponding design variables are identified from the history of accepted design costs. There is no guarantee that the resulting design is an optimum in the mathematical sense, in fact it most likely is not. However, if the cooling schedule is set appropriately for the given problem the resulting design is a good design and in many cases is better than a locally optimal design resulting from a gradient search. 2.5 Performance Tailoring Results The results of the PT optimizations run on the SCI development model are presented in the following section. Three optimization algorithms, SQP, MC SQP and SA are compared for performance and efficiency. The SQP algorithm is part of the MATLAB optimization toolbox, and the SA algorithm is implemented as outlined in Figure 2-4. At each optimization or search iteration a new SCI state-space model is built 60
and the RMS OPD is computed using Equations 2.13 and 2.14. The SCI model is built in MATLAB as indicated in the algorithm outlined in Figure 2-5. The input to the model generation code is a vector of tailoring design variables, x. These design variables are used to build the appropriate finite element model that includes the truss elements, optics, design masses and ACS model. The eigenvalue problem is solved, resulting in the natural frequencies and mode shapes for this particular design. The rigid body modes are removed and a state-space model, SY S (x) is assembled. Data: tailoring parameters, �x Result: SCI system model, SY S(�x) begin Build finite element model including: truss, optics, design masses and ACS → K(�x),M(�x) Solve eigenvalue problem → Ω(K, M), Φ(K, M) Remove rigid body modes Build state-space model → A(Ω,Z),B(Φ),C(Φ) Build data structure and return model end Figure 2-5: Development model generation algorithm. 2.5.1 Comparison of Optimization Algorithms The model performance is evaluated in the nominal configuration to provide a base- line for comparison with the optimized designs. Building the model with the nominal tailoring parameters listed in Table 2.3 and running the disturbance analysis results in a nominal RMS OPD of 471 µm. The cost, optimal tailoring parameters, compu- tation time and number of function evaluations for each optimized design are listed in Table 2.4. Table 2.4: PT optimization results, J0 = 471 µm. Tailoring Parameters Algorithm J [µm] d1[m] d2[m] m1[kg] m2[kg] time [min] nfun SA 101.57 0.03 0.03 0.3614 2.520 8.74 1619 SQP 100.53 0.03 0.03 0.0 1.934 1.28 29 MC SQP 100.53 0.03 0.03 1.941 0.0 18.1 447 61
- Page 9 and 10: List of Figures 1-1 Timeline of Ori
- Page 11 and 12: List of Tables 1.1 Effect of simula
- Page 13 and 14: Nomenclature Abbreviations ACS atti
- Page 15: dk optimization search direction f
- Page 18 and 19: 1.1 Space-Based Interferometry NASA
- Page 20 and 21: unfettered by the Earth’s atmosph
- Page 22 and 23: the SCI, both the size and flexibil
- Page 24 and 25: maybethatitbecomescertain that the
- Page 26 and 27: Table 1.1: Effect of simulation res
- Page 28 and 29: Table 1.2: Effect of simulation res
- Page 30 and 31: has been found that structural desi
- Page 32 and 33: precision telescope structure for m
- Page 34 and 35: attractive, and more conservative a
- Page 36 and 37: to solve the performance tailoring
- Page 39 and 40: Chapter 2 Performance Tailoring A c
- Page 41 and 42: ometer (SCI). In the following sect
- Page 43 and 44: The equations of motion of the unda
- Page 45 and 46: The frequency response functions fr
- Page 47 and 48: the output covariance matrix, Σz,
- Page 49 and 50: where the subscript indicates the i
- Page 51 and 52: 2.3.3 Design Variables The choice o
- Page 53 and 54: and then, by inspection, the inerti
- Page 55 and 56: algorithms begin at an initial gues
- Page 57 and 58: at least locally optimal, and the s
- Page 59: initial design variable state, x =
- Page 63 and 64: # Designs 25 20 15 10 5 Accepted, b
- Page 65 and 66: does not provide information on why
- Page 67 and 68: energy is distributed almost evenly
- Page 69 and 70: also symmetric as seen in the figur
- Page 71 and 72: Chapter 3 Robust Performance Tailor
- Page 73 and 74: through careful and experienced mod
- Page 75 and 76: described above. However, one can r
- Page 77 and 78: ic, σz(�x, �p), that is depend
- Page 79 and 80: Magnitude, OPD/F x [µm/N] Magnitud
- Page 81 and 82: % Energy 100 90 80 70 60 50 40 30 2
- Page 83 and 84: metric to the cost function. Note,
- Page 85 and 86: tion: ∂hi (z,�x, �pi) ∂�x
- Page 87 and 88: values are chosen from their statis
- Page 89 and 90: Table 3.3: Algorithm performance: a
- Page 91 and 92: Statistical Robustness The statisti
- Page 93 and 94: Performance [µm] 1400 1200 1000 80
- Page 95 and 96: (Figure 3-6(b)). The nominal perfor
- Page 97 and 98: Norm. Cum. Var. [µm 2 ] PSD [µm 2
- Page 99 and 100: energy by mode for easy comparison.
- Page 101 and 102: Y−coordinate [m] Y−coordinate [
- Page 103 and 104: RMS performance, [µm] 400 350 300
- Page 105 and 106: The requirement chosen here is some
- Page 107 and 108: Chapter 4 Dynamic Tuning Robust Per
- Page 109 and 110: on a physical truss. Since tailorin
defined by the user and indicate the number of designs that must be accepted or<br />
rejected, respectively, at a given temperature. When either A = kA or R = kR the<br />
temperature is reduced according to a user-defined cooling constant, kC:<br />
Tk = Tk−1 (1 − kC) (2.45)<br />
and the A and R counters are reset to zero. The search continues in this manner until<br />
the termination conditions are met. The termination conditions are chosen by the<br />
user based on the specific problem statement. In the PT SCI search the algorithm<br />
terminates when kR designs are rejected at three temperatures in a row.<br />
Equation 2.44 is called the Boltzman probability, and is the critical feature of<br />
the SA algorithm. At high temperatures the probability of accepting a design that<br />
increases the cost is high. By randomly accepting less favorable designs the algorithm<br />
allows itself to climb out of areas of local minimum and more thoroughly search the<br />
space. As the temperature decreases the acceptance probability, pA, becomes lower,<br />
and the search eventually converges in the neighborhood of a favorable design. The<br />
minimum cost and corresponding design variables are identified from the history of<br />
accepted design costs. There is no guarantee that the resulting design is an optimum<br />
in the mathematical sense, in fact it most likely is not. However, if the cooling<br />
schedule is set appropriately for the given problem the resulting design is a good<br />
design and in many cases is better than a locally optimal design resulting from a<br />
gradient search.<br />
2.5 <strong>Performance</strong> <strong>Tailoring</strong> Results<br />
The results of the PT optimizations run on the SCI development model are presented<br />
in the following section. Three optimization algorithms, SQP, MC SQP and SA are<br />
compared for performance and efficiency. The SQP algorithm is part of the MATLAB<br />
optimization toolbox, and the SA algorithm is implemented as outlined in Figure 2-4.<br />
At each optimization or search iteration a new SCI state-space model is built<br />
60