Next Article in Journal
Myoepithelioma of the Hand: A Systematic Review
Previous Article in Journal
Nonlinear Analysis and Closed-Form Solution for Overhead Line Magnetic Energy Harvester Behavior
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Modeling and Vibration Characteristic Analysis of Fiber Woven Composite Shaft–Disk Rotor with Weight-Reducing Holes

by
Haibiao Zhang
1,2,3,
Mengyu Shen
1,2,
Tao Liu
4,*,
Zhen Li
1,2 and
Qingshan Wang
1,2,*
1
College of Mechanical and Electrical Engineering, Central South University, Changsha 410083, China
2
State Key Laboratory of Precision Manufacturing for Extreme Service Performance, Central South University, Changsha 410083, China
3
Hunan Aviation Powerplant Research Institute, AECC, Zhuzhou 412002, China
4
School of Civil Engineering, Central South University, Changsha 410075, China
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2024, 14(19), 9148; https://doi.org/10.3390/app14199148
Submission received: 29 August 2024 / Revised: 27 September 2024 / Accepted: 4 October 2024 / Published: 9 October 2024

Abstract

:
In order to achieve the goal of a lightweight shaft–disk rotor, this paper applies the fiber woven composite material to the disk structure, and at the same time considers the design of the weight-reducing holes on the porous disk. It introduces the domain decomposition and coordinate mapping technology for this, and then establishes the dynamics model of the fiber woven composite material shaft–disk rotor. The model is based on the differential quadrature finite element method, which is suitable for fiber woven composite rotors with arbitrary complex hole patterns. The validity of the model is verified by comparing the results with the literature, finite element simulation, and experiments, and the mechanism of the influence of the material parameters and pore parameters on the vibration characteristics of the system is investigated, which provides the data support and theoretical basis for the analysis of the dynamics of the fiber woven composite rotor.

1. Introduction

In practical engineering applications, in order to reduce the weight of the shaft–disk rotor, improve the vibration noise phenomenon, eliminate internal stress, and other purposes, the rotor disk is often used to open the weight reduction holes, as well as for the replacement of rotor materials and other methods. At this time, the performance enhancement effect brought about by the material change and the hole parameters on the vibration performance of the system are important aspects of rotor dynamics research.
Among the existing emerging composites, fiber woven composites, in addition to possessing many advantages such as low density, excellent physical and mechanical properties, good heat resistance, etc., also have better designability, processability, and lower manufacturing cost compared to ordinary fiber-reinforced, functional gradient, and other composites, which have broad application prospects in the fields of aerospace and rail transportation. Therefore, this material is one of the ideal materials for shaft–disk rotor preparation.
The current research on fiber woven composites is mainly based on simulation and experimentation, and generally adopts the multi-scale modeling method [1] to achieve the modeling of fiber woven composites from the microscopic scale to the mesoscopic scale (representing the volumetric unit) and then to the macroscopic scale. The research includes the physical and mechanical properties [2] and the damage behaviors [3,4] of the fiber woven composites and their structures. In terms of dynamics research, single structures such as fiber woven composite plates/shells are mostly taken as the analysis objects. Poojary and Hegde [5] carried out an experimental study on the dynamic properties of intra-ply woven carbon fiber/epoxy hybrid composites using fiber orientation as a comparative variable and found that the structural intrinsic frequency decreases with the increase in the orientation angle. Damijan et al. [6] fabricated woven carbon fiber-reinforced composites gears by using a hot-press curing method, and found that, after high week fatigue testing, woven carbon fiber-reinforced composite gears have longer fatigue lives and better mechanical properties than ordinary composite gears. The above studies are based on experiments, which are highly referable and relevant, but this also results in the singularity of the research object. Huang et al. [7] applied the three-dimensional finite element method to model the carbon fiber plain weave material and used a two-step homogenization method to obtain the mechanical properties of the material, then investigated the free vibration characteristics of the carbon fiber plain weave cylindrical shells and came up with the conclusion that increasing the total fiber volume fraction leads to a linear increase in the intrinsic frequency of the structure It can be seen that the related dynamics research also mostly adopts finite element software, which has the disadvantages of difficult modeling and being time-consuming when facing complex research objects.
For rotating members, open holes, while affecting the internal mechanical properties of the structure [8], also change their dynamic properties, such as natural frequency and critical speed. Therefore, the problem of modeling the dynamics of open-hole structures and their dynamic properties has attracted the attention of many scholars. In addition to numerical analytical methods to characterize the hole boundaries, existing studies also treat the holes as virtual structures. Wang et al. [9,10] treated irregular holes as virtual plates with zero density and Young’s modulus, treated a perforated plate as an equivalent rectangular plate with inhomogeneous thickness, and solved the free vibration characteristics of the structure by using a spectral geometry method. Furthermore, subdomain decomposition of a porous plate/shell, thus converting the open-hole problem into a coupled problem of multiple non-hole plates/shells, is also a common approach to porous structures. Such methods are not limited to specific hole types and specific plate and shell structures, and are highly applicable and operable. Kwak [11] took laminated closed cones, cylindrical shells, and annular plates containing regular holes as the object of study and proposed a new, meshless method to solve their free vibration characteristics, in which the open-hole structure is decomposed into multiple non-hole unclosed shells, which are mainly aimed at simple collocated structures with regular hole types and plates/shells. Song [12] combined the artificially connected spring technique and the domain partition integration method to develop a free vibration analytical model for arbitrarily shaped plates with different cutouts, and the analysis found that the variation in the hole size would have an irregular effect on the natural frequency of the plate. Sun et al. [13], based on the proportional boundary finite element method, the exact integration method, and the degree of freedom transformation technique, performed free vibration analysis on single and laminated single and stacked composite plates with multiple cutouts, and the results were summarized in the following table. The free vibration analysis of composite panels with multiple cut-outs is summarized by the fact that the boundary conditions, thickness-to-length ratio, and hole size have significant effects on the distribution of the intrinsic frequency of the panels with the cut-outs. Zhong et al. [14,15] divided the open-pore plywood into multiple quadrilateral/triangular units, solved the unit energy equations by means of coordinate transformation, coupled the units with the help of virtual springs, and then deduced the uniform in-plane vibration solutions for irregular plywood with complex holes under different boundaries.
Rotor systems are widely used in aviation, shipbuilding, power generation, and other fields. The related dynamics analysis is an important part of its dynamics design. The thin-disk rotor is a typical structure in the rotor system, which often works at a high rotational speed and in a high load-bearing environment. The traditional mass-point equivalent simulation of the thin disk modeling method is no longer applicable, and there is a need to take into account the flexibility of the thin-disk structure, which puts more demand on the dynamics of the rotor-modeling methodology and the breadth of the analysis. This puts forward higher requirements for rotor dynamics modeling methods and analysis. For the shaft–disk rotor, the existing research focuses on investigating its dynamic characteristics, such as vibration [16], response [17] and nonlinearity [18,19]. Heydari and Khorram [20] treated both the shaft and disk as flexible bodies and used the assumed modal method to develop a vibration model of the rotor system. They pointed out through analysis that the position of the flexible disk had a significant effect on the system’s intrinsic frequency and critical speed. Tuzzi et al. [21], on the other hand, investigated the coupling behavior of the flexible shaft–disk assembly under asymmetric support conditions. Zhao et al. [22] proposed a unified modeling method for rotating flexible shaft–disk systems based on the Kirchhoff and Euler–Bernoulli beam theories, as well as a unified modeling approach for rotating flexible shaft–disk systems. Furthermore, the dynamic modeling and characterization of complex systems, such as the shaft–disc–vane system [23,24,25], shaft–disc–drum system [26,27], and cracked rotor [28,29,30], are also directions of rotor research that are constantly being more deeply investigated.
With the rise of new materials with excellent mechanical properties, such as composites, scholars have begun to consider changing the materials of the rotor structure to meet the requirements of the system, which include light weight and high speed. Venkatachalam et al. [31,32] applied orthotropic anisotropic materials to the shaft–disk rotor structure, and, with the help of experiments and finite element methods, found that the use of composites with a fiber-laying direction of 90° could significantly reduce the vibration of the rotor. Wang et al. [33,34] proposed and validated a dynamic model of a flexible functional gradient rotor system with a variable thickness disc and comprehensively analyzed the effects of geometrical parameters, material parameters, and the shape of the disc cross-section on its vibration characteristics. They pointed out that the material parameters would affect the overall vibration characteristics of the system to a certain extent. Zhao et al. [35] proposed a modeling method for a graphene nanosheet-reinforced blade–disk system and found that it was possible to use an orthotropic anisotropic material for the rotor structure with the help of experimental and finite element methods. They investigated the impact of graphene nanosheet parameters on the system and concluded that incorporating graphene sheets can enhance the mechanical properties of the system. Furthermore, thinner and larger graphene sheets are observed to increase the structural stiffness. Cai et al. [36] also applied graphene nanosheets and porous copper–metal foams to a rotating shaft–disk blade system and performed resonance analyses and parameter interval uncertainty analyses. The effects of porosity and material parameters were comprehensively explored.
In summary, in the application of fiber woven composites, most of the simple components, such as plates, shells, beams, etc., are the macroscopic carriers of the materials, and their kinetic properties are either studied by experimental means or only with a single component as the object of study. There is still a lack of in-depth exploration of the numerical modeling and analysis of the dynamics of the complex coupled structure of the fiber woven materials; in the case of the open-pore structure, the object of study is generally a simple plate and shell structure. There is no focus on the application of simple structures to complex systems. On the whole, the existing studies seldom consider fiber woven composites, weight-reducing holes and shaft–disk rotors in an integrated way, and the modeling method of the fiber woven composite shaft–disk rotor system with weight-reducing holes has yet to be perfected.
In view of the above research deficiencies, this paper takes the fiber woven composite material shaft–disk rotor as the research object, adopts the domain decomposition and coordinate mapping technology for the porous disk, and establishes a unified dynamic model of the fiber woven composite shaft–disk system through simple beam and plate theory and the differential quadrature finite element method (DQFEM). The model has been verified by the literature, finite element simulation, and experimental comparison, and has been proven to have good convergence, accuracy, and practicability. The corresponding kinetic study was carried out to investigate the effects of material parameters and hole parameters on the rotor’s natural frequency and critical speed, and the conclusions obtained will be helpful to guide the lightweight and vibration-damping design of the shaft–disk rotor.

