Method and device for measuring mechanical parameters of multilayer composite thin film structure

Abstract

A method for measuring mechanical parameters of a multilayer composite thin film structure and belongs to the technical field of online tests of micro-electro-mechanical system (MEMS for short) material parameters. Equivalent Young modulus and equivalent residual stress of each layer of the multilayer composite thin film structure can be obtained in one step by means of solving an equation set on the basis of a relationship between first-order resonance frequency of multilayer composite fixed-fixed beams and multilayer composite cantilever beams and parameters such as material characteristics and structure size, the online test of multilayer thin film materials can be realized, the test structure and calculating method are simple, and the accuracy is higher. The present invention further discloses a device for measuring mechanical parameters of the multilayer composite thin film structure.

Claims

1. A method for measuring mechanical parameters of a multilayer composite thin film structure, the multilayer composite thin film structure being an initially flat and straight N-layer composite thin film structure, where N is an integer greater than or equal to 2, wherein the method comprises the following steps: step 1: preparing a group of at least N initially flat and straight N-layer composite cantilever beam test structures and a group of at least N initially flat and straight N-layer composite fixed-fixed beam test structures by adopting a process for preparing the multilayer composite thin film structure, wherein the material parameters except the geometric size of each test structure are the same; and for each group of test structures, if a vector consisting of the width and length of each layer of each test structure is used as the dimension vector of the test structure, the dimension vector group of the group of test structures is linearly independent; step 2: measuring the first-order resonance frequency of each test structure; step 3: respectively substituting the first-order resonance frequency of the N-layer composite cantilever beam test structures into a first-order resonance frequency expression of N-layer composite cantilever beams to obtain an equation set consisting of at least N linear equations, and solving the equation set to obtain the equivalent Young modulus of each layer of the multilayer composite thin film structure, wherein the expression of the first-order resonance frequency f.sub.1cf of the N-layer composite cantilever beams is specifically as follow: $f_{1\mathrm{cf}}=\frac{{1.875}^2}{2\pi l^2}\sqrt{\frac{\underset{i=1}{\overset{N}{.Math.}}{\tilde{E}}_i{w}_i.Math.\left[{\left(z_i-z_c\right)}^3-{\left(z_{i-1}-z_c\right)}^3\right]}{3\underset{i=1}{\overset{N}{.Math.}}{\rho }_i{w}_i{h}_i}}$ where l is the length of the N-layer composite cantilever beam test structures; the width of the i.sup.th layer thin film material is w.sub.i, the thickness is h.sub.i, the density of the material is ρ.sub.i, and the equivalent Young modulus is {tilde over (E)}.sub.i; z.sub.i is the position of the top surface of the i.sup.th layer thin film material on a z axis, and z.sub.c is the height of the neutral surface; and step 4: respectively substituting the first-order resonance frequency of the N-layer composite fixed-fixed beam test structures and the equivalent Young modulus of each layer of the multilayer composite thin film structure obtained in step 3 into the first-order resonance frequency expression of the N-layer composite fixed-fixed beams to obtain an equation set consisting of at least N equations, and solving the equation set to obtain the equivalent residual stress of each layer of the multilayer composite thin film structure, wherein the expression of the first-order resonance frequency f.sub.1s of the N-layer composite fixed-fixed beams is specifically as follow: $f_{1s}=\frac{{4.730}^2}{2\pi l^2}\sqrt{\frac{\underset{i=1}{\overset{N}{.Math.}}{\tilde{E}}_i{w}_i\left[{\left(\underset{i=0}{\overset{i}{.Math.}}{h}_i-z_c\right)}^3-{\left(\underset{i=0}{\overset{i-1}{.Math.}}{h}_i-z_c\right)}^3\right]}{3\underset{i=1}{\overset{N}{.Math.}}{\rho }_i{w}_i{h}_i}+\frac{l^2}{4{\pi }^2}\frac{\underset{i=1}{\overset{N}{.Math.}}{\tilde{\sigma }}_i{w}_i{h}_i}{\underset{i=1}{\overset{N}{.Math.}}{\rho }_i{w}_i{h}_i}}$ where l is the length of the N-layer composite fixed-fixed beam test structures; the width of the i.sup.th layer thin film material is w.sub.i, the thickness is h.sub.i, the density of the material is ρ.sub.i, the equivalent Young modulus is {tilde over (E)}.sub.i, and the equivalent residual stress is {tilde over (σ)}.sub.i; and z.sub.c is the height of the neutral surface.

