Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Solving, analyzing, manufacturing, and experimental testing of thickness distribution for a cycloid-like variable curvature mirror

Open Access Open Access

Abstract

A cycloid-like variable curvature mirror (VCM) for zoom-imaging systems was investigated. An analytical-deformation solution to a thin-elastic plate with a cycloid-like thickness distribution and simply supported boundary condition under uniform pressure was found using a small parameter method. The finite-element analysis of the thin-elastic plate and designed VCM showed a good correlation with the analytical solution. The VCM was manufactured and polished to the initial shape with a root mean square (RMS) of 1/80λ. Finally, with air-pressure-based actuation testing under 0.07 MPa, the VCM deforms approximately 36.89 µm and maintains the RMS surface performance of 1/10λ, 1/40λ with and without spherical aberrations, respectively.

© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Variable curvature mirrors (VCMs) are advanced active optic components whose radius of curvature can be adjusted as needed by an external force to realize zoom imaging without moving elements and thereby replace a conventional mechanical zoom system [16], thermal lens compensation [710], and interferometric telescopes for pupil size stabilization upstream of the beam recombination [1113]. The earliest VCM that can be searched dates back to 1972, which was designed and fabricated at the Hebrew University of Jerusalem using a solid-state laser [14]. In addition, the former Soviet Union [8,15,16], France [1113], Germany [7,17], the US [1820], Japan [21,22], and China [10,2325] have also conducted extensive research on VCMs for different applications, forms, and actuation methods. Initially, VCMs were designed for use in solid-state lasers to improve the quality of the output beam. With such development, more and more applications have gradually been seen. During its development, various forms and actuation methods have been invented for various purposes and conditions. Most of them are of uniform thickness with a uniform actuation. With the gradual deepening of research, the benefits of a VCM with a variable thickness distribution have been discovered, and this type of VCM has been increasingly studied, such as by Kvantovaya Electron Co. [16], Lawrence Berkeley Laboratory [18], Sumitomo Electric Industries [22], and Laboratoire d’Astrophysique de Marseille-LAM, whose study on VCMs with such distributions has been considered as the most systematic cases [1113]. In G. R. Lemaitre’s book Astronomical Optics and Elasticity Theory, several thickness distributions are presented which can achieve an ideal curvature radius variation during actuation, such as in exponential and cycloid-like types [26]. In our previous paper [27], we worked out and verified a simplified analytical deformation solution of an exponential thin elastic plate under a uniform load with a simply supported edge condition using a small-parameter method based on the Theory of plates and shells of Timoshenko [28]. As in [26], the cycloid-like thickness distribution, which is in the form of $h = {h_0}{\left( {1 - \frac{{{r^2}}}{{{a^2}}}} \right)^{1/3}}$, where h is the thickness varying with the radius, ${h_0}$ is the central thickness, a is the radius of the plate, r is the radius, is also a common form of a VCM. In this research, a systematic study of an elastic plate and a VCM with such form thickness was conducted. This VCM is expected to be used in zoom imaging systems in the future. Because the VCM under study has the thickness of a cycloid-like form, to analyze the deformation characteristics of cycloid-like VCMs, the first step is to obtain the analytical solution of the corresponding thin elastic plate, as in our previous study. In this study, we developed a simplified analytical solution of a cycloid-like elastic plate by applying the small-parameter method [29]. In addition, manufacturing, optical polishing, and air pressure procedures have also been conducted in later studies.

2. Analytical study

According to the theory of plates and shells [28], the curvature variation process can be physically modeled as a thin circular plate deformation caused by a uniform load with variable thickness and under a simply supported boundary condition, as shown in Fig. 1.

 figure: Fig. 1.

Fig. 1. Physical model of simply supported elastic thin plate with variable thickness

Download Full Size | PDF

Whereas the thin circular plate is under a simply supported boundary condition with a uniform load, the stress condition and moment equilibrium of the plate are as shown in Fig. 2 [29].

 figure: Fig. 2.

Fig. 2. Cross-section stress analysis of physical model of simply supported elastic thin plate

Download Full Size | PDF

The corresponding moment equilibrium equation is as follows:

$${M_r} + \frac{{d{M_r}}}{{dr}}r - {M_t} + Qr = 0,$$
where r represents the radial coordinate, $dr$ is the infinitesimal of r, ${M_r}$ is the radial bending moment per unit length, ${\textrm{M}_t}$ is the tangential bending moment per unit length, and Q is the shearing force per unit length.

The equations of ${M_r}$, ${M_t}$, and Q are