2. Dynamic Modeling Process of the Fiber Woven Composite Shaft–Disk Rotor

As shown in Figure 1a, the fiber woven composite shaft–disk rotor consists of a perforated ring plate and several shafts with different radii. The rotor system is placed in the absolute coordinate system C0, and a rotational coordinate system C1 is used to describe the rotational motion of the rotor, which rotates uniformly about the X1 axis with an angular velocity of Ω. A body-fixed coordinate system C2 is used to describe the displacement deformation of the perforated ring plate. The parameters of the rotor model geometry are shown in Figure 1b.

2.1. Application of Fiber Woven Composite

In this paper, fiber woven composites are applied to gear webs. Fiber woven composites are composites with fiber braids as reinforcements and resin, metal, ceramics, etc., as matrix materials, and the two are cured and molded by a curing process. The weaving process can be divided into plain, twill, satin, and more complex three-dimensional woven composites. The following study selects orthogonal plain-weave carbon fiber-reinforced resin matrix composites as the ring plate material, as shown in Figure 2. Meanwhile, the focus of the study is not on solving the mechanical parameters of the material, but on exploring the effect of its use in the rotor model. Therefore, reference [7], the Young’s modulus, shear modulus, Poisson’s ratio, and density for different total fiber volume fractions (Vf) are known quantities and are given in Figure 3.

2.2. Domain Decomposition and Coordinate Mapping Technique for Perforated Webs

In order to reduce weight, facilitate machining and handling, and enable the system’s critical speed to avoid the operating speed, holes are often cut in the rotor web. These are commonly round holes, sector holes, etc. The web structure with holes, unlike a simple ring plate, is not easy to cut in the circumferential and radial directions for dividing regular sector cells when modeling. Therefore, a domain decomposition method for the web with arbitrary complex hole patterns is given below, and a coordinate mapping technique is introduced for subsequent dynamic modeling of the rotor.
As an example, a ring plate with N round holes with inner/outer diameters of Ri/Ro can be uniformly decomposed into N sector plates with holes in the center, as shown in Figure 4. Since the modeling of the ring plate is usually performed in a rotating coordinate system, the sector plate is also divided into sub-domains in the r-θ coordinate system to ensure that the segmentation lines are easy-to-describe straight lines in the r-θ coordinate system. As in Figure 4, the ith sector plate is decomposed into four irregular quadrilateral domains. It should be noted that the domain decomposition should avoid the tangent point of the hole-shaped boundary and the parallel line of the r/θ coordinate axis; otherwise, it will affect the accuracy of the subsequent coordinate mapping, which will lead to the computational results of the model with the problems of non-convergence and large error.
After the subdomain decomposition of the open ring plate, each subdomain can be directly regarded as a plate cell for solving the energy matrix of the cell based on the DQFEM. However, the unit is an irregular quadrilateral, which needs to be converted to a regular rectangular unit by means of the coordinate mapping technique. As shown in Figure 5, for any irregular quadrilateral cell in the r-θ coordinate system, assuming that it has p × q mapping nodes, the nodes here are only used to reflect the mapping relationship, which is different from the nodes used in the differential product in Section 2.3.
By mapping Equation (1) [37], irregular quadrilaterals in the r-θ coordinate system can be converted into square cells with −1 ≤ ξ ≤ 1 and −1 ≤ η ≤ 1 in the ξ-η coordinate system.
r = R q p ψ q p , θ = Θ q p ψ q p , R q p = r 11 , r 12 , , r i j , , r q p , Θ q p = θ 11 , θ 12 , , θ i j , , θ q p , ψ q p = φ 11 , φ 12 , , φ i j , , φ q p T
(rij, θij) denotes the coordinates of the mapping points distributed on the edges or inside of the quadrilateral cell. φij denotes the mapping function corresponding to the ij mapping point (rij, θij), and the specific expression [15] is
φ i j ξ , η = s = 1 , s i p ξ ξ s ξ i j ξ s l = 1 , l j q η η l η i j η l
After the ring plate irregular unit has been coordinate-mapped, the two-dimensional differential product weight coefficient matrix for modeling using the DQFEM needs to be changed accordingly. The existing coordinates (r, θ) are functions with respect to (ξ, η), and the displacement function f of the structure is a composite function of (ξ, η). The basic two-dimensional differential product weight coefficient matrices A02, B02 and the integral weight coefficient matrix C02 are derived in the ξ-η coordinate system. And for the irregular cells in the r-θ coordinate system, the weight coefficient matrices need to be corrected by the partial derivatives of r and θ with respect to ξ and η. The correction formula is derived as follows. The partial derivatives of the displacement function f with respect to ξ and η can be expressed as:
f / ξ f / η = r ξ θ ξ r η θ η f / r f / θ = J f / r f / θ
where J is the Jacobi matrix, and the partial derivatives of the displacement function f with respect to r and θ can be obtained from the inverse of J as follows:
f / r f / θ = J 1 f / ξ f / η = θ η θ ξ r η r ξ / f / ξ f / η r ξ θ η θ ξ r η
f r = 1 J θ η f ξ θ ξ f η , f θ = 1 J r η f ξ + r ξ f η
According to Equation (1), the partial derivative operator in Equations (4) and (5) can be obtained as:
r ξ = i j = 11 p × q r i j φ i j ξ , η ξ , r η = i j = 11 p × q r i j φ i j ξ , η η , θ ξ = i j = 11 p × q θ i j φ i j ξ , η ξ , θ η = i j = 11 p × q θ i j φ i j ξ , η η
The original weighting coefficient matrix expressions are:
f ξ , η ξ i j = m = 1 M A i m f m j = A 02 1 f , f ξ , η η i j = n = 1 N B j n f i n = B 02 1 f
Here and in the following equation, i and j are the number of radial and circumferential nodes in the quadrilateral cell divided based on DQFEM, respectively. i = 1, 2, …, M, and j = 1, 2, …, N. The actual weight coefficients’ matrix expression is corrected to:
f r i j = 1 J i j θ η i j A 02 θ ξ i j B 02 , f θ i j = 1 J i j r η i j A 02 + r ξ i j B 02
For ease of presentation, the above equation can be rewritten as:
A 2 = 1 Y . θ η . A 02 θ ξ . B 02 , B 2 = 1 Y . r η . A 02 + r ξ . B 02 , a α = a α 11 a α M 1 a α 12 a α M N T , a = r , θ ; α = ξ , η , Y = J 11 J M 1 J 12 J M N T
where the symbol aα denotes the partial derivative of a with respect to α and “.*” denotes the Hadamard product. Similarly, the modified matrix of integral weight coefficients is:
C 2 = Y . C 02
In order to facilitate the parametric study, four types of holes are given: (1) no holes; (2) round holes; (3) sector holes; and (4) curved holes. The shape and geometry of the various hole types are shown in Figure 6, and the geometric parameters of the different hole types should be set to ensure that the area of a single hole is the same in the subsequent parameterization study.

