Next Article in Journal
Influence of the Bias Voltage on Effective Electron Velocity in AlGaN/GaN High Electron Mobility Transistors
Previous Article in Journal
Wettability Behaviour of Metal Surfaces after Sequential Nanosecond and Picosecond Laser Texturing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Analysis of Vibration Characteristics of Spatial Non-Uniform Tensioned Thin-Film Structures Based on the Absolute Nodal Coordinate Formulation

State Key Laboratory of Electromechanical Integrated Manufacturing of High-Performance Electronic Equipments, Xidian University, No. 2 South Taibai Road, Xi’an 710071, China
*
Authors to whom correspondence should be addressed.
Micromachines 2024, 15(9), 1147; https://doi.org/10.3390/mi15091147
Submission received: 19 June 2024 / Revised: 1 September 2024 / Accepted: 11 September 2024 / Published: 12 September 2024
(This article belongs to the Section E:Engineering and Technology)

Abstract

:
Due to their lightweight characteristics, spatial thin-film structures can generate vibrations far exceeding their film thickness when subjected to external loads, which has become a key factor limiting their performance. This study examines the vibration characteristics of tensioned membrane structures with non-uniform elements subjected to impacts in air, leveraging the Absolute Nodal Coordinate Formulation (ANCF). This model takes into account the wrinkling deformation of thin films under pre-tension and incorporates it into the dynamic equation derived using the absolute node coordinate method. A detailed discussion was conducted on the influence of non-uniform elements, situated at different locations and side lengths, on the vibration characteristics of the thin film. The analytical results obtained from the vibration model were compared with the experimental results, validating the effectiveness of the vibration model. This provides a theoretical foundation for the subsequent vibration control of thin films.

1. Introduction

Large-scale thin-film structures have received widespread international attention due to their lightweight and foldable characteristics [1,2,3,4]. The introduction of non-uniform structures into these thin-film structures has generated significant interest, such as in thin-film antenna structures and solar sail structures equipped with actuators in aerospace engineering [5,6,7]. Thin films are highly susceptible to vibrations far exceeding their thickness when subjected to impact loads and can develop wrinkles under preloads. The introduction of non-uniform structures not only affects the static characteristics of the thin film, such as stress distribution and wrinkling behavior, but also alters its vibration characteristics under impact [8,9,10]. Therefore, studying the vibration characteristics of non-uniform thin-film structures holds profound theoretical value for designing high-performance spatial thin-film structures.
Currently, scholars have conducted extensive research on the vibration characteristics of thin-film structures [11,12,13]. However, the current research is limited to the free vibration of uniform thin-film structures under no external force, and relatively little attention has been paid to the impact vibration characteristics of non-uniform tensioned thin-film structures. Li established a free vibration model for rectangular thin films considering the added mass of air based on the stress superposition method and discussed several characteristics related to the fundamental frequency and mode shapes of rectangular thin films [14]. Liu combined stability theory to establish a wrinkling prediction model for thin films under tension and further developed a dynamic model considering thin film wrinkling, describing the impact of wrinkling on the vibration characteristics of the thin-film structures [15]. Li proposed an eigenvalue method capable of analyzing the vibration characteristics of wrinkled and deformed thin films, based on the stress distribution pattern in the wrinkled region, and a dynamic model for rectangular sheared wrinkled thin films was established. The research results showed that the vibration frequency of the thin-film structure increased significantly with the vibration of wave peaks [16]. Liu studied the coupling effect between wrinkling and vibration and found that wrinkling has a significant impact on the dynamics of thin-film structures under asymmetric tensile loads [17]. Lu established a nonlinear dynamic model for tensioned thin-film antennas during motion, describing the influence of the rigid–flexible coupling effect on structural stiffness and damping characteristics [18]. Wang studied the free transverse vibration of annular wrinkled thin films, established a dimensionless Hamiltonian equation of motion for annular wrinkled thin films, and solved it using the finite difference method to obtain their vibration frequencies and mode shapes [19]. Fan established an analytical model for post-buckling membrane vibrations on a flexible substrate and systematically studied the analytical prediction of the system’s natural frequency and mode shape [20]. Wang established a linear self-vibration equation for flexible substrate thin films and studied the effects of substrate stiffness and damping on the natural frequency and vibration mode of the substrate thin films [21].
Previous studies have primarily focused on analyzing the vibration characteristics of the thin film itself, neglecting the vibration of the cables during the tensioning process. Therefore, they could not accurately capture the vibration characteristics of the entire thin-film structure. In the ANCF, similar to the classical finite element method, the object is discretized using finite elements, but the gradients of absolute position vectors are used as node coordinates instead of rotations. This description allows for a constant mass matrix with zero Coriolis and centrifugal forces. Conversely, the stiffness matrix is a highly nonlinear function of the element node coordinates [22,23].
Skrinjar used the absolute node coordinate method to model the prestressed rigid–flexible composite body in rigid–flexible assembly and studied the influence of different parameters on its dynamic response [24]. Ren proposed an adaptive triangular element based on the absolute node coordinate method, established a fine theory of thin plates and thin films, and analyzed the dynamic response of thin plates and the wrinkling and buckling behavior of thin films through an adaptive algorithm based on the stress state [25]. Wang conducted an unfolding simulation of foldable thin-film structures based on the absolute node coordinate method and integrated the nonlinear unfolding of foldable thin films into multi-body systems using the principle of virtual work [26]. Liu established a thin-film model using reduced plate elements and studied the wrinkling and buckling behavior of tensioned thin films and the dynamic response of rotating solar sails [27]. Yuan conducted a dynamic modeling of a folded-space thin-film structure subjected to contact impact during the unfolding process and discretized the thin film using triangular elements of the absolute node coordinate method to determine the stress–strain relationship of the thin film [28]. Zhao et al. conducted a modal analysis on rotating thin plates using the absolute nodal method and studied the effect of rotational angular velocity on natural frequency [29]. Wang proposed a multi-node flat plate element based on the absolute node coordinate method to calculate elastic and inertial functions and defined the displacement field using the bidirectional univariate Lagrange interpolation method, verifying the effectiveness of the ANCF element in static and dynamic analyses [30]. Olshevskiy proposed a fully parameterized ANCF flat plate element using only the transverse slope to simulate large displacements and rotations of thick plates [31]. Shababa defined a displacement rotation field based on the ANCF and used a first-order difference method to avoid the redundancy problem encountered in formulas with large rotation vectors [32]. Mikkola defined the relationship between B-spline surface control points and ANCF positions and gradient vectors through linear transformation, resulting in a unique displacement and rotation field that transforms B-spline surfaces into ANCF thin-plate finite elements [33]. Pappalardo developed a rational absolute nodal coordinate formula for thin-plate elements, which can be used for precise geometric modeling and the analysis of flexible continuum structures with complex geometries [34]. Grossi studied the high-frequency modes generated by the application of ANCF in multi-body systems and effectively solved rigid differentiation and memory groups by filtering and suppressing ANCF high-frequency modes, thereby reducing the negative impact of high-frequency modes on computational efficiency [35]. Pappalardo and Zhang developed a new fully parameterized ANCF triangular plate/shell element that defines the inertia and elasticity of the structure through element transformation between volume and Cartesian parameterization [36].
To obtain the dynamic characteristics of a multibody wrinkle thin-film system with cables, thin films, and non-uniform elements under impact, this study utilized reduced beam and thin-film elements to establish a multibody thin-film-structure dynamic model. Based on the non-uniform thin-film-wrinkling prediction model, the out-of-plane wrinkling deformation of thin films was obtained and introduced as an initial boundary into the dynamic model. Impact vibration analysis was performed on the tensioned wrinkled thin films with cables, and the vibration curves at the center of the thin film and the non-uniform element were extracted and compared with the experimental results to validate the effectiveness of the dynamic model.