$${M_r} ={-} D\left( {\frac{{{d^2}w}}{{d{r^2}}} + \frac{\mu }{r}\frac{{dw}}{{dr}}} \right),$$
$${M_t} ={-} \textrm{D}\left( {\frac{1}{r}\frac{{dw}}{{dr}} + \mu \frac{{{d^2}w}}{{d{r^2}}}} \right),$$
$$Q = \frac{1}{{2\pi r}}\mathop \smallint \nolimits_0^r q2\pi rdr,$$
where w is the vertical deflection, $D = \frac{{E{h^3}}}{{12({1 - {\mu^2}} )}}$ is the rigidity of the plate, h is the thickness distribution, $\mu $ is the Poisson’s ratio, and q is the load per unit area. Applying Eqs. (2)–(4) into Eq. (1), we obtain Eq. (5):
$$D\frac{d}{{dr}}\frac{1}{r}\frac{d}{{dr}}\left( {r\frac{{dw}}{{dr}}} \right) + \frac{{dD}}{{dr}}\left( {\frac{{{d^2}w}}{{d{r^2}}} + \frac{\mu }{r}\frac{{dw}}{{dr}}} \right) = \frac{1}{r}\mathop \smallint \nolimits_0^r qrdr,$$
Because the load is uniform, Eq. (5) can also be expressed as following:
$$D\frac{d}{{dr}}\frac{1}{r}\frac{d}{{dr}}\left( {r\frac{{dw}}{{dr}}} \right) + \frac{{dD}}{{dr}}\left( {\frac{{{d^2}w}}{{d{r^2}}} + \frac{\mu }{r}\frac{{dw}}{{dr}}} \right) = \frac{1}{2}qr.$$
Equation (6) represents the universal form of a thin elastic plate under a uniform load. When the thickness distribution h is varied with r, the rigidity of the plate is also changed. In this study, to make the analytical solution more suitable for a greater thickness distribution, the thickness distribution of the plate is expressed in the form of $h = {h_0}{\left( {1 - k\frac{{{r^2}}}{{{a^2}}}} \right)^{1/3}}$[26]. Different values of k represent different thickness distribution modes. To solve the deflection w of the plate, we firstly introduce the following substitutions:
$$x = {({r/a} )^2},\; \; y = \frac{w}{{{h_0}}},\; \; Q = \frac{{3{a^4}({1 - {\mu^2}} )}}{{4Eh_0^4}}q.$$
By combining Eq. (7) into Eq. (6) we obtain
$$\frac{{{d^2}}}{{d{x^2}}}x\frac{{dy}}{{dx}} - \frac{k}{{2({1 - kx} )}}\left[ {2x\frac{{{d^2}y}}{{d{x^2}}} + ({1 + \mu } )\frac{{dy}}{{dx}}} \right] = \frac{Q}{{1 - kx}}.$$
In the case of a plate that is simply supported along the edge, the boundary condition is
$${(w )_{x = 1}} = 0,\; \; {({{M_r}} )_{x = 1}} = 0.$$
Equation (9) can be expressed as $y,\; x$ in Eq. (10) in the following step by solving the differential equation to obtain the corresponding deflection results w of the plate.
$$\left\{ {\begin{array}{c} {y = 0}\\ {2\frac{{{d^2}y}}{{d{x^2}}} + ({1 + \mu } )\frac{{dy}}{{dx}} = 0} \end{array}} \right.\; \; \; \; \; for\; x = 1.$$
We chose a small-parameter method to solve Eq. (8) [29].

First, we set $\varepsilon = k$, and y is developed by the power series for $\varepsilon $ as in Eq. (11),

$$y = {y_0}(x )+ {y_1}(x )\varepsilon + {y_2}(x ){\varepsilon ^2} + {y_3}(x ){\varepsilon ^3} + {y_4}(x ){\varepsilon ^4} + \ldots ,$$
where ${y_0},{y_1},{y_2} \ldots $ are functions of x, respectively. In addition, the right part of the equation, $h = {h_0}{({1 - kx} )^{1/3}}$, is developed by the power series in order of $\varepsilon $,
$$h = {h_0}{({1 - kx} )^{1/3}} = {h_0}\left( {1 - \frac{1}{3}\varepsilon x - \frac{1}{9}{\varepsilon^2}{x^2} - \frac{5}{{81}}{\varepsilon^3}{x^3} - \frac{{10}}{{243}}{\varepsilon^4}{x^4} - \ldots } \right).$$
Second, by solving the differential equations corresponding to the order of $\varepsilon $, we obtain the analytical expressions for ${y_0},{y_1},{y_2},{y_3},{y_4}$. Thus we present five exact solutions in Eqs. (15), (18), (21), (24), and (27).

0) For ${y_0}$, the zero-order equation and the boundary conditions are

$$\frac{{{d^2}}}{{d{x^2}}}x\frac{{d{y_0}}}{{dx}} = Q$$
$${y_0}(1 )= 0,\; \; 2\frac{{{d^2}{y_0}(1 )}}{{d{x^2}}} + ({1 + \mu } )\frac{{d{y_0}(1 )}}{{dx}} = 0$$
The solution is,
$${y_0} = \frac{Q}{4}\left[ {{x^2} - \frac{{2({3 + \mu } )}}{{({1 + \mu } )}}x + \frac{{({5 + \mu } )}}{{({1 + \mu } )}}} \right].$$

1) For ${y_1}$, the first-order equation and the boundary conditions are

$$\frac{{{d^2}}}{{d{x^2}}}x\frac{{d{y_1}}}{{dx}} = \left[ {x\frac{{{d^2}{y_0}}}{{d{x^2}}} + \frac{1}{2}({1 + \mu } )\frac{{d{y_0}}}{{dx}} + Qx} \right],$$
$${y_1}(1 )= 0,\; \; 2\frac{{{d^2}{y_1}(1 )}}{{d{x^2}}} + ({1 + \mu } )\frac{{d{y_1}(1 )}}{{dx}} = 0.$$
The solution is
$${y_1} = \frac{Q}{{144}}\left[ {2({7 + \mu } ){x^3} - 9({3 + \mu } ){x^2} + \frac{{12({ - 4 + 3\mu + {\mu^2}} )}}{{({1 + \mu } )}}x + \frac{{({61 - 16\mu - 5{\mu^2}} )}}{{({1 + \mu } )}}} \right].$$

2) For ${y_2}$, the second-order equation and the boundary conditions are