2.3. Energy Derivation of System Structure Units

As shown in Figure 7, the nodes on the shaft in the rotor have six degrees of freedom (uS, vS, wS, φS,x, φS,y, φS,z), and the nodes on the disk have five degrees of freedom (uD, vD, wD, φD,r, φD,θ).

2.3.1. Potential and Kinetic Energy of the Rotating Shaft Units

Based on the Timoshenko beam theory, the shaft in the rotor can be equated to a six-degree-of-freedom beam, and the elastic deformation displacement at a point (xS, yS, zS) on the beam is:
w ¯ S = w S x S y S φ S , z x S + z S φ S , y x S , u ¯ S = u S x S z S φ S , x x S , v ¯ S = v S x S + y S φ S , x x S
where yS = zS = 0. Based on elastic mechanics, the strain at this point on the axis is deduced to be:
ε S , x = w S x S y S φ S , z x S + z S φ S , y x S , γ S , x y = u S x S z S φ S , x x S φ S , z , γ S , z x = v S x S + y S φ S , x x S + φ S , y
In turn, the potential energy of the shaft is:
V S = 1 2 E S ε S , x 2 + G S γ S , x y 2 + G S γ S , z x 2 d y S d z S d x S = 1 2 0 L E S A S w S x S 2 + E S I S φ S , z x S 2 + E S I S φ S , y x S 2 + κ S G S A S u S x S φ S , z 2 + κ S G S A S v S x S + φ S , y 2 + 2 G S I S φ S , x x S 2 d x S
where ES, κS (= 6(1 + μS)2/(7 + 12μS + 4μS2)) are the modulus of elasticity and shear modification factor of the shaft, respectively, and μS is Poisson’s ratio. Its cross-sectional area AS = πr2, moment of inertia IS = πr4/4, and shear modulus GS = ES/2/(1 + μS). L and r are the length and radius of the shaft section, respectively.
Referring to my existing research [38], the kinetic energy of the shaft is:
T S = 1 2 ρ S 0 L S A S w ˙ S 2 + u ˙ S 2 + v ˙ S 2 + I S 2 φ ˙ S , x 2 + φ ˙ S , y 2 + φ ˙ S , z 2 + 2 Ω A S u S v ˙ S u ˙ S v S 2 Ω I S φ ˙ S , y φ S , z + φ ˙ S , z φ S , y + Ω 2 A S u S 2 + v S 2 Ω 2 I S φ S , y 2 + φ S , z 2 d x S
Ω denotes the angular velocity of the rotor, and ρS denotes the density of the shaft.

2.3.2. Potential and Kinetic Energy of the Disk Units

According to the first-order shear deformation theory (FSDT), the displacement of a point (xD, rD, θD) on the disk is expressed as:
u ¯ D ( r D , θ D , x D , t ) = u D ( r D , θ D , t ) + x D φ D , θ ( r D , θ D , t ) , v ¯ D ( r D , θ D , x D , t ) = v D ( r D , θ D , t ) x D φ D , r ( r D , θ D , t ) , w ¯ D ( r D , θ D , x D , t ) = w D ( r D , θ D , t )
The corresponding strain is:
ε D , r = ε D , r 0 + x D χ D , r 0 = u D r D + x D φ D , θ r D , ε D , θ = ε D , θ 0 + x D χ D , θ 0 = 1 r D v D θ D + u D r D + x D 1 r D φ D , r θ D + φ D , θ r D , γ D , r θ = γ D , r θ 0 + x D χ D , r θ 0 = 1 r D u D θ D + v D r D v D r D + x D 1 r D φ D , θ θ D φ D , r r D + φ D , r r D , γ D , r x = w D r D + φ D , θ , γ D , θ x = 1 r D w D θ D φ D , r
According to Hooke’s law, the stresses of orthogonal plain-weave carbon fiber-reinforced resin matrix composite disks are:
σ D , r σ D , θ τ D , r θ τ D , r x τ D , θ x = Q 11 Q 12 0 0 0 Q 21 Q 22 0 0 0 0 0 Q 66 0 0 0 0 0 Q 44 0 0 0 0 0 Q 55 ε D , r ε D , θ γ D , r θ γ D , r x γ D , θ x = E 11 1 μ 12 μ 21 μ 21 E 11 1 μ 12 μ 21 0 0 0 μ 12 E 22 1 μ 12 μ 21 E 22 1 μ 12 μ 21 0 0 0 0 0 G 12 0 0 0 0 0 G 23 0 0 0 0 0 G 13 ε D , r ε D , θ γ D , r θ γ D , r x γ D , θ x
Since the first-order shear deformation theory is applied to the ring plate, only E11, E22, μ12, and μ21 are introduced for orthotropic anisotropic materials, and they satisfy E11μ21 = E22μ12, i.e., Q12 = Q21. In summary, the potential energy of the ring plate can be found to be:
V D = 1 2 σ D , r ε D , r + σ D , θ ε D , θ + τ D , r θ γ D , r θ + κ D τ D , r x γ D , r x + κ D τ D , θ x γ D , θ x r D d r D d θ D d x D = 1 2 Q 11 r D u D r D + x D φ D , θ r D 2 + Q 22 r D v D θ D + u D x D φ D , r θ D φ D , θ 2 + 2 Q 12 u D r D + x D φ D , θ r D v D θ D + u D x D φ D , r θ D φ D , θ + Q 66 r D 1 r D u D θ D + v D r D v D r D + x D φ D , θ θ D r D φ D , r r D + φ D , r 2 + κ D Q 44 r D w D r D + φ D , θ 2 + κ D Q 55 r D 1 r D w D θ D φ D , r 2 d r D d θ D d x D
where the shear correction factor κD = 5/6.
For a multi-component coupled geared rotor, a multi-body dynamics approach is used to calculate the kinetic energy of the disk in the system. In the absolute coordinate system C0, the position vector of a point (xD, rD, θD) on the disk is shown in Equation (19).
r D = R 0 + R D + U D A 1 A Ω
where R0 (= [X, Y, Z]) is the origin’s position of the rotating coordinate system C1 in the absolute coordinate system C0. RD denotes the position vector of the coupling nodes of the shaft and the disk, whose value is RD = [wSC + xSC, uSC, vSC], and the subscript “SC” denotes the coupling nodes of the shaft and the disk. A1 and AΩ are the coordinate transformation relation matrices from C2 to C1 and C1 to C0, and their expressions are detailed in the literature [38]. The symbol UD in Equation (19) denotes the position vector of a point (xD, rD, θD) on the disk in the follower coordinate system with the following expression:
U D = R DP + U DP = x D r D cos θ D r D sin θ D T + w D u D + x D φ D , θ cos θ D v D x D φ D , r sin θ D u D + x D φ D , θ sin θ D + v D x D φ D , r cos θ D T
Combining the above equations, the velocity of the node on the disk can be obtained as:
r ˙ D = R D + U D A 1 A ˙ Ω + R ˙ D + U ˙ D A 1 + U D A ˙ 1 A Ω
This, in turn, leads to an expression for the kinetic energy of the disk as:
T D = 1 2 ρ D r ˙ D r ˙ D T r D d r D d θ D d x D
where ρ denotes the density. The kinetic energy formula for the disk is too lengthy to show specifically.