2. ANCF Elements with Insufficient Gradients

2.1. Reducing Beam Elements

The full parameter beam element can be used to describe the deformation of the cross-section of Timoshenko beams. The cable structure in this article can be regarded as a Euler–Bernoulli beam, and the deformation of the cross-section is usually not considered. Therefore, a reduced beam element can be used for description. In the membrane structure under impact, the cables are mainly subjected to tensile deformation along their axis and vibration deformation along the normal direction of the membrane. Reduced beam elements can be used to model the cables. Divide the lasso into m elements, each containing two nodes and 12 degrees of freedom. Take any element within the lasso and establish a global coordinate system, as shown in Figure 1. Among them, r b A and r b B are the loss paths of node A and node B in the cable relative to the global coordinate origin O, respectively, γ is the inclination angle of the tangent at any point within the unit relative to the X-axis and γ . The gradient with respect to the position coordinate, denoted as r / x , should satisfy the following relationship:
cos γ = r x / x r x / x 2 + r y / x 2 , sin γ = r y / x r x / x 2 + r y / x 2 .
The shape function and node vector coordinates of the reduced beam element are defined as follows:
r b = S b x e b ,
where S b x is a reduced beam element form function matrix represented by spatial coordinate x, which can be expressed as follows:
S b x = S 1 b I 2 S 2 b I 2 S 3 b I 2 S 4 b I 2 ,
where the elements of the shape function can be expressed as follows:
S 1 b = 1 ξ b 2 + 2 ξ b 3 S 2 b = L ξ b 2 ξ b 2 + ξ b 3 S 3 b = 3 ξ b 2 2 ξ b 3 S 4 b = L ξ b 3 ξ b 2 ,
where ξ b = x / L , L is the length of the element and e b is the node coordinate vector composed of node positions and slope coordinates, which are absolute variables defined in the global coordinate system. For node A, the node coordinate vector is defined as follows:
e b A = r b A T r b A x T T = r b x A r b y A r b z A r b x A x r b y A x r b z A x T .
As shown in Figure 1, the node coordinate vector of the element can be represented as follows:
e b = e b A T e b B T T .
The strain energy of the reduced beam element is composed of two parts: axial deformation and bending deformation, and the specific expressions are as follows:
U b = 1 2 0 l E b A b ε b 2 + E b I b κ b 2 d x ,
where E b is the elastic modulus of the cable, A b is the cross-sectional area of the cable, and I b is the moment of inertia of the cable. The axial strain of cable ε b and curvature κ b can be expressed as follows:
ε b = 1 2 r b , x T r b , x 1 ,
κ b = r b , x × r b , x x r b , x 3 .
where r b , x is represented as r b taking a partial derivative of x and r b , x x is represented as r b taking the double partial derivative of x; this representation is universal throughout the entire text.
The elastic force of the reduced beam element can be expressed as follows:
Q e b = U b e .

2.2. Thin-Film Elements