$$\frac{{{d^2}}}{{d{x^2}}}x\frac{{d{y_2}}}{{dx}} = \left[ {x\frac{{{d^2}{y_1}}}{{d{x^2}}} + {x^2}\frac{{{d^2}{y_0}}}{{d{x^2}}} + \frac{1}{2}({1 + \mu } )\left( {\frac{{d{y_1}}}{{dx}} + \frac{{d{y_0}}}{{dx}}x} \right) + Q{x^2}} \right],$$
$${y_2}(1 )= 0,\; \; 2\frac{{{d^2}{y_2}(1 )}}{{d{x^2}}} + ({1 + \mu } )\frac{{d{y_2}(1 )}}{{dx}} = 0.$$
The solution is
$$\begin{aligned}{y_2} = \frac{Q}{{2304}}\left[ \vphantom{\frac{Q}{{2304}}}({119 + 24\mu + {\mu^2}} ){x^4} + 8({ - 21 - 10\mu - {\mu^2}} ){x^3} + 24({ - 4 + 3\mu + {\mu^2}} ){x^2}\right.\\ \left.- \frac{{4({59 - 79\mu + 13{\mu^2} + 7{\mu^3}} )}}{{({1 + \mu } )}}x + \frac{{({381 - 187\mu + 19{\mu^2} + 11{\mu^3}} )}}{{({1 + \mu } )}} \right]\end{aligned}$$

3) For ${y_3}$, the third-order equation and the boundary conditions are

$$\frac{{{d^2}}}{{d{x^2}}}x\frac{{d{y_3}}}{{dx}} = \left[ {x\frac{{{d^2}{y_2}}}{{d{x^2}}} + {x^2}\frac{{{d^2}{y_1}}}{{d{x^2}}} + {x^3}\frac{{{d^2}{y_0}}}{{d{x^2}}} + \frac{1}{2}({1 + \mu } )\left( {\frac{{d{y_2}}}{{dx}} + \frac{{d{y_1}}}{{dx}}x + \frac{{d{y_0}}}{{dx}}{x^2}} \right) + Q{x^3}} \right],$$
$${y_3}(1 )= 0,\; \; 2\frac{{{d^2}{y_3}(1 )}}{{d{x^2}}} + ({1 + \mu } )\frac{{d{y_3}(1 )}}{{dx}} = 0$$
The solution is
$$\begin{aligned} {y_3} = \frac{Q}{{691200}}\left[\vphantom{\frac{Q}{{691200}}} 6({3689 + 863\mu + 55{u^2} + {\mu^3}} ){x^5} + 75({ - 357 - 191\mu - 27{\mu^2} - {\mu^3}} ){x^4} +\right. \\ 400({ - 28 + 17\mu + 10{\mu^2} + {\mu^3}} ){x^3} - 150({59 - 79\mu + 13{\mu^2} + 7{\mu^3}} ){x^2} - \\ \left. \frac{{30({841 - 1424\mu + 638{\mu^2} - 16{\mu^3} - 39{\mu^4}} )}}{{({1 + \mu } )}}x - \frac{{({ - 49921 + 27532\mu - 9282{\mu^2} + 116{\mu^3} + 451{\mu^4}} )}}{{({1 + \mu } )}} \right] \end{aligned}$$

4) For ${y_4}$, the fourth-order equation and the boundary conditions are

$$\begin{aligned}\frac{{{d^2}}}{{d{x^2}}}x\frac{{d{y_4}}}{{dx}} = \left[ x\frac{{{d^2}{y_3}}}{{d{x^2}}} + {x^2}\frac{{{d^2}{y_2}}}{{d{x^2}}} + {x^3}\frac{{{d^2}{y_1}}}{{d{x^2}}} + {x^4}\frac{{{d^2}{y_0}}}{{d{x^2}}} + \frac{1}{2}({1 + \mu } )\left( {\frac{{d{y_3}}}{{dx}} + \frac{{d{y_2}}}{{dx}}x} \right.\right.\\ \left.\left.+ {\frac{{d{y_1}}}{{dx}}{x^2} + \frac{{d{y_3}}}{{dx}}{x^3}} \right) + Q{x^4}\right]\end{aligned}$$
$${y_4}(1 )= 0,\; \; 2\frac{{{d^2}{y_4}(1 )}}{{d{x^2}}} + ({1 + \mu } )\frac{{d{y_4}(1 )}}{{dx}} = 0$$
The solution is
$$\begin{aligned} {y_4} = \frac{Q}{8294400}\left[ (180761 + 45976\mu + 3558{u^2} + 104{\mu^3} + {\mu^4} ){x^6} - 18(11067 + 6278\mu + \right.\\ 1028{\mu^2} + 58{\mu^3} + {\mu^4}){x^5} + 150({ - 476 + 261\mu + 187{\mu^2} + 27{\mu^3} + {\mu^4}} ){x^4}+ \\ 100({ - 413 + 494\mu - 12{\mu^2} - 62{\mu^3} - 7{\mu^4}} ){x^3} - 45({841 - 1424\mu + 638{\mu^2} - 16{\mu^3} - 39{\mu^4}} ){x^2}\\ - \frac{{18({6677 - 12541\mu + 7838{\mu^2} - 1966{\mu^3} - 115{\mu^4} + 107{\mu^5}} )}}{{({1 + \mu } )}}x - \\ \left.\frac{{6({ - 48196 + 23725\mu - 12048{\mu^2} + 2702{\mu^3} + 148{\mu^4} - 123{\mu^5}} )}}{{({1 + \mu } )}} \right] \end{aligned}$$
In [19], it is shown that with an increase in the order of y, the influence decreases, and when it is higher than 4th order, the influence can be neglected as it is being too small. Therefore, we solve this only up to the 4th order. The overall deformation result is given by Eq. (11) after incorporating Eqs. (27), (24), (21), and (18).

Then, the central deflection (x = 0), which is calculated in Eq. (11), can be developed as follows:

$$\begin{aligned} y(0 ) = Q\left [ \frac{1}{4}\frac{{({5 + \mu } )}}{{({1 + \mu } )}} + \frac{1}{{144}}\frac{{({61 - 16\mu - 5{\mu^2}} )}}{{({1 + \mu } )}}\varepsilon + \frac{1}{{2304}}\frac{{({381 - 187\mu + 19{\mu^2} + 11{\mu^3}} )}}{{({1 + \mu } )}}{\varepsilon^2}\right. \\ - \frac{1}{{691200}}\frac{{({ - 49921 + 27532\mu - 9282{\mu^2} + 116{\mu^3} + 451{\mu^4}} )}}{{({1 + \mu } )}}{\varepsilon^3} - \\ \left.\frac{1}{{8294400}}\frac{{6({ - 48196 + 23725\mu - 12048{\mu^2} + 2702{\mu^3} + 148{\mu^4} - 123{\mu^5}} )}}{{({1 + \mu } )}}{\varepsilon^4}\; \; \right ] \end{aligned}$$

3. Finite element analysis (FEA)

To verify the accuracy of the analytical solution in analyzing the deformation of the elastic plate under different pressures, we built an elastic plate model by CATIA, conducted a finite element analysis (FEA) using MSC Patran/Nastran and compared the FEA results of the elastic plate model to the value calculated using Eq. (11). First, we conduct an FEA of the elastic thin plate model with a simply supported edge condition and pressure on the other side, as shown in Fig. 1, and extract its deformation values along the radial direction. We then compared it to the values calculated by using our analytical solution. In [27], we chose the verification parameters of the elastic plate according to the designed VCM actuated using a uniform air pressure with a closed cavity structure. Therefore, in this study, during the verification process, we use the parameters in [27] as a reference to verify the correctness of the solution, and we obtain a thickness distribution of $h = {h_0}{({1 - kx} )^{1/3}}$. Therefore, the diameter of an elastic plate is $2a = 134.5\; \textrm{mm}$, the central thickness is ${h_0} = 8.5\; \textrm{mm}$, and the material is ultralumin ($E = 72\textrm{GPa},\; \mu = 0.33$). Based on optimization, the parameters of the VCM are as follows: In addition, we choose $k = 0.9$ as the lower surface thickness distribution coefficient, and the thickness expression is $h = 8.5{({1 - 0.9x} )^{1/3}}$. The structure of the elastic plate model with several of the main parameters is shown in Fig. 3, and the VCM structure is shown in Fig. 4, with a collar thickness of $c = 0.62\; mm$.

 figure: Fig. 3.

Fig. 3. 1/2 cross section of elastic thin plate structure parameter

Download Full Size | PDF

 figure: Fig. 4.

Fig. 4. VCM profile with main parameters

Download Full Size | PDF

4. Deformation comparison of FEA and analytical model

During the FEA process, we set the elastic thin plate under a simply supported edge condition and apply 0.1 MPa at its surface. We then obtain the deflection values, as shown in Table 1, and plot the deflection curve along its radius. A comparison between the FEA results of the elastic thin plate and the analytical solution is shown in Fig. 5. From Fig. 5, the difference between the two curves is extremely small, and the deformations show a good consistency. By subtracting the normalized deflection between the values of the two curves, we extract the difference between them, as shown in Fig. 6. The maximum difference is 0.36 µm at the central and edge parts, and the normalized residues value is almost 0.

 figure: Fig. 5.

Fig. 5. FEA and analytical deflection curves of elastic thin plate at 0.1 MPa

Download Full Size | PDF

 figure: Fig. 6.

Fig. 6. FEA and analytical curves of subtraction of normalized deflection

Download Full Size | PDF

Tables Icon

Table 1. Comparison of FEA and analytical deformation values of elastic plate along the radius under 0.1 MPa

Next, compare the central deflection values obtained through the FEA of the designed VCM and our analytical solution under a pressure of 0.01 to 0.1 MPa, as shown in Fig. 7. From Fig. 7, we can see that the overall central deflection values of the mirror are smaller than those of the thin elastic plate, which does not have a collar structure. In addition, we abstracted the deflection values of the FEA and analytical solution along the radius direction under 0.03, 0.06, and 0.09 MPa and plot them in Fig. 8.

 figure: Fig. 7.

Fig. 7. Comparison of central deflection values given by Eq. (28) and finite element analysis

Download Full Size | PDF

 figure: Fig. 8.

Fig. 8. Comparison of surface deformation along radius direction between values given by Eq. (11) and values obtained through FEA

Download Full Size | PDF

From the comparison between the FEA and the analytical results shown in Figs. 7 and 8, we can see that the deflection curves obtained by these two methods are consistent, although the two values are not exactly the same. By adding a thin collar to the VCM structure, a deflection can be achieved by using an air pressure actuation. Meanwhile, the existence of a thin collar adds a semi-clamp effect to the mirror boundary condition instead of the simply supported edge condition [11]. This semi-clamp boundary condition reduces the deformation of the VCM compared to the elastic plate model, which has also been discussed in our previous study [27].

5. Air pressure actuation experiment

5.1 Mirror manufacturing and actuation system building

The prototype mirror is manufactured using ultralumin LC4 with the parameters referred to in the FEA part, and there is a difference between the FEA and the analytical part in that the curvature radius of the mirror is 2807 mm. The overall manufacturing steps are as follows: a) obtain a blank using a CNC machining center, b) use single-point diamond turning to improve the accuracy of the mirror surface figure by better than λ/3, c) after the mechanical process, apply several thermal treatments to release the stress introduced by the mechanical process, d) improve the mirror’s surface by up to λ/50 through classical polishing, e) use ion-beam figuring (IBF) to improve the accuracy of the surface figure by up to λ/80. The mirror after applying IBF and the surface testing results are shown in Fig. 9.

 figure: Fig. 9.

Fig. 9. (a) VCM after ion beam figuring. (b) Surface figure accuracy after IBF.

Download Full Size | PDF

For the actuation system, we have a pressure sensor, proportional valve, switch valve, micro pump, and integrated processor. The maximum pressure provided by the micropump is 0.1 MPa, and the control accuracy is 0.001 MPa. The air pressure actuation system structure is shown in Fig. 10, and the major equipment, pressure sensor, integrated processor, and micro pump are shown in Fig. 11.

 figure: Fig. 10.