2.4. Overall Model of the Rotor System

For the overall structure, the energy expressions (Equations (13), (14), (18) and (22)) for each structural element are converted to matrix form based on the DQFEM:
V j e = 1 2 q j e T K j e q j e , j = S , D q S e = w S u S v S φ S , x φ S , y φ S , z T ; q D e = w D u D v D φ D , r φ D , θ T
T j e = 1 2 q ˙ j e T M j e q ˙ j e + 1 2 q ˙ j e T G j e q j e + 1 2 q j e T C j e q j e , j = S , D j = S , q j e = q S e T ; j = D , q j e = q SC e q D e T ; q SC e = w SC w SC v SC φ SC , x φ SC , y φ SC , z T
where Ke, Me, Ge, and Ce denote the cell stiffness matrix, mass matrix, gyro matrix, and damping matrix, respectively.
The coupling relationship between the shaft segments and the disk is included in the above kinetic energy expression. For the whole geared rotor system, it is also necessary to add the coupling potentials and boundary potentials of the individual shaft segments, shafts, and disks to the total potential energy, and the connections between the substructures and the classical boundary conditions are simulated using the artificial spring technique. For the connection between n shaft segments, three sets of translational springs and three sets of rotational springs are defined to simulate the connection, with spring stiffnesses kSS,u, kSS,v, kSS,w, kSS,x, kSS,y, and kSS,z (all of which have values of 1 × 1013). The corresponding potential energy VSS is:
V SS = 1 2 i = 1 n 1 k SS , m m S , r i + 1 m S , l i 2 , ( m = w , u , v , φ x , φ y , φ z )
where the subscripts ri+1 and li denote the rightmost node of the (I + 1)st shaft segment and the leftmost node of the i-th shaft segment, respectively. n denotes the total number of shaft segments. For the simulation of the coupling point on the shaft connected to the nodes of the inner ring plate, three sets of translational springs and two sets of rotational springs with spring stiffnesses kS,u, kS,v, kS,w, kS,x, kS,φr, and kS,φθ (all of which have values of 1 × 1013) are defined. The corresponding potential energy VC is:
V C = 1 2 0 2 π k SD , m m D , a 2 | a = R i r D d θ D , ( m = w , u , v , φ r , φ θ )
For the simulation of the boundary conditions, three sets of translational springs and three sets of rotational springs are defined with spring stiffnesses kS,u, kS,v, kS,w, kS,x, kS,y, and kS,z. The boundary potential energy is:
V B = 1 2 k S , m , a m S , a 2 + k S , m , b m S , b 2 | a = 0 , b = L , ( m = w , u , v , φ x , φ y , φ z )
We define three boundaries: free (F), simply supported (S), and solidly supported (C) with the following spring stiffness values:
F : k S , u = k S , v = k S , w = 0 , φ S , x = φ S , y = φ S , z = 0 S : k S , u = k S , v = k S , w = 1 × 10 13 , φ S , x = φ S , y = φ S , z = 0 C : k S , u = k S , v = k S , w = 1 × 10 13 , φ S , x = φ S , y = φ S , z = 1 × 10 13
In summary, the total potential and kinetic energies of the system are:
V = V S + V D + V SS + V C + V B
T = T S + T D
Based on the DQFEM, the differential equations of motion for a single gear rotor system can be derived by substituting the matrix form of the total potential and kinetic energies of the system into the Lagrangian equation:
M T + M q ¨ + G T G q ˙ + K T + K C T + C q = 0 , q = q S q D
where M, G, C, and K are the overall mass matrix, gyro matrix, damping matrix, and stiffness matrix assembled from the individual unit matrices, boundary potential matrix, and coupling potential matrix between the structures. By solving the eigenvalues and eigenvectors of Equation (31), the vibration characteristics of the system, such as the natural frequency and vibration mode, can be obtained.

3. Numerical Analysis

Based on the proposed DQFEM, this study utilizes MATLAB programming to calculate the vibration characteristics of a fiber woven composite shaft–disk rotor with weight-reducing holes. Initially, a convergence analysis of the proposed DQFEM model is conducted; subsequently, the computational results are compared with those obtained from the commercial finite element software ANSYS and experimental modal testing to verify the accuracy of the established model. Finally, a parametric study on the hole parameters of the disk is carried out.

3.1. Convergence Study

From Section 2.2, it can be seen that the disk has a definite number of units, so the convergence analysis mainly focuses on the number of units of the shaft section (taking the number of unit nodes of the shaft section as 3) and the number of unit nodes of the ring plate. The convergence analysis is carried out by taking the rotor made of fiber woven composite with round holes as an example, and the structure of the rotor is as follows: the length of the shaft L = 0.2 m, the radius r1 = r2 = 0.02 m; the thickness of the disk hD = 0.008 m, the position Ld = 0.1 m, the inner diameter Ri = 0.02 m, the outer diameter Ro = 0.1 m, the number of round holes N = 4, the radius of the round holes rh = 0.025 m, and the eccentricity of the round holes Rh = 0.06 m. The number of shaft segment units not used as variables is 5, and the number of ring plate unit nodes is 10. The boundaries of the rotor are set at the two ends of the rotor shaft, both of which are simply supported boundaries and are denoted as S–S.
Unless otherwise stated, the shaft in the following calculations is made of steel with Young’s modulus E = 210 GPa, Poisson’s ratio μ = 0.3, and density ρ = 7800 kg/m3, and the disk material is an orthotropic plain-woven carbon fiber-reinforced resin matrix composite when the total fiber volume fraction is equal to 0.28. The mechanical properties of the materials with different total fiber volume fractions are all given in Figure 3.
The first six orders of the natural frequencies of the fiber woven composite rotor are given in Table 1 and Table 2 for different numbers of shaft segment units and different numbers of ring plate nodes, respectively. From the tables, it can be found that the natural frequency of the rotor converges rapidly with the increase in the number of shaft segment units and the number of ring plate nodes, and stabilizes when the number of shaft units is greater than or equal to 3 and the number of ring plate nodes is greater than or equal to 10. Therefore, the number of shaft segment units is set to 3 and the number of ring plate nodes is set to 10 in the following examples.