Thin-plate elements can be used to describe the cross-sectional deformation of thin plates, while for thin films and non-uniform elements, their thickness is much smaller than the length in the other two directions, and their compression deformation in the thickness direction is generally not considered, so the slope vector in the thickness direction can be ignored. Taking the film as an example, the film element is divided into k 1 × k 2 elements, each containing four nodes and 36 degrees of freedom. Any element within the film is taken and a global coordinate system is established, as shown in Figure 2.
Among them, r m A , r m B , r m C , and r m D are the loss of diameter of nodes A, B, C, and D within the thin film element relative to the global coordinate origin, respectively. n is the normal vector of the middle plane of the thin film. Assuming that the normal vector perpendicular to the middle plane of the element before deformation remains perpendicular to the middle plane after deformation, the normal vector n can be expressed as follows:
n = r m , x × r m , y , n = r m , x × r m , y T r m , x × r m , y .
The shape function and node vector coordinates of the thin-film element can be represented as follows:
r m = S m x , y e m ,
where S m x , y is a thin-film element shape function matrix represented by spatial coordinates x and y, which can be expressed as follows:
S m x , y = S 1 m I 3 S 2 m I 3 S 12 m I 3 ,
where the elements of the shape function can be expressed as follows:
S 1 m = ξ m 1 η m 1 2 ξ m 2 + 2 η m 2 η m ξ m 1 , S 2 m = b ξ m ξ m 1 2 η m 1 , S 3 m = b η m ξ m 1 η m 1 2 , S 4 m = ξ m η m 1 2 ξ m 2 + 2 η m 2 η m 3 ξ m , S 5 m = b ξ m 2 ξ m 1 η m 1 , S 6 m = c ξ m η m η m 1 2 , S 7 m = η m ξ m 1 2 ξ m 2 + 2 η m 2 3 η m ξ m , S 8 m = b ξ m η m ξ m 1 2 , S 9 m = b η m 2 ξ m 1 η m 1 , S 10 m = ξ m η m 1 3 ξ m 3 η m + 2 ξ m 2 + 2 η m 2 , S 11 m = b ξ m 2 η m ξ m 1 2 , S 12 m = c ξ m η m 2 η m 1 ,
where ξ m = k 1 / b , η m = k 2 / c , b, and c are the length and width of the element, respectively.
For thin-film elements, the node coordinate vectors of node A and thin-film elements can be represented as follows:
e m A = r m A T r m A x T r m A y T r m A z T T   = r m x A r m y A r m z A r m x A x r m x A y r m x A z r m y A x r m y A y r m y A z r m z A x r m z A y r m z A z T ,
e m = e m A T e m B T e m C T e m D T T .
The elastic potential energy in vibration can be expressed as follows:
U 1 m = Ω E m t m 24 1 ν m 2 κ x x m 2 + κ y y m 2 + ν m E m t m 12 1 ν m 2 κ x x m κ y y m + E m t m 12 1 ν m 2 κ x y 2 d x d y .
The tensile strain potential energy is as follows:
U 2 m = E m t m 2 1 ν m 2 Ω ε x m 2 + ε y m 2 + 2 ν m ε x m ε y m 2 + 1 ν m 2 ε x y m 2 d x d y ,
where E m , ν m , and t m are the elastic modulus, Poisson’s ratio, and thickness of the thin film, respectively. According to continuous medium mechanics, the elastic force vector of an element is the partial derivative of the element’s strain energy with respect to the generalized coordinates. Therefore, the elastic force vector of the thin-film element can be represented as follows:
Q e m = U 1 m e m + U 2 m e m = Q 1 m + Q 2 m .
The strain and curvature vector can be expressed as follows:
ε i j m = 1 2 e m T S , i T S , j T e m δ i j , κ i j m = r m , i j T n n 3 .

2.3. Element Quality Matrix and External Force Matrix

In the ANCF, the element kinetic energy in global coordinates can be expressed as follows:
F e = 1 2 V e ρ r ˙ T r ˙ d V e ,
where ρ is the density of the element, and V e is the volume of the element. Substituting Equation (10) into Equation (19) gives the following:
F e = 1 2 e ˙ T M e e ˙ .
The element quality matrix can be expressed as follows:
M e = V e ρ S T S d V e .
Taking thin-film elements as an example, when an external force T acts on a thin-film element, the virtual work of that force can be expressed as follows:
δ W = T T δ r = T T S δ e = F e T δ e .
At this time, the external force matrix of the thin-film element can be represented as follows:
F e T = S x , y T T .

2.4. Element Rigid Body Motion

Taking the thin-film element as an example, the element node coordinate vector of the thin-film element during motion can be represented as follows:
e m A = r x A r y A r z A N 11 N 21 N 31 N 12 N 22 N 32 T , e m B = r x A + N 11 b r y A + N 21 b r z A + N 31 b N 11 N 21 N 31 N 12 N 22 N 32 , e m C = r x A + N 11 b + N 12 c r y A + N 21 b + N 22 c r z A + N 31 b + N 32 c N 11 N 21 N 31 N 12 N 22 N 32 , e m D = r x A + N 12 c r y A + N 22 c r z A + N 32 c N 11 N 21 N 31 N 12 N 22 N 32 ,
where N i j represents the i-th row and j-th column of the 3 × 3 rigid body rotation matrix. r x A , r y A , and r z A are the global translation position vectors of the middle surface of the thin-film element, respectively. Euler angles ϕ , η , and ψ represent the rotation angles of the thin-film element about the coordinate axes. Therefore, the rotation matrix N can be expressed as follows:
N = cos ψ cos ϕ cos η sin ϕ sin ψ sin ψ cos ϕ cos η sin ϕ cos ψ 0 cos ψ sin ϕ + cos η cos ϕ sin ψ sin ψ sin ϕ + cos η cos ϕ cos ψ 0 sin η sin ψ sin η cos ψ 0 .

3. Dynamic Equations of Thin-Film Structures Considering Wrinkles

3.1. Wrinkling Deformation under Pre-Tension

