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

Research on a bimorph piezoelectric deformable mirror for adaptive optics in optical telescope

Open Access Open Access

Abstract

We have proposed a discrete-layout bimorph piezoelectric deformable mirror (DBPDM) and developed its realistic electromechanical model. Compared with the conventional piezoelectric deformable mirror (CPDM) and the bimorph piezoelectric deformable mirror (BPDM), the DBPDM has both a larger stroke and a higher resonance frequency by integrating the strengths of the CPDM and the BPDM. To verify the advancement, a 21-elements DBPDM is studied in this paper. The results have suggested that the stroke of the DBPDM is larger than 10 microns and its resonance frequency is 53.3 kHz. Furthermore, numerical simulation is conducted on the deformation of the mirror using the realistic electromechanical model, and the dependence of the influence function upon the size of the radius of push pad is analyzed.

© 2017 Optical Society of America

1. Introduction

Stars that can be observed by astronomers using a telescope are usually thousands to millions of light-years away. After star light travels such a long distance to our eyes, a ground-based astronomy telescope would cause serious wavefront errors, mainly because of the distorts due to atmosphere turbulence [1]. Adaptive optics (AO) is a means for real time compensation of the image degradation. The technique was first proposed by Babcock [2]. It consists of using a deformable mirror to correct the instantaneous wave-front distortions [3, 4]. The deformable mirror is usually supported by an array of active elements, which adjusts the mirror surface shape to compensate the light’s wavefront error caused by the atmospheric distortions [5–7]. The band-pass of atmospheric turbulence is more than103Hz [8]. The effective bandwidth of this precision adaptive structure is typically limited by the controller’s signal processing speed and the resonance frequency of the deformable mirror assembly [9]. Recently, the aperture becomes larger, the cost of adaptive optics telescopes becomes much more expensive [10]. Thus, the system requires its actuators to have large stroke, high speed, high accuracy, high resonance frequency, and low cost. There are two kinds of deformable mirrors have been widely used in adaptive optics [11]. One kind of those deformable mirrors, known as conventional piezoelectric deformable mirrors (CPDM), is usually supported by an array of cylindrical piezoelectric actuators, coupled to a continuous mirrored face sheet [12–15]. Another kind, called bimorph piezoelectric deformable mirror (BPDM), is made of a glass sheet bonded to a piezoceramics layer with an appropriate electrode pattern on the back [16–20]. However, each of them has its advantages: the CPDM has a higher temporal response [21], while the BPDM has lower cost and larger stroke. A discrete-layout bimorph piezoelectric deformable mirror, which integrates the strengths of the CPDM and the BPDM, is developed to make the best of the advantages and minimize the disadvantages.

Several kinds of bending actuators had been analyzed in our previous works, the results have shown that the bending actuators, of a simple structure and at a low cost, have the large strokes and the high resonance frequencies [22–24]. As an application, here we propose a discrete-layout bimorph piezoelectric deformable mirror (DBPDM), which consists of a continuous mirrored face sheet and an array of one kind of these bending actuators—the bimorph beam piezoelectric bending actuator (see Fig. 1). The mirror face sheet connects to the mid-points of beams by push pads. The cross section of the push pad is a square with the side length being wpad[see Fig. 2(a)and Fig. 2(c)]. The cross-section area of a push pad is Apad = (wpad)2. As shown in Fig. 2(a), the structure of a simply supported bimorph beam piezoelectric bending actuator (BBPBA) consists of identical piezoelectric ceramic twin layers polarized in the thickness direction and a metallic layer in the middle. The metallic layer strengthens the bending stiffness of the BBPBA, and this configuration is useful to achieve a high resonance frequency of whole deformable mirror.

 figure: Fig. 1

Fig. 1 The configuration of a discrete-layout bimorph piezoelectric deformable mirror.

Download Full Size | PDF

 figure: Fig. 2

Fig. 2 (a) A simply-support BBPBA. (b) The cross section of the BBPBA. (c) The dome height of BBPBA and loading a concentrated force. (d)The cross section of the push pad.

Download Full Size | PDF