3.2. Verification

3.2.1. Validation of Literature and FEM

The model of the shaft–disk system is validated, and Table 3 demonstrates the literature comparison results of the shaft–disk system with the S–S boundary, whose structural parameters are L = 0.35 m, Ld = 0.37 L, r = Ri = 0.01 m, Ro = 0.1 m, and hd = 0.002 m and whose material parameters are E = 200 GPa, μ = 0.266, and ρ = 7860 kg/m3. The comparative values of the intrinsic frequencies are very similar, which indicates that the constructed shaft–ring plate model is equally reliable.
Table 4 compares the coupled natural frequencies and modes of the shaft–disk rotor of the fiber woven material, and the coupled natural frequencies of the steel rotor are also given in the table. The parameter settings of the disk are the same as those in Table 4, except that Ld = 0.09 m and hd = 0.006 m. The shaft parameters L1 = 0.08 m, L2 = 0.02 m, L3 = 0.04 m, and L4 = 0.1 m correspond to the radii r1 = r4 = 0.01 m, r2 = 0.02 m, and r3 = 0.015 m. It can be seen that: the 1st- and 4th-order coupled modes are dominated by the displacement deformation of the disk, and the deformation of each structure is generated in the 2nd- and 3rd-order modes, among which the displacements of the disk in the 2nd order are mainly determined by the coupling points on the shaft. The frequency deviation does not exceed 2.12%, while the modes are basically consistent, which directly reflects the accuracy of the present model. In addition, comparing the steel rotor, it is found that the use of fiber woven material leads to a reduction in the rotor’s 1st coupling frequency, and at the same time, it can reduce the mass of the ring plate by 82.6%, which helps to achieve the goal of light weight. Considering that there is almost no vibration displacement of the shaft in the 1st-order coupled mode of the system, its coupling characteristics are not obvious, so the subsequent parametric study is dominated by the 2nd- and 3rd-order coupled intrinsic frequencies.
Finally, the results of the Campbell diagram comparison of the 2nd- and 3rd-order coupled natural frequencies of the rotor system with fiber woven material are given in Figure 8. The curves are in excellent agreement. Taken together, it seems that the present model is equally applicable to rotating rotor systems.

3.2.2. Experimental Verification

In this section, modal experiments are carried out for the disk as well as the shaft–disk rotor, four structures (round-hole disk, sector-hole disk, shaft–disk (round-hole) rotor, and shaft–disk (sector-hole) rotor) are selected as experimental objects, and the experimental pieces are all prepared from 45 steel. The parameters of the disk in the shaft–disk rotor are the same as those of the corresponding single plate, as follows: Ri = 0.02 m, Ro = 0.2 m, hd = 0.008 m, and N = 4, of which the round hole Rh = 0.11 m, rh = 0.05 m, the sector hole R1 = 0.07 m, R2 = 0.15 m, and dh = 0.02 m. The rotors’ shafts are the same, and the disks are located in the middle of the 2nd section of the shaft. The geometric parameters of shaft are: L1 = 0.3 m, L2 = 0.1 m, L3 = 0.2 m, r1 = 0.015 m, r2 = 0.02 m, and r3 = 0.015 m.
In the experiment, we used the hammering method to test the natural frequency and modal vibration pattern of the structure, and the required equipment came from the DongHua DHDAS dynamic signal test and analysis system. The experimental schematic and field diagrams are shown in Figure 9 and Figure 10.
The overall experimental process is divided into three major parts:
(1)
Experimental pre-processing. The experimental parts are processed, and according to the size of the experimental parts in the modal analysis software, the model is drawn and the measurement point numbering is completed. Then, the corresponding measurement points are drawn. The experimental parts are set as the free boundary, and the elastic rope suspension is used to simulate the boundary conditions. The simulation of the experimental parts is completed and the nodes with large deformation in each order of the modal state are selected as the response test points in order to paste the sensor during the experiment.
(2)
Experimental measurement. The experiment of disk structure adopts the method of single-point vibration pickup; the experiment of shaft–disk structure adopts the method of single-point excitation and multi-point response, so the modal test of the disk only needs to place a single acceleration sensor on the response test point, while the shaft–disk system needs to select a measurement point on the shaft and the disk and paste acceleration sensors. Afterwards, the device connection is completed and the device parameters are debugged. Finally, the force hammer is used to tap each measurement point in turn to obtain and store the response data.
(3)
Post-experiment processing. After completing the experiment, the obtained data are processed through the modal analysis software, the intrinsic frequency is calculated, and the modal vibration pattern is plotted.
Table 5 give the experimental results and their comparisons with the model, and Difference 1 and Difference 2 in the tables represent the deviations of this method from the experimental and finite element results, respectively. The natural frequencies of the two experimental parts are in good agreement, the comparison deviation of the shaft–disk rotor is not more than 2.5%. Considering the processing error of the experimental parts, the measurement error, and the imperfect simulation of the boundary conditions, in addition to many other factors, it can be assumed that the above deviations are in the acceptable range. Taking into account that the number of measurement points is much smaller than the number of nodes in the model, the comparison of the experimental modal results is not good, but in general, the two modal patterns are basically the same. The above results further illustrate the correctness of the modeling methodology of open-plate and rotor dynamics in this chapter.

3.3. Parametric Studies

Based on the constructed rotor model, the vibration characteristics of a rotating fiber woven composite rotor with a simply supported boundary are investigated. Unless otherwise stated, the system geometrical parameters and boundary settings in the following examples are the same as those in Table 4.
Firstly, the effect of the hole types on the rotor’s Campbell diagram is explored. In order to avoid an influence of the hole’s area on the structure, the number of holes N is set as 4 uniformly, and at the same time, the shape is set to make the area of the single holes as similar as possible, with the specific dimensions as follows: round holes: Rh = 0.06 m, rh = 0.025 m; sector-shaped holes: R1 = 0.045 m, R2 = 0.075 m, dh = 0.015 m; curved holes: R1 = 0.043 m, R2 = 0.077 m, φ = π/6. The Campbell diagrams of the 2nd and 3rd coupled natural frequencies of the system are shown in Figure 11. It is clear that the hole type has an effect on the 2nd forward (solid line) and backward (dashed line) coupled frequency values of the system, especially for the latter. This is due to the fact that this order mode is dominated by the shaft and there is a huge difference between the modes of the disk with and without holes at the same time. The 3rd natural frequency is affected by the vibration of the shaft and the disk at the same time, and the overall frequency shows the law of “No hole > Round hole > Curved hole > Sector hole”. The corresponding 3rd critical speed has the same trend.
Figure 12 exemplifies the effect of the total fiber volume fraction Vf on the vibration characteristics of the system, with the system pore parameters set as above. In conjunction with Figure 3, the Young’s modulus, shear modulus, and density of the fiber woven material itself continue to increase as Vf grows from 0.28 to 0.34, with a small decrease in Poisson’s ratio. The joint change of these parameters makes the overall stiffness of the rotor systems with different open hole types increase, so that the 3rd coupling frequency shows a similar linear reduction law: the growth rate of the 3rd coupling forward frequency gradually slows down, and the deceleration of the backward frequency gradually grows. The comparison between the subplots reveals that rotors with different hole types have different coupling frequency reduction rates, among which the reduction rate of the rotor with sector holes is the largest.
Taking the circular hole pattern with Rh = 0.06 m and rh = 0.02 m as an example, Figure 13, Figure 14 and Figure 15 demonstrate the effects of the radius rh, the number of holes N, and the hole eccentricity Rh on the Campbell diagrams of the 2nd and 3rd coupling frequencies. Figure 13 shows the effect of the round hole radius rh on the natural frequency of the system when N = 4 and Rh = 0.06 m. The change in rh can also be regarded as the change in the hole area, which is reflected as the change in stiffness in the overall structure. From the figure, it is found that when rh is raised from 0 to 0.03 m, the 2nd frequency increases gradually and the 3rd decreases continuously. The backward frequency deceleration of the former mutates to 0 at different rotational speeds, and the mutation rotational speed increases with the increase in rh. The mutation phenomenon is related to the modal change, referring to Figure 15 and Figure 16, for the system with small rh. The 1st and 2nd modal diagrams are similar, and during the increase in the rotational speed, these two natural frequencies are shifted, which leads to the mutation of the frequency.
Figure 14 displays the Campbell diagrams for fiber woven composite rotor systems with different hole numbers N. From Figure 14a, it is found that the 2nd natural frequency of the system improves with the increase in N as a whole, but then the decreasing amplitude of the backward frequency produces abrupt changes at different rotational speeds, while the smaller N is the lower the speed at which the abrupt change occurs. Figure 14b shows that the 3rd order coupling frequency of the system grows with increasing N, and the decrement amplitude increases.
Figure 15 and Figure 17 present the Campbell plots and modes of the coupling frequency of the rotor of the fiber woven composite for different hole eccentricities Rh at N = 4 and rh = 0.015 m. With the increase in Rh, the 2nd natural frequency of the system increases, and the 3rd natural frequency changes are not obvious, while the regularity is not strong. Taken together, the effect of the hole type and hole radius on the coupling frequency of the fiber woven material rotor system is more significant.