The schematic diagram of the non-uniform film structure is shown in Figure 3. We found through experimental research that the wrinkling morphology of non-uniform films is related to the size of non-uniform elements; we discussed this phenomenon in detail in reference [37]. For the thin-film structures of non-uniform elements with periodic distribution on the diagonal, when the size of the non-uniform elements a > 2 λ and a 2 λ , the wrinkling models of non-thin-films under uniform tension are as follows [37]:
w 1 = A 11 sin π ρ w R w 1 sin π θ λ + A 12 sin π ρ w l + R w 2 2 R w 2 sin π θ β λ , A 11 = e 3 t m x 2 A R w 1 2 , 0 , A 12 = e 3 t m ( x d + λ 2 2 ) 2 A l , β ,
w 2 = A 21 sin π ρ w R w cos π θ λ + A 22 sin π ρ w l R w 2 cos π θ β λ , A 21 = e 3 t m x λ 2 2 2 A R w 1 2 , 0 , A 22 = e 3 t m ( x 2 λ 1 + λ 2 2 ) 2 A l , β ,
where R w is the folding radius, which can be expressed as follows:
R w = exp ν m R , R w 1 = R c / sin α , R w 2 = R w R w 1 ,
where A i j is the amplitude of wrinkles, which can be expressed as follows:
A ρ w , θ = 2 λ 4 T ν m cos θ π + 2 π 2 E m ρ w t m 4 T cos θ π + 2 π 2 E m ρ w t m y d y ,
where λ is wrinkle radius can be expressed as follows:
λ = E m ρ w π 2 π + 2 t m 3 R w 2 48 T 1 ν m 2 cos θ 1 4 , λ 1 = R c sin α + d 2 tan λ R w , λ 2 = R c sin α + d 2 tan α n w + β d 2 ,
where n w is the number of wrinkles.
According to Equations (28)–(32), the wrinkling deformation of non-uniform thin films under uniform tension can be obtained, the wrinkled prediction results are shown in Figure 4, where reddish and bluish colors indicate the fold deformation of the film in the positive and negative directions, respectively. Finally, combined with Equations (26) and (27), it can be substituted as the initial form into the absolute node coordinate method.

3.2. The Handling of Constraint Conditions

There are motion constraints between thin films and cables, as well as between thin films and non-uniform components. The coordinates of two different element nodes at the end of the cable and the corner vertex of the film have the same spatial node displacement and different cross-sectional slopes. The coordinates of element nodes in contact with the film and non-uniform elements have the same spatial node displacement and cross-sectional slope. For the convenience of handling constraint conditions, the elastic force of the element can be expressed as the product of the stiffness matrix and node coordinates:
Q e = K e e .
Taking the connection between the cable and the thin film as an example, in order to illustrate the constraints between two nodes, each node only takes three vectors representing displacement for derivation, and assuming that they have the same y-direction displacement, the constrained equation is as follows:
k 11 0 k 13 k 14 k 15 + k 12 k 16 0 1 0 0 1 0 k 31 0 k 33 k 34 k 35 + k 32 k 36 k 41 0 k 43 k 44 k 45 + k 42 k 46 k 51 + k 12 1 k 53 + k 32 k 54 + k 42 k 55 + 2 k 52 + k 22 + 1 k 56 + k 62 k 61 0 k 63 k 64 k 65 + k 62 k 66 e 1 e 2 e 3 e 4 e 5 e 6 = f 1 0 f 3 f 4 f 5 f 2 f 6 .
The constraint conditions between thin films and non-uniform elements can be processed using the above methods for computer programming in the future.

3.3. Dynamic Equations

The elements are connected at the nodes, and the global stiffness matrix, global external force matrix, and global mass matrix of the thin-film structure can be expressed by using the Boolean matrix H between the element node coordinates and the generalized coordinates:
K = n H T K e e , F = n H T F e , M = n H T M e H .
According to the Lagrange equation, the dynamic equation of the thin-film structure can be obtained as follows:
M e ¨ + K e = F .
In Equations (35) and (36), e ¨ is the acceleration vector and M is a constant matrix used to describe the mass distribution inside the structure, which can be calculated and stored in advance.
Due to the fluid structure coupling effect, the mass matrix of the thin film structure will increase when the thin film vibrates in the air. Therefore, the Minami additional mass ratio formula was used to modify the dynamic equation of the thin film. The modified dynamic equation is as follows [38]:
1 + α m M e ¨ + K e = F ,
where α m is the additional mass ratio representing the ratio of air mass to the mass of the vibrating film in the air, regardless of vibration frequency and amplitude. The equivalent air layer height corresponding to the uniformly distributed additional mass in the film structure is approximately 68% of the length of the film structure, which can be expressed as follows:
α m = 0.68 ρ a l m / m m ,
where ρ a is the mass density of the air, m m is the mass of the thin-film structure, and l m is the length of the thin film.
When the film center is subjected to an impact force of 0.2 N, a dynamic model is established for the film structure with a rigid element size of 20 × 20 × 1   mm , a film size of 300 × 300 × 0.025   mm , and the cable is a steel cable with a diameter of 1 mm and a length of 100 mm. The parameters of the film and cable remain unchanged in the following text. The Newmark method was used to solve the dynamic equation [33], and the node deformation of the film system at different times is obtained as shown in Figure 5. The red-boxed portion of Figure 5 represents the displacement cloud diagram in the direction normal to the film, where reddish and bluish colors indicate the displacement of the film in the positive and negative directions, respectively. From the red box on the right side of Figure 5, it can be seen that the non-uniform film exhibited wrinkles during vibration.

4. Experimental Verification and Vibration Characteristic Analysis

4.1. Experimental Measurement System

In order to verify the effectiveness of the absolute node coordinate method used in this article, an experimental system was designed as shown in Figure 6. The laser rangefinder had a resolution of 0.002 mm, a sampling frequency of 500 Hz, and an impact effect was applied using a small, inverted pendulum.

4.2. The Influence of Non-Uniform Element Position on the Vibration Characteristics of Thin Films