2. A device for measuring mechanical parameters of a multilayer composite thin film structure, the multilayer composite thin film structure being an initially flat and straight N-layer composite thin film structure, where N is an integer greater than or equal to 2, wherein the device comprises: a group of at least N initially flat and straight N-layer composite cantilever beam test structures prepared by adopting a process for preparing the multilayer composite thin film structure, wherein the material parameters except the geometric size of each test structure are the same; and if a vector consisting of the width and length of each layer of each test structure is used as the dimension vector of the test structure, the dimension vector group of the group of test structures is linearly independent; a group of at least N initially flat and straight N-layer composite fixed-fixed beam test structures prepared by adopting the process for preparing the multilayer composite thin film structure, wherein the material parameters except the geometric size of each test structure are the same; and if a vector consisting of the width and length of each layer of each test structure is used as the dimension vector of the test structure, the dimension vector group of the group of test structures is linearly independent; a resonance frequency measuring unit used for measuring the first-order resonance frequency of each test structure; a calculation unit used for respectively substituting the first-order resonance frequency of the N-layer composite cantilever beam test structures into a first-order resonance frequency expression of N-layer composite cantilever beams to obtain an equation set consisting of at least N linear equations, and solving the equation set to obtain the equivalent Young modulus of each layer of the multilayer composite thin film structure; and then substituting the first-order resonance frequency of the N-layer composite fixed-fixed beam test structures and the equivalent Young modulus of each layer of the multilayer composite thin film structure into the first-order resonance frequency expression of the N-layer composite fixed-fixed beams to obtain an equation set consisting of at least N equations, and solving the equation set to obtain the equivalent residual stress of each layer of the multilayer composite thin film structure, wherein, the expression of the first-order resonance frequency f.sub.1cf of the N-layer composite cantilever beams is specifically as follow: $f_{1\mathrm{cf}}=\frac{{1.875}^2}{2\pi l^2}\sqrt{\frac{\underset{i=1}{\overset{N}{.Math.}}{\tilde{E}}_i{w}_i.Math.\left[{\left(z_i-z_c\right)}^3-{\left(z_{i-1}-z_c\right)}^3\right]}{3\underset{i=1}{\overset{N}{.Math.}}{\rho }_i{w}_i{h}_i}}$ where l is the length of the N-layer composite cantilever beam test structures; the width of the i.sup.th layer thin film material is w.sub.i, the thickness is h.sub.i, the density of the material is ρ.sub.i, and the equivalent Young modulus is {tilde over (E)}.sub.i; z.sub.i is the position of the top surface of the i.sup.th layer thin film material on a z axis, and z.sub.c is the height of the neutral surface; and the expression of the first-order resonance frequency f.sub.1s of the N-layer composite fixed-fixed beams is specifically as follow: $f_{1s}=\frac{{4.730}^2}{2\pi l^2}\sqrt{\frac{\underset{i=1}{\overset{N}{.Math.}}{\tilde{E}}_i{w}_i\left[{\left(\underset{i=0}{\overset{i}{.Math.}}{h}_i-z_c\right)}^3-{\left(\underset{i=0}{\overset{i-1}{.Math.}}{h}_i-z_c\right)}^3\right]}{3\underset{i=1}{\overset{N}{.Math.}}{\rho }_i{w}_i{h}_i}+\frac{l^2}{4{\pi }^2}\frac{\underset{i=1}{\overset{N}{.Math.}}{\tilde{\sigma }}_i{w}_i{h}_i}{\underset{i=1}{\overset{N}{.Math.}}{\rho }_i{w}_i{h}_i}}$ where l is the length of the N-layer composite fixed-fixed beam test structures; the width of the i.sup.th layer thin film material is w.sub.i, the thickness is h.sub.i, the density of the material is ρ.sub.i, the equivalent Young modulus is {tilde over (E)}.sub.i, and the equivalent residual stress is {tilde over (σ)}.sub.i; and z.sub.c is the height of the neutral surface.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1a is a top view of a two-layer composite cantilever beam test structure.

(2) FIG. 1b is a front view of a two-layer composite cantilever beam test structure.

(3) FIG. 2a is a top view of a two-layer composite fixed-fixed beam test structure.

(4) FIG. 2b is a front view of a two-layer composite fixed-fixed beam test structure.