4. Conclusions

In this paper, the application of a fiber woven composite material to the disk of a shaft–disk rotor is proposed, and the domain decomposition method and coordinate mapping technique are introduced for the porous disk. Then, a dynamic modeling method applicable to the rotor with any complex hole types is given. Based on Timoshenko beam theory and FSDT, the energy expression of the shaft–disk rotor model is derived, and a rotor dynamics model of fiber woven material applicable to arbitrary hole shapes is established by DQFEM. The constructed model is verified to have good convergence, accuracy, and practicability after comparative tests in the literature, finite element simulation, and an experiment. A subsequent parametric study based on the proposed model was carried out, and the following conclusions were drawn:
(1)
The existing fiber woven composite rotor model has good convergence and accuracy, while the disk’s mass is reduced by 82.6% and the first-order coupled natural frequency is slightly reduced compared to the steel rotor.
(2)
The 1st and 4th coupled modes of the rotor are dominated by the displacement of the disk, while the 2nd- and 3rd-order modes contain the deformation of each structure, in which the disk displacement in the 2nd mode is mainly determined by the displacement of the coupling point on the shaft.
(3)
The hole type on the disk directly affect the 3rd critical speed of the rotor, which is ranked in descending order as “no hole > round hole > curved hole > sector hole”.
(4)
Taking the rotor with circular holes as an example, the increase in the total fiber volume fraction Vf, the hole eccentricity Rh and the hole number N, as well as the decrease in the hole radius rh, all enhance the 3rd coupled natural frequency and the critical speed of the fiver woven composite rotor system. Among them, the effect of the hole radius is the most significant.

Author Contributions

Conceptualization, T.L. and Q.W.; methodology, H.Z.; software, M.S.; validation, H.Z., M.S. and Z.L.; formal analysis, H.Z.; investigation, H.Z.; resources, H.Z.; data curation, Z.L.; writing—original draft preparation, H.Z.; writing—review and editing, T.L.; visualization, M.S.; supervision, Q.W.; funding acquisition, Q.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China grant number [52075554], the Natural Science Foundation of Hunan Province of China grant number [2021JJ30841 and 2022JJ20070], Central South University Innovation-Driven Research Program, China grant number [2023CXQD049].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

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

Acknowledgments