The effect of non-uniform cell position on the dynamic characteristics of the thin-film structure was investigated such that the size parameter of the non-uniform element was 10 mm, the pre-tension was 15 N, and an impact force of 0.2 N was applied. The experimental results of the trend of the vibration frequency variation with the non-uniform element position parameter when the non-uniform unit position parameter was 80, 85, 90, 95, 100, 105, and 110 are shown in Figure 7. Then, the dynamic analysis of the film structure with non-uniform cell position parameters h of 80 mm, 90 mm, 100 mm, and 110 mm was carried out to extract the vibration curves of the film center and non-uniform cells, which were compared with the experimental results, as shown in Figure 8.
From Figure 7 and Figure 8, it can be seen that as the size parameter h of the film gradually increased, the vibration frequency of the film structure gradually increased. During vibration, the frequency of vibration between the transients of the thin-film structure changed slightly under the influence of the introduction of non-uniform elements as the positional parameter of the non-uniform elements, h, gradually increased, compared to when the position parameter of the non-uniform element was 80 mm; then, the vibration frequencies of non-uniform component position parameters decreased by 1.6%, 2.5%, and 5.1% at 90, 100, and 110, respectively. This was because when the position of the non-uniform element changed, the stress distribution of the thin-film structure changed accordingly, but the total mass of the thin-film structure remained unchanged. Therefore, the vibration frequency and amplitude of the thin-film structure underwent slight changes under the influence of the changed stress distribution.
To study the influence of non-uniform elements on the vibration characteristics of thin-film structures, the size and position parameters of non-uniform elements were kept constant. We compared the vibration curves of the center point of the non-uniform and uniform thin-film structures through experiments and compared the vibration curves of non-uniform element centers and uniform thin films at corresponding positions; the results are shown in Figure 9. According to the experimental results, the maximum and minimum amplitude differences of the last vibration cycle of the uniform film and non-uniform film at the center point of the film were 0.366 mm and 0.513 mm, respectively; The maximum and minimum amplitude differences of the last vibration cycle corresponding to the center point of the non-uniform element were 0.409 mm and 0.539 mm, respectively. It can be seen that non-uniform elements had a significant impact on the amplitude of the thin-film vibration. At the same time, the vibration frequency decreased by 6.82% compared to the non-uniform thin-film structure.

4.3. The Influence of Non-Uniform Element Size on the Vibration Characteristics of Thin Films

To investigate the effect of non-uniform cell size on the vibration characteristics of the film, the non-uniform element position parameter h was made to be 80 mm, the pre-tension was 15 N, and an impact force of 0.2 N was applied. The experimental results of the trend of the vibration frequency with the non-uniform element size parameter are shown in Figure 10 when the non-uniform element size parameters were 10, 12.5, 15, 17.5, and 20. Then, the dynamic analysis of the thin-film structure containing non-uniform elements of different sizes was carried out, respectively, and the vibration curves of the thin film and the center of the non-uniform unit were extracted when the non-uniform unit size parameters were 10, 15, and 20, respectively, and were compared with the experimental results, as shown in Figure 11. By comparing the analytical and experimental results, it can be seen that the vibration amplitude of the non-uniform film structure increased slowly with the increase of the non-uniform film size parameter a. Meanwhile, the vibration frequency of the film decreased by 4.88% and 12.2% when the non-uniform element size parameters were 15 and 20 mm, respectively.

4.4. The Influence of Pre-Tension on the Vibration Characteristics of Non-Uniform Thin Films

To study the effect of non-uniform cell size on the vibration characteristics of the film, an impact force of 0.2 N was applied such that the position parameter h was 80 mm, and the cell size parameter was 10 mm. When the applied pre-tension was 10 N, 12.5 N, 15 N, 17.5 N, and 20 N, the experimental results of the vibration frequency with the trend of the non-uniform film size parameter are shown in Figure 12, and the experimental results of the maximum vibration amplitude with the trend of the non-uniform cell size parameter are shown in Figure 13. Then, the dynamic analysis of the film structure was carried out and the vibration curves of the film and the center of the non-uniform unit at pre-tensions of 10 N, 15 N, and 20 N were extracted and compared with the experimental results, which are shown in Figure 14. It can be concluded that when the pre-tension was set to 15 N and 20 N, the vibration frequency increased by 20.58% and 38.24%, respectively, compared to the non-uniform film with a pre-tension of 10 N. Meanwhile, the vibration amplitude of the non-uniform film structure decreased significantly when the pre-tension was increased, and the maximum vibration amplitude at the center of the film decreased by 14.5% and 35.7%. The maximum amplitude at the center of the non-uniform element decreased by 23.6% and 40.8%, respectively, when the pre-tension was up to 15 N and 20 N, as compared to non-uniform films with a pre-tension of 10 N. Therefore, the selection of effective pre-tensioning is essential to improve the performance of the films.

5. Conclusions

This article establishes a dynamic model of non-uniform tensioned thin-film structures based on the absolute node coordinate method, taking into account the wrinkling deformation of the thin films under pre-tension and the additional air mass when the thin film vibrates in the air. The effectiveness of the model was verified through experiments. The influence of non-uniform element position parameters, size parameters, and film pre-tension on the vibration characteristics of thin-film structures was studied. According to the analysis results, the position and size of non-uniform elements had relatively little effect on the vibration frequency of the thin film. When the position parameter of the thin film increased from 80 mm to 90 mm, 100 mm, and 110 mm, the vibration frequency of the thin film decreased by 1.6%, 2.5%, and 5.1%, respectively. When the size parameter of the film increased from 10 mm to 15 mm and 20 mm, the vibration frequency of the film decreased by 4.8% and 12.2%, respectively. The introduction of non-uniform elements had a significant impact on the vibration amplitude of the thin film, but a relatively small impact on the vibration frequency of the thin film. When the position parameters, size parameters, and pre-tension of the non-uniform elements were 80 mm, 15 mm, and 15 N, respectively. Compared to the uniform thin film structure, the introduction of non-uniform elements increased the vibration frequency of the thin film by 7%, and the vibration amplitudes at the center of the thin film and the center of the non-uniform element increased by 40.2% and 47.3%, respectively. The pre-tension of the film had a significant impact on the vibration amplitude and frequency of the non-uniform films. When the position and size parameters of the non-uniform element were 80 mm and 15 mm, respectively, and the pre-tension of the film increased from 10 N to 15 N and 20 N, the vibration frequency of the film increased by 20.58% and 38.24%, respectively. The maximum vibration amplitude of the film at the center point decreased by 14.5% and 35.7%, respectively, and the maximum vibration amplitude at the center point of the non-uniform element decreased by 23.6% and 40.8%, respectively. Therefore, when designing thin-film structures, it is necessary to choose reasonable non-uniform element sizes and thin-film pre-tensions to reduce the impact on the performance of the thin film structure.

