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

    Meter-Scale Thin-Walled Structure with Lattice Infill for Fuel Tank Supporting Component of Satellite:Multiscale Design and Experimental Verification

    2024-02-19 12:00:36XiaoyuZhangHuizhongZengShaohuiZhangYanZhangMiXiaoLipingLiuHaoZhouHongyouChaiandLiangGao

    Xiaoyu Zhang,Huizhong Zeng,Shaohui Zhang,Yan Zhang,Mi Xiao,Liping Liu,Hao Zhou,★,Hongyou Chai and Liang Gao

    1Beijing Key Laboratory of Lightweight Multi-Functional Composite Materials and Structures,Beijing Institute of Technology,Beijing,100081,China

    2Beijing Key Laboratory of Intelligent Space Robotic Systems Technology and Applications,Beijing Institute of SpacecraftSystem Engineering,China Academy of Space Technology,Beijing,100094,China

    3Key Laboratory of Metallurgical Equipment and Control Technology,Ministry of Education,Wuhan University of Science and Technology,Wuhan,430081,China

    4State Key Laboratory of Digital Manufacturing Equipment and Technology,Huazhong University of Science and Technology,Wuhan,430074,China

    ABSTRACT

    Lightweight thin-walled structures with lattice infill are widely desired in satellite for their high stiffness-to-weight ratio and superior buckling strength resulting from the sandwich effect.Such structures can be fabricated by metallic additive manufacturing technique,such as selective laser melting(SLM).However,the maximum dimensions of actual structures are usually in a sub-meter scale,which results in restrictions on their appliance in aerospace and other fields.In this work,a meter-scale thin-walled structure with lattice infill is designed for the fuel tank supporting component of the satellite by integrating a self-supporting lattice into the thickness optimization of the thin-wall.The designed structure is fabricated by SLM of AlSi10Mg and cold metal transfer welding technique.Quasi-static mechanical tests and vibration tests are both conducted to verify the mechanical strength of the designed large-scale lattice thin-walled structure.The experimental results indicate that the meter-scale thin-walled structure with lattice infill could meet the dimension and lightweight requirements of most spacecrafts.

    KEYWORDS

    Thin-walled structure;lattice infill;supporting component;selective laser melting;satellite

    1 Introduction

    Lightweight and high-efficiency structures are highly desired in the fields of aircraft and spacecraft for improving system performance and conserving propulsion energy during launch or in flight[1,2].In particular,thin-walled structures with lattice infills have attracted numerous attention as an innovative potential structural pattern characterized by improved performance and reduced weights[3,4].This kind of structures,also termed as infill structures,can effectively combine the advantages of the Lightweight lattice and high load-bearing shell structure.With the remarkable development of additive manufacturing (AM) technology,complex lattice infill structures can be rapidly fabricated.Besides,the infilled lattice can facilitate the AM process by providing support for the thin-walled shell.Recently,lightweight metal lattice structures have been successfully fabricated through the selective laser melting(SLM)technique[1,5–8].

    Considering the variation in geometric shape,dimensions and size dependent mechanical properties,the elastic deformation,mechanical strength,energy absorption and vibration isolation properties of lattice structures fabricated through SLM have been studied adequately by means of experiments,calculations and mechanical theory[9–15].To achieve extreme lightweight property,lattice structures are usually employed in the space domain.On the basis of lattice structures,payloads supporting structures and thermal controllers of satellites have been designed and verified through experiments[2,16].These secondary structures or thermal functional structures of satellites usually bear relatively small mechanical loads,and their dimensions are usually in the range of approximately 100~300 mm.For example,the satellite main structure made of lattice sandwich structures was designed,verified [1],and launched into space with the QIANSHENG-1 satellite on August 17,2019.The dimensions of the prototype structure and the actual structure of the satellite launched into space are 400 mm3×400 mm3×400 mm3and 525 mm3×480 mm3×425 mm3,respectively.However,according to several different mission requirements for the satellites,the dimensions of the fuel tank supporting components are usually in the range of approximately 1~3 m.Their structure patterns are generally assumed to be a thin-wall stiffened structure,and the most applied manufacturing method is mainly machining or casting.Due to the high stiffness required,the tie rods and other external structures are usually designed on the fuel tank supporting component to further improve its stiffness,which leads to the complex structure and insufficient lightweight.Therefore,the meter-scale thin-walled structure with lattice infill needs to be urgently developed to meet the dimension requirements and lightweight property of most spacecrafts.

    2 Related Works

    In general,a lattice infill structure consists of a solid shell and lattice infill,in which the solid shell acts as an exterior thin surface layer to surround these lattices.It is worth noting that the kind of structure has typical two-scale characteristics.Thus,a multiscale design method,with flexible designability for both macroscopic solid shell and microscopic infilled lattice,would be a highly effective tool to create such a thin-walled structure with lattice infill.Recently,a significant advance in multiscale design methods enables scholars and engineers to pursue lightweight and high-performance infill structures.Wu et al.[17] designed the infill porous structures by introducing local volume constraint into a density-based topology optimization method.Clausen et al.[18]presented a coating approach to topology optimization to design the macroscale layout for both a solid shell and a porous triangular infill.Wadbro et al.[19]proposed a multiscale topology optimization method for an infill structure with a thin skin and uniform optimized lattices.Then,they further extended this method to suppress the undemand vibration of such an infill structure[4].Zhou et al.[20]presented a concurrent topology optimization method of shells with self-supporting infills,and a precise overhang control was tailored to satisfy the additive manufacturing constraint.Liu et al.[21]developed a hybrid MMCMMV approach for designing additive manufacturing-oriented shell structures with graded lattice infills.Recently,the authors[22,23]proposed a multiscale topology optimization method for designing geometrically asymmetric sandwich structures with graded or uniform microstructure infills.Later,this method was further extended to design the graded lattice sandwich structure with lightweight lattices as self-supporting infills [24,25].It can be noted that the multiscale design method for thinwalled structures with lattice infills has been a hot topic.Unexpectedly,its extraordinary potential has not yet been extensively exploited to study satellite components.Therefore,an effective multiscale design method is still in great demand to enhance the innovative design ability of satellite components for better performance and lighter weight.

    In presented work,the multiscale design method is further explored for a meter-scale thin-walled structure with lightweight lattice infill as fuel tank supporting component of a satellite.To achieve the goal of better performance and lighter weight while meeting all design requirements,the design procedure is performed to optimize thin-walled structure with lattice infill by integrating the selfsupporting lattice into the thickness optimization of thin-wall structure.The maximum dimension of the satellite component is 1.255 m in diameter,and it is fabricated through SLM of AlSi10Mg and cold metal transfer welding technique.To the best of our knowledge,this metal lattice thin-walled structure used as the fuel tank supporting component has the largest dimension.Both quasi-static mechanical tests and vibration tests are conducted to verify the mechanical strength of the designed large-scale lattice thin-walled structure.

    3 Infill Lattice Design

    3.1 Design Constraints for AM

    AM techniques,such as SLM,have become the first choice of fabricating 3D lattice structures with complex geometry due to its dramatic manufacturing freedoms.However,there are still some geometrical constraints of AM that need to be carefully considered during the design process [26–30],mainly manifesting in the following four aspects: the absence of enclosed voids,the absence of an islet,minimum/maximum dimensions,and an overhang angle that is measured with respect to the horizontal plane.It is worth noting that only the overhang angle affects the configuration design of 3D lattices composed of slender structures.Specifically,the dimensional error and surface roughness of 3D lattices fabricated by SLM increase as the angle between the lattices and the horizontal direction decreases [10,31–33].Each part of the 3D lattice with an overhang angle of less than 35°must be sufficiently supported to prevent structural material from being distorted or even collapsed during fabrication process.To avoid the supporting material in the fabrication process,self-supporting lattices with overhang angles larger than 35°are designed in this work.Different from traditional lattices composed of struts at any angle,the self-supporting lattices only consist of inclined struts,except for horizontal struts.

    3.2 Lattice Geometry Description

    Inspired by the multi-fold rotational symmetry of crystallography in solid physics,several classes of self-supporting lattices for AM can be designed by rotating a slant strut around the vertical direction at an angle distance of 120°,90°,or 60°[34].In this work,the lattice with four-fold rotational symmetry(i.e.,body-centered cubic cell) is employed for its low density and high-efficiency load-bearing,as shown in Fig.1.For the sake of high-accuracy and free-form geometry design,a topology description function(TDF)based on a level set is used to describe the explicit geometry of rod-shaped parts within lattice unit cell[24,35,36],as expressed in Eq.(1).

    whereDis a fixed Eulerian domain andΩis a subset ofDoccupied by a rod-shaped part.?Ωis the structural boundary and x denotes the point coordinates inD.The rod-shaped part consists of a cylinder and two spheres,and their TDFs areφpart,c(x,y,z),φpart,s1(x,y,z) andφpart,s2(x,y,z),respectively.The two spheres are respectively centered at(x1,y1,z1)and(x2,y2,z2),which are also the center of the ends of the cylinder.Namely,the diameters of the ends of the cylinder and two spheres are equal,which are denoted bytpart.The center point coordinate of the cylinder is(x0,y0,z0),and its length is denoted byLpart,which can be calculated asLpart=

    Figure 1 :Schematic illustration of a lattice unit cell

    Andθpartis the slant angle of a rod-shaped part with respect to the horizontal plane,which can be calculated by using the point coordinates of both ends(x1,y1,z1)and(x2,y2,z2),as follows:

    whereLd=dx2=x2-x0,dy2=y2-y0,dz2=z2-z0,dx=x-x0,dy=y-y0,dz=z-z0.

    It is observed from Fig.1b that the lattice unit cell can be viewed as a combination of four rodshaped parts.Thus,its TDF of a lattice unit cell can be expressed by

    whereφS(x)=max(φi),andφi=φpart,i(x),i=1,2,3,4.φpart,iis the TDF of a rod-shaped part in the lattice unit cell.ΩSdenotes the region occupied by the lattice unit cell,and obviously?ΩSis the structural boundary of the lattice unit cell.

    3.3 Lattice Property Estimation

    Benefitting from the explicit geometry description based on a level set function,its effective density of a lattice unit cell can be calculated directly in a numerical integration way[25,37],as follows:

    whereΩSand|ΩS|represent the design domain and the volume of a lattice unit cell,respectively.H(·)is the Heaviside function to indicate different parts of the design domain.

    With the help of the numerical homogenization method as an effective means of estimating the effective mechanical properties of periodic materials,its effective elasticity tensors of a lattice unit cell can be calculated as follows[24,38,39]:

    wherei,j,kandlare equal to 1,2,...,d,anddis the spatial dimension.Dpqrsdenotes the elasticity tensor of the solid material.The displacement field u(ij)corresponds to strain fieldwhich can be calculated by solving the following equation at microscale with the given macroscopic strain:

    where v(kl)represents the virtual displacement field,and(ΩS)denotes the kinematically admissible displacement field satisfying the periodic conditions for the lattice unit cell.

    Based on the effective elasticity tensors of a lattice unit cell,its macroscopic stiffness matrix k can be calculated by

    where B is the strain-displacement matric.It is noted that the calculated macroscopic stiffness matrix k will be used as the elasticity equivalent stiffness of lattice infill for the subsequent thickness optimization of thin-wall.

    4 Thickness Optimization of Thin-Walled Structure

    The workflow of designing thin-walled structure with lattice infill is shown in Fig.2.In the structural requirement analysis phase,the design constraints and design objective are determined in accordance with the satellite requirements.In the infill lattice design phase,the structural geometry of the infill lattice is represented by the TDF in accordance with its design parameters.And the equivalent stiffness of infill lattice is calculated through the numerical homogenization method.In the thin-wall thickness optimization phase,the thickness of thin-wall is optimized for the minimum mass with the constraint of sufficient stiffness and fundamental frequency.In the geometry reconstruction phase,the apparent thin-wall is built,and the joint structures for screws are built as solids to improve the local mechanical strength.And the lattice infill is filled in the internal enclosed space of the thin-wall.Finally,the 3D model is outputted to be fabricated through AM.

    The connection relationship among the main supporting component,fuel tanks,and other part of satellite is shown in Fig.3.The main supporting component is firmly connected to the other part of satellite by 16 screws,that is,a cylindrical shell.The diameter of main supporting component is 1255 mm.The weights of both fuel tanks supported by the main supporting component are 71.5 kg.They are fixed to the main supporting component by 16 screws.

    In the thickness optimization model of thin-wall,the design constraints mainly include the dimension constraints and mechanical properties constrains.The dimension constraint is that the height of the fuel tank supporting component is less than 151 mm,which is limited by the spatial layout of various equipment inside the satellite.What’s more,the fuel tank supporting component is designed with six parts,as shown in Fig.4.This is mainly because the SLM equipment with the limited fabrication dimension of 600 mm3×600 mm3×600 mm3is used in this work.

    Figure 2 :Workflow of designing thin-walled structure with lattice infill

    Figure 3 : Schematic of connection relationship among main supporting component,fuel tanks,and other part of satellite

    The mechanical properties constraint is that the fundamental frequency of the fuel tank supporting component is larger than 80 Hz,and the maximum stress is less than 100 MPa.The design objective is to minimize the mass of the fuel tank supporting component.The macroscopic shape of the thinwalled structure is a spherical crown with the height of 151 mm.A widely used body-centered cubic cell with dimensions of 5 mm3×5 mm3×5 mm3,and struts diameter of 0.5 mm is employed as the lattice infill.The finite element model of the thin-walled structure with lattice infill to be optimized is built,as shown in Fig.4.The thin-wall is simulated with shell elements,and the lattice infill is simulated with homogenized solid elements.The boundary condition is the zero displacement constraint of the 16 screws at the edge of the fuel tank supporting component.Each fuel tank is simplified as a mass element with mass of 71.5 kg and inertia moment of 14.1 kg·m2in three orthogonal directions.The mass element is connected to the fuel tank supporting component via the multipoint constraint element at the 16 screw positions.And thus,the optimization model of thin-walled structure can be formulated as follows:

    where t is then-dimensional vector of the thickness variables of shell elements.m(t)is the objective function,i.e.,the mass of the fuel tank supporting component.f(t) is the fundamental frequency,andf0is the minimum natural frequency required.σmax(t) is the maximum stress,andσmis the allowable maximum stress.Tloweris the lower limit of the thickness variableti,and its value of 0.5 mm adopted here,is the minimum dimension where most SLM equipment can steadily achieve without much degradation of the mechanical property.The thickness of the thin-wall is extensively optimized everywhere with the objective of minimum mass using Optistruct software.

    The optimized results show that the variable thickness varies from 0.5 to 1.8 mm,and the mass of the fuel tank supporting component is 12.83 kg.The fundamental frequency is 86.1 Hz,and the maximum stress is 80 MPa,as shown in Figs.5b and 5c.To sufficiently assess the effectiveness of the proposed design method,the comparison is conducted between the thin-walled structure with lattice infill and its traditional frame-type one.As shown in Figs.6a and 6b,the mass of the frametype structure is 15.4 kg.The fundamental frequency is 52.9 Hz,and the maximum stress is 70 MPa.It can be viewed that although the maximum stress of thin-walled structure with lattice infill is larger than that of the frame-type structure,its mass is reduced by 16.7%,and its fundamental frequency is increased nearly over half.The mentioned-above discussions and analysis demonstrate the effectiveness of the proposed method and the ultimate performance of designed thin-walled structure with lattice infill.

    Figure 5 :Mechanical analysis of thin-walled lattice structure for the fuel tank supporting component

    Figure 6 : Mechanical analysis of traditional frame-type structure for the fuel tank supporting component

    The apparent thin-wall is built with Pro/E software in accordance with the optimized results.The joint structures for screws are built as solids to improve the local mechanical strength.The lattice infill is filled in the internal enclosed space of the thin-wall.The resulting model to be fabricated using AM technique is built by the merging the geometric models obtained in the former two processes as shown in Fig.7.

    Figure 7 :Geometric modeling of fuel tank supporting component

    5 Fabrication of Fuel Tank Supporting Component

    AlSi10Mg alloy powder,with an average particle diameter of 30 μm,is used to fabricate the bracket via the SLM equipment(BLT-S600)supplied by BLT Company.The main process parameters are as follows:laser power 200–300 W,scan speed 1000 mm/s–2000 mm/s,hatch spacing 0.15 mm–0.20 mm,and layer thickness 30 μm.The SLM-fabricated structure is annealed at 200°C–300°C for 3 h to release residual stress and enhance ductility properties.Before welding,the region to be welded is machined to remove the oxide and impurities,in order to improve the welding quality and dimensional precision.Cold metal transfer welding is used in this work because of its superiority in welding aluminum alloy thin-wall structures with excessively small heat input.The width of welding lines ranges from 5 mm to 8 mm,and no sag in welding lines or other visible defects are found by X inspection.After the welding,heat treatment is conducted to release the residual stress.Finally,the integrated structure is machined at the interfaces between the fuel tank supporting component and fuel tanks or cylindrical shell of the satellite to ensure the installation accuracy being superior to 0.1 mm.The SLM fabricated parts before and after welding is shown in Fig.8.

    Figure 8 :Photos of the SLM fabricated parts before and after welding

    6 Quasi-Static Mechanical and Vibration Experiments

    Quasi-static mechanical tests are conducted to examine the mechanical strength of the fuel tank supporting component.The schematic and photo of the test system are illustrated in Figs.9a and 9b.The fuel tank supporting component is connected to a cylindrical shell,which is an imitation of the principal cylinder of the satellite and is fixed to the ground.Two circular steel disks with diameter of 507 mm are used as the imitation of two fuel tanks.Each circular disk is connected to the fuel tank supporting component by 16 screws.Two loading beams are connected to the circular disks with a hinge joint.Vertical(Z direction)and horizontal forces(X or Y direction)are simultaneously applied to each disk by adjusting the angle between the loading beam and disk.The resultant force applied to each beam is measured using a force sensor.The test procedure is shown in Table 1.The resultant force of components in X and Z directions is loaded from 0 to 8791 N in 10 min.Then,the peak force is held for 1 min.The force is unloaded from 8791 to 0 N in 10 min.The test procedure for the resultant force of components in Y and Z directions is the same as that in X and Z directions.The vertical displacement at the loading points and the internal edge of the fuel tank supporting component are measured with four displacement sensors,as shown in Fig.9c.The strain at 34 points,especially the positions close to the welding lines and screws,are measured using strain rosettes.The positions of the strain rosettes are shown in Fig.9d.The strain rosettes on the top surface and that on the lower surface of the fuel tank supporting component are 17.

    Table 1 : Quasi-static mechanical test conditions

    Besides the quasi-static mechanical tests,vibration experiments are also implemented to verify the rationality of structure design and manufacturing process quality of the SLM and cold metal transfer welding.The fuel tank supporting component is assembled to the principle cylinder of satellite by screws.The entire satellite is exposed to the vibration loads in X (horizontal),Y (horizontal),and Z directions (vertical),respectively.The photo of the vibration test system is illustrated in Fig.10a.The acceleration of fuel tank supporting component is measured with three acceleration sensors.The positions of three acceleration sensors are shown in Fig.10b.The conditions of sinusoidal vibration experiments are listed in Table 2.

    Table 2 : Vibration tests conditions

    Figure 10 :Photo of the installed fuel tank supporting component(a),and schematic of acceleration monitoring positions for vibration experiments(b)

    The frequency of the sinusoidal vibration is increased monotonically from 5 to 100 Hz.The maximum magnitudes of the loads are 0.9 g(from 8 to 100 Hz)in the X or Y direction and 1.2 g(from 8 to 100 Hz)in the Z direction.However,the magnitude is decreased to 0.65 g from 54 to 87 Hz in the X direction and from 71 to 84 Hz in the Y direction to avoid the strength failure of other equipment on the satellite caused by local structural resonance,which is larger than the vibration magnitude during real launching.The sweep rate of vibration is 2 Octave per minute,i.e.,the frequency of the sinusoidal vibration increased by quadruple in every minute.

    7 Results and Discussion

    The vertical displacements at the monitoring positions are plotted in Fig.11.It can be seen that the displacement increases as the loading force increases.The maximum displacement is 2.65 mm.The displacement decreases as the unloading force decreases.The residual displacement is 0 mm~0.15 mm,which is small enough to be neglected for large structures with assembling stress.It should be noted that no “instant displacement/force drop” phenomenon occurs during the loading and unloading process,indicating that no fracture occurs in the fuel tank supporting component.

    Figure 11 :Displacement results of quasi-static mechanical tests

    The strains measured at different positions of the fuel tank supporting component are plotted in Fig.12.The strain increases when the loading force increases,and decreases when the unloading force decreases.It agrees well with the variation tendency of the displacement.The maximum strain is 390 με,and the residual strain is 0~45 με.The maximum stress and residual stress can be estimated by the strain multiplied by the elastic modulus of the AlSi10Mg thin-walled structure,i.e.,about 60 GPa.The maximum stress and residual stress are estimated to be 23.4 and 2.7 MPa,respectively.The stresses are far below the yield strength of AlSi10Mg,which is approximately 160 MPa.The strain measurement results of quasi-static mechanical tests verify the elastic deformation of the fuel tank supporting component.That is to say,the fuel tank supporting component has not been damaged during the quasi-static mechanical tests.

    Figure 12 :Strain results of quasi-static mechanical tests

    The acceleration response monitored in the X,Y and Z direction are plotted in Fig.13.It can be seen that the Z directional sinusoidal vibration loads stimulate the first order structural resonance of the fuel tank supporting component at the frequency of 80.8 Hz.The value is smaller than the calculated value 86.1 Hz.This finding is mainly because the real boundary condition of the fuel tank supporting component is unideal clamped boundary due to the flexibility of the cylindrical shell,fuel tanks and other parts of satellite.The maximum acceleration responses in the X,Y and Z direction are 3.5,3.8,and 5.4 g,respectively.

    Figure 13 :Acceleration response monitored in X direction(a),Y direction(b)and Z direction(c)

    Low-magnitude vibration excitation is applied to the satellite before and after vibration experiments to examine whether the resonant frequencies are altered for the structural degradation or not [1,40],considering that the resonant frequencies can be changed as a result of structural failure or damage [41].The low-magnitude vibration excitation is 0.1 g,which increases from 5 to 100 Hz monotonously at a sweep rate of 4 Oct/min.The acceleration responses monitored by the sensors are plotted in Fig.14.The acceleration response stimulated by low-magnitude vibration excitation before and after vibration experiments almost overlap with each other,in spite of small discrepancy that usually resulted from the variation in assembly stress of the assembly containing structural parts and connections.Therefore,it is justifiable to conclude that the fuel tank supporting component fabricated through SLM has not been damaged during the high-magnitude vibration tests,and there is no structural fracture or plastic deformation found.That is to say,the fuel tank supporting component made of thin-wall and lattice infill manufactured by the SLM technique could sustain the vibration during the rocket launching.

    Figure 14 : Acceleration response comparison before and after vibration experiments in X direction(a),Y direction(b),and Z direction(c)in response to low-magnitude vibration excitation

    8 Conclusions

    To summarize,the largest known thin-walled structure with lattice infill with the diameter of 1.255 m is designed and fabricated through SLM technique and cold metal transfer welding.The structure is used as the fuel tank supporting component of the satellite.Quasi-static mechanical tests and vibration tests verify the strength of the fuel tank supporting component.The recovery deformation in quasi-static mechanical tests and unaltered resonance frequencies in vibration tests confirm that no structural failure occurs during the mechanical experiments.The design and verification approach presented in this work can benefit the application of large scale thin-walled structures with lattice infill used as fuel tank supporting components in the aerospace field and others.

    Funding Statement: The authors are grateful for the support by National Key Research and Development Program of China(2021YFF0500300,2020YFB1708300)and the National Natural Science Foundation of China(52205280,12172041).

    Author Contributions: The authors confirm contribution to the paper as follows: study conception and design:Xiaoyu Zhang,Hao Zhou,Liang Gao;data collection:Huizhong Zeng,Shaohui Zhang,Liping Liu,Hao Zhou; analysis and interpretation of results: Xiaoyu Zhang,Yan Zhang.Mi Xiao,Hongyou Chai,Liang Gao;draft manuscript preparation:Xiaoyu Zhang,Hao Zhou,Yan Zhang.Mi Xiao.All authors reviewed the results and approved the final version of the manuscript.

    Availability of Data and Materials:The data that support the findings of this study are available from the first and corresponding authors upon reasonable request.

    Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

    99国产精品免费福利视频| 黑人猛操日本美女一级片| av在线app专区| 久热这里只有精品99| 午夜福利免费观看在线| 大香蕉久久成人网| 成人18禁高潮啪啪吃奶动态图| xxxhd国产人妻xxx| 精品高清国产在线一区| 国产精品一区二区在线不卡| 高清黄色对白视频在线免费看| 一区二区日韩欧美中文字幕| 男人爽女人下面视频在线观看| 丰满迷人的少妇在线观看| 首页视频小说图片口味搜索 | 日韩一区二区三区影片| 美女午夜性视频免费| 午夜福利在线免费观看网站| 欧美精品一区二区大全| 亚洲五月婷婷丁香| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美精品亚洲一区二区| 又大又黄又爽视频免费| 国产伦人伦偷精品视频| 国产在线一区二区三区精| a 毛片基地| 少妇粗大呻吟视频| 19禁男女啪啪无遮挡网站| 天天躁夜夜躁狠狠久久av| 黄色片一级片一级黄色片| 国产精品一区二区免费欧美 | 自线自在国产av| 人人妻人人爽人人添夜夜欢视频| 麻豆国产av国片精品| www.熟女人妻精品国产| 国产亚洲av高清不卡| 青春草视频在线免费观看| 99国产综合亚洲精品| 国产黄色免费在线视频| 欧美日韩视频精品一区| 午夜福利免费观看在线| videos熟女内射| 在线观看免费高清a一片| 在线观看免费午夜福利视频| 亚洲中文av在线| 大陆偷拍与自拍| 欧美人与性动交α欧美精品济南到| 欧美精品av麻豆av| 丰满饥渴人妻一区二区三| 成人午夜精彩视频在线观看| 最新的欧美精品一区二区| 亚洲av男天堂| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品一区二区在线观看99| 成年动漫av网址| 在线观看免费视频网站a站| 国产高清视频在线播放一区 | 人人妻人人澡人人看| 一级片免费观看大全| 国产色视频综合| 国产片内射在线| 下体分泌物呈黄色| 欧美亚洲日本最大视频资源| 亚洲av在线观看美女高潮| 国产日韩一区二区三区精品不卡| 亚洲欧洲国产日韩| 天天操日日干夜夜撸| 老汉色av国产亚洲站长工具| 99国产精品免费福利视频| 女人久久www免费人成看片| 国产黄色免费在线视频| 成人18禁高潮啪啪吃奶动态图| 国产亚洲av高清不卡| 男人添女人高潮全过程视频| 国产精品九九99| 国产精品av久久久久免费| 免费看不卡的av| 一级毛片我不卡| 亚洲中文字幕日韩| 国产精品一区二区在线不卡| 日日摸夜夜添夜夜爱| 激情视频va一区二区三区| 亚洲久久久国产精品| 欧美成狂野欧美在线观看| 国产成人av激情在线播放| 亚洲精品第二区| 新久久久久国产一级毛片| 18禁国产床啪视频网站| 亚洲三区欧美一区| 99九九在线精品视频| 婷婷丁香在线五月| 婷婷色麻豆天堂久久| 99精国产麻豆久久婷婷| 桃花免费在线播放| 国产精品99久久99久久久不卡| av在线老鸭窝| 精品少妇黑人巨大在线播放| 国产真人三级小视频在线观看| 国产亚洲一区二区精品| av国产久精品久网站免费入址| 国产激情久久老熟女| 精品欧美一区二区三区在线| 国产高清videossex| 亚洲欧美中文字幕日韩二区| 男女午夜视频在线观看| 欧美变态另类bdsm刘玥| 亚洲一区中文字幕在线| 久久 成人 亚洲| 国产精品一区二区免费欧美 | 18在线观看网站| 美女扒开内裤让男人捅视频| 热99国产精品久久久久久7| 亚洲精品国产av蜜桃| 精品人妻一区二区三区麻豆| 久久久久久人人人人人| 国产老妇伦熟女老妇高清| 国产日韩一区二区三区精品不卡| 欧美精品人与动牲交sv欧美| 久久精品成人免费网站| 如日韩欧美国产精品一区二区三区| 日本一区二区免费在线视频| 一级,二级,三级黄色视频| 一二三四社区在线视频社区8| 国产精品香港三级国产av潘金莲 | 亚洲欧美日韩另类电影网站| 考比视频在线观看| 精品国产一区二区久久| 精品一区在线观看国产| 夜夜骑夜夜射夜夜干| 一边摸一边抽搐一进一出视频| 国产一区亚洲一区在线观看| 在线精品无人区一区二区三| 久久毛片免费看一区二区三区| 久久人人97超碰香蕉20202| 国产精品久久久久久精品古装| 美女中出高潮动态图| 一级,二级,三级黄色视频| 2018国产大陆天天弄谢| 欧美日韩av久久| 免费看不卡的av| 精品福利永久在线观看| 国产色视频综合| 伦理电影免费视频| 精品久久久精品久久久| 久久久久国产精品人妻一区二区| 97精品久久久久久久久久精品| 咕卡用的链子| 2018国产大陆天天弄谢| 亚洲激情五月婷婷啪啪| 纵有疾风起免费观看全集完整版| 看十八女毛片水多多多| 国产成人一区二区三区免费视频网站 | 18禁国产床啪视频网站| a级片在线免费高清观看视频| 午夜免费男女啪啪视频观看| 麻豆国产av国片精品| 中国国产av一级| 曰老女人黄片| avwww免费| 欧美变态另类bdsm刘玥| 狂野欧美激情性bbbbbb| 黄片播放在线免费| 日本色播在线视频| 99国产综合亚洲精品| av网站免费在线观看视频| 少妇人妻久久综合中文| 亚洲成国产人片在线观看| 国产亚洲欧美精品永久| 这个男人来自地球电影免费观看| 成人免费观看视频高清| 成年动漫av网址| 亚洲国产看品久久| 80岁老熟妇乱子伦牲交| 日韩一卡2卡3卡4卡2021年| 亚洲国产最新在线播放| 日韩电影二区| 亚洲欧美成人综合另类久久久| 国产亚洲欧美精品永久| 69精品国产乱码久久久| 日韩熟女老妇一区二区性免费视频| kizo精华| 亚洲精品美女久久av网站| 在线观看一区二区三区激情| 在线观看人妻少妇| 91麻豆精品激情在线观看国产 | 国产精品一区二区在线不卡| 一级毛片电影观看| 首页视频小说图片口味搜索 | 9热在线视频观看99| 夜夜骑夜夜射夜夜干| 亚洲国产成人一精品久久久| 视频区欧美日本亚洲| 1024视频免费在线观看| √禁漫天堂资源中文www| 天堂8中文在线网| 国产高清国产精品国产三级| 国产高清videossex| 国产亚洲午夜精品一区二区久久| 日本vs欧美在线观看视频| 黄片播放在线免费| 欧美精品亚洲一区二区| av在线app专区| 欧美97在线视频| 久久久久久久大尺度免费视频| 亚洲欧洲日产国产| 国产亚洲欧美精品永久| 亚洲欧美日韩高清在线视频 | 如日韩欧美国产精品一区二区三区| 精品久久蜜臀av无| 国产日韩一区二区三区精品不卡| 男男h啪啪无遮挡| 看免费av毛片| 咕卡用的链子| 亚洲精品日本国产第一区| 超碰成人久久| 欧美+亚洲+日韩+国产| 中文字幕人妻丝袜制服| 人人妻人人添人人爽欧美一区卜| 天天躁狠狠躁夜夜躁狠狠躁| 欧美日韩一级在线毛片| 免费在线观看日本一区| 久久久久国产精品人妻一区二区| 免费少妇av软件| 免费少妇av软件| 亚洲伊人色综图| 在线观看www视频免费| 激情视频va一区二区三区| 三上悠亚av全集在线观看| 国产一卡二卡三卡精品| 免费女性裸体啪啪无遮挡网站| 国产一区二区三区综合在线观看| 一级黄片播放器| 男人爽女人下面视频在线观看| 精品国产超薄肉色丝袜足j| 9191精品国产免费久久| 欧美日韩视频精品一区| 亚洲激情五月婷婷啪啪| 亚洲熟女精品中文字幕| 老司机午夜十八禁免费视频| 欧美中文综合在线视频| 两个人看的免费小视频| 国产一区二区三区av在线| 免费观看a级毛片全部| 国产成人精品久久久久久| 国产色视频综合| 人人澡人人妻人| 成年av动漫网址| 90打野战视频偷拍视频| 麻豆国产av国片精品| 国产男人的电影天堂91| 免费久久久久久久精品成人欧美视频| 精品欧美一区二区三区在线| 首页视频小说图片口味搜索 | 成人国语在线视频| 欧美精品av麻豆av| 日日夜夜操网爽| 欧美乱码精品一区二区三区| 国产日韩欧美视频二区| 国产淫语在线视频| 十八禁高潮呻吟视频| 成年动漫av网址| 国产精品二区激情视频| 欧美黄色片欧美黄色片| 国产三级黄色录像| 最新的欧美精品一区二区| 国产片内射在线| 丝袜脚勾引网站| 国产xxxxx性猛交| 午夜免费观看性视频| 美女高潮到喷水免费观看| 最近最新中文字幕大全免费视频 | 爱豆传媒免费全集在线观看| av片东京热男人的天堂| netflix在线观看网站| 黄片小视频在线播放| 男的添女的下面高潮视频| 美女扒开内裤让男人捅视频| 免费日韩欧美在线观看| 婷婷色综合www| 一边摸一边抽搐一进一出视频| 大型av网站在线播放| 日韩伦理黄色片| 午夜久久久在线观看| 欧美乱码精品一区二区三区| 欧美+亚洲+日韩+国产| 久久人妻福利社区极品人妻图片 | 午夜精品国产一区二区电影| 你懂的网址亚洲精品在线观看| 久久久精品国产亚洲av高清涩受| 狂野欧美激情性xxxx| 9色porny在线观看| 啦啦啦视频在线资源免费观看| 手机成人av网站| 国产精品久久久久久精品古装| 无遮挡黄片免费观看| 欧美日本中文国产一区发布| 最黄视频免费看| 精品少妇久久久久久888优播| 久久人人爽人人片av| 日本91视频免费播放| 91成人精品电影| 国产精品香港三级国产av潘金莲 | 免费女性裸体啪啪无遮挡网站| 久久免费观看电影| av网站在线播放免费| 一级毛片我不卡| 黄网站色视频无遮挡免费观看| 亚洲精品在线美女| 在线观看人妻少妇| 亚洲第一青青草原| 最近中文字幕2019免费版| 亚洲熟女精品中文字幕| 校园人妻丝袜中文字幕| 另类精品久久| 天天影视国产精品| videosex国产| 又大又爽又粗| 婷婷色综合大香蕉| 亚洲av国产av综合av卡| 最新在线观看一区二区三区 | 国产无遮挡羞羞视频在线观看| 色94色欧美一区二区| 欧美日韩视频高清一区二区三区二| av在线app专区| 一区二区av电影网| 一区二区三区激情视频| 伊人亚洲综合成人网| 91精品国产国语对白视频| 欧美人与善性xxx| av片东京热男人的天堂| 在线 av 中文字幕| av在线app专区| 精品一品国产午夜福利视频| 精品人妻熟女毛片av久久网站| 高清黄色对白视频在线免费看| 国产一区二区在线观看av| 99九九在线精品视频| 中文字幕制服av| 亚洲男人天堂网一区| 大香蕉久久成人网| 99热全是精品| 永久免费av网站大全| 高清欧美精品videossex| 一二三四社区在线视频社区8| 久久人人爽av亚洲精品天堂| 亚洲欧美精品综合一区二区三区| 国产主播在线观看一区二区 | 看免费成人av毛片| 国产女主播在线喷水免费视频网站| 中文字幕人妻熟女乱码| 校园人妻丝袜中文字幕| 新久久久久国产一级毛片| 岛国毛片在线播放| 亚洲免费av在线视频| 五月天丁香电影| 亚洲成人国产一区在线观看 | 日本91视频免费播放| 又大又爽又粗| 亚洲精品久久午夜乱码| 一级黄色大片毛片| 纯流量卡能插随身wifi吗| 女人高潮潮喷娇喘18禁视频| 久久天堂一区二区三区四区| 尾随美女入室| 婷婷色综合大香蕉| 黄色毛片三级朝国网站| 女人高潮潮喷娇喘18禁视频| 色精品久久人妻99蜜桃| 亚洲自偷自拍图片 自拍| 黄色视频不卡| 久久这里只有精品19| 高清欧美精品videossex| 涩涩av久久男人的天堂| 国产精品.久久久| 色精品久久人妻99蜜桃| 少妇人妻 视频| 国产激情久久老熟女| 国产日韩一区二区三区精品不卡| 亚洲精品一卡2卡三卡4卡5卡 | 欧美精品人与动牲交sv欧美| 免费一级毛片在线播放高清视频 | 夜夜骑夜夜射夜夜干| 九草在线视频观看| 91精品国产国语对白视频| a级片在线免费高清观看视频| 在线av久久热| 女性被躁到高潮视频| 亚洲男人天堂网一区| 亚洲人成电影观看| √禁漫天堂资源中文www| 国产成人精品久久二区二区免费| 国产一区二区 视频在线| 午夜福利视频精品| 不卡av一区二区三区| 999精品在线视频| 99国产综合亚洲精品| 欧美另类一区| 久久亚洲精品不卡| 一级毛片 在线播放| 无遮挡黄片免费观看| 操美女的视频在线观看| 男人舔女人的私密视频| 成人国产av品久久久| bbb黄色大片| 亚洲专区中文字幕在线| 永久免费av网站大全| 2018国产大陆天天弄谢| tube8黄色片| 七月丁香在线播放| 三上悠亚av全集在线观看| 免费看十八禁软件| 黄网站色视频无遮挡免费观看| 国产一区二区三区av在线| 午夜视频精品福利| 亚洲免费av在线视频| 男人舔女人的私密视频| 午夜av观看不卡| 午夜两性在线视频| 国产精品九九99| 婷婷色综合大香蕉| 国产亚洲午夜精品一区二区久久| 91九色精品人成在线观看| 欧美日韩黄片免| 国产在线免费精品| 一区在线观看完整版| 又黄又粗又硬又大视频| 久久人人97超碰香蕉20202| 老汉色av国产亚洲站长工具| 男女边吃奶边做爰视频| 国产日韩欧美在线精品| 日韩一卡2卡3卡4卡2021年| 天天躁狠狠躁夜夜躁狠狠躁| 性少妇av在线| 日韩精品免费视频一区二区三区| 国产一区二区三区av在线| 亚洲精品一卡2卡三卡4卡5卡 | 9色porny在线观看| av线在线观看网站| 亚洲av国产av综合av卡| 欧美性长视频在线观看| 久久精品熟女亚洲av麻豆精品| 黑丝袜美女国产一区| 成年美女黄网站色视频大全免费| 久久亚洲国产成人精品v| 十八禁高潮呻吟视频| 久久精品国产亚洲av高清一级| 久久99热这里只频精品6学生| 亚洲国产毛片av蜜桃av| 黄色片一级片一级黄色片| 欧美日韩黄片免| 男女之事视频高清在线观看 | 日韩大码丰满熟妇| 一级片'在线观看视频| 人妻人人澡人人爽人人| 97在线人人人人妻| 亚洲一区二区三区欧美精品| 免费在线观看日本一区| 亚洲图色成人| 国产精品成人在线| 国产老妇伦熟女老妇高清| 国产高清不卡午夜福利| 久久久久网色| 男女免费视频国产| 亚洲熟女毛片儿| 十八禁人妻一区二区| 日韩一区二区三区影片| 美女脱内裤让男人舔精品视频| 91精品国产国语对白视频| 日本欧美视频一区| 欧美精品一区二区免费开放| 精品第一国产精品| 亚洲av综合色区一区| 看十八女毛片水多多多| 777米奇影视久久| 在线 av 中文字幕| 亚洲av日韩精品久久久久久密 | 啦啦啦在线免费观看视频4| 亚洲久久久国产精品| 久久精品国产综合久久久| 国产欧美亚洲国产| 在现免费观看毛片| 国产女主播在线喷水免费视频网站| 国产精品一二三区在线看| 免费看十八禁软件| 国产精品久久久久久精品古装| 久久99一区二区三区| 一本大道久久a久久精品| 久久精品aⅴ一区二区三区四区| 亚洲av国产av综合av卡| 日韩av免费高清视频| 欧美日韩亚洲高清精品| 大话2 男鬼变身卡| 久久综合国产亚洲精品| 亚洲伊人色综图| 中文字幕人妻熟女乱码| 国产老妇伦熟女老妇高清| 国产免费现黄频在线看| 欧美 亚洲 国产 日韩一| 精品免费久久久久久久清纯 | 人人妻,人人澡人人爽秒播 | 久久免费观看电影| 日日摸夜夜添夜夜爱| 叶爱在线成人免费视频播放| 一本综合久久免费| 欧美黑人欧美精品刺激| 亚洲免费av在线视频| 国产成人91sexporn| 国产免费视频播放在线视频| 国产亚洲av片在线观看秒播厂| 嫁个100分男人电影在线观看 | av有码第一页| 少妇猛男粗大的猛烈进出视频| 人成视频在线观看免费观看| 97人妻天天添夜夜摸| 久9热在线精品视频| 国产爽快片一区二区三区| 日韩av免费高清视频| 人妻 亚洲 视频| 免费一级毛片在线播放高清视频 | 高清不卡的av网站| 亚洲欧美精品自产自拍| 亚洲一区中文字幕在线| av福利片在线| 午夜精品国产一区二区电影| 亚洲精品国产一区二区精华液| 久久精品熟女亚洲av麻豆精品| 校园人妻丝袜中文字幕| 亚洲中文av在线| 丝袜脚勾引网站| 一级片免费观看大全| 欧美+亚洲+日韩+国产| 菩萨蛮人人尽说江南好唐韦庄| 在线观看免费视频网站a站| 亚洲精品国产色婷婷电影| 国产精品.久久久| 亚洲欧美中文字幕日韩二区| 免费人妻精品一区二区三区视频| 国产精品99久久99久久久不卡| 只有这里有精品99| av在线老鸭窝| 黑人猛操日本美女一级片| 曰老女人黄片| 99国产精品一区二区蜜桃av | 性高湖久久久久久久久免费观看| 中文乱码字字幕精品一区二区三区| 亚洲国产欧美一区二区综合| 亚洲av成人不卡在线观看播放网 | 高潮久久久久久久久久久不卡| 久久中文字幕一级| 热re99久久精品国产66热6| 亚洲成国产人片在线观看| 男女床上黄色一级片免费看| 在现免费观看毛片| 精品熟女少妇八av免费久了| 91成人精品电影| 黄色a级毛片大全视频| 国产成人精品久久久久久| 国产高清视频在线播放一区 | 欧美日韩国产mv在线观看视频| 国产成人精品无人区| 18禁观看日本| 成人国语在线视频| 捣出白浆h1v1| 国语对白做爰xxxⅹ性视频网站| 欧美精品一区二区免费开放| 国产精品一二三区在线看| 成年女人毛片免费观看观看9 | 中文字幕人妻丝袜制服| 精品国产超薄肉色丝袜足j| 麻豆国产av国片精品| 久久国产精品男人的天堂亚洲| xxxhd国产人妻xxx| 国产亚洲精品久久久久5区| 悠悠久久av| 亚洲免费av在线视频| 亚洲精品美女久久久久99蜜臀 | 亚洲一码二码三码区别大吗| 亚洲精品第二区| 一本色道久久久久久精品综合| 久久狼人影院| 精品国产国语对白av| 中文欧美无线码| 亚洲图色成人| 人人澡人人妻人| 老司机影院毛片| 另类精品久久| 精品国产乱码久久久久久小说| 久久毛片免费看一区二区三区| 亚洲欧美一区二区三区久久| 亚洲激情五月婷婷啪啪| 波多野结衣av一区二区av| 亚洲精品国产色婷婷电影| 天堂中文最新版在线下载| 一级黄色大片毛片| 成年女人毛片免费观看观看9 | 中国美女看黄片| 亚洲国产欧美日韩在线播放| 亚洲欧洲日产国产| 人人妻人人爽人人添夜夜欢视频| 亚洲国产精品成人久久小说| 欧美亚洲日本最大视频资源| 黄片小视频在线播放| 久9热在线精品视频| 啦啦啦啦在线视频资源| 亚洲天堂av无毛| 黑人巨大精品欧美一区二区蜜桃| 69精品国产乱码久久久| 999久久久国产精品视频| 十分钟在线观看高清视频www| 精品人妻1区二区| 十分钟在线观看高清视频www| 男女无遮挡免费网站观看| 免费av中文字幕在线| 国产成人一区二区在线| 亚洲国产精品一区三区|