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

    A closed-form nonlinear model for spatial Timoshenko beam flexure hinge with circular cross-section

    2019-12-28 07:54:48DanWANGJieZHANGJiangzhenGUORuiFAN
    CHINESE JOURNAL OF AERONAUTICS 2019年11期

    Dan WANG, Jie ZHANG, Jiangzhen GUO, Rui FAN

    School of Mechanical Engineering and Automation, Beihang University, Beijing 100083, China

    KEYWORDS

    Abstract Beam flexure hinges can achieve accurate motion and force control through the elastic deformation. This paper presents a nonlinear model for uniform and circular cross-section spatial beam flexure hinges which are commonly employed in compliant parallel mechanisms. The proposed beam model takes shear deformations into consideration and hence is applicable to both slender and thick beam flexure hinges.Starting from the first principles,the nonlinear strain measure is derived using beam kinematics and expressed in terms of translational displacements and rotational angles. Second-order approximation is employed in order to make the nonlinear strain within acceptable accuracy. The natural boundary conditions and nonlinear governing equations are derived in terms of rotational Euler angles and subsequently solved for combined end loads. The resulting end load-displacement model,which is compact and closed-form,is proved to be accurate for both slender and thick beam flexure using nonlinear finite element analysis. This beam model can provide designers with more design insight of the spatial beam flexure and thus will benefit the structural design and optimization of compliant manipulators.

    1. Introduction

    Compliant mechanisms, on account of their inherent advantages such as high precision,zero backlash,long operation life and design simplicity, have been widely used in a variety of applications, e.g. micro-alignment assembly,1micropositioning stage,2-5and force/acceleration sensor.6,7One of the most commonly used flexure elements is a spatial beam flexure hinge (SBFH), as illustrated in Fig.1(b), which is usually mounted between two rigid components of the host mechanism shown in Fig.1(a).

    Since a beam flexure hinge transmits motions and forces through its elastic deformation, a large body of investigations have focused on the load-displacement model of flexure beam.These studies provided the research community with several applicable methods, such as the finite-element-method(FEM) based solution,8,9the Pseudo-Rigid-Body model(PRBM),10compliance-based matrix method11,12and the constraint based spatial-beam constraint model (SBCM).13-16

    Fig.1 A compliant mechanism and its flexure element.

    The FEM based model is probably the most accurate computational method to estimate the performance of a SBFH,but it is usually denounced for its low efficiency and lack of parametric design insight. The recently reported PRBM reduces the spatial beam to a general spatial spring model which is accurate for deriving the displacement of a beam over a large range of deformation. However, this PRBM is only valid for a certain loading case and may be invalid when a six-component general load is concerned. Moreover, the PRBM fails to capture the compliance contributions of the bending displacements to the axial and twisting displacements.For example, Hao et al.15provided an analytical formulation for the spatial beam flexure,for which load-displacement relation for the two bending directions were modeled in decoupled bending planes,respectively.Therefore,this model is incapable of acquiring the intercoupling between the two bending directions.

    In Awtar’s recent work,13,14a closed-form nonlinear spatial beam model SBCM is given in matrix form. The SBCM is developed based on the Euler’s deformation assumption and can capture the coupling between the two bending directions,the load stiffening effect, the trapeze effect coupling and the kinematic and elastokinematic nonlinearities, which all result from the axial force and twisting moment. Furthermore, the SBCM is proved to be accurate when the transverse displacement is within 10% of the total beam length. In the previous work, authors17applied the SBCM to the design of a 6-DOF compliant parallel manipulator, obtaining accurate kinematic solutions as compared to the nonlinear finite element analysis results.

    However,most of the existing beam formulations and models are based on Euler’s deformation assumption valid only for long, slender spatial beams. As for the behavior of thick beams, Timoshenko beam theory is more effective since the shear distortion of the beam’s cross-section cannot be ignored anymore.For example,in the application of load bearings,the host compliant mechanism employing SBFHs needs to be designed in high stiffness while maintaining the product’s scale.Hence the related beam flexures should also be designed with higher stiffness and may be no more slender.In this case,shear distortion of the beam’s cross-section cannot be ignored anymore and most of the existing beam models will be unavailable.Timoshenko beam is a classic but not dated theory,leading to many interesting research topics: Selig and Ding18derived a screw form of Hooke’s law for Timoshenko beams;Rahmani and Pedram19studied the vibration properties of a functionally graded nano-beam with nonlocal Timoshenko beam theory; in a study of fully compliant bistable mechanisms, Chen and Ma20developed an improved planar beam constraint model based on Timoshenko beam theory, which took the shear effects into consideration and significantly enhanced the prediction accuracy over the beam model.

    The above analysis highlights the need for a general spatial beam model that can be used for a general spatial beam flexure both thick and slender. This paper starts from the basis of mechanics, and derives a general nonlinear model applicable for a thick beam(5≤slenderness ratio <20)within a small displacement range. It is noteworthy that the proposed model is accurate when the beam displacement is within 10% of beam length and rotary angle less than 0.1 rad.Moreover,the model proposed here can be utilized as more accurate deflected elements in chained spatial beam model to describe the large spatial deflection of flexible beams,as described in Chen and Bai’s work.21A further research on this application will be discussed in the future works.

    Accordingly, the remainder of this paper is organized as follows. Section 2 presents the second-order approximation of the nonlinear strain which is subsequently used in Section 3 to derive the generalized natural beam boundary conditions.In Section 4, the end load-displacement relation for a slender SBFH is derived. In Section 5, a novel end loaddisplacement model that can be used for thick SBFH is deduced based on the model in Section 4. Then the proposed models are validated and compared with finite element analysis(FEA) method in Section 6. Finally, conclusions are drawn in Section 7.

    2. Nonlinear strain formulation

    For a spatial beam suffering a general end-loading, the induced deformation can be completely defined with respect to the related cross-section that is perpendicular to the centroidal axis prior to deformation. In general, the beam deformation may be classified into five groups, including translation and rotation of the cross-section, rotation of the deformed cross-section due to the shear force, in-plane distortion of the cross-section, in-plane dilation/contraction of the cross-section and the out-of-plane warping of the crosssection.

    The analysis in this paper is performed far below the elastic limit of the beam material and no external load is applied on the lateral surfaces of the beam, and thus in-plane dilation/-contraction of the cross-section arises only due to the Poisson’s effect and is usually negligible. The SBFH concerned in this analysis is a cylindrically constructed beam with circular cross-section. Hence the out-of-plane warping of the crosssection can also be safely ignored. Translation and rotation of the cross-section, according to the Euler’s assumption, are vital to bending deformation and hence are surely kept in this analysis. The shear deformation that occurs in the bending plane of the SBFH can be safely ignored when the beam body is long, uniform and slender.14However, the SBFH in this paper is not always slender and therefore the shear deformation is taken into consideration according to Timoshenko and Goodier’s beam theory.22The in-plane distortion induced by beam torsion can be nonzero for a spatial beam flexure with combined end-loads and is also included in this analysis.

    In Fig.2,a SBFH of circular cross-section subjected to general end-loading is illustrated. Translation, rotation and the shear deformation of the cross-section are shown; however,the in-plane distortion of the cross-section, which is included in this analysis, is not shown in Fig.2 for the purpose of simplification.

    Fig.2 Spatial kinematics of SBFH.

    As shown in Fig.2, coordinate frame XYZ is the fixed global Cartesian coordinate system and point Q denotes the center of a cross-section at distance x0from the fixed end of the SBFH before deformation. When a general end-loading [FXL, FYL, FZL, MXL, MYL, MZL] is applied at the free end of the SBFH, point Q translates to Qdwith mutually perpendicular translations [u, v, w] and rotations[θd, β, α]. Point Qdlies on the deformed centroid axis of the SBFH. The three Euler angles [θd, β, α] can be used to describe the relationship between the fixed frame XYZ and the local frame XdYdZdwith a transformation matrix T.The deformed axis Xdis located at point Qdand tangent to the deformed centroid axis of the SBFH while the two perpendicular axes Ydand Zdare defined within the deformed cross-sectional plane at point Qd. For an Euler beam which is long and slender, shear deformation at cross-section Qdcan be ignored because of the equivalency of two related terms, -β(α) and dv/dx (dw/dx). However, when the SBFH is thick, according to Timoshenko’s beam theory, the shear distortion can no longer be neglected, as shown in Fig.2.In this paper, YZX-typed Euler angles are adopted for calculation simplification of the shear deformation.

    Since the load equilibrium should be established in the deformed beam configuration while all the translational deformations and end loads are completely defined in the fixed coordinate frame, the transformation matrix is employed in this analysis to relate the local coordinate frame XdYdZdand fixed frame XYZ.With the Euler angles expressed in Fig.2,the final transformation matrix T can be written as where TX,TZand TYindicate rotations around X-, Z- and Yaxis respectively.

    As the transformation matrix T varies along the deformed centroid axis of the beam, the curvature matrix of the beam can be obtained by calculating the derivative of matrix T with respect to the deformed X-coordinate variable Rnin the following form:

    A useful relation,which will benefit the simplification of the strain formulation, can be derived from the curvature matrix given in Eq. (2) and can be described as

    Obviously, nonlinear curvatures in Eq. (3) are directly derived from the transformation matrix without any approximation, and thus they are theoretically more accurate than their counterpart in previous works.13,14,23

    The position vector of an arbitrary point [x, y, z] on the undeformed SBFH can then be translated to the deformed beam configuration using the transformation matrix T in the following manner:

    where Rdand R0are the same point on the deformed and undeformed SBFH, respectively.

    With the beam deformation completely defined,the Green’s strain measure at any point within the deformed beam configuration can be retrieved as

    where du, dv and dw are small increments along the three undeformed coordinate axes X, Y and Z.

    Using Eqs. (1)-(6), the complete form expression for nonlinear strain with respect to the concerned deformations can be obtained.The full form of the nonlinear strain formulation is too complicated for both of computation and strain analysis.Hence appropriate simplifications will be applied in this analysis.The beam deformation concerned in this analysis is small,i.e.the transverse displacements are confined to be within 10%of the beam length L and the rotational displacements are less than 0.1 rad. Over this displacement range, the derivative u′ is of the order 0.01 and the v′ and w′ are of the order 0.1. Therefore,the results derived using Eqs.(1)-(6)can be simplified by truncating the high-order terms in order to approximate the result to the second order. On the other hand, the following approximation can also be applied to the results considering the small rotations:

    Thus,the final expression for nonlinear strain of the SBFH can be retrieved in a second-order form as

    The strain formulation given in Eqs.(7)-(9)is referred to as the second-order approximation in the rest of this paper. The first three terms in the axial strain, given in Eq. (7), capture the elastic extension under bending moments. The remainder in Eq.(7)reveals the influence of bending and torsion on axial strain, by the aid of beam curvatures given in Eq. (3). Shear strains given in Eq. (8) and Eq. (9) each consists of two parts:shear strain arises from pure bending moment and shear strain related to shear distortion.It is evident that shear strains given in Eqs. (8) and (9) are more accurate in solving the issue of a Timoshenko beam.Strain components εXXand εZZare associated with in-plane contraction and arise only due to the Poisson’s effect in the absence of lateral forces. When the out-ofplane warping is ignored, they will equal zeros and thus not appear in the final expression of the nonlinear strains. Moreover, the strain component εYZalso equals zeros along with the neglect of the cross-section warping.

    Consequently, the final expression for nonlinear strain,approximated to the second order, will subsequently be used to derive the general closed-form nonlinear model of the SBFH.

    3. Load equilibrium and nonlinear analysis of SBFH

    In the view of load equilibrium, the moments at any crosssection of the beam can be obtained by integrating the force acting on an infinitesimal area of the cross-section multiplied by its perpendicular distance from the centroid axis. Since the in-plane distortion and the out-of-plane warping of the cross-section are ignored, the bending moments related to the target cross-section can be simply derived as

    where E,A and I are the elastic modulus,the area and the area moment of the circular cross-section, respectively; MYdand MZdare bending moments about the deformed axes Ydand Zd, respectively.

    The twist in the deformed beam, which occurs in the deformed circular cross-section plane, can be obtained by

    where G is the shear modulus and the torsion constant J can be expressed as

    Eqs. (10)-(12) can be expressed in a matrix form as

    The moments in the deformed beam configurations can be transformed to the fixed frame using the transformation matrix T. Substituting matrix T into Eq. (13), and along with Eq. (3), the relationship between Euler angles and moments can be expressed in the fixed frame XYZ as

    Then the following relations can be obtained:

    It can be seen that Eqs.(16)and(17)are independent of θd,and hence can be solved in the absence of the first equation.Compared with Eq. (17),Eq. (16) has an extra nonlinear term which may largely complicate the solving process.However,it is evident that the nonlinear term in Eq. (16) is two orders of magnitude smaller than the remaining terms and can be dropped in this analysis. The simplified equations are listed below:

    The two expressions on the right side of Eqs. (18) and (19)are simply the effective bending moments in the Y and Z directions approximated to the second order. Furthermore, Eqs.(18) and (19) actually reveal the cross coupling effect between the two bending directions. However, it should be noted that the three equations that are already derived, i.e. Eqs. (15),(18) and (19), are still not enough to formulate the loaddisplacement relations of the SBFH which is generally a sixinput-and-six-output nonlinear model. As a result, three more equations are required to provide additional relations between the translational displacements and the general end-loading.

    By applying load equilibrium in the deformed beam configuration, moments at any point Q can be represented as

    where MX, MYand MZdenote the moments at point Q that has a distance×from the fixed end of the SBFH. The general end-loading [FXL, FYL, FZL, MXL, MYL, MZL], as depicted in Fig.2, will keep constant during the static analysis.

    Substituting Eq.(20)into Eqs.(18)and(19)and truncating the high-order terms, we can derive the following equations:

    Redundant constant values in Eq.(21)can be eliminated by differentiating Eq.(21)with respect to×and the results can be written in matrix form:

    where α′′and β′′are the second derivatives of α and β respectively.

    Eq. (22) has a pair of coupled but linear equations. The coupling arises from the axial twisting mXl, which may contribute additional loads in two bending directions.In addition,this equation has two desired unknowns (β and α) and two undesired unknowns(w′and v′),which make it mathematically unsolvable. Two supernumerary relations are still needed to relate the desired unknowns and undesired unknowns. On the one hand, when the spatial beam is long and slender,desired unknowns (β and α) are approximately equal to the two undesired ones (w′ and v′) within a small displacement range.14,18On the other hand, when the spatial beam is short and thick, a certain relationship can also be found between the desired and undesired unknowns according to Timoshenko’s beam theory.17,24

    It is noteworthy that Eq.(22)provides a constraint relation which can be extended to governing equations for both long slender beams and short thick beams within a small displacement range.

    4. End load-displacement relations of a slender SBFH

    An Euler beam refers to the long,slender and uniform beam in this analysis. Within a small displacement range and the second-order approximation, the differences between the two Euler angles, β and α, and the two derivatives, -w′ and v′,are quite small, as shown in Fig.3. Therefore, the two derivatives can be safely replaced by the two Euler angles based on Eq. (22), which can be rewritten as

    Thus the beam governing equations can be obtained as

    Fig.3 Spatial kinematics of beam deformation: slender beam.

    The general solutions for the two Euler angles β and α can be derived by solving Eq. (24) and given by

    Then,the transverse displacements w and v can be obtained by integrating Eq. (25) and expressed in a simplified form:

    where ci(i=1-8) are the constants of integration.

    The geometric boundary conditions to Eq. (24) can be described as

    Referring to Eqs. (20), (23), (25) and (26), end loaddisplacement relations at the free end of the slender SBFH can be determined in a compact form as shown in Eq. (28).

    In this 4×4 matrix K, components are actually transcendental expressions and too complex to be given in full form.Therefore, Taylor series expansion is applied here to matrix K in terms of mXLand fXLin order to retrieve a concise form.The final expression of the transcendental components of matrix K, dropping third- and higher-order terms, can then be given by

    The results illustrated in Eqs.(28)and(29)capture essential stiffness characteristics of the SBFH under general end-loading and within the displacement range of interest. The stiffness matrix K will become the linear elastic stiffness associated with the four transverse bending displacements when the twisting moment mXLand axial force fXLare zeros. Other inconstant terms in Eq.(29)capture the load-stiffening in these directions in the presence of mXLand fXL. Furthermore, the approximated stiffness matrix K can produce an acceptable accuracy of less than 1% error while making a notable decrement on computation burden and providing more design insight.

    In order to derive the load-displacement relation along X direction, the displacement component ulshould be expressed in terms of bending displacements, i.e., β, α, w and v. As shown in Fig.3, the axial force at point Q along the deformed centroid axis equals

    Since the translational and rotational displacement in this analysis are both small, axial force in Eq. (30) has ignorable difference with the force along X-axis of the fixed frame.Thus the normalized axial displacement in the fixed coordinate frame XYZ can then be obtained by integrating Eq. (30) with respect to x

    where kSLis the normalized elastic stiffness and has the following form:

    By substituting Eqs.(3),(25)and(26)into Eq.(31),the displacement in X direction can be expressed in transcendental form. For the sake of brevity, the derived transcendental expression is also expanded in terms of mXLand fXL.By truncating the second-and higher-order terms,the final expression for total axial displacement of the beam’s free end can be stated in the following compact form:

    where

    The first term in Eq. (32) captures the normalized elastic extension of the spatial beam while kSLis the normalized elastic stiffness in the X-direction.Constant terms of series expansion in Eq. (33) capture the purely kinematic contributions of the transverse forces and bending moments to the axial stretching,while the remaining inconstant terms reveal the elastokinematic coupling between the transverse displacements and the axial extension.

    The twisting displacement in the deformed beam configuration is expressed in a derivative form in Eqs. (15)-(17). The beam twisting θlin the fixed frame XYZ can be expressed as25

    By substituting Eqs. (15)-(17) and Eq. (26) into Eq. (34),the twist angle at the free end of the SBFH can be expressed in transcendental functions of mXLand fXL. To allow further design insight and computational simplification, all transcendental functions derived should be expanded in terms of mXLand fXL. In the load and displacement ranges of interest,acceptable approximations, which will incur an error of less than 1%, can be derived by truncating the second- and higher-order terms. The resulting simplified expression for the twist angle of the SBFH can be given in the following form:

    where

    The first term in Eq. (35) is only dependent on the twisting moment and therefore a purely elastic component of θl. Constant terms in c41, c32, c23and c14are free of any end load and only depend on the transverse displacements and bending angles, thus capturing the purely kinematic contributions of the transverse loads to the twist angle. The remaining inconstant terms with respect to mXLand fXLare load-dependent and will incur the elastokinematic coupling between the transverse displacements and the beam twisting.

    In Eqs. (28), (32) and (35), the end load-displacement relations for a long, slender spatial beam with combined endloading are presented, where the shear deformation of the cross-section are ignored according to Euler’s assumptions.However, slender SBFHs will cause superfluous compliances which are not always expected in flexure hinge design. For instance, a beam hinge designed for higher load bearing may have a slenderness ratio less than 10 in order to obtain extra stiffness enhancement, in which case the shear deformation,as captured by Timoshenko and Goodier,22will become prominent and cannot be ignored. Accordingly, loaddisplacement relations for a stubby beam flexure hinge will be subsequently developed based on the analysis above.

    5. End load-displacement relations of a thick SBFH

    Fig.4 Spatial kinematics of beam deformation: thick beam.

    In the deformed configuration of a stubby SBFH as shown in Fig.4,shear deformations that are caused by transverse forces are non-ignorable.Applying load equilibrium at the free end of the deformed beam configuration, shear distortions can be approximately related to transverse forces in the fixed coordinate frame XYZ by

    where μ is the shear coefficient for a solid circular crosssectional spatial beam, as expressed by26

    with ν being the Poisson’s ratio.

    By substituting Eq. (37) into Eq. (22), natural boundary conditions for the thick SBFH can be rewritten as

    Since terms on the right side of Eq. (38) are all constant,governing equations for the thick SBFH are identical with those given in Eq. (24), and meanwhile the general solutions for the two Euler angles β and α are also the same as those shown in Eq. (26).

    Using Eq.(37)and Eq.(25),the transverse displacements w and v are given by

    Referring to Eqs. (20), (26), (27) and (39), the end loaddisplacement relations are obtained at the free end of a thick SBFH, as shown in Eq. (40).

    where

    The expressions in Eq. (41), which are derived by Taylor series expansion, indeed give more insight into the effects of the axial load and twisting moment on the four transverse displacements. It can be seen that terms relating to mXLand fXLbecome zeros and the resulting stiffness matrix turns to the purely linear stiffness matrix in the absence of axial load and twisting moment. In each bending direction, there exists a prominent load-stiffening component associated with the first order of axial load fXLwhile components with the secondorder fXLterms cast relative small influences on transverse displacements. Furthermore, the stiffness matrix components relating to third- or higher-order mXLand fXLterms always have very small coefficient and thus are removed from Eq.(41).

    Next, the displacement along X-axis can be divided into two parts. The first part indicates the purely elastic deformation caused by the elastic stretching of the SBFH arc-length due to the axial force. The second part is induced by the two transverse displacement, indicating the geometric constraint of beam-arc length conservation. Substituting Eqs. (25) and(29) into Eq. (31), we can state the total axial displacement at the free end of SBFH in the following compact form:

    where

    Note that matrix components given in Eq. (43) can also be divided into two groups. Terms in the first group accord with their corresponding ones in Eq. (33), showing the kinematic and elastokinematic contributions of the axial load and twisting moment to the displacement in the X-direction. It can be seen that the coefficients of terms fXLand mXL, except kSL,are usually small, which indicates that the two transverse and the two rotational displacements exert relative small effect on the axial displacement.Other terms in Eq.(33),which relate to constant ψ, provide essential corrections to compensate the displacements caused by shear distortions.Generally speaking,constant ψ is quite small for a slender Euler beam and negligible in the slender beam model.However,while the beam is no longer slender,i.e.the slenderness ratio is less than 20,the constant ψ can no more be neglected and the terms with respect to constant ψ will become prominent.

    Purely kinematic terms in Eqs.(33)and(43)generally dominate the axial displacements in terms of magnitude, but they make no contribution to the axial stiffness. On the contrary,the elastokinematic terms, although with small magnitude,have prominent influences on the axial compliance. The correction terms in Eq.(43),which are usually small in magnitude,play a vital role in keeping Eq. (42) accurate while solving the axial displacement for a thick SBFH. But higher-order terms with respect to mXL, fXLand ψ are dropped due to their insignificant contribution.

    Finally,solutions given in Eqs.(3),(25)and(39)are substituted into Eq. (31), which upon integration provides the solution for the twisting of the stubby SBFH. Similarly, the solution contains a 4×4 matrix and each nonzero element is a transcendental function of the axial load fXLand twisting moment mXL. After applying the Tayler series expansions to the matrix element and truncating the high-order terms, the final expression for the twisting displacement of the thick SBFH can be written as

    where

    The first term in Eq. (44) is the purely elastic component.The kinematic and elastokinematic terms in Eq. (44) are exactly the same as their counterparts in Eq. (36), while the correction terms are complicated in form and usually have relative small magnitude. It can be seen in Eq. (45) that elastokinematic terms with respect to mXLbecome the most important contributors while the purely kinematic terms are not as important as those in Eq. (43). Elastokinematic terms with respect to axial load fXLcan also be found in Eq. (45),which reflect the effect of axial forces on torsion. In general,this term has very small effect on beam twisting, but can become significant under high axial load and therefore is not ignorable.

    Consequently,a closed-form load-displacement model for a thick SBFH can be established using Eqs. (40), (42) and (44).Compared to the load-displacement model for slender SBFH in Section 4, this model is improved by taking into consideration of the shear distortion which is usually concerned in the Timoshenko’s beam theory. The two load-displacement models proposed above are similar in the expression form and identical in execution. In addition, the model developed for the thick SBFH can be easily transformed into the model for the slender SBFH by setting constant ψ to zero.In other words,load-displacement model established here can be used as a general beam model for all circular cross-sectional spatial beam flexure hinges,whether it is slender or thick.This model is subsequently referred to as the spatial Timoshenko beam model(STBM).

    6. Model validation via finite element analysis

    The results derived in the previous sections are validated via nonlinear FEA software ABAQUS. Simulations are performed in two groups with designated beam parameters given in Table 1.In these simulations,established load-displacement solutions to both slender and thick SBFHs are verified by applying designed combined loads to the free end of the beams. In these simulations, the 2-node cubic beam elementB33 is adopted to model the slender SBFH, while the 2-node linear beam elements B31 are chosen for simulations of the thick SBFH. All simulations are performed in ABAQUS software with large deformation switch turned on.

    Table 1 Geometric and material parameters of SBFHs.

    Table 2 Simulation load spectrums.

    In order to keep the analysis within the displacement range of interest, combined loads applied to the end of tested SBFH are carefully designed in magnitude as well as load history.Since there are large differences existing between structural compliance of the two sorts of beam flexures, end loads are also quite different in magnitude, as shown in Table 2.

    For each of the two loading histories,end displacements of the tested SBFHs are firstly obtained numerically by solving the STBM and then compared with the results that are obtained by ABAQUS software. Comparison results are plotted in Figs.5 and 6 with respect to simulation time, where the prefix ‘‘Num” and ‘‘FEA” denote results derived from the STBM and simulation respectively in this paper.

    The comparison results given in Figs.5 and 6 illustrate that the presented load-displacement models are good approximations to the FEA finite element model within the interested displacement range. Another noteworthy fact is that the slender SBFH and the thick SBFH exhibit similar behavior under the given load histories.

    As shown in Figs. 5 and 6, increments of the transverse forces FYLand FZLhave little effect on the axial displacement,along which direction the SBFH has the maximum structural stiffness. This is also accordance with the small coefficients that are given in Eqs. (33) and (43). The twisting moment,which arises at the fourth second of the load history,has a relative small impact on the transverse displacements. In the last second, the bending moment MZLleads to a prominent increment on the displacement v, which is mostly caused by the purely kinematic terms of large magnitude listed in Eqs. (29)and (41).

    The accuracy of the proposed models can be proved by error distribution plotted in Fig.7 and the maximum relative discrepancies are listed in Tables 3 and 4.

    For the slender SBFH, absolute transverse errors along Yand Z-axis have obvious change during the load history,while the corresponding discrepancy of other four displacements are all small in magnitude and fluctuating around zero. Fig.7(a)shows that: from the first second, the error along Y-axis (v)begins to climb; the error along Z-axis (w) also starts to increase from the next second;these two errors reach a plateau from the fourth to fifth second and become maximum at the last second, but the corresponding relative errors are still within 1.2%,as shown in Table 3 and Fig.7(a).The maximum relative error of the other four displacements lies in the rotational displacement about the Y-axis,where the extreme value is 4.3%.

    Fig.5 Simulation results of slender SBFH.

    Fig.6 Simulation results of thick SBFH.

    Fig.7 Discrepancy between numerical and FEA results.

    Table 3 Maximum relative discrepancies between FEA and STBM: slender beam.

    Table 4 Maximum relative discrepancies between FEA and STBM: thick beam.

    For the thick SBFH, magnitudes of the errors along two transverse axes and the rotational error about Y-axis are relatively large and cannot be neglected. The following trends can be observed from Fig.7(b): the error along Y-axis responds from the first second and reaches the local maximum value (in magnitude) at second two; other discrepancies start to enhance during second two and second three; all the displacement errors turn into a stable state from the fourth and fifth second and the transverse errors along Y- and Z-axis reduce to a low level; all the errors experience dramatic increase and reach the maximum in the last load step. As shown in Table 4, the two maximum transverse discrepancies are 0.04 mm and 0.033 mm, while the corresponding relative discrepancies are both 1.2%. The relative discrepancy about Y-axis is 4.6%, which is much larger than the other five components; however it is still accurate to be within 5% with respect to FEA and is acceptable in practice.

    It can be seen from Fig.7 that large deviations occur when the two transverse forces start to act and plateau at the end of the third second. In physical terms, it can be deduced that the large deviations mostly stem from approximation errors of the coupling effect between the two transverse axes.With the same reason,the two transverse displacement errors start to increase with the introduction of the bending moment about Z-axis in the sixth second.

    7. Conclusions

    In this paper, a closed-form load-displacement model for circular cross-sectional spatial beam is developed. Compared with the existing beam models, the proposed one takes shear distortion of the beam cross-section into consideration and hence can be used to solve problems related to thick beams(5≤slenderness ratio <20).

    Starting from the first principles, the nonlinear strain measure is derived using beam kinematics and expressed in terms of translational displacements and rotational angles. Natural boundary conditions and governing equations for a constraint spatial beam are then formulated. Based on the Euler’s deformation assumption,the load-displacement model for a slender SBFH is subsequently derived in terms of axial load and twisting moment. The proposed slender-beam model can capture most of the geometric nonlinearities of a SBFH, including load-stiffening in the bending directions, coupling between the bending directions, the kinematic and elastokinematic component in the axial displacement and twisting angle, and the coupling between the axial and torsional directions. However, the slender-beam model is only available to slender beams for which shear distortion in the beam cross-section can be neglected.

    For the sake of generality,we develop a novel model STBM by introducing a shear distortion term into the slender-beam model. The proposed model STBM is developed based on Timoshenko’s beam theory and featured by the appropriate engineering approximations. This makes the final STBM suitable for both the thick beam and the slender beam while reducing the computational load to a manageable level.Finally,the STBM is used to solve the load-displacement relations for a slender beam and a thick beam. Results from the STBM are fully compared with those derived using nonlinear FEA and the comparison shows that the proposed STBM is always accurate and within 5% in any displacement direction of the two beams.

    The main contribution of this paper exists in the proposition of the STBM which will facilitate the structural design and optimization of compliant mechanisms that employ SBFHs. The STBM can capture most of the main geometric nonlinearities that affect the mechanical behavior of a spatial beam that is either slender or thick, and hence the analytical model will benefit the designer with more parametric design insight. In addition, the derivation process of the STBM can be extended and used in modeling a beam with any other bisymmetrical cross-section with two identical principal bending moments of area.

    Acknowledgement

    This study was supported by the National Natural Science Foundation of China (No. 51305013).

    伊人亚洲综合成人网| av福利片在线| 青春草亚洲视频在线观看| 午夜激情av网站| 国产一区二区三区av在线| 国产精品久久久久久精品古装| 满18在线观看网站| 韩国高清视频一区二区三区| 国产成人免费无遮挡视频| 法律面前人人平等表现在哪些方面 | 美国免费a级毛片| 日日夜夜操网爽| 看免费av毛片| 日本av手机在线免费观看| 丰满人妻熟妇乱又伦精品不卡| 久久久久视频综合| 国产国语露脸激情在线看| 成人手机av| 啪啪无遮挡十八禁网站| 国产一卡二卡三卡精品| 纯流量卡能插随身wifi吗| 淫妇啪啪啪对白视频 | 视频在线观看一区二区三区| 啦啦啦视频在线资源免费观看| 亚洲一码二码三码区别大吗| 国产不卡av网站在线观看| 成年人午夜在线观看视频| 精品少妇久久久久久888优播| 亚洲第一欧美日韩一区二区三区 | 国产xxxxx性猛交| 午夜福利视频精品| av有码第一页| 色婷婷久久久亚洲欧美| 精品久久久久久电影网| 超色免费av| 欧美日本中文国产一区发布| 久久亚洲精品不卡| 女人被躁到高潮嗷嗷叫费观| 日本欧美视频一区| 欧美在线一区亚洲| av一本久久久久| 9热在线视频观看99| 十八禁高潮呻吟视频| 国产在线免费精品| 亚洲国产av影院在线观看| 中文字幕av电影在线播放| 热re99久久精品国产66热6| a级片在线免费高清观看视频| 久久久久久久久久久久大奶| 搡老乐熟女国产| av免费在线观看网站| 欧美97在线视频| 久久99一区二区三区| 一本色道久久久久久精品综合| 王馨瑶露胸无遮挡在线观看| av在线app专区| 美女脱内裤让男人舔精品视频| av片东京热男人的天堂| av电影中文网址| 国产在线一区二区三区精| 日韩熟女老妇一区二区性免费视频| 成年人黄色毛片网站| 国产成人av激情在线播放| 777久久人妻少妇嫩草av网站| 日韩 欧美 亚洲 中文字幕| 亚洲欧美一区二区三区黑人| 国产精品 欧美亚洲| 美女国产高潮福利片在线看| 久久亚洲精品不卡| 日本av免费视频播放| 黄片播放在线免费| 性色av一级| 国产av一区二区精品久久| 国产亚洲精品第一综合不卡| 亚洲国产av新网站| 叶爱在线成人免费视频播放| 老汉色av国产亚洲站长工具| 黄色 视频免费看| 两个人免费观看高清视频| 国产伦人伦偷精品视频| 精品一区二区三区四区五区乱码| 99国产精品99久久久久| 不卡一级毛片| 另类亚洲欧美激情| 国产精品九九99| 亚洲一区二区三区欧美精品| 夜夜夜夜夜久久久久| 中文字幕高清在线视频| 狂野欧美激情性bbbbbb| 国产又爽黄色视频| 91麻豆精品激情在线观看国产 | 啦啦啦 在线观看视频| 高清黄色对白视频在线免费看| 成人手机av| 欧美另类亚洲清纯唯美| 真人做人爱边吃奶动态| av免费在线观看网站| 亚洲成人免费av在线播放| 男女床上黄色一级片免费看| 男女边摸边吃奶| 一级,二级,三级黄色视频| 大型av网站在线播放| 午夜福利在线观看吧| 国产区一区二久久| 亚洲av电影在线观看一区二区三区| 国产免费一区二区三区四区乱码| 国产日韩欧美视频二区| 国产精品麻豆人妻色哟哟久久| 国产在线视频一区二区| 大型av网站在线播放| 亚洲国产av新网站| 天堂8中文在线网| 国产精品99久久99久久久不卡| 嫩草影视91久久| 国产91精品成人一区二区三区 | 一级片免费观看大全| 国产极品粉嫩免费观看在线| 丰满人妻熟妇乱又伦精品不卡| 韩国高清视频一区二区三区| 天天添夜夜摸| 高清av免费在线| 国产精品99久久99久久久不卡| 黄片小视频在线播放| 国产欧美日韩一区二区精品| 中文字幕高清在线视频| 精品人妻在线不人妻| 午夜福利在线观看吧| 成人国语在线视频| 日本猛色少妇xxxxx猛交久久| 久久av网站| 不卡一级毛片| 国产精品久久久久久精品古装| 亚洲欧美色中文字幕在线| 欧美日本中文国产一区发布| 国产精品久久久久久人妻精品电影 | 999精品在线视频| 男女午夜视频在线观看| 国产一卡二卡三卡精品| 亚洲精品中文字幕一二三四区 | 亚洲成国产人片在线观看| 一级毛片女人18水好多| 欧美日韩精品网址| 午夜福利一区二区在线看| 国产av又大| 91大片在线观看| 久久综合国产亚洲精品| 亚洲久久久国产精品| 热99久久久久精品小说推荐| 日韩视频在线欧美| www.自偷自拍.com| 悠悠久久av| 精品国产一区二区三区四区第35| 美女高潮到喷水免费观看| 十八禁人妻一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 看免费av毛片| 亚洲欧美激情在线| 日韩精品免费视频一区二区三区| 99国产综合亚洲精品| 午夜精品久久久久久毛片777| 伦理电影免费视频| 午夜激情av网站| 天天影视国产精品| 日韩视频一区二区在线观看| 十分钟在线观看高清视频www| 不卡av一区二区三区| 少妇被粗大的猛进出69影院| 后天国语完整版免费观看| 99久久综合免费| 日韩视频在线欧美| 国产精品秋霞免费鲁丝片| 无限看片的www在线观看| 他把我摸到了高潮在线观看 | 亚洲av男天堂| 国产日韩欧美视频二区| 久久女婷五月综合色啪小说| 国产亚洲午夜精品一区二区久久| 99国产极品粉嫩在线观看| 国产成人精品无人区| 久久国产精品人妻蜜桃| 天天影视国产精品| av福利片在线| 黑人巨大精品欧美一区二区mp4| 捣出白浆h1v1| 欧美一级毛片孕妇| 久久久久久人人人人人| 亚洲精品日韩在线中文字幕| 91精品伊人久久大香线蕉| 久久精品熟女亚洲av麻豆精品| 一区二区三区精品91| 在线观看免费高清a一片| 久久国产精品男人的天堂亚洲| 久久精品成人免费网站| 美国免费a级毛片| a级片在线免费高清观看视频| 菩萨蛮人人尽说江南好唐韦庄| 777久久人妻少妇嫩草av网站| 19禁男女啪啪无遮挡网站| 一本色道久久久久久精品综合| 亚洲三区欧美一区| 亚洲欧美一区二区三区黑人| 中文字幕高清在线视频| 91麻豆精品激情在线观看国产 | 久久精品aⅴ一区二区三区四区| 色精品久久人妻99蜜桃| 精品国产国语对白av| 日韩欧美一区二区三区在线观看 | 国产精品.久久久| 亚洲精品国产一区二区精华液| 亚洲精品国产区一区二| 99国产精品一区二区蜜桃av | av一本久久久久| www.自偷自拍.com| 中文欧美无线码| 女人爽到高潮嗷嗷叫在线视频| 少妇被粗大的猛进出69影院| av又黄又爽大尺度在线免费看| 国产精品久久久久久精品古装| 日本a在线网址| 久久国产精品人妻蜜桃| avwww免费| 99香蕉大伊视频| 女性被躁到高潮视频| 久久久久精品人妻al黑| 成人三级做爰电影| 极品少妇高潮喷水抽搐| 麻豆国产av国片精品| 欧美国产精品va在线观看不卡| 深夜精品福利| 欧美中文综合在线视频| 亚洲精品一区蜜桃| 午夜精品久久久久久毛片777| 亚洲国产精品成人久久小说| 人人妻人人澡人人爽人人夜夜| 久久久精品94久久精品| 天天躁狠狠躁夜夜躁狠狠躁| 99热国产这里只有精品6| 亚洲精品一卡2卡三卡4卡5卡 | 黄色a级毛片大全视频| 少妇精品久久久久久久| 国产欧美日韩综合在线一区二区| 午夜两性在线视频| 永久免费av网站大全| 亚洲精品av麻豆狂野| 欧美精品高潮呻吟av久久| 国产熟女午夜一区二区三区| 亚洲精品第二区| 肉色欧美久久久久久久蜜桃| 国产免费一区二区三区四区乱码| 19禁男女啪啪无遮挡网站| √禁漫天堂资源中文www| 亚洲 国产 在线| 十八禁人妻一区二区| 不卡av一区二区三区| 五月开心婷婷网| 国产亚洲午夜精品一区二区久久| 国产xxxxx性猛交| 午夜精品久久久久久毛片777| 不卡av一区二区三区| 国产主播在线观看一区二区| 欧美在线黄色| 女性被躁到高潮视频| 免费日韩欧美在线观看| 免费在线观看影片大全网站| 精品人妻熟女毛片av久久网站| 欧美日韩福利视频一区二区| 日日夜夜操网爽| 亚洲av成人不卡在线观看播放网 | 免费久久久久久久精品成人欧美视频| a级毛片黄视频| 欧美午夜高清在线| 欧美成狂野欧美在线观看| 午夜福利视频精品| 一个人免费看片子| 高清黄色对白视频在线免费看| 国产人伦9x9x在线观看| 国产高清视频在线播放一区 | 精品一区二区三区av网在线观看 | 超碰97精品在线观看| 国产不卡av网站在线观看| 免费高清在线观看视频在线观看| 久久青草综合色| 亚洲三区欧美一区| 啦啦啦免费观看视频1| 91av网站免费观看| 亚洲精品乱久久久久久| 国产精品 欧美亚洲| 超碰成人久久| 日韩视频一区二区在线观看| 日日摸夜夜添夜夜添小说| 他把我摸到了高潮在线观看 | 免费在线观看视频国产中文字幕亚洲 | 国产一区二区在线观看av| 韩国高清视频一区二区三区| 汤姆久久久久久久影院中文字幕| 水蜜桃什么品种好| 母亲3免费完整高清在线观看| 精品卡一卡二卡四卡免费| 一二三四社区在线视频社区8| 亚洲国产中文字幕在线视频| 国产在线观看jvid| 亚洲免费av在线视频| 少妇粗大呻吟视频| 99香蕉大伊视频| 韩国高清视频一区二区三区| 免费人妻精品一区二区三区视频| 侵犯人妻中文字幕一二三四区| 91成年电影在线观看| 91麻豆av在线| 在线观看www视频免费| 日韩视频在线欧美| 精品亚洲成a人片在线观看| 亚洲视频免费观看视频| 日韩一卡2卡3卡4卡2021年| 久久久久网色| 国产一区二区激情短视频 | av在线app专区| 国产精品偷伦视频观看了| av电影中文网址| 久久精品熟女亚洲av麻豆精品| 亚洲精品日韩在线中文字幕| 窝窝影院91人妻| av在线老鸭窝| 欧美成狂野欧美在线观看| 又黄又粗又硬又大视频| 国产三级黄色录像| 少妇的丰满在线观看| 麻豆乱淫一区二区| 视频区欧美日本亚洲| 中文字幕最新亚洲高清| 天天躁狠狠躁夜夜躁狠狠躁| 最黄视频免费看| 国产精品1区2区在线观看. | 啦啦啦视频在线资源免费观看| 成年人午夜在线观看视频| 欧美激情久久久久久爽电影 | 动漫黄色视频在线观看| 国产精品一区二区在线观看99| 欧美精品av麻豆av| 亚洲全国av大片| 天天添夜夜摸| 黑人猛操日本美女一级片| 成人国产一区最新在线观看| 国产精品 国内视频| 青青草视频在线视频观看| 嫩草影视91久久| 啪啪无遮挡十八禁网站| 午夜激情av网站| 欧美xxⅹ黑人| 欧美精品啪啪一区二区三区 | av在线播放精品| www.999成人在线观看| 久久久久久久久免费视频了| 深夜精品福利| 啦啦啦中文免费视频观看日本| 亚洲国产精品999| 中国国产av一级| 久久精品国产亚洲av香蕉五月 | 国产麻豆69| 中文字幕制服av| 欧美日韩福利视频一区二区| 国产深夜福利视频在线观看| 久久这里只有精品19| 韩国精品一区二区三区| 久久这里只有精品19| 韩国精品一区二区三区| 夜夜夜夜夜久久久久| 极品少妇高潮喷水抽搐| 一本—道久久a久久精品蜜桃钙片| 老熟女久久久| 亚洲精品乱久久久久久| 男女无遮挡免费网站观看| av又黄又爽大尺度在线免费看| 午夜影院在线不卡| 亚洲欧洲日产国产| 97在线人人人人妻| 国内毛片毛片毛片毛片毛片| 岛国毛片在线播放| 精品一区二区三区四区五区乱码| 黑人操中国人逼视频| 久久狼人影院| 满18在线观看网站| 亚洲精品日韩在线中文字幕| 亚洲欧美一区二区三区久久| 美女高潮喷水抽搐中文字幕| 中文精品一卡2卡3卡4更新| 国产一级毛片在线| 久久久久国内视频| 色播在线永久视频| 亚洲av美国av| 中文字幕色久视频| 妹子高潮喷水视频| av网站免费在线观看视频| 亚洲av片天天在线观看| 免费在线观看完整版高清| 一本—道久久a久久精品蜜桃钙片| 侵犯人妻中文字幕一二三四区| 亚洲av美国av| 菩萨蛮人人尽说江南好唐韦庄| 成年人午夜在线观看视频| 久热爱精品视频在线9| videosex国产| 欧美中文综合在线视频| 狠狠精品人妻久久久久久综合| 久久人妻熟女aⅴ| 熟女少妇亚洲综合色aaa.| 黄色 视频免费看| 欧美在线一区亚洲| 精品高清国产在线一区| 亚洲精品一二三| 自拍欧美九色日韩亚洲蝌蚪91| 成人国产一区最新在线观看| 亚洲精品日韩在线中文字幕| 又黄又粗又硬又大视频| 美女高潮到喷水免费观看| 久久久精品94久久精品| 99久久人妻综合| 亚洲精品久久久久久婷婷小说| 亚洲欧美清纯卡通| 精品少妇黑人巨大在线播放| 18禁黄网站禁片午夜丰满| 97人妻天天添夜夜摸| 天天添夜夜摸| 啦啦啦视频在线资源免费观看| 国产欧美日韩一区二区三 | 中文欧美无线码| 亚洲精品国产区一区二| 午夜老司机福利片| 女人高潮潮喷娇喘18禁视频| 少妇的丰满在线观看| 成年动漫av网址| www.熟女人妻精品国产| 国产一区二区三区综合在线观看| 国产av国产精品国产| 久久国产精品大桥未久av| av电影中文网址| 亚洲黑人精品在线| 亚洲精品国产区一区二| 久久久精品国产亚洲av高清涩受| 啦啦啦免费观看视频1| 丰满人妻熟妇乱又伦精品不卡| 一本一本久久a久久精品综合妖精| 久久国产精品男人的天堂亚洲| 亚洲欧美成人综合另类久久久| 亚洲精品一区蜜桃| 大码成人一级视频| 另类亚洲欧美激情| 看免费av毛片| 免费观看人在逋| 日韩 亚洲 欧美在线| 亚洲一区二区三区欧美精品| 老熟妇仑乱视频hdxx| 69精品国产乱码久久久| 亚洲国产成人一精品久久久| 视频区图区小说| 欧美精品高潮呻吟av久久| 久久久精品94久久精品| 岛国毛片在线播放| 亚洲va日本ⅴa欧美va伊人久久 | 色婷婷久久久亚洲欧美| 久久久久国内视频| 午夜精品国产一区二区电影| 国产成人免费观看mmmm| 考比视频在线观看| 黄色毛片三级朝国网站| 男女免费视频国产| 一边摸一边做爽爽视频免费| 欧美日韩福利视频一区二区| 亚洲第一av免费看| 黄片播放在线免费| 丝袜美腿诱惑在线| 日韩 亚洲 欧美在线| 一区在线观看完整版| videos熟女内射| 一级片'在线观看视频| 另类亚洲欧美激情| 99国产精品99久久久久| 国产成人精品久久二区二区91| 亚洲三区欧美一区| 老司机午夜十八禁免费视频| 久久狼人影院| 日本av手机在线免费观看| 精品国产一区二区三区久久久樱花| 五月天丁香电影| 日韩大码丰满熟妇| 乱人伦中国视频| 婷婷丁香在线五月| www.999成人在线观看| 丰满少妇做爰视频| a级毛片在线看网站| 欧美 日韩 精品 国产| av片东京热男人的天堂| 久久久欧美国产精品| av在线播放精品| 黄频高清免费视频| 一边摸一边做爽爽视频免费| 国产一区二区在线观看av| 久久亚洲精品不卡| 久久青草综合色| 国产伦人伦偷精品视频| 在线观看www视频免费| 亚洲av男天堂| 亚洲精品一二三| 亚洲av片天天在线观看| 一级片免费观看大全| 久久久久久人人人人人| 久久久国产一区二区| 黄片小视频在线播放| 99九九在线精品视频| 午夜福利在线观看吧| 91字幕亚洲| 777米奇影视久久| 蜜桃在线观看..| 精品卡一卡二卡四卡免费| 国产精品 国内视频| 亚洲精品国产av蜜桃| 国产亚洲欧美精品永久| 精品少妇黑人巨大在线播放| 99国产精品一区二区三区| 18在线观看网站| 亚洲精品一区蜜桃| 中文字幕最新亚洲高清| 俄罗斯特黄特色一大片| 狂野欧美激情性bbbbbb| 国产精品秋霞免费鲁丝片| 日韩欧美一区视频在线观看| 国产av国产精品国产| 男女下面插进去视频免费观看| 老汉色∧v一级毛片| 欧美黄色淫秽网站| 宅男免费午夜| 夜夜骑夜夜射夜夜干| 最新的欧美精品一区二区| 日本五十路高清| 中文字幕最新亚洲高清| 亚洲九九香蕉| 另类精品久久| 爱豆传媒免费全集在线观看| 十八禁网站免费在线| 色婷婷久久久亚洲欧美| 久久香蕉激情| 啦啦啦免费观看视频1| 国产99久久九九免费精品| 天堂中文最新版在线下载| 韩国高清视频一区二区三区| 中文字幕色久视频| 久久久水蜜桃国产精品网| 色老头精品视频在线观看| 欧美日韩一级在线毛片| 美女中出高潮动态图| 黄色毛片三级朝国网站| 午夜福利在线观看吧| av天堂在线播放| 久久影院123| 午夜福利免费观看在线| 国产精品免费视频内射| av片东京热男人的天堂| 欧美黄色片欧美黄色片| 免费日韩欧美在线观看| 免费在线观看黄色视频的| 国产av一区二区精品久久| 精品国产国语对白av| 在线观看免费日韩欧美大片| 天天影视国产精品| 啦啦啦视频在线资源免费观看| 一级片免费观看大全| 亚洲伊人色综图| 9色porny在线观看| 亚洲国产精品一区二区三区在线| 免费在线观看影片大全网站| 欧美中文综合在线视频| 老司机影院毛片| 精品少妇一区二区三区视频日本电影| 国产97色在线日韩免费| 99国产精品免费福利视频| 老司机深夜福利视频在线观看 | 午夜影院在线不卡| 国精品久久久久久国模美| av超薄肉色丝袜交足视频| 97在线人人人人妻| 色视频在线一区二区三区| cao死你这个sao货| 男人舔女人的私密视频| 欧美大码av| 久久人妻福利社区极品人妻图片| 亚洲中文字幕日韩| 爱豆传媒免费全集在线观看| 国产在线免费精品| 99热全是精品| 婷婷丁香在线五月| 90打野战视频偷拍视频| 一边摸一边做爽爽视频免费| 成年女人毛片免费观看观看9 | 成年女人毛片免费观看观看9 | 免费高清在线观看日韩| 中文字幕精品免费在线观看视频| 永久免费av网站大全| 国产精品一二三区在线看| 国产欧美日韩综合在线一区二区| 精品亚洲成a人片在线观看| 亚洲国产欧美在线一区| 日韩欧美一区视频在线观看| 午夜免费鲁丝| 俄罗斯特黄特色一大片| 狂野欧美激情性bbbbbb| 免费观看a级毛片全部| 丁香六月天网| 男女国产视频网站| 黄色视频不卡| 午夜福利影视在线免费观看| 一级片免费观看大全| 性色av一级| 99久久精品国产亚洲精品| 精品人妻在线不人妻| 啦啦啦中文免费视频观看日本|