Author Contributions

Conceptualization, P.S.; methodology, P.S.; validation, P.S.; formal analysis, P.S.; investigation, P.S.; resources, J.H., J.Z., F.M. and P.Z.; data curation, P.S.; writing—original draft preparation, P.S.; writing—review and editing, P.S. and J.H.; visualization, P.S.; supervision, J.H.; project administration, J.H.; funding acquisition, J.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China under Grant 52035010, in part by the Shaanxi Innovation Team Project under Grant 2018TD-012, in part by the Shaanxi Key Industry Chain Project under Grant 2020ZDLGY14-08, in part by the National 111 Project under Grant B14042, in part by the Shaanxi Provincial Fund under Grant 2022JQ-366, in part by the National Natural Science Foundation of China under Grant 52205411, in part by the National Key Laboratory Fund under Grant 2022-JCJQ-LB-018, and in part by the National Natural Science Foundation of China under Grant 52275372.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding authors.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Liu, M.J.; Huang, J.; Wang, Y. Analysis of wrinkled membrane structures based on a wrinkle-wave model. AIP Adv. 2017, 7, 1–7. [Google Scholar] [CrossRef]
  2. Chandra, M.; Kumar, S.; Chattopadhyaya, S.; Chatterjee, S.; Kumar, P. A review on developments of deployable membrane-based reflector antennas. Adv. Space Res. 2021, 68, 3749–3764. [Google Scholar] [CrossRef]
  3. Dai, X.J.; Yuan, T.Y.; Zheng, Z.; Shao, X.X.; Li, L.J.; Cheng, X.M.; Zhou, J.L.; Yang, F.J.; He, X.Y. Experimental study of wrinkling behavior of membrane structures via visual method. Thin Wall Struct. 2019, 149, 106537. [Google Scholar] [CrossRef]
  4. Attipou, K.; Hu, H.; Potier-Ferry, M.; Mohri, F.; Potier-Ferry, M.; Belouettar, S. Thermal wrinkling of thin membranes using a Fourier-related double scale approach. Thin Wall Struct. 2015, 94, 532–544. [Google Scholar] [CrossRef]
  5. Luo, Y.J.; Xing, J.; Kang, Z.; Zhan, J.J.; Li, M. Uncertainty of membrane wrinkling behaviors considering initial thickness imperfections. Int. J. Solids Struct. 2020, 191–192, 264–277. [Google Scholar] [CrossRef]
  6. Huang, J.; Lou, A.; Feria, A.; Kim, Y. An inflatable L-band microstrip SAR array. In Proceedings of the IEEE Antennas and Propagation Society International Symposium, Atlanta, GA, USA, 21–26 June 1998. [Google Scholar]
  7. Fleurent-Wilson, E.; Pollock, T.E.; Su, W.J.; Warrier, D.; Salehian, A. Wrinkle localization in membrane structures patched with macro-fiber composite actuators: Inflatable space antenna applications. J. Intel. Mat. Syst. Str. 2014, 25, 1978–2009. [Google Scholar] [CrossRef]
  8. Wang, C.G.; Tan, H.F.; Du, X.W.; Wan, Z.M. Wrinkling prediction of rectangular shell-membrane under transverse in-plane displacement. Int. J. Solids Struct. 2007, 50, 6507–6516. [Google Scholar] [CrossRef]
  9. Sun, P.; Huang, J.; Zhang, J.Y.; Meng, F.B. Wrinkling patterns and stress analysis of tensile membrane with rigid elements. Appl. Sci. 2022, 12, 6630. [Google Scholar] [CrossRef]
  10. Sun, P.; Huang, J.; Zhang, J.Y.; Meng, F.B. Wrinkling suppression in thin film using position distribution of microstructures. Int. J. Mech. Mater. Des. 2024, 20, 3–13. [Google Scholar] [CrossRef]
  11. Shao, Q.; Lu, Y.F.; Shi, C.; Yue, H.H.; Lv, L.L.; Liu, R.Q.; Fang, G.Q. Hybrid vibration control of the membrane antenna experiencing maneuver with cable actuators. Thin Wall Struct. 2023, 192, 111149. [Google Scholar] [CrossRef]
  12. Zhang, J.H.; Wu, N.; Tong, A.; Liu, Y.H. Structural dynamic responses of a stripped solar sail subjected to solar radiation pressure. Chin. J. Aeronaut. 2020, 33, 2204–2211. [Google Scholar] [CrossRef]
  13. Zhang, F.; Gong, S.P.; Gong, H.R.; Baoyin, H. Solar sail attitude control using shape variation of booms. Chin. J. Aeronaut. 2022, 35, 326–336. [Google Scholar] [CrossRef]
  14. Zhang, X.X.; Zhou, C.Y. Dynamic analysis of spinning solar sails at deployment process. Chin. J. Aeronaut. 2017, 30, 1719–1728. [Google Scholar] [CrossRef]
  15. Li, B.Y.; Liu, R.Q.; Qiang, C.; Guo, H.W.; Lin, Q.H. Stress Superposition Method and free vibration of corner tensioned rectangular thin membranes. Thin Wall Struct. 2021, 159, 107201. [Google Scholar] [CrossRef]
  16. Liu, M.J.; Huang, J.; Liu, M.Y. Multiple-model switching control for vibration suppression of planar membrane structures. Adv. Mech. Eng. 2019, 11, 168781401988377. [Google Scholar] [CrossRef]
  17. Li, Y.L.; Lu, M.Y.; Tan, H.F.; Tan, Y.Q. A study on wrinkling characteristics and dynamic mechanical behavior of membrane. Acta Mech. Sin. 2012, 28, 201–210. [Google Scholar] [CrossRef]
  18. Liu, X.; Cai, G.P. Nonlinear dynamic analysis of wrinkled membrane structure. Eng. Comput. 2023, 40, 41–46. [Google Scholar] [CrossRef]
  19. Lu, Y.F.; Shao, Q.; Lv, L.L.; Yue, H.H. Nonlinear Dynamics of a Space Tensioned Membrane Antenna during Orbital Maneuvering. Aerospace 2022, 9, 794. [Google Scholar] [CrossRef]
  20. Wang, C.G.; Liu, Y.P.; Lan, L.; Tan, H.F. Free transverse vibration of a wrinkled annular thin film by using finite difference method. J. Sound. Vib. 2016, 363, 272–284. [Google Scholar] [CrossRef]
  21. Fan, X.; Wang, Y.; Li, Y.; Fu, H. Vibration Analysis of Post-Buckled Thin Film on Compliant Substrates. Sensors 2020, 20, 5425. [Google Scholar] [CrossRef]
  22. Wang, Y.; Cui, X.B.; Fu, H.R.; Zhao, Q.; Li, Y.H. Dynamic Behaviors of Postbuckled Thin Film on Flexible Substrates Considering Viscoelastic Effects. J. Appl. Mech-T Asme 2021, 88, 041007. [Google Scholar] [CrossRef]
  23. Dufva, K.; Shabana, A.A. Analysis of thin plate structures using the absolute nodal coordinate formulation. Proc. Inst. Mech. Eng. Part K J. Multi-Body Dyn. 2005, 219, 345–355. [Google Scholar] [CrossRef]
  24. García-Vallejo, D.; Mayo, J.; Escalona, J.L.; Escalona, J. Efficient Evaluation of the Elastic Forces and the Jacobian in the Absolute Nodal Coordinate Formulation. Nonlinear Dyn. 2004, 35, 313–329. [Google Scholar] [CrossRef]
  25. Skrinjar, L.; Slavic, J.; Boltezar, M. Absolute Nodal Coordinate Formulation in a Pre-Stressed Large-Displacements Dynamical System. Stroj. Vestn. J. Mech. Eng. 2017, 63, 417–425. [Google Scholar] [CrossRef]
  26. Ren, H.; Fan, W. An adaptive triangular element of absolute nodal coordinate formulation for thin plates and membranes. Thin Wall Struct. 2023, 182, 110257. [Google Scholar] [CrossRef]
  27. Wang, T.F.; Wu, Z.Y.; Wang, J.; Lan, P.; Xu, M.L. Simulation of membrane deployment accounting for the nonlinear crease effect based on absolute nodal coordinate formulation. Nonlinear Dyn. 2023, 111, 2521–2535. [Google Scholar] [CrossRef]
  28. Liu, C.; Tian, Q.; Yan, D.; Hu, H.Y. Dynamic analysis of membrane systems undergoing overall motions, large deformations and wrinkles via thin shell elements of ANCF. Comput. Method. Appl. M. 2013, 258, 81–95. [Google Scholar] [CrossRef]
  29. Zhao, J.; Tian, Q.; Hu, H.Y. Modal Analysis of a Rotating Thin Plate via Absolute Nodal Coordinate Formulation. J. Comput. Nonlin Dyn. 2011, 6, 041013. [Google Scholar] [CrossRef]
  30. Wang, J.; Wang, T.F.; Zhang, Y.L.; Bian, H.Y.; Liu, W.J. A novel plate element based on absolute nodal coordinate formulation with collocation strategy. Acta Mech. 2023, 234, 6677–6690. [Google Scholar] [CrossRef]
  31. Olshevskiy, A.; Olshevskiy, O.; Dai, M.D.; Kim, C.W. The simplest 3-, 6- and 8-noded fully-parameterized ANCF plate elements using only transverse slopes. Multibody Syst. Dyn. 2015, 34, 23–51. [Google Scholar] [CrossRef]
  32. Mikkola, A.; Shabana, A.A.; Sanchez-Rebollo, C.; Jimenez-Octavio, J.R. Comparison between ANCF and B-spline surfaces. Multibody Syst. Dyn. 2013, 30, 119–138. [Google Scholar] [CrossRef]
  33. Shabana, A.A. ANCF Consistent Rotation-Based Finite Element Formulation. J. Comput. Nonlin. Dyn. 2016, 11, 014502. [Google Scholar] [CrossRef]
  34. Pappalardo, C.M.; Yu, Z.Q.; Zhang, X.S.; Shabana, A.A. Rational ANCF Thin Plate Finite Element. J. Comput. Nonlin. Dyn. 2016, 11, 051009. [Google Scholar] [CrossRef]
  35. Grossi, E.; Shabana, A.A. Analysis of high-frequency ANCF modes: Navier–Stokes physical damping and implicit numerical integration. Acta Mech. 2019, 230, 2581–2605. [Google Scholar] [CrossRef]
  36. Pappalardo, C.M.; Zhang, Z.; Shabana, A.A. Use of independent volume parameters in the development of new large displacement ANCF triangular plate/shell elements. Nonlinear Dyn. 2018, 91, 2171–2202. [Google Scholar] [CrossRef]
  37. Sun, P.; Huang, J.; Zhang, J.Y.; Meng, F.B.; Zhao, P.B. Prediction of wrinkle patterns in tensioned thin-film structures containing rigid elements. Aeronaut. J. 2024, 1–11. [Google Scholar] [CrossRef]
  38. Minami, H. Added mass of a membrane vibrating at finite amplitude. J. Fluid. Struct. 1998, 12, 919–932. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of reduced beam elements.
