• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Structural design and modal behaviors analysis of a new swept baffled inflatable wing

    2023-07-04 08:02:08NuoMaLiLiuFanminMengJunhuiMeng
    Defence Technology 2023年6期

    Nuo Ma, Li Liu, Fanmin Meng, Junhui Meng

    School of Aerospace Engineering, Beijing Institute of Technology, Beijing 100081, PR China

    Keywords:Inflatable wing Structural design Wet mode Flutter analysis

    ABSTRACT

    1. Introduction

    By maintaining a holding pattern in the air for a long time and rapidly attacking ability,the loitering munition can attack non-lineof-sight targets with exact time, attitude, and direction, which is considered to be rapidly changing the battlefield[1].It is difficult to balance the suppression and damage performance of loitering munitions due to the limited size, weight, and cost [2], because of the contradiction between the fixed configuration of loitering munitions and different tasks. Specifically, loitering munitions for reconnaissance and electronic warfare sacrifice speed for longer flight time, while loitering munitions for attack need higher dynamic performance to achieve accurate strikes against time sensitive targets[2].Therefore,realizing more diversified performances(especially speed) and functions with limited volume will be the main development direction of loitering munitions in the future.

    The inflatable wing is a flexible inflatable membrane structure applied to the aerodynamic surface of aircraft [3—5]. It is suitable for loitering munitions because of their advantages of lightweight,cabinet storage volume,and rapid deployment.For example,the I-2000 UAV achieved a 42.7% weight reduction by replacing rigid wings with inflatable wings [3]. Due to the high flexibility, the inflatable wing can achieve adaptive change in span length and even rapid discard by pneumatic actuation.Such structural features make the inflatable winged loitering munitions have the potential to achieve wider range of available speeds and multi-mission capacity in flight, as shown in Fig.1. Obviously, the performance of inflatable structures is faced with more severe challenges due to various mission conditions. Considering the limitation of the traditional scheme,the structure of the inflatable wing needs to be redesigned.

    Fig.1. Flight profile of warpable inflatable winged loitering munitions.

    Fig.2. Distorted section of a closed membrane structure with a flexible baffle restraint.

    Different from the traditional wing supported by the skeleton with adamant material, the inflatable wing is composed of inflatable beams formed by high strength membranes.In particular,the airfoil can be taken shape by proper arrangement and restraint of the inflatable beams’ layout [6]. These features result in a prominent difference in design procedures compared to the traditional rigid wing,which cannot be ignored in the application.To be compatible with the bubble effect of the inflatable membrane,the target airfoil of the inflatable wing is approximated by a series of inscribed circles [7—9]. Therefore, the designed inflatable airfoil will form a rippled surface, resulting in deterioration of the aerodynamic characteristics,and more baffles should be used to remedy the airfoil, especially at the trailing edge [10]. Thus, the center of gravity is moved toward the trailing edge, which will further deteriorate the flutter performance of the inflatable wing. The bearing capacity of the inflatable structure is directly related to its diameter and inflation pressure [10], and the stiffness distribution of the inflatable wing is jointly determined by the inflatable beams.Under the constraint of airfoil,the stiffness distribution of inflatable wing is most independent of the number and distribution of inflatable beams, which leads to the worse designability of aeroelastic performance. With the increasing of the flight speed, the aeroelastic performance of inflatable winged loitering munitions is put forward a higher requirement. Therefore, it is an important issue to improve the aeroelastic performance of the inflatable wing to ensure its stable flight.

    For different requirements, scholars have developed and redesigned the structure of inflatable wings, and certain results in corresponding indicators have been achieved [11—13]. However,most of the research on new configurations of the inflatable wing are aimed to improve the bearing capacity or reduce structure mass, while the aeroelastic performance is neglected. To improve the aeroelastic properties of the wing structure, the aeroelastic tailoring technology has been developed by scholars [14,15]. By designing the order, number, and fabricating direction of the composite material layers laid on the wing structure, a targeted design of the stiffness distribution and the desired response under aerodynamic loads can be realized [16—18]. As the basic bearing membrane structural element of the inflatable wing,the inflatable beams shown anisotropic mechanical properties,and the inflatable wing with better aeroelastic performance is excepted be realized by changing the baffles sweep angle of the inflatable beams.However,the designed stiffness distribution is still difficult to be implemented easily due to the complexity of geometry.

    The inflatable beams will remain cylindrical after inflation according to the bubble principle [9]. Under the assumption of the swept arrangement, the inflatable beams have continuously changing chordwise positions, and the uniform cylindrical beam will not be able to meet the requirement of continuous thickness variation. Therefore, the geometry realization remains the main difficulty in implementing the designable stiffness distribution design based on the swept arrangement of inflatable beams. To overcome the difficulties mentioned above and improve the aeroelastic performance of the inflatable wing,an inflatable wing with a new structural configuration, referred to as the swept baffled inflatable wing is proposed,and the design method is discussed in this paper.The proposed method is used to design prototypes with typical parameters. According to the results of experiment and simulation, the elastic axis is customized with the sweep arrangement of the inflatable beams, indicate that the aeroelastic performance can be changed as well.

    The structural dynamic parameters are indispensable input in traditional aeroelastic analysis. Different from traditional structures, the structural mass of the inflatable wing is relatively light and the overall density is low, which is close to air for magnitude.Therefore, when the inflatable wing vibrates in the atmospheric environment,the air layer associated with the structure wall cannot be ignored. In this case, the structural dynamic characteristics are named as wet mode, and some targeted analysis methods are developed. Considering the characteristics of inflatable structures,the available wet mode analysis methods can be divided into fluid structure interaction(FSI)method[19,20]and added mass method[21—24], and have obtained certain applications respectively. The FSI method equates the influence of the flow field to an acoustic body and carries out modal analysis of the whole system. The FSI method has a more perfect mechanism explanation,but also leads to the difficulty of balancing its accuracy and cost. The additional mass method only retains the contribution of fluid to the equivalent mass of the system,greatly improving the efficiency of solution.In addition, some wet modal analysis methods for thin plate structures in fluid have been developed as well[25,26].However,due to significant geometric differences,it is obviously difficult to apply to inflatable structures. Therefore, there is still room for further improvement of a universal wet mode analysis method.

    Focus on the shortcomings of the existing methods, this paper proposes an improved wet mode analysis method by considering the spatial distribution and equivalent stiffness of the disturbed fluid,namely the added mass-stiffness method.Effectiveness of the proposed method is verified by ground vibration tests. On this basis, with the obtained modal parameters as input, the modal behavior of the swept baffled inflatable wing is analyzed and compared in detail. The remainder of this paper is organized as follows. In Section 2, the structural design method and aerodynamic performance of the swept baffled inflatable wing is analyzed, and the shape preserving performance of the swept baffled type is verified by experiment.In Section 3,the added mass stiffness method is proposed, and the example of an inflatable beam is selected to verify the accuracy of the established model.The wet mode of the inflatable wing is analyzed and the effects of some main factors are discussed in Section 4. The flutter performance is analyzed based on the frequency domain method in Section 5.

    2. Structural design of swept baffled inflatable wing

    2.1. Design of the traditional baffled type

    According to the bubble principle [9], an unconstrained 2D closed membrane structure tends to be circular under fully inflated conditions. With a flexible baffle arranged inside the membrane structure, two cavities are formed under inflation pressure, as shown in Fig. 2. When the inflation pressure in the cavities is consistent and stable, only the tensile force acts as a traction constraint on the baffle, and the membrane structure between nodes and tends to a new arc with as the common chord.Therefore,by arranging the specific length of the baffles at different positions inside the membrane structure. The membranes between the two nodes will form arcs with different radiuses, and the target airfoil can be approximated by the contour under the condition of the same pressure in the entire structure[6,9].Baffled inflatable wings with a limited span can be obtained by stretching and closing the designed 2D structure,as shown in Fig. 3.

    Fig. 3. Baffled inflatable wing.

    For thenth cavity of the inflatable wing shown in Fig. 4, equilibrium equations can be established according to the inflated static equilibrium as follows:

    Fig. 4. Force equilibrium of the n th cavity of the inflatable wing.

    wherepis the inflation pressure;Lnis the length of thenth baffle;are the normal stresses in the direction of the leading edge and trailing edge on the arc of thenth cavity,respectively;is the angle between the membrane and horizontal axis at the node in the direction of the leading edge and trailing edge of thenth cavity.

    According to the nodal force equilibrium [7]:

    For a 2D airfoil withNcavities under aerodynamic loads, the following equilibrium is satisfied in the direction of lift.

    whereLis the lift;α is the angle of attack(AOA);are the nodal force components at the upper and lower surface,respectively.

    For a single-layer plate-shaped inflatable membrane structure resembling an inflatable wing with a shape preserving ability, the existence of sections in the main bearing direction whose shape satisfies the bubble principle should be ensured.Thus,at least one section that is composed of straight baffles and surface arcs are required to ensure the shape preserving ability.

    2.2. Design of the swept baffled type

    Abstracting from the structural geometry of the baffled inflatable wing,the inflatable beam can be regarded as the basic bearing element of the inflatable wing [10], which is arranged and connected according to certain rules to satisfy the restriction of an airfoil. Therefore, the stiffness distribution and geometric characteristics of the inflatable wing can be changed by rearranging the inflatable beams.An effect closed to an aeroelastic tailoring design will be formed under the rearrangement,and the aerodynamic and aeroelastic characteristics of the inflatable wing will be changed accordingly.For a swept baffled rectangular inflatable wing,whose geometrical parameters and airfoil have been determined, the design procedure is shown in Fig. 5 and can be divided into the following steps.

    Fig. 5. Design procedure of the swept baffled inflatable wing.

    Fig. 6. Sections of the swept baffled inflatable wing with θ = 24?.

    (1) Initial selection of design variables.The main variables of the swept baffled inflatable wing are the baffles sweep angle of the inflatable beams θ and the spanwise distance between two adjacent beamsxd.

    (2) Formation of the inflatable beam. The chord length of the target airfoil can be stretched as a ratio ofrb, which can be determined as follows.

    After the trailing edge of the stretched airfoil is rounded, a complete inflatable beam can be formed by enveloping the normal circular section of the middle arc.

    (3) Formation of the effective wing segment. According to the previously selectedxdvalue,the inflatable beams are linearly arrayed in the spanwise direction. The effective wing segment can be intercepted, and the surface of the wing segment can be extracted as the skin.

    (4) Formation of the tip and baffles. After projecting the intersecting line between the inflatable beams to the wing plane,a baffle can be formed by filling in the region between the projected and intersecting lines. Then, the wingtip can be formed by filling in the region between the skin and baffle edges.

    According to the above design procedure, a series of swept inflatable wings with 13 cross-section inflatable beams,800 mm in span and 400 mm in chord,are designed.The range of θ is limited by the primary load bearing capability of the inflatable wing.Considering the integrity of the bearing element, at least one inflatable beam is required to be ensured to extend from root to tip.In conjunction with other geometric parameters of the prototypes,the boundary of θ is approximately±24?.The rounding radiusRris selected as 12 mm,andxis selected as 59.8 mm to ensure that the rounding of the trailing edge is tangent in the prototypes. These prototypes with critical θ will be used as the basis of subsequent research in this study. A geometric comparison shows that the inflatable wings with symmetrical or asymmetrical airfoils can be effectively designed by the method. Because of the swept arrangement of the inflatable beams,sections satisfying the bubble principle will form an angle with the chord direction, as shown section B-B in Fig.6.Sections of the chordwise inflatable beams are chamfered into ellipses,and the ripples are smoothed,as shown in Fig. 7. Therefore, the swept baffled type achieves a better approximation effect to the target airfoil than that of the baffle type under the same number of inflatable beams, which can effectively improve the aerodynamic performance.

    Fig. 7. Airfoil comparison of the inflatable wings.Note that two traditional configurations of inflatable wings can also be designed by the above design procedure. When θ = 0?, the inflatable beam is a uniform cylinder with an infinite length,which corresponds to a traditional baffled type.When θ =90?,a chordwise baffled inflatable wing will be obtained,as shown in Fig.8.The chordwise baffled type has a smoother chordwise profile and a predictable ideal aerodynamic performance before stalling. Because of the complete loss of spanwise support, the aerodynamic load is balanced by the tension provided by the distributed suspension cables.Therefore,the chordwise baffled type is mainly used for aircraft with low speed and load requirements,such as parasols or paragliders, which are not considered in present study.

    Fig. 8. Designed inflatable wings with θ = 90?.

    2.3. Shape preserving performance and aerodynamic performance

    To verify the shape preserving ability of the designed swept baffled inflatable wing, prototypes of the swept baffled type and traditional baffled type are manufactured with the same membrane material and design parameters as the simulation model. The VSTARS dual camera system based on the digital image correlation(DIC) method is used for testing. The process is shown in Fig. 9.Orientation of the test system is completed by photographing and calculating the 3D-coordinate information of spatially coded mark points at first. The measurement mark points are pasted on inflatable wings to capture and process the 3D-coordinate information of feature mark points.The V-STARS dual camera system is used to process the 3D-coordinate information of the measurement mark point set. A comparison between the marking points on the inflatable wings and target airfoil is shown in Fig.10.The geometric root-mean-square errors (RMES) between the inflatable wing and target airfoil are compared, and the results show that the swept baffled inflatable wing has better shape preserving ability.

    Fig. 9. The process of shape preserving test.

    Fig.10. Scatter plot comparison of inflatable wings.

    Studies on the aerodynamic characteristics of baffled inflatable wings show that the rippled surface will increase drag with delayed separation at a low Reynolds number [27]. Compared with the baffled type, ripples of the swept baffled type are smooth and extended spanwise obliquely. To verify the aerodynamic characteristics of the swept baffled inflatable wing, the steady flow field with different AOAs was solved by the CFD method based on the Reynolds averaged N—S (RANS) equation. The computational simulations are performed by ANSYS/Fluent.To ensure the reliability of the solution,the turbulence model is selected as SSTk—ω[27],and the upwind scheme is set to the second order based on reference to previous research.The size of the flow field is set to 12.5 times the characteristic size of inflatable wings.An unstructured mesh of the order of 5 million is adopted and locally densified near the wing surface, as shown in Fig.11. The boundary condition is set as the pressure far field, and the velocity of the incoming flow is set as 100 km/h (Ma0.082).

    Fig.11. Flow field model of inflatable wing.

    Fig.12. Comparation of the lift coefficient.

    Fig.13. Comparation of the drag coefficient.

    Fig.14. Comparation of the lift drag ratio.

    Fig.15. Comparation of the A-curve.

    Comparisons of the aerodynamic coefficient,lift drag ratio,and A-curve of inflatable wings under different angles of attack are shown from Figure 12—15, and listed in Table 1. For the lift drag ratio, the swept baffled types achieve an approximate increase before stalling compared with the baffled type. For the backward swept baffled type with θ= 24?, the trend of the lift coefficient is consistent with the baffled type, while a higher lift drag ratio is achieved after stalling.For the forward swept baffled type with θ=- 24?, there is a more severe lift loss and drag cost after stalling.

    Table 1Comparation of lift drag ratio.

    3. Wet modal modeling

    The structural dynamic parameters of inflatable wings need to be accurately analyzed as the input of aeroelastic analysis.Focus on the structural characteristics of inflatable wings and the shortcomings of existing methods, an additional mass-stiffness method is developed in this paper. Based on the added mass method, the proposed method fully considers the spatial distribution and stiffness contribution of fluid, and achieves both accuracy and efficiency. The basic principle is briefly introduced as follows.

    3.1. Pre-stressed membrane theory

    For inflatable structures,the membrane balances the out-plane load with the in-plane tension, which is mainly provided by the pressure difference between the internal pressure of the inflation and the external pressure of the atmosphere.Assume that the prestressed membrane is isotropic, and the thickness ish, the equilibrium equation of the pre-stressed membrane element is shown as follows:

    whereTx,Tyis the in-plane tension of membrane element in the direction ofxandy,pis the pressure.Tx,Tycan be expressed as follows:

    whereEmis the Youngs' modulus of the membrane material.

    The statics equilibrium equation of the pre-stressed membrane structure can be obtained with Eqs. (10) and (9).

    When the pre-stressed membrane structure vibrates freely,the shape of the membrane with inflation and deformation is taken as the reference configuration. Assume that the plane namedxoyis consistent with the plane where the membrane deforms, and the vibration displacement , which is vertical to thexoyplane, is considered as minim, thus the tension change caused by vibration can be ignored.As shown in Fig.16,the motion differential equation of pre-stressed membrane structure when vibrates freely can be expressed as follows:

    Fig.16. Vibration of the pre-stressed membrane structure [28].

    where ρ is the surface density of the pre-stressed membrane element dxdy,TxandTyis the membrane tension of membrane element.

    3.2. Wet modal modeling based on added mass-stiffness method

    Because of the lightweight of the inflatable wing structure,ambient air occupies an indispensable proportion of mass in the vibration system, and the added mass effect should be taken into consideration [29,30]. Therefore, a dynamic characteristic analysis should be conducted considering the influence of the flow field on the inflatable wing. This is referred to as a wet modal analysis.Existing wet modal analysis methods mainly include the FSI method and added mass method[31].In FSI method,the flow field is treated as the acoustic field,and an asymmetric solver is adopted[20]. Under the acoustic pressure energy diffusion effect, the accuracy of the FSI method is determined by geometry scale of flow field,which leads to an imbalance between cost and accuracy[23].The added mass method simplifies the flow field to mass nodes distributed on the surface of the structure, which has the advantages of a high efficiency and accuracy [31]. However, the contribution of the flow field to the system stiffness is ignored,resulting in an insufficient description of the wet mode mechanism.

    Focusing on the inflatable wing, a wet modal analysis method based on the added mass method and further considering the contribution of the flow field to the system stiffness is adapted,which is named as the added mass stiffness method here. The external added mass of slender structurecan be expressed as follows [32]:

    whereBis the width of the structure,dis the average draft,andMVis displacing. For the inflatable membrane structure, the internal pressured gas vibrates with the structure,and masscan be expressed as follows:

    where ρflu(p)is the density function about the pressure of the fluid;Vis the volume of the pressured gas. The fluid added mass of inflatable membrane structureMaddedcan be expressed as follows.

    The shape function of the flow field is determined according to the principle of the average. The volume of flow fieldVfcan be expressed as the following.

    According to the principle of average, the flow field is determined based on the structural geometry or boundary conditions.The shape function of flow fieldNfcan be expressed as follows:

    whereNsis the shape function of the inflatable structures,n is the normal vector of the inflatable structure shape function and Δxis the distance between the inner and outer boundary of the flow field.

    The shape function of the flow field is the same as that of the inflatable structure at the inner boundary.Therefore,the flow field shares nodes with the inflatable structure at the inner boundary.The shape function of the flow field at the outer boundary is determined according to the distance between the inner and outer boundary evenly. For the constrained boundary,the flow field and inflatable structure are fixed because of the absence of vibration.

    For homogeneous fluids, the form of the governing equation in micelle is similar to that of a solid [33]. Air is able to resist stress owing to compressibility and viscosity, while vibrating accompanied by the inflatable structure [34]. According to the relationship between the elastic coefficients of isotropic materials,an equivalent modulus of air under a small disturbance can be expressed as follows [35]:

    whereKis the bulk modulus.The equivalent model of the flow field can be obtained based on Eqs.(15)—(18)and the contribution of the flow field to system stiffnessKaddedcan be expressed as follows:whereIfis the inertia of the flow field.

    3.3. The validation of an inflatable tube

    In order to verify the accuracy of the established model, the inflatable ethylene-terafluoroethlene (EFTE) tube is selected as an example to be verified,and compared with the experiment results.The geometric, material parameters and boundary conditions are consistent with previous studies[23].The finite element(FE)model established by added mass stiffness method in this paper is shown in Fig.17,the comparison with other wet modal analysis method is listed in Table 2, and the modal shapes are shown in Fig.18.

    Table 2Natural frequencies comparison of the inflatable beam.

    Fig.17. The FE model of the inflatable beam based on added mass-stiffness method.

    Fig.18. Mode shapes of the inflatable beam.

    According to the results shown in Table 2, there exist a significant difference between dry mode and wet mode,which proves the necessity of wet modal analysis.Because there is no need to adopt the asymmetric solver, compared with the FSI method, the efficiency of added mass stiffness method is significantly improved.While compared with the traditional added mass method[36],the contribution of the flow field to the system stiffness is well considered, which improves the estimation performance for the natural frequencies of inflatable structures. For the validation example in this paper, the relative error of first two natural frequencies between experiment results and the simulation results based on the added-mass stiffness method is 3.2%and 0.8%,which is able to satisfy the accuracy need for the wet modal analysis of inflatable wings.

    4. Wet modal analysis and experimental verification

    4.1. Modeling

    To further verify the accuracy of the proposed added massstiffness method, the modal experiments of the inflatable wings are conducted. The modal experiment system is set up shown as Fig.19, including 4 parts, which are the inflatable wing structure,the force and acceleration sensors, the measurement and analysis system,and the pedestal.The modal experiment image is shown as Fig. 20.

    Fig.19. Vibration experimental system of the inflatable wing.

    Fig. 20. Modal experimental image of the inflatable wing.

    Fig. 21. Modal experimental results of the inflatable wing with θ = 0?.

    As shown in Fig.22,the baffles sweep angles θ of three inflatable wings tested in this paper is 0?in Fig.22(a),-24?in Fig.22(c),and 24?in Fig. 22(d) respectively. The geometry parameters of the inflatable wings are the same as the ones designed in subsection 2.2, where the airfoil is NACA0015, chord length is 400 mm, span length is 800 mm, and the number of cross-section inflatable beams is 13. All three inflatable wings are made of high strength composite flexible membrane material, and the material parameters are shown in Table 3. The flexible fiber composite material is orthotropic material, and the mechanical property difference between warp and weft is not remarkable. For the sake of engineering estimating, the membrane is considered isotropic[37—39]. Besides, the inflation pressure is set as 30 KPa. The pedestal structure is made of the aluminum profile and sandbags,where the inflatable wing is able to be fixed at one end while free at the other. The environment parameter of the modal experiment is(25 ± 3)?C, and the relative humidity is 50%.

    Table 3Material parameters of the inflatable membrane.

    Fig. 22. Modal experimental results of the inflatable wing with θ = - 24?.

    Fig. 23. Modal experimental results of the inflatable wing with θ = 24?.

    As shown in Fig.20(b),the hammer chosen in this paper is PCB 086C03,and the air capsule hammerhead is adapted to improve the effect of excitation for flexible inflatable wings [40]. The force and acceleration sensors are used to measure the excitation and response signal of the tested structure.As shown in Fig.20(a),nine accelerometer named PCBTM 333B30 are placed evenly on and under the inflatable wing surface, and the excitation point is selected as the node near the trailing edge at the tip of the inflatable wing, which is node 9 for the inflatable wings with θ equaling 0?and 24?, and node 7 when θ = - 24?.The measurement and analysis system is LMS Test. Lab SCADS III, which is used to define the sampling frequency and sampling time, and process the force and acceleration signal measured from sensors. The sampling frequency is set as 512 Hz,and the sampling time is 8 s. To eliminate the measurement error as much as possible, each test is repeated and averaged for 5 times under certain pressure. The averaged frequency response function(FRF)is identified with poly-reference least squares complex frequency method (PolyMAX). The tests results are shown as Fig. 21, Figs. 22 and 23 which provides the reference for simulation of inflatable wings.

    For the inflatable wing with θ = 0?, also known as the traditional baffled inflatable wing, the first five modes with distinct modal shapes are got under 100 Hz,which include the 1st bending,1st torsion, 2nd bending,1st chordwise bending, and 2nd torsion.Among the five modes,the 1st chordwise bending mode is a unique modal shape of the inflatable wings compared with the rigid wings,which will participate in the flutter [28]. Meanwhile, for the inflatable wings with θ=±24?under the same bandwidth, there are only four modes and the eliminated one is the 1st chordwise bending mode, which will help with the flutter performance.

    4.2. Wet mode analysis and discussion of factors of the inflatable wings

    Based on the added mass stiffness method, a simulation of inflatable wings with different θ are performed by a comparative analysis. The FE model is shown as Fig. 24. Considering the influence of external and internal flow field, the internal flow field is discounted to the external flow field to simply the modeling process, and the volume of external flow field is obtained by Eq. (16).The inflatable wing is modeled with Shell181 element and is fixed at one end.Different from the FSI method,the added-mass stiffness method is able to save the computing cost for only considering the local flow field participate in vibration. Compared with the traditional added-mass method, the proposed method shows the advantage in efficiency and accuracy of estimating the dynamic characteristic of the inflatable structures. The simulation results identified with Lanczos method are listed in Table 4 and shown in Figs.25-27.The inflation pressure is set as 30 kPa,which is the same as the one in modal experiment and suitable enough to make sure the membrane is completely tensioned.

    Table 4Natural frequencies comparison of the inflatable wings.

    Fig. 24. The FE model of the inflatable wing based on added mass-stiffness method.

    Fig. 25. Simulation results of the inflatable wing with θ= 0?.

    Fig. 26. Simulation results of the inflatable wing with θ= - 24?.

    Fig. 28. Natural frequencies of the inflatable wing under various internal pressure.

    The simulation results of the inflatable wing with θ=0?reflect five modes,the 1st bending,1st torsion,2nd bending,1st chordwise bending,and 2nd torsion,which shows the same tendency with the modal experiment results.For the swept inflatable wings with θ =±24?,there are some changes in natural frequencies resulted from the deflection of stiffness axis, and the 1st chordwise bending is also eliminated,due to the exist of the baffles sweep angle θ,which means the baffles inside the inflatable wings are able to provide the structural support both spanwise and chordwise.The average error shown in Table 4 indicated the consistency of the natural frequency between the simulation and experiment results, which proves the effectiveness of the proposed method. The averaged error of the inflatable wings with θ=0?and θ=-24?is 8.3% and 5.3%, which are all under 10%and proves the accuracy of the proposed method.For the inflatable wing with θ =24?,the natural frequency is lower due to the process limitation of the tested inflatable wing. The manometer is hung at the free end,and despite the slinging of the elastic rope, the mass of vibration system increases inevitably.Generally, the proposed method is precise enough to estimate the dynamic characteristics of the inflatable wings.The ratio of the first bending frequency to the torsional frequency of the swept baffled type is higher than that of the baffled type, which is beneficial for increasing the flutter critical speed according to Pines’theory[41].Based on the conclusions above,the swept configuration may affect the flutter characteristics of the inflatable wings.

    The swept baffled inflatable wings should be fully inflated to maintain the system stiffness and resist bulking. Therefore, the influence of internal pressure should be analyzed. The specific effect of baffles sweep angle on dynamic characteristics of inflatable wings is still unclear and worth discussing.Besides,for aircraft with twin-boom configuration, the horizontal tail is fixed at both ends and has special boundary conditions. Therefore, the influence of boundary conditions of inflatable wing on the structural dynamic characteristics should be deeply analyzed as well.

    4.2.1. Effect of internal pressure

    For the models with various internal pressure, the inflatable wing with θ=0?is selected and the pressure is set from 4 to 28 kPa.The inflatable wing shows five modes under 100 Hz according to the simulation results are shown as Fig.28,which is consistent with the experiment and simulation results above. All natural frequencies show an upward tendency due to the increasing structure stiffness as the internal pressure rises,and the growth tendency is also slowing down, which is more obvious for the high order modes.For inflatable wings,the stiffness distribution of structure is settled when the internal pressure is enough to fully tension the membrane structure. There will not be distinctive changes of the natural frequencies when the internal pressure increases further.Based on the experiment and simulation results, the reference pressure for the inflatable wing in this paper is 30 kPa, which is proper for both simulation and experiments.

    4.2.2. Effect of baffles sweep angle

    The swept angel θ is assumed positive for the sweptback wings.According to the definition, it's the traditional baffled inflatable wing when θ = 0°. In the progress of modeling, at least one inflatable beam should extend form the root to tip, which guarantees the ability for inflatable wings to resist bulking. Hence,there is an upper limit for the baffles sweep angle according to the aspect ratio and number of cross-section inflatable beams. The simulation for inflatable wings with various baffles sweep angles is performed and the results are shown as Fig. 29. The internal pressure is set as 30 kPa,and the geometry and material parameters are same as the ones above.

    Fig. 29. Natural frequencies of the inflatable wing with various baffles sweep angle θ.

    Fig. 30. Simulation results of the inflatable wing with θ= 0?.

    Fig. 31. Simulation results of the inflatable wing with θ = ±24?.

    Under 100Hz, five modes are obtained for the inflatable wing with θ = 0°, while only four modes for other configurations with nonzero baffles sweep angle θ.The eliminated mode is 1st bending chordwise, which is consistent with the experiment results above.Generally,the natural frequencies show an upward tendency as the absolute value of baffles sweep angle θ decreases. The length of main bearing inflatable beam for the configuration with nonzero baffles sweep angle is larger than the one with θ = 0°. Due to the exist of swept angel θ,the inflatable beams have been participated in circumferential bearing. The membrane structure is barely able to withstand the out-plane load,and thus the stiffness distribution is changed and the equivalent stiffness is also decreased, and results in the decrease of the natural frequencies. Specially, the natural frequency of Mode3 for inflatable wing with θ=-24°is anomalous, and the average error for Mode3 is unsatisfactory.Compared with other order modes. In the modeling process for swept inflatable wings, the intercepted is performed based on the length of the wingspan,there is still the uncertainty of the relative location to cut off, which considers as the reason of unsatisfactory averaged error for Mode3, the 2nd bending.

    4.2.3. Effect of boundary condition

    The internal pressure is set as 30 kPa, and the inflatable wings with θ=0 and θ=±24°are selected. Based on the one end fixed,the boundary condition of two end fixed is also performed,and the results are shown as Figs. 30 and 31.

    Figs. 30 and 31 shows the first five order modal parameters of the inflatable wing with θ=0°and θ=±24°when fixed at two ends. Generally, the natural frequencies of the models with two ends fixed are larger than the ones with one end fixed. For the inflatable wing with θ = 0°, the modal shape is symmetric due to the symmetric configuration. And compared with the inflatable wing with one end fixed,which can be considered to be a cantilever plate, the bending deformation chordwise is obtained for Mode1 and Mode4 in addition to the bending and torsion deformation spanwise.For the inflatable wing with θ =±24°,the modal shape is unsymmetrical because of the swept angle.

    5. Analysis of the aeroelasticity performance

    5.1. Flutter modeling

    For the segment of the inflatable wing shown in Fig. 32, the governing equations can be expressed as follows [36]:

    Fig. 32. Motion model of the inflatable wing segment.

    wheremis the mass of the segment;his the deflection at the elastic axis of the segment;Sα is the static moment of the segment;KhandKαare the bending and torsional stiffness of the segment, respectively;Iα is the inertia of the elastic axis section;Mis the aerodynamic moment.The wing is assumed to vibrate harmonically at the flutter critical speed. The generalized displacement can be expressed as follows:

    For inflatable wings suitable for low-speed flight, the flutter determinant can be obtained simultaneously with Eqs. (22) and(23).

    where ωαand ωhare the natural frequencies corresponding to the bending and torsion mode, respectively;xα is the length between the center of gravity and elastic axis;rαis the radius of gyration of the elastic axis;ais the length between the elastic axis and half chord;xα andacan be expressed as follows:

    wherexGis the chordwise coordinate of the center of gravity.

    All elements excepta11will be affected by the center of rigidity distribution in the flutter determinant. Because the flutter critical velocity is determined by eigenvalues of the corresponding coefficient matrix of the flutter determinant, the center of rigidity distribution will affect the velocity critical flutter.Based on the theory mentioned above,the flutter critical velocity of the inflatable wing is analyzed by the V-g method.The structural damping coefficientgis introduced into the motion equation of the wing in the V-g method,and the sum of the structural damping force and resilience can be expressed as follows:

    Substituting Eq. (26) into Eq. (22), the elements of the flutter determinant are transformed into the following:

    An intermediate variableZis introduced,as shown in Eq.(28);gcan be expressed in the form of Eq.(29).Flutter frequency ω can be obtained by solving the flutter determinant wheng=0 [41].

    5.2. Discussion of the baffles sweep angle for aeroelastic performance

    By definition,loading on the elastic axis will not twist the wing.A simple static simulation is repeated to determine the location of the elastic axis, as shown in Fig. 33. The center of gravity can be directly obtained through computer aided design (CAD) software and the center of pressure is assumed to be a quarter of the chord length from the leading edge in low-speed flight. The relevant parameters of the inflatable wings are listed in Table 5,according to previous research on the flutter analysis of an inflatable wing with similar geometric parameters [36]. The Theodorsen theory is selected as an unsteady aerodynamic model, while natural frequencies and modal shapes functions of the wet modal are used as the structural discipline input;gis taken as 0 to obtain a conservative value[42].The flutter critical speed and V-g graphs are listed in Table 6 and shown in Fig. 34, respectively.

    Table 5Flutter model parameters of the inflatable wing.

    Table 6Flutter critical speed comparison of the inflatable wing.

    Fig. 33. Calculation process of elastic axis distribution.

    Fig. 34. Comparison V-g graphs of the inflatable wings.

    Fig. 35. Elastic axis distribution of the inflatable wing.

    Fig. 36. Comparison V-ω graphs of the inflatable wings.

    The flutter critical velocities of the inflatable wing increase with the increase of the θ,as shown in Table 6 and shown in Fig.34.The main reason is that the swept arrangement of the baffles’ leads to the change of the elastic axis distribution of the inflatable wing.As shown in Fig.35,the elastic axis of the inflatable wing changes with the direction of the baffle,and even cross the center of pressure axis and the center of gravity axis.This will lead to large changes in the aerodynamic term and inertia term in different segments. For the 24?backward swept configuration, the elastic axis is closer to the trailing edge at the wingtip than the center of gravity, thus inhibiting the classical flutter caused by the coupling of low-order modes[42]. For the - 24?forward swept configuration, the elastic axis is closer to the leading edge at the wing tip, which leads to the reduction of the flutter critical speed. At the same angle, both the forward swept and the backward swept configurations will lead to significant changes in the distribution of the elastic axis over the entire wing. Since the mode shape functions increase monotonically in the spanwise, the effect of the sweep angle on the flutter speed can be estimated mainly by the relative position of the elastic axis at the wing tip.

    At the same time, the frequencies of torsional modes’ decrease with the increase in speed,and the most obvious decrease is in the 24?swept back configuration, which indicates that it may have a lower divergence speed,as shown in Fig.36.This may also be due to the distribution of the elastic axis, which results in the maximum torque at the wingtip of the swept back configuration under similar aerodynamic loads.Note that with the aerodynamic model selected in this paper, all cases do not have torsional divergence, but this does not mean that inflatable wings are safe at the speed before flutter. For the nonlinear behavior of inflatable wing failure, previous research have given more accurate and specific methods[43].The focus of this paper is to reveal the influence trend of sweep angle on the aeroelastic behavior of inflatable wings, rather than accurately predict the divergence speed.

    In general, the contribution of swept back configuration and swept forward configuration to aeroelastic characteristics of inflatable wing is opposite.Flutter and divergence can be inhibited by designing different sweep angles, respectively. This indicates that the proposed structure may be enlightening for the aeroelastic tailoring design of inflatable wings. Considering that the swept baffled type has better aerodynamic performance, the inflatable winged loitering munitions with this structural scheme is expected to achieve higher useable speed.

    6. Conclusions

    (1) Compared with the traditional baffled inflatable wing, the swept baffled inflatable wing can achieve better approximation and aerodynamic performance to target the airfoil after inflation.The proposed method can effectively provide the geometry of inflatable wings with different baffles sweep angles.

    (2) The stiffness of the surrounding fluid will affect the modal parameters of the structure. Therefore, compared with the traditional added mass method, the proposed added mass stiffness method can predict the wet modal parameters of inflatable structures more accurately.

    (3) Take the internal pressure, the baffles sweep angle and the boundary condition as effect factors into consideration. The natural frequencies are rising as the internal pressure is increasing, and the absolute value of baffles sweep angle is decreasing. The larger natural frequencies and more chordwise bending modal shapes are obtained for the inflatable wings with both ends fixed.

    (4) The aeroelastic characteristics of the inflatable wing vary with the baffles sweep angle. With the increase in baffles sweep angle,the flutter performance of the inflatable wing is improved and the divergence performance is degraded.This conclusion is considered having potential value for aeroelastic tailoring design of inflatable wing.

    Declaration of competing interest

    The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Junhui Meng reports financial support was provided by National Natural Science Foundation of China.

    Acknowledgements

    This research was supported by National Natural Science Foundation of China (Grant No.11902029). The authors thank all the people involved in the past and present progress of the experiment. The authors also are grateful to the reviewer and the executive editor for their precious suggestions about this paper.

    黄色毛片三级朝国网站| 国产精品一区二区精品视频观看| 亚洲国产高清在线一区二区三| 不卡一级毛片| 婷婷六月久久综合丁香| 露出奶头的视频| 这个男人来自地球电影免费观看| 午夜激情福利司机影院| 天堂影院成人在线观看| 欧美不卡视频在线免费观看 | а√天堂www在线а√下载| 丝袜人妻中文字幕| 老司机福利观看| 两人在一起打扑克的视频| 午夜免费激情av| 国产v大片淫在线免费观看| 午夜精品一区二区三区免费看| 国产午夜精品论理片| 亚洲人成网站在线播放欧美日韩| 久久婷婷成人综合色麻豆| 久久精品国产综合久久久| 99国产极品粉嫩在线观看| 中文在线观看免费www的网站 | 亚洲人成网站高清观看| 欧美日韩国产亚洲二区| 大型黄色视频在线免费观看| 免费高清视频大片| 色精品久久人妻99蜜桃| 黄色a级毛片大全视频| 人人妻人人看人人澡| 亚洲专区国产一区二区| 丁香欧美五月| 淫秽高清视频在线观看| 欧美乱妇无乱码| 日本撒尿小便嘘嘘汇集6| 日本一本二区三区精品| 亚洲色图av天堂| 亚洲av成人不卡在线观看播放网| 1024香蕉在线观看| 天天添夜夜摸| 1024视频免费在线观看| 少妇的丰满在线观看| 国内精品久久久久精免费| 成人高潮视频无遮挡免费网站| 国内少妇人妻偷人精品xxx网站 | 国产高清有码在线观看视频 | 两个人看的免费小视频| 午夜激情av网站| 在线十欧美十亚洲十日本专区| 欧美日韩黄片免| av在线播放免费不卡| 90打野战视频偷拍视频| 久久午夜亚洲精品久久| 黑人欧美特级aaaaaa片| 欧美日韩国产亚洲二区| 日韩欧美国产一区二区入口| 视频区欧美日本亚洲| 亚洲五月婷婷丁香| 欧美日韩亚洲国产一区二区在线观看| 成人高潮视频无遮挡免费网站| 人人妻,人人澡人人爽秒播| 亚洲国产欧洲综合997久久,| 亚洲精品久久国产高清桃花| 亚洲精品国产一区二区精华液| netflix在线观看网站| 精品国产超薄肉色丝袜足j| 精品久久久久久,| 天堂√8在线中文| 欧美黑人巨大hd| 亚洲一区高清亚洲精品| 国产免费男女视频| 亚洲国产精品sss在线观看| 国产一区二区在线av高清观看| 啪啪无遮挡十八禁网站| 精品午夜福利视频在线观看一区| 国产成+人综合+亚洲专区| 999久久久精品免费观看国产| 国产不卡一卡二| 精品久久久久久久毛片微露脸| a在线观看视频网站| 欧美日韩亚洲国产一区二区在线观看| 国产成+人综合+亚洲专区| 日本免费一区二区三区高清不卡| 国产三级黄色录像| 禁无遮挡网站| 国产成人av激情在线播放| 亚洲色图av天堂| 亚洲美女黄片视频| 视频区欧美日本亚洲| 俺也久久电影网| 法律面前人人平等表现在哪些方面| 欧美日本视频| www.自偷自拍.com| 国产精品一区二区三区四区久久| 久久午夜亚洲精品久久| 亚洲 欧美一区二区三区| 欧美日韩黄片免| 国产精品,欧美在线| 母亲3免费完整高清在线观看| 又爽又黄无遮挡网站| 欧美日本视频| 午夜福利在线在线| 免费观看人在逋| 在线观看免费午夜福利视频| 欧美久久黑人一区二区| 亚洲精品在线美女| 亚洲熟妇熟女久久| 嫁个100分男人电影在线观看| 夜夜夜夜夜久久久久| 亚洲精品av麻豆狂野| 亚洲人成网站在线播放欧美日韩| 亚洲美女视频黄频| 午夜福利在线观看吧| 桃色一区二区三区在线观看| 日韩欧美在线乱码| 天天躁夜夜躁狠狠躁躁| 九九热线精品视视频播放| 久久久久九九精品影院| 精品午夜福利视频在线观看一区| www日本在线高清视频| 99国产极品粉嫩在线观看| 亚洲中文av在线| 在线观看免费日韩欧美大片| 高潮久久久久久久久久久不卡| 国产亚洲精品久久久久久毛片| 男女做爰动态图高潮gif福利片| 狠狠狠狠99中文字幕| 最新在线观看一区二区三区| 99国产精品一区二区三区| 久久欧美精品欧美久久欧美| 91成年电影在线观看| 亚洲人成电影免费在线| 亚洲av片天天在线观看| 色在线成人网| 欧美zozozo另类| 性欧美人与动物交配| 亚洲午夜理论影院| 人人妻,人人澡人人爽秒播| 九色成人免费人妻av| 制服诱惑二区| 啦啦啦韩国在线观看视频| 中文亚洲av片在线观看爽| 欧美午夜高清在线| 俺也久久电影网| 久久午夜综合久久蜜桃| 成人av一区二区三区在线看| a级毛片a级免费在线| 女警被强在线播放| 国产成+人综合+亚洲专区| 最好的美女福利视频网| 国产精品自产拍在线观看55亚洲| 狂野欧美白嫩少妇大欣赏| 日韩欧美免费精品| 中文字幕精品亚洲无线码一区| 人妻丰满熟妇av一区二区三区| 亚洲中文字幕日韩| 久久精品国产清高在天天线| 久久久久九九精品影院| 免费观看精品视频网站| 亚洲国产高清在线一区二区三| 丰满的人妻完整版| 18禁观看日本| 床上黄色一级片| 1024手机看黄色片| 久久精品国产综合久久久| 日韩欧美国产一区二区入口| 19禁男女啪啪无遮挡网站| 高清毛片免费观看视频网站| 这个男人来自地球电影免费观看| 听说在线观看完整版免费高清| 哪里可以看免费的av片| 长腿黑丝高跟| 一区二区三区高清视频在线| 亚洲专区国产一区二区| 日韩精品免费视频一区二区三区| 久久天躁狠狠躁夜夜2o2o| 久久伊人香网站| 激情在线观看视频在线高清| 在线观看一区二区三区| 久久久久九九精品影院| 日日摸夜夜添夜夜添小说| 18禁观看日本| 国内毛片毛片毛片毛片毛片| 色av中文字幕| 小说图片视频综合网站| 日韩三级视频一区二区三区| 热99re8久久精品国产| 香蕉久久夜色| 色噜噜av男人的天堂激情| 一个人观看的视频www高清免费观看 | 国产私拍福利视频在线观看| 国产欧美日韩一区二区精品| 国产午夜精品论理片| 五月伊人婷婷丁香| 又大又爽又粗| 国产精品野战在线观看| 美女免费视频网站| 99国产精品一区二区蜜桃av| 制服人妻中文乱码| 伊人久久大香线蕉亚洲五| 亚洲九九香蕉| 国产成人啪精品午夜网站| 国产精品免费视频内射| АⅤ资源中文在线天堂| 99久久精品国产亚洲精品| 国产亚洲精品av在线| 人人妻人人澡欧美一区二区| 国产成人精品久久二区二区91| 成人亚洲精品av一区二区| 白带黄色成豆腐渣| 国产高清视频在线观看网站| 夜夜夜夜夜久久久久| 久久久久国产精品人妻aⅴ院| av国产免费在线观看| 亚洲精品中文字幕一二三四区| 国产激情久久老熟女| 一级黄色大片毛片| 男女下面进入的视频免费午夜| 欧美黑人巨大hd| 国内精品久久久久久久电影| 搞女人的毛片| 岛国视频午夜一区免费看| 一进一出抽搐动态| 一本大道久久a久久精品| 亚洲国产看品久久| 男女视频在线观看网站免费 | avwww免费| 老鸭窝网址在线观看| 精品久久久久久久久久久久久| 亚洲 欧美 日韩 在线 免费| 露出奶头的视频| 国产成人欧美在线观看| 高潮久久久久久久久久久不卡| 久久欧美精品欧美久久欧美| 欧美成人一区二区免费高清观看 | 好看av亚洲va欧美ⅴa在| 亚洲欧美日韩无卡精品| 又大又爽又粗| 亚洲av成人一区二区三| 日韩精品免费视频一区二区三区| 青草久久国产| 人人妻人人看人人澡| 日韩成人在线观看一区二区三区| 久久精品成人免费网站| 国产熟女xx| 午夜影院日韩av| 91av网站免费观看| 久久午夜综合久久蜜桃| 久久精品人妻少妇| 日本 欧美在线| 啦啦啦韩国在线观看视频| 2021天堂中文幕一二区在线观| 国产成人精品久久二区二区91| 日韩欧美国产在线观看| 少妇粗大呻吟视频| 狠狠狠狠99中文字幕| 中文字幕最新亚洲高清| 国产爱豆传媒在线观看 | 国产精品久久电影中文字幕| 日本成人三级电影网站| 久久精品影院6| 老司机午夜十八禁免费视频| 国产精品一区二区三区四区免费观看 | 亚洲中文字幕日韩| 亚洲一码二码三码区别大吗| 丁香六月欧美| 88av欧美| 亚洲欧美激情综合另类| 999精品在线视频| 国产精品精品国产色婷婷| 成年女人毛片免费观看观看9| 麻豆成人av在线观看| av有码第一页| 欧美中文综合在线视频| 91字幕亚洲| 性色av乱码一区二区三区2| 亚洲男人的天堂狠狠| 老汉色∧v一级毛片| 天天添夜夜摸| 国产av一区在线观看免费| 91字幕亚洲| 亚洲男人的天堂狠狠| 久热爱精品视频在线9| or卡值多少钱| 九色国产91popny在线| 久久精品人妻少妇| 这个男人来自地球电影免费观看| 久久精品亚洲精品国产色婷小说| 一进一出好大好爽视频| 亚洲一区高清亚洲精品| 国产精品一区二区三区四区免费观看 | 黄色视频不卡| 欧洲精品卡2卡3卡4卡5卡区| 亚洲全国av大片| 久久精品国产亚洲av高清一级| 91老司机精品| 99re在线观看精品视频| 别揉我奶头~嗯~啊~动态视频| 亚洲一区二区三区不卡视频| 搡老岳熟女国产| 男插女下体视频免费在线播放| 99久久99久久久精品蜜桃| 国产精品香港三级国产av潘金莲| 18禁黄网站禁片午夜丰满| 岛国在线观看网站| 88av欧美| 欧美丝袜亚洲另类 | 亚洲在线自拍视频| 99国产极品粉嫩在线观看| 国产精品久久视频播放| 久久久久免费精品人妻一区二区| 88av欧美| 很黄的视频免费| av福利片在线观看| 在线观看美女被高潮喷水网站 | 三级男女做爰猛烈吃奶摸视频| 九色成人免费人妻av| 成年免费大片在线观看| 动漫黄色视频在线观看| 久久亚洲真实| 欧美日韩国产亚洲二区| 香蕉av资源在线| xxxwww97欧美| 国产一级毛片七仙女欲春2| 久久香蕉国产精品| 一卡2卡三卡四卡精品乱码亚洲| 亚洲成人久久性| 亚洲黑人精品在线| 亚洲av中文字字幕乱码综合| 黄片小视频在线播放| 国产日本99.免费观看| 精品国产亚洲在线| 69av精品久久久久久| 国产av一区在线观看免费| 国产一级毛片七仙女欲春2| 国产欧美日韩精品亚洲av| 久久草成人影院| 久久伊人香网站| 757午夜福利合集在线观看| 岛国视频午夜一区免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲午夜精品一区,二区,三区| 天堂动漫精品| 婷婷丁香在线五月| 欧美av亚洲av综合av国产av| 亚洲人成电影免费在线| 亚洲中文字幕日韩| 亚洲黑人精品在线| 久久精品亚洲精品国产色婷小说| 99在线人妻在线中文字幕| 国产精品久久电影中文字幕| 久久人人精品亚洲av| 一本久久中文字幕| 色老头精品视频在线观看| 色综合站精品国产| 在线播放国产精品三级| 日韩欧美精品v在线| 成人手机av| 国内精品一区二区在线观看| 亚洲精品粉嫩美女一区| 久久国产精品人妻蜜桃| 亚洲国产精品999在线| 两性夫妻黄色片| 老司机在亚洲福利影院| 国产亚洲精品久久久久5区| 日本一本二区三区精品| 久久99热这里只有精品18| 国产97色在线日韩免费| 两个人的视频大全免费| 麻豆av在线久日| 国产精品亚洲美女久久久| 亚洲熟妇熟女久久| 亚洲人成网站高清观看| 99久久精品热视频| a在线观看视频网站| 欧美性长视频在线观看| 天天添夜夜摸| 国产激情偷乱视频一区二区| 亚洲片人在线观看| 亚洲人与动物交配视频| 日韩中文字幕欧美一区二区| 99久久无色码亚洲精品果冻| 欧美日本视频| 变态另类丝袜制服| 久久久久精品国产欧美久久久| 99国产极品粉嫩在线观看| 一边摸一边抽搐一进一小说| 国产激情偷乱视频一区二区| 十八禁网站免费在线| 国产成人一区二区三区免费视频网站| 老司机靠b影院| 在线观看午夜福利视频| 欧美精品亚洲一区二区| 国产成人精品无人区| 国产av麻豆久久久久久久| 国产91精品成人一区二区三区| 久久精品国产99精品国产亚洲性色| 午夜视频精品福利| a在线观看视频网站| 岛国在线观看网站| 性色av乱码一区二区三区2| 精品高清国产在线一区| 国产黄a三级三级三级人| 观看免费一级毛片| 国产又色又爽无遮挡免费看| 天堂影院成人在线观看| 一二三四社区在线视频社区8| 国产探花在线观看一区二区| 舔av片在线| 国产免费男女视频| a级毛片在线看网站| 小说图片视频综合网站| 99在线视频只有这里精品首页| 精品电影一区二区在线| 99热6这里只有精品| 精品电影一区二区在线| 久99久视频精品免费| 欧美成人午夜精品| 黄色视频,在线免费观看| 国产av麻豆久久久久久久| 国产精品一区二区三区四区久久| 成人18禁高潮啪啪吃奶动态图| 一区福利在线观看| 成年女人毛片免费观看观看9| 人人妻,人人澡人人爽秒播| 在线观看日韩欧美| 999精品在线视频| 禁无遮挡网站| 一级黄色大片毛片| cao死你这个sao货| 国内精品一区二区在线观看| 一夜夜www| 午夜日韩欧美国产| 国产av一区二区精品久久| 最近最新中文字幕大全免费视频| 高潮久久久久久久久久久不卡| 精品乱码久久久久久99久播| 看黄色毛片网站| 一进一出好大好爽视频| 国产av一区在线观看免费| 在线观看一区二区三区| www.精华液| 男女下面进入的视频免费午夜| 99久久综合精品五月天人人| 国产成人系列免费观看| 午夜激情福利司机影院| 天天躁夜夜躁狠狠躁躁| 久久香蕉激情| 身体一侧抽搐| 欧美乱妇无乱码| 国产精品免费视频内射| 亚洲熟妇中文字幕五十中出| 熟女少妇亚洲综合色aaa.| 五月伊人婷婷丁香| 日本撒尿小便嘘嘘汇集6| 日本五十路高清| 久久久精品国产亚洲av高清涩受| 日韩免费av在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 两个人的视频大全免费| 在线观看日韩欧美| 亚洲全国av大片| 狂野欧美白嫩少妇大欣赏| 黄色 视频免费看| 久久久精品欧美日韩精品| 黑人操中国人逼视频| 亚洲国产欧美人成| 日日爽夜夜爽网站| 亚洲熟妇熟女久久| 啦啦啦观看免费观看视频高清| 88av欧美| 国产激情欧美一区二区| 琪琪午夜伦伦电影理论片6080| 国产精品电影一区二区三区| 一a级毛片在线观看| 亚洲国产中文字幕在线视频| 草草在线视频免费看| 成熟少妇高潮喷水视频| 国产一区二区在线观看日韩 | 97人妻精品一区二区三区麻豆| 又爽又黄无遮挡网站| 97人妻精品一区二区三区麻豆| 欧美成人免费av一区二区三区| 久久亚洲精品不卡| 中文在线观看免费www的网站 | 首页视频小说图片口味搜索| 无遮挡黄片免费观看| www国产在线视频色| 搡老岳熟女国产| 午夜福利成人在线免费观看| 国产免费av片在线观看野外av| 欧洲精品卡2卡3卡4卡5卡区| 免费看a级黄色片| 欧美日韩一级在线毛片| 女人高潮潮喷娇喘18禁视频| 国产v大片淫在线免费观看| 在线观看一区二区三区| 亚洲欧洲精品一区二区精品久久久| 十八禁网站免费在线| 99国产精品99久久久久| 欧美性猛交╳xxx乱大交人| 熟妇人妻久久中文字幕3abv| 亚洲国产精品成人综合色| 夜夜看夜夜爽夜夜摸| 麻豆成人午夜福利视频| 欧美+亚洲+日韩+国产| 桃色一区二区三区在线观看| 巨乳人妻的诱惑在线观看| 欧美中文日本在线观看视频| 国产欧美日韩一区二区三| 日韩大尺度精品在线看网址| 90打野战视频偷拍视频| 免费在线观看视频国产中文字幕亚洲| 在线观看午夜福利视频| 12—13女人毛片做爰片一| 国产麻豆成人av免费视频| av天堂在线播放| 成年女人毛片免费观看观看9| 巨乳人妻的诱惑在线观看| 一卡2卡三卡四卡精品乱码亚洲| 男女做爰动态图高潮gif福利片| 最近最新中文字幕大全电影3| 国产高清视频在线观看网站| 亚洲成av人片在线播放无| 欧美日韩乱码在线| 久久国产精品影院| 99精品久久久久人妻精品| 亚洲av成人不卡在线观看播放网| xxxwww97欧美| 18禁黄网站禁片免费观看直播| 精品久久久久久成人av| 久久天堂一区二区三区四区| 两个人看的免费小视频| 琪琪午夜伦伦电影理论片6080| 女人爽到高潮嗷嗷叫在线视频| 亚洲七黄色美女视频| 国产精品免费视频内射| 毛片女人毛片| 精品久久久久久,| 欧洲精品卡2卡3卡4卡5卡区| 五月玫瑰六月丁香| 两个人看的免费小视频| 丰满的人妻完整版| 国产av麻豆久久久久久久| 亚洲九九香蕉| 国产精品99久久99久久久不卡| 1024手机看黄色片| 男女午夜视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| 97人妻精品一区二区三区麻豆| 高潮久久久久久久久久久不卡| or卡值多少钱| 国产精品一区二区三区四区免费观看 | 色尼玛亚洲综合影院| 国产av在哪里看| 可以免费在线观看a视频的电影网站| 男女下面进入的视频免费午夜| 观看免费一级毛片| 国产精品,欧美在线| 精品高清国产在线一区| 国产高清有码在线观看视频 | 亚洲人成77777在线视频| 欧美绝顶高潮抽搐喷水| 狠狠狠狠99中文字幕| ponron亚洲| 狂野欧美激情性xxxx| 黄色女人牲交| 国产99久久九九免费精品| 国产精品自产拍在线观看55亚洲| 亚洲人成网站在线播放欧美日韩| 午夜老司机福利片| 国产不卡一卡二| 精品第一国产精品| 宅男免费午夜| 99久久99久久久精品蜜桃| 亚洲黑人精品在线| 老汉色∧v一级毛片| 久久精品国产99精品国产亚洲性色| 国产精品一及| 国产激情偷乱视频一区二区| 亚洲七黄色美女视频| 国产亚洲精品久久久久久毛片| 最近在线观看免费完整版| 午夜福利欧美成人| 日韩大尺度精品在线看网址| 丝袜美腿诱惑在线| 亚洲欧美精品综合一区二区三区| 亚洲欧美激情综合另类| 久久久久精品国产欧美久久久| 中文字幕人妻丝袜一区二区| 无限看片的www在线观看| 午夜精品久久久久久毛片777| 国产精品香港三级国产av潘金莲| 亚洲国产欧美网| 两性夫妻黄色片| 久久精品夜夜夜夜夜久久蜜豆 | 国产精品久久久人人做人人爽| 国产午夜精品论理片| 99在线视频只有这里精品首页| 久久亚洲精品不卡| 久久久久免费精品人妻一区二区| 麻豆国产av国片精品| 日本 av在线| 国产伦在线观看视频一区| 国产真实乱freesex| 亚洲激情在线av| www日本在线高清视频| 中文字幕精品亚洲无线码一区| 88av欧美| 国产精品99久久99久久久不卡| 最新在线观看一区二区三区| 床上黄色一级片| 熟女电影av网| 色综合欧美亚洲国产小说| 91成年电影在线观看|