Fig. 10. Air pressure actuation system

Download Full Size | PDF

 figure: Fig. 11.

Fig. 11. Photographs of pressure sensor, integrated processor, and micro air pump

Download Full Size | PDF

5.2 Analysis of actuation experiment and actuation results

After all of the apparatuses are assembled and connected together, the experimental testing system is as shown in Fig. 12. During the testing process, the surface performance of the VCM was measured using a ZYGO interferometer with a 6-inch aperture.

 figure: Fig. 12.

Fig. 12. Assembled performance testing system of VCM

Download Full Size | PDF

During the actuation process, we tested the performance of a mirror with a 120-mm diameter surface under pressure values of 0.02, 0.036, 0.053, and 0.070 MPa, and the accuracies of the corresponding surface figure are as shown in Fig. 13.

 figure: Fig. 13.

Fig. 13. Accuracy of VCM surface figure under pressure: (a) 0.02, (b) 0.036, (c) 0.053, and (d) 0.07 MPa (for a demonstration of the actuation effect, see Visualization 1, where the pressure used is 0.04 MPa)

Download Full Size | PDF

Through a decomposition of the performance of the mirror surface, the variation in the Zernike coefficients along with the pressure increase, particularly the Astigmatism 0/90, Astigmatism ±45, X coma, Y Coma, and Sphere 3 are as shown in Fig. 14. From Fig. 14, the coefficient values of Astigmatism 0/90, Astigmatism ±45, and Y Coma have nearly no change during the pressurization process. X Coma showed a slight increase. Only Sphere 3 increases sharply with an increase in the pressure value. Furthermore, the overall surface aberration was analyzed, as shown in Fig. 15. It is clear that the spherical aberration is also the main part after each actuation, when the others only change slightly during this process.

 figure: Fig. 14.

Fig. 14. Zernike coefficient variation with pressure

Download Full Size | PDF

 figure: Fig. 15.

Fig. 15. RMS evolution with pressure

Download Full Size | PDF

Through the analysis of the actuation testing results, the surface figure accuracy degrades to approximately λ/10 at 0.07 MPa, and the main aberration of this degradation is a spherical aberration, as values shown in Table 2. Based on this, the spherical aberration was removed from each test result. We found that although the surface performance slightly decreases with an increase in the pressure value, the surface performance of the VCM is still as high as λ/40 under a pressure of 0.07 MPa, as shown in Fig. 16. Therefore, a spherical aberration is the main aspect that needs to be removed in the following process of mirror research.

 figure: Fig. 16.

Fig. 16. Accuracy of VCM surface figure under (a) 0.02, (b) 0.036, (c) 0.053, (d) 0.070 MPa after removing the spherical aberration

Download Full Size | PDF

Tables Icon

Table 2. Experimental results of VCM surface figure accuracy and corresponding deformation

After the air pressure actuation, the deformation values were also analyzed. We found that the actuation deformation of the VCM under each air pressure value is larger than that of the VCM in Section 4. Based on the analysis, we found that this is because the tested VCM has an original curvature radius on its reflected surface of 2807 mm. By contrast, the analytical and FEA models are flat, and the other side obeys the thickness distribution. As a result, the thickness from the central part to the edge of the VCM is smaller than that of the flat surface model. It is then thinner than the overall model, except for the central thickness. Under the same air pressure value, the deformation is larger than that of the analytical model. Accordingly, a VCM model with an original curvature radius was built and analyzed using an FEA.

From the comparison, although the deformation results between the FEA model and the real fabrication model, which have the same theoretical structural parameters, are extremely close to each other, they are still different. The testing deformation of the mirror under 0.07 MPa is 36.89 µm, which is equal to that of a VCM model with a central thickness of 8.2-mm, as shown in Fig. 17. As is well known, the most influential reasons for a mirror deformation are the central thickness and the thickness of the collar, and thus there must be some operation that can influence these parameters during the manufacturing and polishing process. After examining the entire manufacturing process, we found that there are three processes that can reduce the thickness. First, the blank manufacturing process was not as accurate as the required parameter. Second, as before the classical polishing process, diamond turning was applied to the blank surface to improve the accuracy of the original surface. Third, the classical polishing process also contributes to its reduction, although it is extremely small. Because the mirror structure is not tested before it is sealed, its thickness cannot be verified. However, the thickness of its surface and collar could be two possible reasons.

 figure: Fig. 17.

Fig. 17. Comparison of VCM deformation between FEA and real air pressure actuation (with original curvature radius)

Download Full Size | PDF

6. Conclusion

In the first part of this paper, an analytical solution to a thin elastic plate with a cycloid-like thickness distribution with simply supported edge conditions and a uniform pressure was described by introducing an approximation of the small parameter method. The analytical solution was in a simple form for analysis. In the second part, the deformation of the analytical results of a thin elastic plate was compared to the FEA results of both the plate and the VCM. The analytical and FEA results for the thin elastic plate are almost the same. Therefore, the accuracy of the analytical solution is verified indirectly. In addition, the analytical results are consistent with the FEA results of the VCM, although the results are not exactly the same because of the difference in the real structures between the VCM and the elastic plate. Furthermore, in the following section, a VCM prototype manufactured and tested through the air pressure actuation is described. Deflection and surface performance values were collected and analyzed. The deflection and surface performance of a VCM under 0.07 MPa are 36.89 µm and 0.099 λ (with a spherical aberration), respectively. In this VCM, the main aberration that contributes to the surface performance deterioration is a spherical aberration. After removing the spherical aberration under each pressure, the surface performance of the VCM is still extremely good even at 0.07 MPa. Its value was approximately 0.024 λ. Therefore, the elimination of the spherical aberration of the VCM with a variable thickness distribution during the actuation process will be the focus of the following study. However, as shown in the experimental figures, the volume of the entire testing system in our lab is quite large. Considering its practicality in image systems, the devices must be further improved in terms of volume and integration to adapt to the system requirements.

Funding

West Light Foundation of the Chinese Academy of Sciences (XAB2015A09); Research Fund for Young Star of Science and Technology in Shaanxi Province (2016KJXX-08); China Scholarship Council.

Acknowledgments

The authors thank Prof. Emmanuel Hugot from Aix Marseille University, CNRS, Laboratoire d’Astrophysique de Marseille (LAM) for providing the useful advice.

Disclosures

The authors declare no conflicts of interest.

Data availability

The data underlying the results presented in this paper are not publicly available at this time, but may be obtained from the authors upon reasonable request.

References

1. S. Kendrew, P. Doel, D. Brooks, A. M. King, C. Dorn, C. Yates, R. M. Dwan, I. Richardson, and G. Evans, “Prototype carbon fiber composite deformable mirror,” Opt. Eng. 46(9), 094003 (2007). [CrossRef]  

2. S. Kendrew, P. Doel, D. Brooks, A. M. King, C. Dorn, R. M. Dwan, C. Yates, G. Dando, I. Richardson, and G. Evans, “Carbon fiber composite deformable mirrors: Developments at UCL,” Proc. SPIE 6272, 62722Q (2006). [CrossRef]  

3. H. Zhao, X. W. Fan, Z. H. Pang, G. Ren, W. Wang, Y. Xie, Z. Ma, Y. Du, Y. Su, and J. Wei, “Carbon-fiber-reinforced polymer variable-curvature mirror used for optical zoom imaging: Prototype design and experimental demonstration,” Opt. Eng. 54(2), 025103 (2015). [CrossRef]  

4. C. C. Wilcox, M. E. L. Jungwirth, D. V. Wick, M. S. Baker, C. G. Hobart, R. C. Romeo, and R. N. Martin, “Closed-loop performance of an actuated deformable carbon fiber reinforced polymer mirror,” Proc. SPIE 8373, 83730S (2012). [CrossRef]  

5. B. E. Bagwell, D. V. Wick, W. D. Cowan, O. B. Spahna, W. C. Sweatt, T. Martinez, S. R. Restaino, J. R. Andrews, C. C. Wilcox, D. M. Payne, and Robert Romeo, “Active zoom imaging for operationally responsive space,” Proc. SPIE 6467, 64670D (2007). [CrossRef]  

6. D. V. Wick, T. Martinez, D. M. Payne, W. C. Sweatt, and S. R. Restaino, “Active optical zoom system,” Proc. SPIE 5798(1), 151–157 (2005). [CrossRef]  

7. U. J. Greiner and H. H. Klingenberg, “Thermal lens correction of a diode-pumped Nd: YAG laser of high TEM00 power by an adjustable-curvature mirror,” Opt. Lett. 19(16), 1207–1209 (1994). [CrossRef]  

8. V. V. Apollonov, G. V. Vdovin, L. M. Ostrovskaya, V. N. Rodin, and S. A. Chetkin, “Active correction of a thermal lens in a solid-state laser. I. Metal mirror with a controlled curvature of the central region of the reflecting surface,” Sov. J. Quantum. Electron. 21(1), 116–118 (1991). [CrossRef]  

9. J. Schwarz, M. Ramsey, D. Headley, P. Rambo, I. Smith, and J. Porter, “Thermal lens compensation by convex deformation of a flat mirror with variable annular force,” Appl. Phys. B: Lasers Opt. 82(2), 275–281 (2006). [CrossRef]  

10. Z. Feng, L. Bai, Z. Zhang, and G. Lin, “Thermal deformation compensation of high-energy laser mirrors,” Opt. and Precis. Eng. 18(8), 1781–1787 (2010). [CrossRef]  

11. M. Ferrari and F. Derie, “Variable curvature mirrors for the VLTI,” Proc. SPIE 3350, 830–838 (1998). [CrossRef]  

12. Z. Challita, E. Hugot, F. Madec, M. Ferrari, D. L. Mignant, S. Vives, and J. G. Cuby, “Active optics: Variable curvature mirrors for ELT laser guide star refocusing systems,” Proc. SPIE 8172, 81721A (2011). [CrossRef]  

13. E. Hugot, M. Ferrari, G. R. Lemaitre, F. Madec, S. Vives, E. Chardin, D. L. Mignant, and J. G. Cuby, “Active optics for high-dynamic variable curvature mirrors,” Opt. Lett. 34(19), 3009–3011 (2009). [CrossRef]  

14. E. Bin Nun and F. Dothan Deutsch, “Mirror with adjustable radius of curvature,” Rev. Sci. Instrum. 44(4), 512–513 (1973). [CrossRef]  

15. T. B. Volyak, I. K. Krasjuk, and P. P. Pashinin, “Adaptive mirror for large-aperture laser beams,” Proc. SPIE 0473, 280–283 (1985). [CrossRef]  

16. P. I. Belomestnov, E. I. Vyazovich, R. I. Soloukhin, and Y. A. Yakobi, “Tunable resonator with a mirror of variable curvature,” Sov. J. Quantum. Electron. 3(4), 347–348 (1974). [CrossRef]  

17. U. K. Jarosch, “Adaptive metal mirror for high-power CO2 lasers,” Proc. SPIE 2774, 457–467 (1996). [CrossRef]  

18. T. Warwick, M. Howells, and M. Shlezinger, “A variable radius mirror for imaging the exit slit of an SGM undulatore beamline at the ALS,” Rev. Sci. Instrum. 66(2), 2270–2272 (1995). [CrossRef]  

19. J. Schwarz, M. Geissel, P. Rambo, J. Porter, D. Headley, and M. Ramsey, “Development of a variable focal length concave mirror for on-shot thermal lens correction in rod amplifiers,” Opt. Express 14(23), 10957–10969 (2006). [CrossRef]  

20. J. H. Clark III, F. E. Penado, M. Heilman, and W. E. Clark, “Experimental verification of compliant mirror wavefront correction using a single actuator,” Proc. SPIE 10747, 1074708 (2018). [CrossRef]  

21. T. Okada, K. Ebata, M. Shiozaki, T. Kyotani, A. Tsuboi, M. Sawada, and H. Fukushima, “Development of adaptive mirror for CO2 laser,” Proc. SPIE 3888, 509–520 (2000). [CrossRef]  

22. T. Okada, Sumitomo Electric Industries, Ltd., Osaka. “Variable-curvature reflecting mirror,” United States, 6021153[P]. 2000-02-01.

23. X. Lu, Y. Zhang, and Z. Cheng, “Flexible focus variable expander,” Proc. SPIE 3550, 419–424 (1998). [CrossRef]  

24. X. Xie, H. Zhao, C. Li, X. Fan, G. Ren, W. Chai, L. Xu, and Y. Su, “Annular force based on variable curvature mirror technique,” Acta. Photon. Sinica. 45(8), 822001 (2016). [CrossRef]  

25. H. Zhao, X. Xie, L. Xu, G. Zou, C. Li, and X. Fan, “Performance comparison between two kinds of variable curvature mirrors: Mathematical analysis, prototype design, experimental demonstration, and application potentials in realizing non-moving element optical zooming,” Proc. SPIE 10837, 108371C (2019). [CrossRef]  

26. G. R. Lemaitre, Astronomical optics and elasticity theory, (Springer, 2009).

27. X. Xie, E. Hugot, C. Li, L. Xu, S. Lemared, M. Liu, X. Fan, and H. Zhao, “Elastic bending of variable curvature mirrors: Validation of a simplified analytical method,” Appl. Opt. 58(26), 7121–7126 (2019). [CrossRef]  

28. S. P. Timoshenko and S. W. Krieger, Theory of plates and shells, 2nd ed, (McGraw-Hill, 1959).

29. K. Yeh, “Bending of a thin elastic plate of variable thickness,” Acta. Phys. Sinica. 11(3), 207–218 (1955). [CrossRef]  

Supplementary Material (1)

NameDescription
Visualization 1       The visualization1 shows how the mirror's surface change when it is actuated. The pressure value used just for demonstrating the actuating process and the value is 0.04MPa.

Data availability

The data underlying the results presented in this paper are not publicly available at this time, but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (17)

Fig. 1.
Fig. 1. Physical model of simply supported elastic thin plate with variable thickness
Fig. 2.
Fig. 2. Cross-section stress analysis of physical model of simply supported elastic thin plate
Fig. 3.
Fig. 3. 1/2 cross section of elastic thin plate structure parameter
Fig. 4.
Fig. 4. VCM profile with main parameters
Fig. 5.
Fig. 5. FEA and analytical deflection curves of elastic thin plate at 0.1 MPa
Fig. 6.
Fig. 6. FEA and analytical curves of subtraction of normalized deflection
Fig. 7.
Fig. 7. Comparison of central deflection values given by Eq. (28) and finite element analysis
Fig. 8.
Fig. 8. Comparison of surface deformation along radius direction between values given by Eq. (11) and values obtained through FEA
Fig. 9.
Fig. 9. (a) VCM after ion beam figuring. (b) Surface figure accuracy after IBF.
Fig. 10.
Fig. 10. Air pressure actuation system
Fig. 11.
Fig. 11. Photographs of pressure sensor, integrated processor, and micro air pump
Fig. 12.
Fig. 12. Assembled performance testing system of VCM
Fig. 13.
Fig. 13. Accuracy of VCM surface figure under pressure: (a) 0.02, (b) 0.036, (c) 0.053, and (d) 0.07 MPa (for a demonstration of the actuation effect, see Visualization 1, where the pressure used is 0.04 MPa)
Fig. 14.
Fig. 14. Zernike coefficient variation with pressure
Fig. 15.
Fig. 15. RMS evolution with pressure
Fig. 16.
Fig. 16. Accuracy of VCM surface figure under (a) 0.02, (b) 0.036, (c) 0.053, (d) 0.070 MPa after removing the spherical aberration
Fig. 17.
Fig. 17. Comparison of VCM deformation between FEA and real air pressure actuation (with original curvature radius)

Tables (2)

Tables Icon

Table 1. Comparison of FEA and analytical deformation values of elastic plate along the radius under 0.1 MPa

Tables Icon

Table 2. Experimental results of VCM surface figure accuracy and corresponding deformation

Equations (28)

Equations on this page are rendered with MathJax. Learn more.

M r + d M r d r r M t + Q r = 0 ,
M r = D ( d 2 w d r 2 + μ r d w d r ) ,
M t = D ( 1 r d w d r + μ d 2 w d r 2 ) ,
Q = 1 2 π r 0 r q 2 π r d r ,
D d d r 1 r d d r ( r d w d r ) + d D d r ( d 2 w d r 2 + μ r d w d r ) = 1 r 0 r q r d r ,
D d d r 1 r d d r ( r d w d r ) + d D d r ( d 2 w d r 2 + μ r d w d r ) = 1 2 q r .
x = ( r / a ) 2 , y = w h 0 , Q = 3 a 4 ( 1 μ 2 ) 4 E h 0 4 q .
d 2 d x 2 x d y d x k 2 ( 1 k x ) [ 2 x d 2 y d x 2 + ( 1 + μ ) d y d x ] = Q 1 k x .
( w ) x = 1 = 0 , ( M r ) x = 1 = 0.
{ y = 0 2 d 2 y d x 2 + ( 1 + μ ) d y d x = 0 f o r x = 1.
y = y 0 ( x ) + y 1 ( x ) ε + y 2 ( x ) ε 2 + y 3 ( x ) ε 3 + y 4 ( x ) ε 4 + ,
h = h 0 ( 1 k x ) 1 / 3 = h 0 ( 1 1 3 ε x 1 9 ε 2 x 2 5 81 ε 3 x 3 10 243 ε 4 x 4 ) .
d 2 d x 2 x d y 0 d x = Q
y 0 ( 1 ) = 0 , 2 d 2 y 0 ( 1 ) d x 2 + ( 1 + μ ) d y 0 ( 1 ) d x = 0
y 0 = Q 4 [ x 2 2 ( 3 + μ ) ( 1 + μ ) x + ( 5 + μ ) ( 1 + μ ) ] .
d 2 d x 2 x d y 1 d x = [ x d 2 y 0 d x 2 + 1 2 ( 1 + μ ) d y 0 d x + Q x ] ,
y 1 ( 1 ) = 0 , 2 d 2 y 1 ( 1 ) d x 2 + ( 1 + μ ) d y 1 ( 1 ) d x = 0.
y 1 = Q 144 [ 2 ( 7 + μ ) x 3 9 ( 3 + μ ) x 2 + 12 ( 4 + 3 μ + μ 2 ) ( 1 + μ ) x + ( 61 16 μ 5 μ 2 ) ( 1 + μ ) ] .
d 2 d x 2 x d y 2 d x = [ x d 2 y 1 d x 2 + x 2 d 2 y 0 d x 2 + 1 2 ( 1 + μ ) ( d y 1 d x + d y 0 d x x ) + Q x 2 ] ,
y 2 ( 1 ) = 0 , 2 d 2 y 2 ( 1 ) d x 2 + ( 1 + μ ) d y 2 ( 1 ) d x = 0.
y 2 = Q 2304 [ Q 2304 ( 119 + 24 μ + μ 2 ) x 4 + 8 ( 21 10 μ μ 2 ) x 3 + 24 ( 4 + 3 μ + μ 2 ) x 2 4 ( 59 79 μ + 13 μ 2 + 7 μ 3 ) ( 1 + μ ) x + ( 381 187 μ + 19 μ 2 + 11 μ 3 ) ( 1 + μ ) ]
d 2 d x 2 x d y 3 d x = [ x d 2 y 2 d x 2 + x 2 d 2 y 1 d x 2 + x 3 d 2 y 0 d x 2 + 1 2 ( 1 + μ ) ( d y 2 d x + d y 1 d x x + d y 0 d x x 2 ) + Q x 3 ] ,
y 3 ( 1 ) = 0 , 2 d 2 y 3 ( 1 ) d x 2 + ( 1 + μ ) d y 3 ( 1 ) d x = 0
y 3 = Q 691200 [ Q 691200 6 ( 3689 + 863 μ + 55 u 2 + μ 3 ) x 5 + 75 ( 357 191 μ 27 μ 2 μ 3 ) x 4 + 400 ( 28 + 17 μ + 10 μ 2 + μ 3 ) x 3 150 ( 59 79 μ + 13 μ 2 + 7 μ 3 ) x 2 30 ( 841 1424 μ + 638 μ 2 16 μ 3 39 μ 4 ) ( 1 + μ ) x ( 49921 + 27532 μ 9282 μ 2 + 116 μ 3 + 451 μ 4 ) ( 1 + μ ) ]
d 2 d x 2 x d y 4 d x = [ x d 2 y 3 d x 2 + x 2 d 2 y 2 d x 2 + x 3 d 2 y 1 d x 2 + x 4 d 2 y 0 d x 2 + 1 2 ( 1 + μ ) ( d y 3 d x + d y 2 d x x + d y 1 d x x 2 + d y 3 d x x 3 ) + Q x 4 ]
y 4 ( 1 ) = 0 , 2 d 2 y 4 ( 1 ) d x 2 + ( 1 + μ ) d y 4 ( 1 ) d x = 0
y 4 = Q 8294400 [ ( 180761 + 45976 μ + 3558 u 2 + 104 μ 3 + μ 4 ) x 6 18 ( 11067 + 6278 μ + 1028 μ 2 + 58 μ 3 + μ 4 ) x 5 + 150 ( 476 + 261 μ + 187 μ 2 + 27 μ 3 + μ 4 ) x 4 + 100 ( 413 + 494 μ 12 μ 2 62 μ 3 7 μ 4 ) x 3 45 ( 841 1424 μ + 638 μ 2 16 μ 3 39 μ 4 ) x 2 18 ( 6677 12541 μ + 7838 μ 2 1966 μ 3 115 μ 4 + 107 μ 5 ) ( 1 + μ ) x 6 ( 48196 + 23725 μ 12048 μ 2 + 2702 μ 3 + 148 μ 4 123 μ 5 ) ( 1 + μ ) ]
y ( 0 ) = Q [ 1 4 ( 5 + μ ) ( 1 + μ ) + 1 144 ( 61 16 μ 5 μ 2 ) ( 1 + μ ) ε + 1 2304 ( 381 187 μ + 19 μ 2 + 11 μ 3 ) ( 1 + μ ) ε 2 1 691200 ( 49921 + 27532 μ 9282 μ 2 + 116 μ 3 + 451 μ 4 ) ( 1 + μ ) ε 3 1 8294400 6 ( 48196 + 23725 μ 12048 μ 2 + 2702 μ 3 + 148 μ 4 123 μ 5 ) ( 1 + μ ) ε 4 ]
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.