Figure 1. Schematic diagram of reduced beam elements.
Micromachines 15 01147 g001
Figure 2. Schematic diagram of thin-film element.
Figure 2. Schematic diagram of thin-film element.
Micromachines 15 01147 g002
Figure 3. Schematic diagram of film and wrinkle parameters.
Figure 3. Schematic diagram of film and wrinkle parameters.
Micromachines 15 01147 g003
Figure 4. Schematic diagram of wrinkle deformation.
Figure 4. Schematic diagram of wrinkle deformation.
Micromachines 15 01147 g004
Figure 5. Node deformation of the thin-film system at different times.
Figure 5. Node deformation of the thin-film system at different times.
Micromachines 15 01147 g005
Figure 6. Measurement system for the thin-film impact vibration experiment: (a) is a non-uniform thin-film structure and a laser rangefinder; (b) is the method of applying the impact force.
Figure 6. Measurement system for the thin-film impact vibration experiment: (a) is a non-uniform thin-film structure and a laser rangefinder; (b) is the method of applying the impact force.
Micromachines 15 01147 g006
Figure 7. The trend of vibration frequency changing with non-uniform element position parameters.
Figure 7. The trend of vibration frequency changing with non-uniform element position parameters.
Micromachines 15 01147 g007
Figure 8. Vibration curves of thin-film structures with non-uniform elements at different positions: (a,c,e,g) are the vibration curves of the center point of the thin film when the non-uniform element position parameter h is 80 mm, 90 mm, 100 mm, and 110 mm, respectively; (b,d,f,h) are the vibration curves of the center point of the film when the non-uniform element position parameter h is 80 mm, 90 mm, 100 mm, and 110 mm, respectively.
Figure 8. Vibration curves of thin-film structures with non-uniform elements at different positions: (a,c,e,g) are the vibration curves of the center point of the thin film when the non-uniform element position parameter h is 80 mm, 90 mm, 100 mm, and 110 mm, respectively; (b,d,f,h) are the vibration curves of the center point of the film when the non-uniform element position parameter h is 80 mm, 90 mm, 100 mm, and 110 mm, respectively.
Micromachines 15 01147 g008
Figure 9. Vibration curves of uniform thin films and non-uniform thin films: (a,b) are the vibration curves of the center of the thin film and the center point of the non-uniform element, respectively.
Figure 9. Vibration curves of uniform thin films and non-uniform thin films: (a,b) are the vibration curves of the center of the thin film and the center point of the non-uniform element, respectively.
Micromachines 15 01147 g009
Figure 10. The trend of vibration frequency changing with non-uniform element size parameters.
Figure 10. The trend of vibration frequency changing with non-uniform element size parameters.
Micromachines 15 01147 g010
Figure 11. Vibration curves of non-uniform elements at different size parameters: (a,c,e) are the vibration curves of the center point of the thin film when the non-uniform element size parameter h is 10 mm, 15 mm, and 20 mm, respectively; (b,d,f) are the vibration curves of the center point of the film when the non-uniform element size parameter h is 10 mm, 15 mm, and 20 mm, respectively.
Figure 11. Vibration curves of non-uniform elements at different size parameters: (a,c,e) are the vibration curves of the center point of the thin film when the non-uniform element size parameter h is 10 mm, 15 mm, and 20 mm, respectively; (b,d,f) are the vibration curves of the center point of the film when the non-uniform element size parameter h is 10 mm, 15 mm, and 20 mm, respectively.
Micromachines 15 01147 g011
Figure 12. The trend of vibration frequency changing with pre-tension of thin film.
Figure 12. The trend of vibration frequency changing with pre-tension of thin film.
Micromachines 15 01147 g012
Figure 13. The trend of maximum vibration amplitude changing with pre-tension of thin film.
Figure 13. The trend of maximum vibration amplitude changing with pre-tension of thin film.
Micromachines 15 01147 g013
Figure 14. Vibration curves of non-uniform thin films under different pre-tensions: (a,c,e) are the vibration curves of the center point of the film under tensions of 10 N, 15 N, and 20 N, respectively; (b,d,f) are the vibration curves of the center point of the film under tensions of 10 N, 15 mm, and 20 N, respectively.
Figure 14. Vibration curves of non-uniform thin films under different pre-tensions: (a,c,e) are the vibration curves of the center point of the film under tensions of 10 N, 15 N, and 20 N, respectively; (b,d,f) are the vibration curves of the center point of the film under tensions of 10 N, 15 mm, and 20 N, respectively.
Micromachines 15 01147 g014
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sun, P.; Huang, J.; Zhang, J.; Meng, F.; Zhao, P. Analysis of Vibration Characteristics of Spatial Non-Uniform Tensioned Thin-Film Structures Based on the Absolute Nodal Coordinate Formulation. Micromachines 2024, 15, 1147. https://doi.org/10.3390/mi15091147

AMA Style

Sun P, Huang J, Zhang J, Meng F, Zhao P. Analysis of Vibration Characteristics of Spatial Non-Uniform Tensioned Thin-Film Structures Based on the Absolute Nodal Coordinate Formulation. Micromachines. 2024; 15(9):1147. https://doi.org/10.3390/mi15091147

Chicago/Turabian Style

Sun, Peng, Jin Huang, Jiaying Zhang, Fanbo Meng, and Pengbing Zhao. 2024. "Analysis of Vibration Characteristics of Spatial Non-Uniform Tensioned Thin-Film Structures Based on the Absolute Nodal Coordinate Formulation" Micromachines 15, no. 9: 1147. https://doi.org/10.3390/mi15091147

APA Style

Sun, P., Huang, J., Zhang, J., Meng, F., & Zhao, P. (2024). Analysis of Vibration Characteristics of Spatial Non-Uniform Tensioned Thin-Film Structures Based on the Absolute Nodal Coordinate Formulation. Micromachines, 15(9), 1147. https://doi.org/10.3390/mi15091147

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop