• <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).

    午夜精品一区二区三区免费看| 我的老师免费观看完整版| 亚洲真实伦在线观看| 777米奇影视久久| 在线观看免费高清a一片| 少妇 在线观看| 联通29元200g的流量卡| 日本免费在线观看一区| 亚洲av成人精品一二三区| 精品少妇久久久久久888优播| 简卡轻食公司| 少妇的逼水好多| 熟女电影av网| 国产成年人精品一区二区| 少妇裸体淫交视频免费看高清| 狂野欧美激情性bbbbbb| 国产成人a区在线观看| 欧美最新免费一区二区三区| 美女国产视频在线观看| 禁无遮挡网站| 男人添女人高潮全过程视频| 国产成人精品久久久久久| 国产精品一及| 国产亚洲精品久久久com| 免费黄网站久久成人精品| 久久精品国产亚洲av涩爱| 国产免费一级a男人的天堂| 中文欧美无线码| 精品国产乱码久久久久久小说| 日韩精品有码人妻一区| 草草在线视频免费看| 亚洲精品久久午夜乱码| 美女视频免费永久观看网站| 成年版毛片免费区| 久久精品久久久久久噜噜老黄| 综合色av麻豆| 久久精品久久久久久久性| 欧美高清成人免费视频www| 日韩欧美 国产精品| 亚洲精品色激情综合| 国模一区二区三区四区视频| 亚洲人成网站高清观看| 欧美xxⅹ黑人| 成人特级av手机在线观看| 人妻少妇偷人精品九色| 欧美高清成人免费视频www| 亚洲精品日韩在线中文字幕| 我要看日韩黄色一级片| 一级av片app| 日韩成人伦理影院| 少妇被粗大猛烈的视频| 国产成人免费无遮挡视频| 干丝袜人妻中文字幕| 麻豆成人av视频| 亚洲一区二区三区欧美精品 | a级毛色黄片| 国产成人午夜福利电影在线观看| 少妇人妻一区二区三区视频| 丝袜美腿在线中文| av线在线观看网站| 国产免费福利视频在线观看| xxx大片免费视频| 日本欧美国产在线视频| 国产高清三级在线| 午夜福利高清视频| 久久综合国产亚洲精品| 特大巨黑吊av在线直播| 99久久九九国产精品国产免费| 看非洲黑人一级黄片| 亚洲国产精品成人综合色| 午夜福利网站1000一区二区三区| 国产亚洲最大av| 51国产日韩欧美| 精品久久久久久久人妻蜜臀av| 亚洲激情五月婷婷啪啪| 国产视频首页在线观看| 国产一级毛片在线| 免费av不卡在线播放| 国产伦精品一区二区三区视频9| av卡一久久| 九九爱精品视频在线观看| .国产精品久久| 亚洲av.av天堂| 能在线免费看毛片的网站| 国产大屁股一区二区在线视频| 亚洲第一区二区三区不卡| 五月玫瑰六月丁香| av女优亚洲男人天堂| 久久这里有精品视频免费| 久久精品人妻少妇| 少妇裸体淫交视频免费看高清| 欧美 日韩 精品 国产| 亚洲色图av天堂| 亚洲怡红院男人天堂| 18禁在线无遮挡免费观看视频| 直男gayav资源| 国产有黄有色有爽视频| 久久久国产一区二区| 久久久久精品性色| 久久久久久久久久成人| 草草在线视频免费看| 亚洲精品国产av成人精品| 国产成人a区在线观看| 午夜福利高清视频| 女人十人毛片免费观看3o分钟| 国产欧美日韩一区二区三区在线 | 国产精品熟女久久久久浪| 欧美一级a爱片免费观看看| 中文天堂在线官网| 精品少妇久久久久久888优播| 亚洲精品国产av蜜桃| 日韩国内少妇激情av| 精品酒店卫生间| 三级国产精品片| 国产大屁股一区二区在线视频| 欧美老熟妇乱子伦牲交| 国产高清不卡午夜福利| 日本免费在线观看一区| 国产精品无大码| 精品午夜福利在线看| 色吧在线观看| 成人漫画全彩无遮挡| 一区二区三区免费毛片| 一级黄片播放器| 亚洲精品乱码久久久v下载方式| 亚洲人成网站在线观看播放| 国产精品一区二区三区四区免费观看| 五月伊人婷婷丁香| 亚洲国产最新在线播放| 国产成人午夜福利电影在线观看| 美女主播在线视频| 日本一二三区视频观看| 伦精品一区二区三区| 26uuu在线亚洲综合色| 最后的刺客免费高清国语| 午夜精品一区二区三区免费看| 亚洲一区二区三区欧美精品 | 另类亚洲欧美激情| 久久久久精品久久久久真实原创| 亚洲av免费高清在线观看| 国产伦精品一区二区三区四那| 黄片无遮挡物在线观看| av卡一久久| 成人美女网站在线观看视频| 亚洲激情五月婷婷啪啪| 久久午夜福利片| 在线观看一区二区三区激情| 自拍偷自拍亚洲精品老妇| 丝袜美腿在线中文| 精品一区二区三卡| 精品国产乱码久久久久久小说| 国产在视频线精品| 精品视频人人做人人爽| 国产一区二区三区av在线| 一区二区三区免费毛片| 18+在线观看网站| 男女下面进入的视频免费午夜| 国语对白做爰xxxⅹ性视频网站| 国产大屁股一区二区在线视频| 亚洲av二区三区四区| 国产精品99久久久久久久久| 久久久a久久爽久久v久久| 国产伦理片在线播放av一区| 永久免费av网站大全| 亚洲天堂av无毛| 亚洲精品国产色婷婷电影| 亚洲欧美清纯卡通| 嘟嘟电影网在线观看| 午夜精品一区二区三区免费看| 成人免费观看视频高清| 亚洲精品aⅴ在线观看| 我的女老师完整版在线观看| 一级二级三级毛片免费看| 国产亚洲av嫩草精品影院| 婷婷色av中文字幕| 国产在线男女| 久久精品国产自在天天线| 亚洲av免费高清在线观看| 中国美白少妇内射xxxbb| 精品一区二区免费观看| 秋霞在线观看毛片| 国产熟女欧美一区二区| 午夜亚洲福利在线播放| 国产色爽女视频免费观看| 国产精品国产三级专区第一集| 国产综合精华液| 狂野欧美白嫩少妇大欣赏| 人人妻人人爽人人添夜夜欢视频 | 亚洲欧美日韩另类电影网站 | 亚洲欧美一区二区三区黑人 | 综合色av麻豆| 又黄又爽又刺激的免费视频.| 男女啪啪激烈高潮av片| 亚洲综合精品二区| 18+在线观看网站| 99久久精品热视频| 亚洲成人一二三区av| 永久网站在线| 久久ye,这里只有精品| 成年免费大片在线观看| 国内精品美女久久久久久| 久久综合国产亚洲精品| 真实男女啪啪啪动态图| 国产真实伦视频高清在线观看| 亚洲欧美清纯卡通| 白带黄色成豆腐渣| 综合色丁香网| 国产成人a区在线观看| 深夜a级毛片| 插阴视频在线观看视频| 国产又色又爽无遮挡免| 色网站视频免费| 欧美极品一区二区三区四区| 国产精品99久久久久久久久| 亚洲精品久久久久久婷婷小说| 日韩不卡一区二区三区视频在线| 亚洲最大成人av| 亚洲欧美清纯卡通| 中文在线观看免费www的网站| www.av在线官网国产| 亚洲av国产av综合av卡| 亚洲精品色激情综合| 丰满少妇做爰视频| 性色av一级| 欧美老熟妇乱子伦牲交| 亚洲精品日韩在线中文字幕| 国产亚洲午夜精品一区二区久久 | 久久久a久久爽久久v久久| 亚洲电影在线观看av| 丰满人妻一区二区三区视频av| 午夜老司机福利剧场| 亚洲图色成人| www.av在线官网国产| 日韩欧美精品v在线| 欧美另类一区| 亚洲av免费高清在线观看| 精品亚洲乱码少妇综合久久| 伦理电影大哥的女人| 欧美精品人与动牲交sv欧美| 五月玫瑰六月丁香| 亚洲av免费高清在线观看| 又爽又黄a免费视频| 三级男女做爰猛烈吃奶摸视频| 国产伦精品一区二区三区视频9| 最后的刺客免费高清国语| 伦精品一区二区三区| 国产亚洲91精品色在线| 精品亚洲乱码少妇综合久久| 美女主播在线视频| 干丝袜人妻中文字幕| 亚洲精品aⅴ在线观看| 性插视频无遮挡在线免费观看| 欧美变态另类bdsm刘玥| 国产男人的电影天堂91| 久久久久久久久久久丰满| 只有这里有精品99| 97精品久久久久久久久久精品| 色网站视频免费| 久久综合国产亚洲精品| 亚洲精华国产精华液的使用体验| 中国三级夫妇交换| 精品人妻一区二区三区麻豆| av国产免费在线观看| 日韩 亚洲 欧美在线| 国产老妇女一区| 免费大片18禁| 18+在线观看网站| 国产一区亚洲一区在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产精品无大码| 亚洲成人久久爱视频| 夫妻性生交免费视频一级片| 丝袜美腿在线中文| 亚洲国产日韩一区二区| 欧美一区二区亚洲| 日日摸夜夜添夜夜添av毛片| a级毛片免费高清观看在线播放| 国产综合懂色| 最新中文字幕久久久久| 超碰97精品在线观看| 日韩av不卡免费在线播放| 婷婷色综合www| 日韩,欧美,国产一区二区三区| 国产一区二区亚洲精品在线观看| 青春草视频在线免费观看| 99久久精品热视频| 高清毛片免费看| 中文欧美无线码| 最近2019中文字幕mv第一页| 交换朋友夫妻互换小说| 男插女下体视频免费在线播放| 国产精品99久久久久久久久| 黑人高潮一二区| 国产人妻一区二区三区在| 国产午夜精品一二区理论片| 亚洲精品日韩在线中文字幕| 观看美女的网站| 香蕉精品网在线| 少妇高潮的动态图| eeuss影院久久| 欧美xxxx黑人xx丫x性爽| 在线观看人妻少妇| 国产男女超爽视频在线观看| 少妇人妻久久综合中文| 精品久久久久久久末码| 中文字幕免费在线视频6| 成人亚洲精品一区在线观看 | 精品少妇黑人巨大在线播放| 久久99热这里只有精品18| 69av精品久久久久久| 蜜桃久久精品国产亚洲av| 欧美三级亚洲精品| 国产成人a∨麻豆精品| 亚洲无线观看免费| 亚洲综合色惰| 亚洲自拍偷在线| 成年免费大片在线观看| 亚洲精品国产av成人精品| 97人妻精品一区二区三区麻豆| 九色成人免费人妻av| 亚洲色图综合在线观看| 国产综合精华液| 成年人午夜在线观看视频| 性色av一级| 亚洲av成人精品一区久久| 国产成人freesex在线| 男女边摸边吃奶| 久久精品综合一区二区三区| 99热全是精品| 18禁动态无遮挡网站| 亚洲色图av天堂| 日本色播在线视频| 欧美少妇被猛烈插入视频| 婷婷色综合www| 中文在线观看免费www的网站| 高清午夜精品一区二区三区| 国产男女内射视频| 国产精品一二三区在线看| 日韩,欧美,国产一区二区三区| 久久99精品国语久久久| 亚洲最大成人av| 免费在线观看成人毛片| 少妇高潮的动态图| 99热网站在线观看| 午夜福利视频1000在线观看| 国产久久久一区二区三区| 一个人看视频在线观看www免费| 三级国产精品欧美在线观看| 看非洲黑人一级黄片| 综合色av麻豆| 国产老妇伦熟女老妇高清| 老师上课跳d突然被开到最大视频| 高清视频免费观看一区二区| 亚洲av日韩在线播放| 久久国产乱子免费精品| 久久久久久久久大av| 有码 亚洲区| 校园人妻丝袜中文字幕| 最近中文字幕高清免费大全6| 国产69精品久久久久777片| 欧美日本视频| 亚洲熟女精品中文字幕| 欧美激情在线99| 91精品国产九色| 黄色一级大片看看| 街头女战士在线观看网站| 激情五月婷婷亚洲| 久久久久国产精品人妻一区二区| 免费av观看视频| 亚洲成人中文字幕在线播放| 国产成人免费无遮挡视频| 99热这里只有精品一区| 岛国毛片在线播放| 五月玫瑰六月丁香| 国语对白做爰xxxⅹ性视频网站| 直男gayav资源| av在线天堂中文字幕| 热re99久久精品国产66热6| 特级一级黄色大片| av在线观看视频网站免费| 噜噜噜噜噜久久久久久91| 国产精品99久久99久久久不卡 | 黄片wwwwww| 日韩成人av中文字幕在线观看| 国产精品久久久久久久久免| 日韩,欧美,国产一区二区三区| 男女下面进入的视频免费午夜| 人妻少妇偷人精品九色| 大香蕉久久网| 日韩强制内射视频| 免费看日本二区| 女人久久www免费人成看片| 国产精品av视频在线免费观看| 亚洲精品影视一区二区三区av| 国产一区有黄有色的免费视频| 亚洲成人精品中文字幕电影| 插阴视频在线观看视频| 久久精品国产亚洲av涩爱| 精品久久久久久久久av| 国产在线一区二区三区精| 蜜臀久久99精品久久宅男| 禁无遮挡网站| 国产精品一区二区三区四区免费观看| av一本久久久久| 亚洲人成网站高清观看| 国产精品秋霞免费鲁丝片| 建设人人有责人人尽责人人享有的 | 黄片无遮挡物在线观看| 欧美日韩国产mv在线观看视频 | 在线观看人妻少妇| 国产精品一区二区性色av| 亚洲精品国产色婷婷电影| 成人鲁丝片一二三区免费| 亚洲在久久综合| 婷婷色综合大香蕉| 一区二区三区免费毛片| 又爽又黄a免费视频| 熟女电影av网| 久久久国产一区二区| 最近的中文字幕免费完整| 亚洲精品一二三| 国产精品福利在线免费观看| 最近中文字幕高清免费大全6| 九九爱精品视频在线观看| 熟女电影av网| 国产精品三级大全| 久久久久性生活片| 十八禁网站网址无遮挡 | 成年免费大片在线观看| 日韩国内少妇激情av| 在线观看一区二区三区| 最近手机中文字幕大全| 一级毛片aaaaaa免费看小| 午夜日本视频在线| 久久综合国产亚洲精品| 欧美高清性xxxxhd video| 你懂的网址亚洲精品在线观看| 日韩av免费高清视频| 亚洲精品成人av观看孕妇| 在线精品无人区一区二区三 | 久久国产乱子免费精品| 美女内射精品一级片tv| av在线老鸭窝| 身体一侧抽搐| 国产高清国产精品国产三级 | 日韩亚洲欧美综合| 午夜免费观看性视频| 中文欧美无线码| 欧美日韩国产mv在线观看视频 | 日韩一区二区视频免费看| 国产欧美另类精品又又久久亚洲欧美| 国产精品三级大全| 日韩免费高清中文字幕av| 联通29元200g的流量卡| 国产精品偷伦视频观看了| 性插视频无遮挡在线免费观看| 国产成人精品一,二区| 日韩av在线免费看完整版不卡| 综合色丁香网| 尤物成人国产欧美一区二区三区| 永久网站在线| 舔av片在线| 丝瓜视频免费看黄片| 亚洲第一区二区三区不卡| 嫩草影院精品99| 99热这里只有精品一区| 永久免费av网站大全| 日本av手机在线免费观看| 久久精品国产鲁丝片午夜精品| 麻豆久久精品国产亚洲av| 偷拍熟女少妇极品色| 直男gayav资源| 亚洲一区二区三区欧美精品 | 国产又色又爽无遮挡免| 日本午夜av视频| 国产片特级美女逼逼视频| 日韩,欧美,国产一区二区三区| av又黄又爽大尺度在线免费看| 国产69精品久久久久777片| av.在线天堂| 在线观看国产h片| 久久精品国产鲁丝片午夜精品| 亚洲激情五月婷婷啪啪| 欧美成人一区二区免费高清观看| 亚洲人成网站在线观看播放| 七月丁香在线播放| 亚洲成人久久爱视频| 久久久久久久久久久免费av| 亚洲精品乱久久久久久| 啦啦啦啦在线视频资源| 美女脱内裤让男人舔精品视频| 国产视频首页在线观看| 亚洲不卡免费看| 男女边摸边吃奶| 国语对白做爰xxxⅹ性视频网站| 国产毛片在线视频| 九色成人免费人妻av| 亚洲美女搞黄在线观看| 国产黄色免费在线视频| 午夜老司机福利剧场| 在线精品无人区一区二区三 | 一区二区av电影网| 久久久久国产精品人妻一区二区| 老司机影院毛片| 美女xxoo啪啪120秒动态图| 久久久a久久爽久久v久久| 免费观看a级毛片全部| 下体分泌物呈黄色| 成人免费观看视频高清| 97在线视频观看| 新久久久久国产一级毛片| 亚洲欧美日韩无卡精品| 毛片女人毛片| 国产黄片美女视频| 美女xxoo啪啪120秒动态图| 男人添女人高潮全过程视频| 搡女人真爽免费视频火全软件| 麻豆精品久久久久久蜜桃| 大话2 男鬼变身卡| 国产免费一区二区三区四区乱码| 在线精品无人区一区二区三 | 少妇人妻久久综合中文| 三级国产精品片| 1000部很黄的大片| 91精品伊人久久大香线蕉| 久久久欧美国产精品| 亚洲av福利一区| 香蕉精品网在线| 欧美亚洲 丝袜 人妻 在线| 久久国内精品自在自线图片| 九草在线视频观看| 一级二级三级毛片免费看| 大又大粗又爽又黄少妇毛片口| 免费播放大片免费观看视频在线观看| 麻豆精品久久久久久蜜桃| 成年女人看的毛片在线观看| 亚洲aⅴ乱码一区二区在线播放| 搞女人的毛片| 亚洲精品国产色婷婷电影| 青春草视频在线免费观看| 国产精品秋霞免费鲁丝片| 特大巨黑吊av在线直播| 午夜日本视频在线| 国产人妻一区二区三区在| 你懂的网址亚洲精品在线观看| 性插视频无遮挡在线免费观看| 色综合色国产| 欧美性猛交╳xxx乱大交人| 成人午夜精彩视频在线观看| 亚洲av福利一区| 观看免费一级毛片| 亚洲美女视频黄频| 国产v大片淫在线免费观看| freevideosex欧美| 水蜜桃什么品种好| 午夜福利网站1000一区二区三区| 18禁裸乳无遮挡免费网站照片| 国产白丝娇喘喷水9色精品| 亚洲精品久久午夜乱码| 久久亚洲国产成人精品v| 国产精品三级大全| 精品久久久噜噜| av又黄又爽大尺度在线免费看| av在线老鸭窝| 最近最新中文字幕免费大全7| av线在线观看网站| 视频区图区小说| 91精品国产九色| 亚洲熟女精品中文字幕| 成人鲁丝片一二三区免费| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美日韩另类电影网站 | 97在线人人人人妻| 国产精品av视频在线免费观看| 久久久久久久大尺度免费视频| 日本与韩国留学比较| 青春草国产在线视频| 亚洲精品久久久久久婷婷小说| 国产中年淑女户外野战色| 精品人妻偷拍中文字幕| 国产人妻一区二区三区在| 制服丝袜香蕉在线| 亚洲国产精品成人久久小说| 1000部很黄的大片| 国产老妇伦熟女老妇高清| 两个人的视频大全免费| 日本欧美国产在线视频| 国产成人精品久久久久久| 婷婷色综合www| av.在线天堂| 国产美女午夜福利| 亚洲欧美日韩无卡精品| 在线亚洲精品国产二区图片欧美 | 欧美日韩视频精品一区| 免费看av在线观看网站| 亚洲熟女精品中文字幕| 亚洲av不卡在线观看| av国产免费在线观看| 国产精品蜜桃在线观看| 精品国产乱码久久久久久小说| 亚洲精品中文字幕在线视频 | 女人十人毛片免费观看3o分钟| 少妇高潮的动态图| 久久午夜福利片| 精品国产露脸久久av麻豆| 亚洲欧美清纯卡通| 永久网站在线| 国产精品一区www在线观看| 22中文网久久字幕| 99热全是精品| 成人黄色视频免费在线看| 成人二区视频| 亚洲成人中文字幕在线播放| 中文字幕制服av| 久久99热这里只有精品18|