(5) Reference signs in the drawings respectively have the following meanings:

(6) 101: anchor area, 102: top layer, 103: bottom layer, and 104: substrate.

DETAILED DESCRIPTION OF THE INVENTION

(7) Aiming at the problem about online test of the mechanical parameters of multilayer composite thin film structures, in the present invention, by using the combined test structure of multilayer cantilever beams and multilayer fixed-fixed beams, based on the relationship between the first-order resonance frequency of the multilayer fixed-fixed beams and the multilayer cantilever beams and the parameters such as material characteristics and structure size, the equivalent Young modulus and the equivalent residual stress of each layer of the multilayer composite thin film structure can be obtained in one step, the calculation process is simple and convenient, and the test accuracy is high.

(8) In order to facilitate the public to understand the technical solution of the present invention, the theoretical principle of the technical solution of the present invention is described in detail at first.

(9) For a composite thin film structure formed by at least two layers of thin film structures through combination, when the width and thickness of the thin film satisfy w.sub.i<5h.sub.i, that is, when the thin film structure is a narrow beam, the equivalent Young modulus {tilde over (E)}.sub.i is the Young modulus itself; and when the width and thickness of the thin film satisfy w.sub.i≥5h.sub.i, that is, when the thin film is a wide beam, the equivalent Young modulus is a relational expression between Young modulus and Poisson ratio, that is

(10) ${\tilde{E}}_i=\frac{E_i}{1-v_i^2}.$
The relationship between the equivalent Young modulus {tilde over (E)}.sub.i of the i.sup.th layer thin film and the thickness and width of the thin film is as follow:

(11) $\begin{array}{cc}{\tilde{E}}_i=\{\begin{array}{ccc}{E}_i& ,& w_i<5h_i\\ \frac{E_i}{1-v_i^2}& ,& w_i\ge 5h_i\end{array}& \left(1\right)\end{array}$

(12) It is assumed that the length of the beam is l before the release process (i.e., when no deformation occurs); the width of the i.sup.th layer thin film material is w.sub.i, the thickness is h.sub.i, the density of the material is ρ.sub.i, and the Young modulus is E.sub.i. If the position of the bottom surface of the multilayer thin film on z axis is z.sub.0, z.sub.0=0, and the position of the top surface of the i.sup.th layer from bottom to top is z.sub.i, then:

(13) $\begin{array}{cc}{z}_i=\underset{i=1}{\overset{i}{.Math.}}{h}_i& \left(2\right)\end{array}$

(14) the height of the neutral surface z.sub.c is:

(15) $\begin{array}{cc}{z}_c=\frac{\underset{i=1}{\overset{N}{.Math.}}{\tilde{E}}_i{w}_i\left(z_i^2-z_{i-1}^2\right)}{2\underset{i=1}{\overset{N}{.Math.}}{\tilde{E}}_i{w}_i\left(z_i-z_{i-1}\right)}& \left(3\right)\end{array}$

(16) the moment of inertia I.sub.i of the cross-sectional area of the i.sup.th layer thin film material with respect to the neutral axis is:
I.sub.i=⅓w.sub.i[(z.sub.i−z.sub.c).sup.3−(z.sub.i-1−z.sub.c).sup.3]  (4)

(17) the bending rigidity, linear density and axial load of the micromechanical thin film beam are defined as EI, ρA and σA respectively, i.e.:

(18) $\begin{array}{cc}\stackrel{\_}{\mathrm{EI}}=\underset{i=1}{\overset{N}{.Math.}}{\tilde{E}}_i{I}_i& \left(5\right)\\ \stackrel{\_}{\rho A}=\underset{i=1}{\overset{N}{.Math.}}{\rho }_i{A}_i& \left(6\right)\\ \stackrel{\_}{\sigma A}=\underset{i=1}{\overset{N}{.Math.}}{\tilde{\sigma }}_i{A}_i& \left(7\right)\end{array}$

(19) where A.sub.i is the cross-sectional area of the i.sup.th layer thin film material, ρ.sub.i is the density of the i.sup.th layer thin film material, and {tilde over (σ)}.sub.i is the equivalent residual stress of the i.sup.th layer thin film material.

(20) For the initially flat and straight multilayer cantilever beam, the differential equation of transverse vibration of the micromechanical multilayer thin film cantilever beam is as follow:

(21) 0$\begin{array}{cc}\stackrel{\_}{\mathrm{EI}}\frac{{\partial }^{4}z\left(x,t\right)}{\partial x^4}=-\stackrel{\_}{\rho A}\frac{{\partial }^{2}z\left(x,t\right)}{\partial t^2}& \left(8\right)\end{array}$

(22) where z(x,t) is the height displacement when the beam transversely vibrates. The first-order resonance frequency formula of the multilayer cantilever beam which is kept flat and straight after release is obtained by solving the equation using the method of separation of variables, taking z(x,t)=Z(x)T(t) and substituting it into the equation (8) for solution:

(23) $\begin{array}{cc}\begin{array}{c}{f}_{1\mathrm{cf}}=\frac{{1.875}^2}{2\pi l^2}\sqrt{\frac{\stackrel{\_}{\mathrm{EI}}}{\stackrel{\_}{\rho A}}}\\ =\frac{{1.875}^2}{2\pi l^2}\sqrt{\frac{\underset{i=1}{\overset{N}{.Math.}}{\tilde{E}}_i{w}_i.Math.\left[{\left(z_i-z_c\right)}^3-{\left(z_{i-1}-z_c\right)}^3\right]}{3\underset{i=1}{\overset{N}{.Math.}}{\rho }_i{w}_i{h}_i}}\end{array}& \left(9\right)\end{array}$

(24) The subscript 1 of the frequency value represents the first-order resonance frequency, and cf represents that the boundary condition is cantilever beam.

(25) Based on the above-mentioned linear equation, the first-order resonance frequency can be respectively and actually measured by using a group of multilayer composite cantilever beam test structures with the same material parameters except the geometric size. According to the measured first-order resonance frequency and the size parameters of the test structures, a linear equation set can be established, and the equivalent Young modulus of each layer of the corresponding multilayer thin film structure can be obtained by solving the equation set.

(26) When the initially flat and straight multilayer thin film fixed-fixed beam vibrates near the equilibrium point, it is linear small-deformation vibration. When the height displacement of the beam is z(x,t) when it transversely vibrates, any cross section which is perpendicular to the length direction inside is not only subjected to the bending moment caused by the internal stress formed by bending and the bending moment generated by the load on the unit length, but also subjected to the bending moment caused by the axial force P to the offset from the balanced position of the cross section. The bending moment is Pz(x,t), and the above-mentioned relationship can be described as:
κEI=M+Pz(x,t)  (10)

(27) Under the condition of small-deformation vibration, the differential equation of the small-angle bending of the micromechanical thin film beam is as follow:

(28) $\begin{array}{cc}\kappa =\frac{d^{2}z\left(x,t\right)}{{\mathrm{dx}}^2}& \left(11\right)\end{array}$

(29) For a beam section with a length of dx and internal acting force and inertial acting force acting on it, the relationship among shear force V, bending moment M and load p(x) of the beam on the unit length is as follow:

(30) $\begin{array}{cc}\frac{d^{2}M}{{\mathrm{dx}}^2}=\frac{\mathrm{dV}}{\mathrm{dx}}=p\left(x\right)& \left(12\right)\end{array}$

(31) For a beam that transversely vibrates, the load of the beam on the unit length is equal to the load that overcomes the inertia force produced by the mass (i.e., linear density ρA) and acceleration of the beam on the unit length, i.e.:

(32) $\begin{array}{cc}p\left(x\right)=-\stackrel{\_}{\rho A}\frac{{\partial }^{2}z\left(x,t\right)}{\partial t^2}& \left(13\right)\end{array}$

(33) Formulas (10), (11) and (13) are substituted into formula (12) to finally obtain:

(34) $\begin{array}{cc}\stackrel{\_}{\mathrm{EI}}\frac{{\partial }^{4}z\left(x,t\right)}{\partial x^4}-\stackrel{\_}{\sigma A}\frac{{\partial }^{2}z\left(x,t\right)}{\partial x^2}=-\stackrel{\_}{\rho A}\frac{{\partial }^{2}z\left(x,t\right)}{\partial t^2}& \left(14\right)\end{array}$

(35) The equation is the differential equation of the transverse vibration of the micromechanical thin film fixed-fixed beam. The equation is solved by the method of separation of variables. It is allowed that z(x,t)=Z(x)T(t), it is substituted into the formula (14), and the general solution of amplitude Z(x) can be obtained as follow:

(36) $\begin{array}{cc}Z\left(x\right)=C_1\sin \left({\lambda }_{1}x\right)+C_1\cos \left({\lambda }_{1}x\right)+C_3\sinh \left({\lambda }_{2}x\right)+C_4\cosh \left({\lambda }_{2}x\right)\mathrm{where}{\lambda }_1=\sqrt{\left(\sqrt{{\beta }_1^2+4{\beta }_2}-{\beta }_1\right)\text{/}2},{\lambda }_2=\sqrt{\left(\sqrt{{\beta }_1^2+4{\beta }_2}+{\beta }_1\right)\text{/}2},{\beta }_1=\frac{\stackrel{\_}{\sigma A}}{\stackrel{\_}{\mathrm{EI}}},{\beta }_2=\frac{\stackrel{\_}{\rho A}}{\stackrel{\_}{\mathrm{EI}}}{\omega }^2.& \left(15\right)\end{array}$

(37) It can be considered that the first-order resonance mode of the thin film fixed-fixed beam is assumed to be the shape of a cosine function, i.e., Z.sub.1(x)=b[1−cos(2□x/l)]. For this form, there is the following relationship:

(38) $\begin{array}{cc}{Z}_1^{\left(4\right)}\left(x\right)=-\frac{4{\pi }^2}{l^2}{Z}_1^{″}\left(x\right)& \left(16\right)\end{array}$

(39) This relational expression is substituted into the equation of amplitude Z(x) to obtain:

(40) $\begin{array}{cc}\left[\stackrel{\_}{\mathrm{EI}}+\frac{l^2}{4{\pi }^2}\stackrel{\_}{\sigma A}\right]Z^{\left(4\right)}\left(x\right)-{\omega }^2\stackrel{\_}{\rho A}Z\left(x\right)=0& \left(17\right)\end{array}$

(41) By solving the formula, the approximate analytical formula of the first-order resonance frequency of the initially flat and straight multilayer thin film fixed-fixed beam with residual stress can be obtained:

(42) $\begin{array}{cc}{f}_{1s}=\frac{{4.730}^2}{2\pi l^2}\sqrt{\frac{\stackrel{\_}{\mathrm{EI}}}{\stackrel{\_}{\rho A}}+\frac{l^2}{4{\pi }^2}\frac{\stackrel{\_}{\sigma A}}{\stackrel{\_}{\rho A}}}& \left(18\right)\end{array}$

(43) The subscript 1 of the frequency value represents the first-order resonance frequency, and the subscript s represents that it is kept flat and straight after release.

(44) It can be seen from the above-mentioned formula that the first-order resonance frequency of the multilayer fixed-fixed beams is related to the material parameters and structural geometric size of the beams. On this basis, the first-order resonance frequency of the multilayer fixed-fixed beam can be regarded as an implicit function of Young modulus and residual stress determined by the above-mentioned formulas:

(45) 0$\begin{array}{cc}{f}_{1s}=\frac{{4.730}^2}{2\pi l^2}\sqrt{\frac{\underset{i=1}{\overset{N}{.Math.}}{\tilde{E}}_i{w}_i\left[{\left(\underset{i=0}{\overset{i}{.Math.}}{h}_i-z_c\right)}^3-{\left(\underset{i=0}{\overset{i-1}{.Math.}}{h}_i-z_c\right)}^3\right]}{3\underset{i=1}{\overset{N}{.Math.}}{\rho }_i{w}_i{h}_i}+\frac{l^2}{4{\pi }^2}\frac{\underset{i=1}{\overset{N}{.Math.}}{\tilde{\sigma }}_i{w}_i{h}_i}{\underset{i=1}{\overset{N}{.Math.}}{\rho }_i{w}_i{h}_i}}& \left(19\right)\end{array}$

(46) where, l is the length of the N-layer composite fixed-fixed beam; the width of the i.sup.th layer thin film material is w.sub.i, the thickness is h.sub.i, the density of the material is ρ.sub.i, the equivalent Young modulus is {tilde over (E)}.sub.i, the equivalent residual stress is {tilde over (σ)}.sub.i; and z.sub.c is the height of the neutral surface.

(47) Based on the implicit function mentioned above, the first-order resonance frequency can be respectively and actually measured by using a group of multilayer composite cantilever beam test structures with the same material parameters except the geometric size. According to the measured first-order resonance frequency and the size parameters of the test structures, a linear equation set can be established, and the equivalent Young modulus of each layer of the multiplayer composite fixed-fixed beam can be obtained by solving the equation set.

(48) However, considering the difficulty in solving the non-linear equation set, in order to reduce the calculation complexity, the equivalent Young modulus of each layer of the multilayer thin film structure can be obtained firstly by using the multilayer composite cantilever beam test structure, and then the equivalent residual stress of each layer of the multilayer thin film structure can be calculated by using the multilayer composite fixed-fixed beam test structure and the obtained equivalent Young modulus of each layer of the multilayer thin film structure. To sum up, the technical solution of the present invention can be obtained as follow:

(49) step 1: preparing a group of at least N initially flat and straight N-layer composite cantilever beam test structures and a group of at least N initially flat and straight N-layer composite fixed-fixed beam test structures by adopting a process for preparing the multilayer composite thin film structure, wherein the material parameters except the geometric size of each test structure are the same; and for each group of test structures, if a vector consisting of the width and length of each layer of each test structure is used as the dimension vector of the test structure, the dimension vector group of the group of test structures is linearly independent;

(50) step 2: measuring the first-order resonance frequency of each test structure;

(51) step 3: respectively substituting the first-order resonance frequency of the N-layer composite cantilever beam test structures into a first-order resonance frequency expression of N-layer composite cantilever beams, i.e., formula (9), to obtain an equation set consisting of at least N linear equations, and solving the equation set to obtain the equivalent Young modulus of each layer of the multilayer composite thin film structure; and

(52) step 4: respectively substituting the first-order resonance frequency of the N-layer composite fixed-fixed beam test structures and the equivalent Young modulus of each layer of the multilayer composite thin film structure obtained in step 3 into the first-order resonance frequency expression of the N-layer composite fixed-fixed beams, i.e., formula (19), to obtain an equation set consisting of at least N equations, and solving the equation set to obtain the equivalent residual stress of each layer of the multilayer composite thin film structure.

(53) In order to further help the public understand, the technical solution of the present invention will be further described below by taking the simplest two-layer thin film structure as an example.

(54) In step 1, a group of at least two initially flat and straight two-layer composite cantilever beam test structures and a group of at least two initially flat and straight two-layer composite fixed-fixed beam test structures are prepared by adopting a process for preparing the multilayer composite thin film structure, wherein the material parameters except the geometric size of each test structure are the same; and for each group of test structures, if a vector consisting of the width and length of each layer of each test structure is used as the dimension vector of the test structure, the dimension vector group of the group of test structures is linearly independent.

(55) The prepared two-layer composite cantilever beam test structure is illustrated in FIGS. 1a and 1b, the prepared two-layer composite fixed-fixed beam test structure is illustrated in FIG. 2a and FIG. 2b, each of which comprises a top layer 102 and a bottom layer 103. In the drawings, 101 represents an anchor area and 104 represents a substrate.

(56) For the two-layer composite cantilever beam consisting of two layers of thin film materials, there are two mechanical parameters to be solved. At this time, at least two two-layer composite cantilever beam test structures with different sizes are needed. If a vector consisting of the width and length of each layer of each test structure is used as the dimension vector of the test structure, in order to enable the established equation to have a solution, the dimension vectors of all test structures must be linearly independent.

(57) Similarly, for the two-layer composite fixed-fixed beam consisting of two layers of thin film materials, since the equivalent Young modulus of each layer is known, the equivalent residual stress of each layer remains unknown. At this time, at least two two-layer composite fixed-fixed beam test structures with different sizes are needed. Similarly, if a vector consisting of the width and length of each layer of each test structure is used as the dimension vector of the test structure, in order to enable the established equation to have a solution, the dimension vectors of all test structures must be linearly independent.

(58) In order to simplify the preparation and calculation process of the test structures, in the present embodiment, the remaining parts except the width of the bottom layers of the two test structures in each group are fully identical.

(59) In step 2, the first-order resonance frequency of each test structure is measured.

(60) The measurement of the first-order resonance frequency can be realized by various existing techniques, such as electrostatic excitation, thermal excitation and acoustic excitation. In the present embodiment, the resonance frequency of the test structure is measured by using the MSV-400-M2 Laser Doppler Vibrometer (LDV) of Polytech Company, Germany. The specific method is as follows: piezoelectric ceramic (PZT) is placed in the LDV test field of view, the sample chip to be tested is pasted on the PZT, and the PZT is driven to vibrate by applying a voltage signal. The chip transmits the mechanical vibration signal to the test structure under the excitation of the PZT, the excitation signal is transmitted to the beam structure through the anchor area of the test structure, and finally the substrate excitation form, which is also called as the inertial force excitation form, in the vibration mechanics is formed. Ideally, the frequency spectrum of the substrate excitation should be as smooth as possible, and the resonant beam will be forced to vibrate under uniform excitation in the whole band. Resonance will be generated near the eigenfrequency due to selective amplification, and a resonance peak will be formed in the vibration spectrum. It is supposed that the first-order resonance frequency of two two-layer composite cantilever beam test structures is respectively f.sub.1, f.sub.2, and the first-order resonance frequency of two two-layer composite fixed-fixed beam test structures are respectively f.sub.3, f.sub.4.

(61) In step 3, the first-order resonance frequency of the two-layer composite cantilever beam test structures is respectively substituted into a first-order resonance frequency expression of two-layer composite cantilever beams, i.e., formula (9), to obtain an equation set consisting of at least two linear equations, and the equation set is solved to obtain the equivalent Young modulus of each layer of the two-layer composite thin film structure.

(62) Assuming that the length of the two test structures is l.sub.1, the thickness of each layer is h, the top layer width is w.sub.1, the bottom layer width is respectively w.sub.2 and w3, and w.sub.2≠w.sub.3. According to formula (8), the following equation set can be obtained:

(63) $\{\begin{array}{c}{f}_{1\mathrm{cf}}\left({\tilde{E}}_1,{\tilde{E}}_2,w_1,w_{2.}{l}_1\right)=f_1\\ f_{1\mathrm{cf}}\left({\tilde{E}}_1,{\tilde{E}}_2,w_1,w_{3.}{l}_1\right)=f_2\end{array}\hspace{1em}$

(64) where {tilde over (E)}.sub.1, {tilde over (E)}.sub.2 are respectively equivalent Young modulus of the top layer and the bottom layer.

(65) Because the above-mentioned formula is a linear analytic formula for the first-order resonance frequency of two-layer cantilever beams, which corresponds to only one frequency value (i.e., the first-order resonance frequency), a finite set of real solutions and complex solutions for Young modulus can be obtained by solving the mechanical parameters of the material according to the known frequency and geometric size, wherein the complex solutions do not comply with the actual situation, but under the general situation, in engineering application, the only set of real solutions which complies with the value range of material parameters in the remaining several sets of real solutions can be identified clearly. The calculated equivalent Young modulus of the top layer and the bottom layer is the equivalent Young modulus of the upper layer and the lower layer of the two-layer composite thin film structures prepared by adopting the same preparation process.

(66) In step 4, the first-order resonance frequency of the two-layer composite fixed-fixed beam test structures and the equivalent Young modulus of each layer of the two-layer composite thin film structure obtained in step 3 are respectively substituted into the first-order resonance frequency expression of the two-layer composite fixed-fixed beams, i.e., formula (19), to obtain an equation set consisting of at least two equations, and the equation set is solved to obtain the equivalent residual stress of each layer of the two-layer composite thin film structure.

(67) Assuming that the length of two two-layer composite fixed-fixed beam test structures is l.sub.2, the top layer width is w.sub.4, and the bottom layer width of the two test structures is respectively w.sub.5 and w.sub.6, and w.sub.5≠w.sub.6, the equivalent Young modulus which has been solved is substituted into the equation set for solving, then the equivalent residual stress of each layer of the test structures can be obtained, and the equation set is as follow:

(68) $\{\begin{array}{c}{f}_{1s}\left({\tilde{E}}_1,{\tilde{E}}_2,{\tilde{\sigma }}_1,{\tilde{\sigma }}_2,w_4,w_5,l_2\right)=f_3\\ f_{1s}\left({\tilde{E}}_1,{\tilde{E}}_2,{\tilde{\sigma }}_1,{\tilde{\sigma }}_2,w_4,w_6,l_2\right)=f_4\end{array}\hspace{1em}$

(69) The calculated equivalent residual stress of the top layer and the bottom layer is the equivalent residual stress of the upper layer and the lower layer of the two-layer composite thin film structures prepared by adopting the same preparation process.