The geometric parameters of the actuator are shown in Fig. 2(a) and 2(b). The total thickness of the BBPBA is 2h, and the thickness of the metallic layer is 2c. The length and the width of the BBPBA are 2L (2L>>wpad) and b, respectively. The top and the bottom electrodes are in parallel and have the same voltage across each one. The metallic middle layer is shared as the third electrode by the upper and lower piezoelectric layers. The BBPBA is simply supported at the boundary at x1 = ± L, where any in-plane displacement is prohibited. When a load voltage V is applied at the BBPBA, one piezoelectric layer will be driven to extend and the other layer to contract, which induces the BBPBA into a state having a bending deformation. As shown in Fig. 2(d), the dome height d is the distance from the horizontal line to the mid-point of the beam plate during the deformation process, and it is also the deflection at x = 0. Distinct from the conventional cylindrical piezoelectric actuators (piezoceramics stacks), whose strokes are obtained by stretching out their length in the thickness direction, the BBPBA has a large stroke due to its bending deformation.

2. Realistic electro-mechanical model of discrete-layout bimorph piezoelectric deformable mirror (DBPDM)

Here, the piezoelectric model of the actuator in this paper is affected by the stress on the cross section of the push pad [see Fig. 2(a)and Fig. 2(d)].

F=ATzdS=T¯zApad=wpad2T¯z.
where Tz and F denote the stress on the cross section and the effectively concentrated force on the cross section, respectively. When the deformable mirror is working, each actuator is acted upon by external forces from other active actuators, and F is the resultant of these external forces.

In the following, we consider the flexural deformation of the BBPBA in the x3 direction by assuming that the BBPBA is slender, and 2L is much larger than the thickness 2h and the width b. We denote the deflection curve by u3(x,t). According to the electrode configurations, the electric field in the BBPBA has the following components: E1 = E2 = 0, E3 = −V/(h-c) at c<z<h, and E3 = V/(h-c) at -h<z<-c. Using the usual one-dimensional stress approximation of beam [25, 26], we have the following stress components: σ1 = σ1(x,t), σ2 = σ3 = σ4 = σ5 = σ6 = 0. From [25] and [27], we can obtain the axial stresses in the upper and lower layers

σ1p=x3s111(2u3/x2)s111d31E3+Fx3(x1+L)/(2Ie),
where Ie = 2Ip + b(h-c)(h + c)2/2 + nIm is the effective moment of inertia of three-layer composite cross-section, and n = Ems11in which Em is Young’s modulus of the metallic middle layer [27]. Here, the middle layer is assumed to be elastic, and its constitutive relation is stated as below:
σ1m=x3Em(2u3/x2)+nFx3(x1+L)/(2Ie).
The bending moment is defined by the following integral over the cross section, which can be integrated explicitly with the expression forσ1 in Eq. (2) and Eq. (3) [25,28]
M=x3σ1dx2dx3=D(2u3/x2)+(h+c)e31pV+ΓF(x1+L).
Where e31p=bd31s111, D = (2/3)[Emc3 + (h3-c3)/s11]b, and Γ = (nc3 + h3-c3)/(3Ie). D is the bending stiffness. The equation of equilibrium in a slender beam can be expressed by the following equation regarding the bending moment:
2M/x12=D(4u3/x14)=0.
Because the BBPBA is simply supported at x1 = ± L, we have
u3(L)=u3(L)=0,M(L)=M(L)=0.
A general solution to Eq. (5) can be written as
u3=H1+H2x1+H3x12+H4x13.
whereH1-H4 are constants to be determined by Eq. (6). Thus, we get the static displacement of the dome height
ds=u3(0)=H1=[2e31pV(c+h)+LΓF]L2/(4D).
where ds denotes the dome height at x1 = 0, which is equal to the extensional displacement of the actuator. From Eq. (8), we can obtain
F=ksds+kvV,
where ks = -2D/(L3Γ) and kv = ep31(c + h)/(LΓ).

The BBPBAs at different sites drive the deformable mirror system. We treat the mirror as linear elastic structure impeding the motion of the BBPBA array. The number of BBPBAs in the mirror is Nact. The force-displacement relationship for the mirror structure in the BBPBA array is [15]

{FN}=[KdmMN]{dsNd0}.
Kdm represents the Nact by Nact stiffness matrix of the mirror, while F and d are vectors with the size being Nact. The displacement dN, is the displacement of the Nth BBPBA normal to the plane of the mirror face sheet, and FN is the corresponding force at the BBPBA location. The superscripts denote the elements of the matrix and the vectors, and the repeated superscripts indicate the summation. d0 is the average displacement of all the BBPBAs in the array. Since the mirror rests only on the BBPBA in the array, the reaction force caused by BBPBA displacements are non-zero relative to d0. Compression occurred in BBPBAs with displacements above the average, while tension occurred in BBPBAs with displacements below the average. When all the BBPBAs had the same voltage and hence the same displacements, the reaction forces on the BBPBAs are zero.

The Kdm can be obtained via an FEA model, which only includes the mirror face sheet and push pads. This method has been discussed in detail in [15]. Using Eqs. (9) and (10), we obtain [29]

[KdmMN]{dNd0}=ks{dN}+kv{VN}.
The average displacement for the BBPBA array is d0 = (1/NactdN. Thus, there is a linear system of Nact + 1 equations for the displacements. This system is solved numerically for the given BBPBA voltages. Finally, the normal displacement w(x,y) in the mirror face sheet is reconstructed from the influence functions fN(x,y)using

w(x,y)=NfN(x,y)dN.

3. Numerical results and discussion

In the following numerical calculations, we take a deformable mirror of 21 elements as an example (see Fig. 3 and Fig. 4). The material of pads and face sheet is single crystal silicon with Young’s modulus Esi = 190 GPa, Poisson’s ratioν = 0.30, mass density is ρ = 2350 kg/m3, and the radius and final thickness of the face sheet are R = 15 mm and h0 = 1 mm, respectively. We use PZT-5H [30] as our piezoelectric material of actuator, whose parameters are listed below:

 figure: Fig. 3

Fig. 3 (a)The static dome height displacement ds versus the electric field for different stresses of Apad. (b) The . resonance frequency of the whole DBPDM is obtained with the mode analysis of whole finite element model.

Download Full Size | PDF

 figure: Fig. 4

Fig. 4 A deformable mirror of 21 elements. (a)The locations of No.1, No.2 and No.3 BBPBAs. (b) FEA model of mirror face sheet and pads.

Download Full Size | PDF

(s11,s12,s13,s44) = (16.5,-4.78,-8.45,20.7,43.5) × 10−12m2/N,(d31,d15,d33) = (−274,741,539) × 10−12C/N, ρ = 7500 kg/m3. We consider a beam flexural-mode piezoelectric actuator with 2L = 5 mm, Δh = h-c = 0.5mm, η = ch = 1/2, b = 2mm, wpad = 1mm, unless stated otherwise. Besides, the metallic layer is an aluminum (Al) alloy with Young’s modulus Em = 70 GPa, Poisson’s ratioν = 0.33, and mass density is ρ = 2700 kg/m3.

The model predictions of ds in response to the electric field are shown in Fig. 3(a) for different stresses. From the figure, ds has a zero offset due to the external force. When T¯z<0, the effectively concentrated force is a compression force, and the zero offset is a negative value. Conversely, the effectively concentrated force is a stretching force, and its zero offset is a positive value. From Eq. (8), it is observed that the zero offset is equal to L4FΓ/(4D) and proportional to F. From Fig. 3(a), the stroke of the actuator is larger than 10 microns at about |E3| = 10MV/m, and it equates to the stroke of the deformable mirror [see Eqs. (10) and (11)]. According to the performance requirements of a deformable mirror with high stroke, its stroke must be at least 10 microns. Thus, the deformable mirror has a large stroke. Here, the finite element model of DBPDM is created with the commercial finite element software ANSYS. With the mode analysis of whole finite element model, we get its resonant frequency of 53.3kHz which is much larger than that of the BPDM. Without the metallic layer, the resonant frequency is 14.3 kHz, which is still very high. The high resonant frequency is caused by the introduction of a metallic layer into the BBPBAs.

With the realistic electro-mechanical model of DBPDM, consisting of Eqs. (1)-(12) and the FEA model of mirror face sheet and pads [see Fig. 4(b)], we can obtain the influence functions of all elements of the deformable mirror. For example, a vector of voltage {VN}will be formed when only the No.1 BBPBA [1 in Fig. 4 (a)] is active with the applied voltage being 1 V, and 0 V are applied to other BBPBAs, respectively. With Eq. (11), we can obtain the displacement vector {dN}, which is used as the input displacement loads to be applied to pads of the FEA model. The resulting displacement field in mirror face sheet defines the influence function fN(x,y). Thus, we can get the influence function for the face sheet surface surrounding No.1 BBPBA [see Fig. 5(a)]. With the same way, we can get the influence functions of the other BBPBAs. Figure 5(b) shows the curves of the influence functions of No.1 BBPBA versus normalized radius for the different side lengths of the cross-section of the pad. It is observed that the top of the peak of the influence function widens with the increasing wpad, and the peak foot of the influence function narrows as wpad increases. Furthermore, the negative values of the influence function decrease with the increasing wpad. These phenomena also appear in the locations of No. 2 and No.3 BBPBAs with an increase of wpad (see Fig. 6). When all sites' influence functions have been obtained, we can use the realistic electromechanical model to simulating wavefront correction, which will be addressed in future.

 figure: Fig. 5

Fig. 5 (a) The 3D influence function for a unit voltage of the No.1 BBPBA. (b)The influence function versus normalized radius for different side lengths of pads at No.1 BBPBA.

Download Full Size | PDF

 figure: Fig. 6

Fig. 6 (a) The influence function versus normalized radius for different side lengths of pads at No.2 BBPBA. (b) The influence function versus normalized radius for different side lengths of pads at No.3 BBPBA.

Download Full Size | PDF

4. Summary

In this paper, a discrete-layout bimorph piezoelectric deformable mirror (DBPDM) and its realistic electromechanical model have been proposed. As an example, a 21-elements DBPDM is studied. The achievement of a high resonance frequency of the whole DBPDM have been proved by the results of the mode analysis of whole finite element model. We have obtained the large stroke of the DBPDM through the analytical analysis of the behaviors of the actuator.

The numerical simulation is performed on the DBPDM using with the realistic electromechanical model. The results show that the size of the radius of the push pad has an impact on the influence function of the deformable mirror. In future, we can use the realistic electromechanical model to simulate wavefront correction errors in the study of the adaptive optics in optical telescope. Besides, the cost of the deformable mirror could be reduced due to the simple structure of the BBPBA. The results in this paper are useful for the design of the deformable mirror.

Funding

National Natural Science Foundation of China (Grant Nos. 10921063, 11190014, 11403109 and 11373073); Natural Science Foundation of Jiangsu Province (Grants No. BK20141042).

References and links

1. P. Y. Lin, H. T. Hsieh, and G. D. J. Su, “Design and fabrication of a large-stroke MEMS deformable mirror for wavefront control,” J. Opt. 13(5), 055404 (2011). [CrossRef]  

2. H. W. Babcock, “The possibility of compensating astronomical seeing,” Publ. Astron. Soc. Pac. 65, 229–236 (1953). [CrossRef]  

3. H. W. Babcock, “The possibility of compensating astronomical seeing,” Publ. Astron. Soc. Pac. 65, 229–236 (1953). [CrossRef]  

4. N. Woolf, “Adaptive optics,” in IAU Colloq. 79: Very Large Telescopes, their Instrumentation and Programs (1984), 221–233.

5. R. Tyson, Principles of Adaptive Optics (CRC Press, 2010).

6. S. Bonora and L. Poletto, “Push-pull membrane mirrors for adaptive optics,” Opt. Express 14(25), 11935–11944 (2006). [CrossRef]   [PubMed]  

7. A. Haber, A. Polo, S. Ravensbergen, H. P. Urbach, and M. Verhaegen, “Identification of a dynamical model of a thermally actuated deformable mirror,” Opt. Lett. 38(16), 3061–3064 (2013). [CrossRef]   [PubMed]  

8. R. Wilson, F. Franza, and L. Noethe, “Active optics: I. a system for optimizing the optical quality and reducing the costs of large telescopes,” J. Mod. Opt. 34(4), 485–509 (1987). [CrossRef]  

9. C. L. Hom, “Simulating electrostrictive deformable mirrors: II. Nonlinear dynamic analysis,” Smart Mater. Struct. 8(5), 700–708 (1999). [CrossRef]  

10. X. D. Lin, C. Xue, X. Y. Liu, J. l. Wang, and P. F. Wei, “Current status and research development of wavefront correctors for adaptive optics,” Chinese Optics 4, 337–351 (2012).

11. P. Y. Madec, “Overview of deformable mirror technologies for adaptive optics and astronomy,” in SPIE Astronomical Telescopes + Instrumentation, 844705 (International Society for Optics and Photonics, 2012), pp. 5–18.

12. H. Hofer, L. Chen, G.-Y. Yoon, B. Singer, Y. Yamauchi, and D. R. Williams, “Improvement in retinal image quality with dynamic correction of the eye’s aberrations,” Opt. Express 8(11), 631–643 (2001). [CrossRef]   [PubMed]  

13. A. Roorda, F. Romero-Borja, W. Donnelly Iii, H. Queener, T. Hebert, and M. Campbell, “Adaptive optics scanning laser ophthalmoscopy,” Opt. Express 10(9), 405–412 (2002). [CrossRef]   [PubMed]  

14. Y. Zhang, J. Rha, R. Jonnal, and D. Miller, “Adaptive optics parallel spectral domain optical coherence tomography for imaging the living retina,” Opt. Express 13(12), 4792–4811 (2005). [CrossRef]   [PubMed]  

15. C. L. Hom, P. D. Dean, and S. R. Winzer, “Simulating electrostrictive deformable mirrors: II. Nonlinear dynamic analysis,” Smart Mater. Struct. 8(5), 691–699 (1999). [CrossRef]  

16. Y. Ning, W. Jiang, N. Ling, and C. Rao, “Response function calculation and sensitivity comparison analysis of various bimorph deformable mirrors,” Opt. Express 15(19), 12030–12038 (2007). [CrossRef]   [PubMed]  

17. R. J. Zawadzki, S. M. Jones, S. S. Olivier, M. Zhao, B. A. Bower, J. A. Izatt, S. Choi, S. Laut, and J. S. Werner, “Adaptive-optics optical coherence tomography for high-resolution and high-speed 3D retinal in vivo imaging,” Opt. Express 13(21), 8532–8546 (2005). [CrossRef]   [PubMed]  

18. J. Zhao, F. Xiao, J. Kang, H. Zhao, Y. Dai, and Y. Zhang, “Quantifying intraocular scatter with near diffraction-limited double-pass point spread function,” Biomed. Opt. Express 7(11), 4595–4604 (2016). [CrossRef]   [PubMed]  

19. J. Ma, Y. Liu, Y. Hu, C. Xu, B. Li, and J. Chu, “Low-cost unimorph deformable mirror with high actuator count for astronomical adaptive optics,” Opt. Eng. 52(1), 016602 (2013). [CrossRef]  

20. P. Wnuk, C. Radzewicz, and J. Krasiński, “Bimorph piezo deformable mirror for femtosecond pulse shaping,” Opt. Express 13(11), 4154–4159 (2005). [CrossRef]   [PubMed]  

21. E. J. Fernández, L. Vabre, B. Hermann, A. Unterhuber, B. Povazay, and W. Drexler, “Adaptive optics with a magnetic deformable mirror: applications in the human eye,” Opt. Express 14(20), 8900–8917 (2006). [CrossRef]   [PubMed]  

22. H. Wang and S. Yang, “Modeling and analysis of the thermal effects of a circular bimorph piezoelectric actuator,” Appl. Opt. 55(4), 873–878 (2016). [CrossRef]   [PubMed]  

23. H. Wang, M. Hu, and Z. Li, “Modelling and analysis of circular bimorph piezoelectric actuator for deformable mirror,” Appl. Math. Mech. 37(5), 639–646 (2016). [CrossRef]  

24. H. Wang, “Analytical analysis of a beam Flexural-Mode Piezoelectric Actuator for Deformable Mirrors,” J. Astron. Telesc. Instrum. Syst. 14, 049001 (2015).

25. S. Jiang, X. Li, S. Guo, Y. Hu, J. Yang, and Q. Jiang, “Performance of a piezoelectric bimorph for scavenging vibration energy,” Smart Mater. Struct. 14(4), 769–774 (2005). [CrossRef]  

26. J. S. Yang and H. Y. Fang, “Analysis of a rotating elastic beam with piezoelectric films as an angular rate sensor,” IEEE Trans. Ultrason. Ferroelectr. Freq. Control 49(6), 798–804 (2002). [CrossRef]   [PubMed]  

27. T. Low and W. Guo, “Modeling of a three-layer piezoelectric bimorph beam with hysteresis,” J. Microelectromech. Syst. 4(4), 230–237 (1995). [CrossRef]  

28. H. Wang, X. Xie, Y. Hu, and J. Wang, “Nonlinear analysis of a 5-layer beam-like piezoelectric transformer near resonance,” Guti Lixue Xuebao 27, 195–201 (2014).

29. A. E. H. Love, A Treatise on the Mathematical Theory of Elasticity (Cambridge University Press, 2013), Vol. 1.

30. B. A. Auld, Acoustic Fields and Waves in Solids (Рипол Классик, 1973).

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 (6)

Fig. 1
Fig. 1 The configuration of a discrete-layout bimorph piezoelectric deformable mirror.
Fig. 2
Fig. 2 (a) A simply-support BBPBA. (b) The cross section of the BBPBA. (c) The dome height of BBPBA and loading a concentrated force. (d)The cross section of the push pad.
Fig. 3
Fig. 3 (a)The static dome height displacement ds versus the electric field for different stresses of Apad. (b) The . resonance frequency of the whole DBPDM is obtained with the mode analysis of whole finite element model.
Fig. 4
Fig. 4 A deformable mirror of 21 elements. (a)The locations of No.1, No.2 and No.3 BBPBAs. (b) FEA model of mirror face sheet and pads.
Fig. 5
Fig. 5 (a) The 3D influence function for a unit voltage of the No.1 BBPBA. (b)The influence function versus normalized radius for different side lengths of pads at No.1 BBPBA.
Fig. 6
Fig. 6 (a) The influence function versus normalized radius for different side lengths of pads at No.2 BBPBA. (b) The influence function versus normalized radius for different side lengths of pads at No.3 BBPBA.

Equations (12)

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

F= A T z dS= T ¯ z A pad = w pad 2 T ¯ z .
σ 1 p = x 3 s 11 1 ( 2 u 3 / x 2 ) s 11 1 d 31 E 3 +F x 3 ( x 1 +L)/(2 I e ),
σ 1 m = x 3 E m ( 2 u 3 / x 2 )+nF x 3 ( x 1 +L)/(2 I e ).
M= x 3 σ 1 d x 2 d x 3 =D ( 2 u 3 / x 2 )+(h+c) e 31 p V+ΓF( x 1 +L).
2 M/ x 1 2 =D( 4 u 3 / x 1 4 )=0.
u 3 (L)= u 3 (L)=0, M(L)=M(L)=0.
u 3 = H 1 + H 2 x 1 + H 3 x 1 2 + H 4 x 1 3 .
d s = u 3 (0)= H 1 =[2 e 31 p V(c+h)+LΓF] L 2 /(4D).
F= k s d s + k v V,
{ F N }=[ K dm M N ]{ d s N d 0 }.
[ K dm M N ]{ d N d 0 }= k s { d N }+kv{ V N }.
w(x,y)= N f N (x,y) d N .
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.