The authors would like to thank the anonymous reviewers for their very valuable comments. The authors also gratefully acknowledge the financial support from the National Natural Science Foundation of China (Grant No. 52075554), the Natural Science Foundation of Hunan Province of China (2021JJ30841 and 2022JJ20070), and Central South University Innovation-Driven Research Program, China (Grant No. 2023CXQD049).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, Y.; Li, H.; Liu, X.; Chen, Y.; Qin, C.; Fang, D. Prediction of mechanical properties of 3D tubular braided composites at different temperatures using a multi-scale modeling framework based on micro-CT. Compos. Sci. Technol. 2024, 245, 110349. [Google Scholar] [CrossRef]
  2. Huang, C.; Norton, M.; Joosten, M.W. Investigating the tensile properties of 3D printed semi-woven intraply hybrid composites reinforced with carbon/glass continuous fibres. Compos. Part A Appl. Sci. Manuf. 2024, 177, 107924. [Google Scholar] [CrossRef]
  3. Kumar, Y.; Rezasefat, M.; Amico, S.C.; Manes, A.; Dolez, P.I.; Hogan, J.D. Comparison of two progressive damage models for predicting low-velocity impact behavior of woven composites. Thin-Walled Struct. 2024, 197, 111611. [Google Scholar] [CrossRef]
  4. Grotz, L.; Kirane, K. Characterizing compressive failure mechanisms and their transitions in woven composites under on and off-axis loading. Compos. Struct. 2024, 330, 117848. [Google Scholar] [CrossRef]
  5. Poojary, U.R.; Hegde, S. Experimental Investigation of Influence of Fibre Orientation on the Dynamic Properties of Carbon Fibre and Intra-Ply Woven Carbon-Kevlar/Epoxy Hybrid Composite. J. Compos. Sci. 2024, 8, 78. [Google Scholar] [CrossRef]
  6. Zorko, D.; Tavčar, J.; Bizjak, M.; Šturm, R.; Bergant, Z. High cycle fatigue behaviour of autoclave-cured woven carbon fibre-reinforced polymer composite gears. Polym. Test. 2021, 102, 107339. [Google Scholar] [CrossRef]
  7. Huang, Q.; Gao, Y.; Hua, F.; Fu, W.; You, Q.; Gao, J.; Zhou, X. Free vibration analysis of carbon-fiber plain woven reinforced composite conical-cylindrical shell under thermal environment with general boundary conditions. Compos. Struct. 2023, 322, 117340. [Google Scholar] [CrossRef]
  8. Reddy, A.B.; Ram, K.S.S. Buckling of functionally graded carbon nanotube reinforced composite cylindrical shell panel with a cutout under uniaxial compression. Mater. Today Proc. 2022, 49, 1865–1869. [Google Scholar] [CrossRef]
  9. Wang, G.; Li, W.; Feng, Z.; Ni, J. A unified approach for predicting the free vibration of an elastically restrained plate with arbitrary holes. Int. J. Mech. Sci. 2019, 159, 267–277. [Google Scholar] [CrossRef]
  10. Wang, G.; Li, W.; Feng, Z.; Ni, J. Sound Transmission between Cavities Via an Elastically Mounted Panel with Apertures, 3rd ed.; Institute of Noise Control Engineering: Wakefield, MA, USA, 2019; Volume 259, pp. 6023–6028. [Google Scholar]
  11. Kwak, S.; Kim, K.; Yun, J.; Kim, S.; Ri, P. Free Vibration Analysis of Laminated Closed Conical, Cylindrical Shells and Annular Plates with a Hole Using a Meshfree Method. Structures 2021, 34, 3070–3086. [Google Scholar] [CrossRef]
  12. Song, Y. Free vibration of arbitrarily shaped plates with complex cutouts. Thin-Walled Struct. 2023, 190, 110979. [Google Scholar] [CrossRef]
  13. Sun, X.; Zhang, P.; Qiao, H.; Lin, K. High-order free vibration analysis of elastic plates with multiple cutouts. Arch. Appl. Mech. 2021, 91, 1837–1858. [Google Scholar] [CrossRef]
  14. Zhong, R.; Wang, Q.; Hu, S.; Qin, B.; Shuai, C. Spectral element modeling and experimental investigations on vibration behaviors of imperfect plate considering irregular hole and curved crack. J. Sound Vib. 2022, 529, 116924. [Google Scholar] [CrossRef]
  15. Zhong, R.; Wang, Q.; Huang, Z.; Chen, L.; Shao, W.; Shuai, C. In-plane dynamic analysis of complex-shaped laminated cracked plates with irregular holes. AIAA J. 2023, 61, 3172–3189. [Google Scholar] [CrossRef]
  16. Cui, Y.S.; Wang, Y.Q. Free vibrations of axially loaded thin-walled shaft-disk rotors subjected to non-uniform temperature field. Thin-Walled Struct. 2024, 196, 111461. [Google Scholar] [CrossRef]
  17. Phadatare, H.P.; Pratiher, B. Dynamic stability and bifurcation phenomena of an axially loaded flexible shaft-disk system supported by flexible bearing. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2020, 234, 2951–2967. [Google Scholar] [CrossRef]
  18. Ri, K.; Han, W.; Pak, C.; Kim, K.; Yun, C. Nonlinear forced vibration analysis of the composite shaft-disk system combined the reduced-order model with the IHB method. Nonlinear Dyn. 2021, 104, 3347–3364. [Google Scholar] [CrossRef]
  19. Zhao, T.Y.; Wang, Y.X.; Pan, H.G.; Yuan, H.Q.; Cai, Y. Nonlinear forced vibration analysis of spinning shaft-disk assemblies under sliding bearing supports. Math. Methods Appl. Sci. 2021, 44, 12283–12301. [Google Scholar] [CrossRef]
  20. Heydari, H.; Khorram, A. Effects of location and aspect ratio of a flexible disk on natural frequencies and critical speeds of a rotating shaft-disk system. Int. J. Mech. Sci. 2019, 152, 596–612. [Google Scholar] [CrossRef]
  21. Tuzzi, G.; Schwingshackl, C.W.; Green, J.S. Study of coupling between shaft bending and disc zero nodal diameter modes in a flexible shaft-disc assembly. J. Sound Vib. 2020, 479, 115362. [Google Scholar] [CrossRef]
  22. Zhao, S.; Zhang, X.; Zhang, S.; Safaei, B.; Qin, Z.; Chu, F. A unified modeling approach for rotating flexible shaft-disk systems with general boundary and coupling conditions. Int. J. Mech. Sci. 2022, 218, 107073. [Google Scholar] [CrossRef]
  23. She, H.; Li, C.; Tang, Q.; Wen, B. Influence mechanism of disk position and flexibility on natural frequencies and critical speeds of a shaft-disk-blade unit. J. Sound Vib. 2020, 469, 115156. [Google Scholar] [CrossRef]
  24. Li, C.; She, H.; Tang, Q.; Wen, B. The coupling vibration characteristics of a flexible shaft-disk-blades system with mistuned features. Appl. Math. Model. 2019, 67, 557–572. [Google Scholar] [CrossRef]
  25. Heydari, H.; Khorram, A.; Afzalipour, L. The influences of stagger and pretwist angles of blades on coupling vibration in shaft-disk-blade systems. J. Vib. Acoust. 2020, 142, 011007. [Google Scholar] [CrossRef]
  26. Yang, T.; Ma, H.; Qin, Z.; Guan, H.; Xiong, Q. Coupling vibration characteristics of the shaft-disk-drum rotor system with bolted joints. Mech. Syst. Signal Process. 2022, 169, 108747. [Google Scholar] [CrossRef]
  27. Zhao, S.; Zhang, L.; Zhu, R.; Han, Q.; Qin, Z.; Chu, F. Modeling approach for flexible shaft-disk-drum rotor systems with elastic connections and supports. Appl. Math. Model. 2022, 106, 402–425. [Google Scholar] [CrossRef]
  28. Li, Z.; Zhao, T.; Kou, H.; Zhang, H.; Yuan, H. Vibration characteristics of multi-stage blade–disk–shaft integrated structure with three-dimensional crack. J. Vib. Eng. Technol. 2021, 9, 597–611. [Google Scholar] [CrossRef]
  29. Yang, L.; Mao, Z.; Chen, X.; Yan, R.; Xie, J.; Hu, H. Dynamic coupling vibration of rotating shaft–disc–blade system—Modeling, mechanism analysis and numerical study. Mech. Mach. Theory 2022, 167, 104542. [Google Scholar] [CrossRef]
  30. Yang, L.-H.; Mao, Z.; Wu, S.-M.; Chen, X.-F.; Yan, R.-Q. Steady-state coupling vibration analysis of shaft–disk–blade system with blade crack. Nonlinear Dyn. 2021, 105, 61–98. [Google Scholar] [CrossRef]
  31. Venkatachalam, R.; Prabu, S.B. Vibration characteristics of orthotropic shaft–disk system with different constrained damping layers: Experimental and numerical study. Int. J. Adv. Manuf. Technol. 2013, 65, 601–610. [Google Scholar] [CrossRef]
  32. Venkatachalam, R.; Balasivanandha Prabu, S. Vibration and damping analysis of orthotropic sandwich shaft-disc system using finite element method. Int. J. Mech. Mater. Des. 2012, 8, 287–296. [Google Scholar] [CrossRef]
  33. Wang, R.; Wang, Q.; Guan, X.; Zhang, Y.; Shao, W. Coupled free vibration analysis of functionally graded shaft-disk system by differential quadrature finite element method. Eur. Phys. J. Plus 2021, 136, 1–27. [Google Scholar] [CrossRef]
  34. Wang, R.; Wang, Q.; Guan, X.; Shao, W. The coupling free vibration characteristics of a rotating functionally graded shaft-disk system in thermal field. Thin-Walled Struct. 2022, 176, 109278. [Google Scholar] [CrossRef]
  35. Zhao, T.Y.; Ma, Y.; Zhang, H.Y.; Pan, H.G.; Cai, Y. Free vibration analysis of a rotating graphene nanoplatelet reinforced pre-twist blade-disk assembly with a setting angle. Appl. Math. Model. 2021, 93, 578–596. [Google Scholar] [CrossRef]
  36. Cai, Y.; Liu, Z.-F.; Zhao, T.-Y.; Yang, J. Parameter interval uncertainty analysis of internal resonance of rotating porous shaft–disk–blade assemblies reinforced by graphene nanoplatelets. Materials 2021, 14, 5033. [Google Scholar] [CrossRef] [PubMed]
  37. Filiz, S.; Bediz, B.; Romero, L.A.; Ozdoganlar, O.B. Three dimensional dynamics of pretwisted beams: A spectral-Tchebychev solution. J. Sound Vib. 2014, 333, 2823–2839. [Google Scholar] [CrossRef]
  38. Shen, M.; Wang, Q.; Wang, R. Investigation on the vibration mechanisms of a rotating FG-GPLRC shaft-disk-shell combined system. Structures 2023, 56, 105049. [Google Scholar] [CrossRef]
Figure 1. Coordinate system representation and geometrical parameters of the fiber woven composite shaft–disk rotor.
Figure 1. Coordinate system representation and geometrical parameters of the fiber woven composite shaft–disk rotor.
Applsci 14 09148 g001
Figure 2. Composition and application of plain fiber woven composites [7].
Figure 2. Composition and application of plain fiber woven composites [7].
Applsci 14 09148 g002
Figure 3. Mechanical properties of orthogonal plain-weave carbon fiber-reinforced resin matrix composites [7].
Figure 3. Mechanical properties of orthogonal plain-weave carbon fiber-reinforced resin matrix composites [7].
Applsci 14 09148 g003
Figure 4. Subdomain decomposition of a perforated ring plate.
Figure 4. Subdomain decomposition of a perforated ring plate.
Applsci 14 09148 g004
Figure 5. Coordinate mapping.
Figure 5. Coordinate mapping.
Applsci 14 09148 g005
Figure 6. Shape and geometry of various hole types.
Figure 6. Shape and geometry of various hole types.
Applsci 14 09148 g006
Figure 7. Degree-of-freedom representation of the nodes on the fiber woven composite shaft–disk rotor.
Figure 7. Degree-of-freedom representation of the nodes on the fiber woven composite shaft–disk rotor.
Applsci 14 09148 g007
Figure 8. Campbell diagram of a fiber woven composite shaft–disk system.
Figure 8. Campbell diagram of a fiber woven composite shaft–disk system.
Applsci 14 09148 g008
Figure 9. Schematic diagram of modal test experiment.
Figure 9. Schematic diagram of modal test experiment.
Applsci 14 09148 g009
Figure 10. Modal test experiment of shaft–disk rotors.
Figure 10. Modal test experiment of shaft–disk rotors.
Applsci 14 09148 g010
Figure 11. Campbell diagrams for fiber woven composite rotor systems with different hole types.
Figure 11. Campbell diagrams for fiber woven composite rotor systems with different hole types.
Applsci 14 09148 g011
Figure 12. Campbell diagrams of the 3rd coupling frequency of the fiber woven composite rotor systems with different total fiber volume fractions Vf.
Figure 12. Campbell diagrams of the 3rd coupling frequency of the fiber woven composite rotor systems with different total fiber volume fractions Vf.
Applsci 14 09148 g012aApplsci 14 09148 g012b
Figure 13. Campbell diagrams of fiber woven composite rotor systems at different apertures rh.
Figure 13. Campbell diagrams of fiber woven composite rotor systems at different apertures rh.
Applsci 14 09148 g013
Figure 14. Campbell diagrams for fiber woven composite rotor systems with different hole numbers N.
Figure 14. Campbell diagrams for fiber woven composite rotor systems with different hole numbers N.
Applsci 14 09148 g014
Figure 15. Campbell diagrams for fiber woven composite rotor systems with different hole eccentricities Rh.
Figure 15. Campbell diagrams for fiber woven composite rotor systems with different hole eccentricities Rh.
Applsci 14 09148 g015
Figure 16. Modal shapes of fiber woven composite rotor systems at different apertures rh.
Figure 16. Modal shapes of fiber woven composite rotor systems at different apertures rh.
Applsci 14 09148 g016
Figure 17. Modal shapes of fiber woven composite rotor systems with different hole eccentricities Rh.
Figure 17. Modal shapes of fiber woven composite rotor systems with different hole eccentricities Rh.
Applsci 14 09148 g017
Table 1. Convergence of the natural frequency (Hz) for fiber woven composite rotors with respect to the number of shaft units.
Table 1. Convergence of the natural frequency (Hz) for fiber woven composite rotors with respect to the number of shaft units.
Natural Frequency OrderNumber of Shaft Units
1234567
1565.25565.05565.03565.03565.03565.03565.03
2646.35646.35646.35646.35646.35646.35646.35
3696.87696.87696.87696.87696.87696.87696.87
4793.65793.65793.65793.65793.65793.65793.65
51588.521588.521588.521588.521588.521588.521588.52
Table 2. Convergence of the natural frequency (Hz) for fiber woven composite rotors with respect to the number of nodes in the disk unit.
Table 2. Convergence of the natural frequency (Hz) for fiber woven composite rotors with respect to the number of nodes in the disk unit.
Natural Frequency OrderNumber of Nodes in the Disk Unit
46810121416
1572.86567.67565.23565.03565.01565.00565.01
2666.23649.99646.77646.35646.28646.28646.28
3811.67702.49697.43696.87696.78696.74696.73
4824.67794.71793.82793.65793.63793.62793.64
51601.521589.621588.521588.361588.321588.311601.52
Table 3. Comparative validation of the natural frequency (Hz) of the shaft–disk system.
Table 3. Comparative validation of the natural frequency (Hz) of the shaft–disk system.
MethodDimensionless Natural Frequency Order
123456
Present163.44231.141109.971323.372744.023692.69
Zhao et al. [22]164.17232.241125.861331.65--
Heydari et al. [20]164.27231.571124.411332.042796.713729.08
FEM163.51230.801107.201324.502727.203700.10
Table 4. Validation of coupled natural frequency (Hz) and mode shape of the fiber woven composite shaft–disk system.
Table 4. Validation of coupled natural frequency (Hz) and mode shape of the fiber woven composite shaft–disk system.
MethodOrder
1234
PresentApplsci 14 09148 i001
431.03605.502455.574604.54
FEMApplsci 14 09148 i002
435.42604.592449.224509.08
Difference−1.01%−0.14%0.26%2.12%
Present
(made of steel)
418.92526.992352.404344.89
Table 5. Comparative results of natural frequency (Hz) and modal experiments of shaft–disk rotors with different hole types.
Table 5. Comparative results of natural frequency (Hz) and modal experiments of shaft–disk rotors with different hole types.
Hole TypeMethodOrder
1234
Round holePresentApplsci 14 09148 i003
166.34229.09359.62537.02
ExperimentApplsci 14 09148 i004
170.46231.5359.33545.44
FEMApplsci 14 09148 i005
166.23228.87359.24536.53
Difference 1−2.42%−1.04%0.08%−1.54%
Difference 20.06%0.10%0.11%0.09%
Sector holePresentApplsci 14 09148 i006
151.61218.34357.90501.74
ExperimentApplsci 14 09148 i007
154.09213.77354.16493.34
FEMApplsci 14 09148 i008
151.12217.85357.74500.78
Difference 1−1.61%2.14%1.06%1.70%
Difference 20.32%0.23%0.05%0.19%
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

Zhang, H.; Shen, M.; Liu, T.; Li, Z.; Wang, Q. Dynamic Modeling and Vibration Characteristic Analysis of Fiber Woven Composite Shaft–Disk Rotor with Weight-Reducing Holes. Appl. Sci. 2024, 14, 9148. https://doi.org/10.3390/app14199148

AMA Style

Zhang H, Shen M, Liu T, Li Z, Wang Q. Dynamic Modeling and Vibration Characteristic Analysis of Fiber Woven Composite Shaft–Disk Rotor with Weight-Reducing Holes. Applied Sciences. 2024; 14(19):9148. https://doi.org/10.3390/app14199148

Chicago/Turabian Style

Zhang, Haibiao, Mengyu Shen, Tao Liu, Zhen Li, and Qingshan Wang. 2024. "Dynamic Modeling and Vibration Characteristic Analysis of Fiber Woven Composite Shaft–Disk Rotor with Weight-Reducing Holes" Applied Sciences 14, no. 19: 9148. https://doi.org/10.3390/app14199148

APA Style

Zhang, H., Shen, M., Liu, T., Li, Z., & Wang, Q. (2024). Dynamic Modeling and Vibration Characteristic Analysis of Fiber Woven Composite Shaft–Disk Rotor with Weight-Reducing Holes. Applied Sciences, 14(19), 9148. https://doi.org/10.3390/app14199148

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