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

    Study on the Structural Strengthening Design under the Ship-ice Collision Load

    2016-05-16 06:16:05,,,
    船舶力學(xué) 2016年6期
    關(guān)鍵詞:外板文心介質(zhì)

    ,,,

    (1.School of Naval Architecture and Ocean Engineering,Jiangsu University of Science and Technology,Zhenjiang 212003, China;2.Department of Naval Architecture,Ocean&Marine Engineering,University of Strathclyde,Glasgow,UK)

    Study on the Structural Strengthening Design under the Ship-ice Collision Load

    ZHANG Jian1,HE Wen-xin1,YUAN Zhi-ming2,CHEN Cong1

    (1.School of Naval Architecture and Ocean Engineering,Jiangsu University of Science and Technology,Zhenjiang 212003, China;2.Department of Naval Architecture,Ocean&Marine Engineering,University of Strathclyde,Glasgow,UK)

    In order to investigate the water media effect on structure response under ship-ice collision load,this paper presents a study on ship-ice collision considering the effect of water media for the first time.The ship-ice collision is simulated in aqueous medium and the coupling behaviors among the ship,water and ice were considered.The calculation model is applied to various working conditions of ship-ice collision calculation.The parametric study on the thickness of shell and diaphragm plates and the spacing of the frames is then conducted.The damage deformations,collision forces,structural energy absorptions varying with thickness of shell plate,thickness of the diaphragm plate and different spacing of frame are discussed.The capability and contribution of each major component in the shoulder of the ship on the resistance of ice loads are analyzed.This study can be used in improving anti-ice load ability of ship structure.

    ship-ice collision in water media;ice model;structure response; strengthening design

    0 Introduction

    In the recent years,with the discovery and exploitation of oil and gas resources in arctic area as well as the development of Northern Sea Route(NSR),more and more attentions have been paid to the study on ship-ice collision in the field of naval architecture and ocean engineering.

    Ship-ice collision has been a topic of interest among researchers especially those from Canada,Norway,America.Gagnon[1]investigated the ship-ice collision through both the experimental and numerical methods.However,the structural damage was not reported.Wang et al[2]used a coupled model to simulate ice-LNG collision.The floating ice model with inducted ice failure model,Liu et al[3-4]systematically studied the external and internal dynamics of a ship collision.They obtained the pattern of the energy dissipation during the collision by using a simplified method.They also simulated the ship-ice collision by a numerical method.Lee etal[5]conducted a collision test between a cylindrical ice and a grillage structure.Special efforts were put on the failure law of the ice model.Zhang et al[6]used both numerical and experimental methods to investigate the collision between the side grillage and ice to validate the feasibility of the ice material model,as well as the numerical method.

    This paper based on finite element software,numerical simulation of ship-ice collision is conducted in the water media,The discussions are highlighted on the reinforcement of the structures in ice area.This study can be used in improving anti-ice load ability of ship structure design.

    1 Ice numerical model

    The accuracy of the numerical simulation depends on two aspects:the material model of both sides during the collision and method of numerical simulation.There are two kinds of materials mainly involving in ship collision:the steel hull and the ice.The material of hull is very mature already.However,The mechanical property of the ice is very complicated.One example is that ice fails in different modes under different strain rates(Liu et al)[7].At a low strain rate,the deformation of ice can be stimulated by the elastic-plastic theory and the nonrecoverable stage of deformation is dominated by the nonlinear properties.At a high strain rate, the ice deformation tends to be smaller.The stress-strain curve is linear to a large degree and the fragility of the ice takes the dominate role,and the ice is fractured,the failure mode(e.g. large scale and partly fragmentation and bending failure)varies.A transitional zone appears between the ductile zone and the brittle zone and is characterized by plastic deformation,many cracks shear failure and partly fragmentation.Schulson et al[8]investigated this strain rate effect,as shown in Fig.1 and suggested that ice stress tends to rise along with the strain rate until brittle failure appears.

    Fig.1 Ice strength plotted against strain rate

    Based on the elastic-plastic theory,the expression of stress-strain relationship of constitutive equation derived from the yield curve formula,the yield curve is usually expressed by the stress tensor.In the case of the ship-ice collision,the ice is in the three-dimensional stressstate.In 1995,Gagnon and Gammon[9]made an ice triaxial experiment,based on the triaxial experiment data,Derradji-Aouat proposed iceberg elliptical yield closed line equation:

    where τ is the Octahedral shear stress,τmaxis the minor axis of the ellipse,η and λ are the center coordinates of the ellipse,p is the hydrostatic pressure,pcis the major axis of the ellipse.

    The iceberg material model,which is based on the‘Tsai-Wu’yield surface:

    Some of the parameters have been customized,taking strain rate into consideration,the ice strain rate is set as field in the material model of ice.The ice material model as shown in Tab.1.In the reference[10],the findings of study had indicated,the ice material model and the numerical method can be used in numerical simulation for ship-ice collision.

    Tab.1 The material parameters of ice

    2 Numerical simulation of ship-ice collision in the water media

    2.1 Ship-water-ice model

    Recently,added water mass method is adopted by many researchers(Liu et al;Myhre; Wang et al)[3-4,11,2].However,the influence from the fluid resistance is neglected and the added mass is imposed to the ship hull mass.Therefore,the accuracy of the predictions is influenced. This section discusses the application of fluid-structure interaction method in the numerical simulation of ship collision.In order to investigate the water media effect on ship motions andthe fluid field after collision,the hybrid models of ship,water and ice have been established to simulate ship-ice collision in the water environment.

    In order to simulate the infinite computational domain,the nonreflecting boundary condition is used.The coupling effects between shell plate and water media need to be defined in the model calculation. A hollow unit has to be established to guarantee that the waves and spatters are formed outside the computational domain and the gap in the stern can be filled immediately.The material property of the hollow unit is defined as fluid property.The finite element model of ship and water media is shown in Fig.2.

    2.2 Flow distribution of water area

    Fig.3 shows resultant velocity induced from the coupling effect of ship-water-ice.Fore body of the ship displaces water which forms waves on all sides.The advancing direction of the vessel could be modified after the ship-ice collision,which occurs at the fore part of the ship. From the top view we can see that the flow in the portside and starboard is asymmetrical and the movement of the ice has an indirect effect on water media.

    Fig.2 FE Model of ship and water

    Fig.3 Resultant velocity under ship-water-ice coupling effect(top view)

    Fig.4 Transverse velocity of water

    Fig.4 shows the horizontal flow field around the vessel.It can be found that the transverse component of velocity is asymmetry in the fore part of the vessel.The distribution of the transverse velocity in the plane of the vessel’s bottom plate could be influenced by the transverse impact from the ice.The flow at starboard can provide resistance,which can keep the vesselfrom rotation.Meanwhile,the flow at the starboard is squeezed by the vessel,which makes the velocity of the fluid on the starboard opposite to that on the portside.It can be concluded that ship-ice collision can influence the flow field around the vessel.

    According to the analysis above,we can find that the fluid-structure interaction method can predict the ship motion in the water and the effect of transverse impact on flow field around vessel.This method can clearly simulate the ship-ice collision in the water media. Therefore,the fluid-structure interaction method can be used for further study on the anti-collision performance of ship structures under the ice loads.

    3 Study on structure reinforcement design in ice area

    According to the previous results on the comparison between the measurement and the numerical simulation,the simulation method for ship-water-ice model can be determined.This method takes the water media into consideration.In this section,this method will be extended to the study on structure reinforcement design.

    3.1 Ship model

    For ship-ice collision,the collision scheme must be determined first.This paper guarantees the total weight of the hull model,moment of inertia and center of gravity location and other physical properties are consistent with the real ship,then,a fine grid is taken to the area which is going to be contact,and the simplified method is adopted for the structure of the ship which is far away from the collision area.The finite element model of the ship is shown in Fig.5 and Fig.6.

    Fig.5 FE model of the ship

    Fig.6 FE model of the bow and shoulder

    The time step integration is the characteristics of the finite element technology,every step of the calculation required to satisfy the Courant criterion,this paper draws lessons from the domestic and international researchers,control the characteristic length of the hull structure model is about 250 mm,the grid size can ensure the precision of calculation,and can make the calculation efficiency to a certain extent.

    The unit size of the calculation has a significant effect on calculation,unit type also has a significant impact on the result and the computation efficiency.Unit of ship finite element model with shell elements is defined by*SECTION_SHELL in LS-DYNA.Hughes-Liu S/R algorithm is chosen as its element algorithm.

    In the process of ship-ice collision,there exist a lot of contact problems,contact bound-ary is constantly changing,so in the calculation we must continue to search the contact boundary,this paper uses contact erosion(Eroding_surface_to_surface).

    3.2 Collision scheme

    According to results by Zhang et al[6]and the recommended ice shape by DNV(DNV Technical Report)[12],three types of ice shape are simulated.It is found that the rhombic ice possesses the most destructive force against ship.In the present study,we select this type of ice to simulate the ship-ice collision problem.The simplified model and its collision zone are shown in Fig.7.

    Fig.7 Ice model and impact location

    Fig.8 The sketch for the ship-ice collision

    This paper focuses on a tanker navigating in ice region and the advancing speed is fixed at 8 m/s.The ice impacts the portside of the ship with the speed of 2 m/s at the angle of 30°. The impact angle here refers to the angle between top surface of cubic ice and base plane. The sketch for the ship-ice collision is shown in Fig.8.

    There are mainly three parameters which could affect the reinforcement design:the thickness of the shell plate and diaphragm plate,and the spacing of the frame.The schemes of the impact are described as follows:

    Plan I:Three cases of the thickness of the shell plate are taken into consideration.The thickness of cases A,B and C is 19 mm,22 mm and 25 mm,respectively.Since the collision region is limited to ice zone only,the shell plate should be only strengthened in ice belt region.The vertical strengthening range is from 2.7 m to 8.8 m.The thickness beyond the ice region is fixed at 12 mm.

    Plan II:Three cases of the thickness of the diaphragm plate are taken into consideration. The thickness of cases A,B and C is 19 mm,22 mm and 25 mm,respectively.

    Plan III:Two cases of the spacing of the frame are taken into consideration.The spacing of cases A and B is 0.8m and 2.4m,respectively.

    3.3 Effects of shell plate thickness on ship-ice collision

    3.3.1 Collision damage deformation

    Through the comparisons in Fig.9,it can be found that the damage deformation of shell plate is mainly caused by tearing and membrane stretching strength.At the beginning,the large stress occurs at the collision area on the shell plate.And then the fracture occurs when the deformation goes beyond failure strain.At the same time,the ship hull has some velocity. The ice corner inserts into the inner structures and cuts the shell plate.The diaphragm plate,platform plate and some side shell frames yields,and then subjected to large deformation,instability and fracture.In terms of the damage level of the shell plate and inner structures,the most serious damage occurs in Case A,followed by Cases B and C.

    Fig.9 Damage deformation of hull structure under three cases in Plan I

    3.3.2 Collision force

    Collision force exhibits nonlinearity and variability as the collision progresses.Qualitatively the process is as follows:(i)Collision force increases as the corner of the ice piece crushes against the hull and the contact area increases,(ii)Process continues,force increases as the shell plates and internal structure deform elastically,then plastically, (iii)Eventually hull structure fails(tearing and buckling)resulting in a force decrease,(iv)As the ship continues to advance,new structure comes into contact with the ice,force increases until the structure fails and the process repeats itself,(v)Further failure of the structure stops.

    Fig.10 shows the time history of the collision force at different shell plate thickness.It can be found that at the initial loading stage,the deformation of shell plate increases gradually.The peak value occurs at the moment when the shell plate is being penetrated through.After that, there is a short period of unloading process.Then,the ice keeps cutting the inner structures in the portside of the hull during the relative motion between the ship and ice.As a result,the fracture of the strong components will cause the increase of the collision force,which can result in a violate fluctuation of the time history of the collision force.It can also be found that the collision forces are significantly influenced by the thickness of the shell plate.As the time is smaller than 0.6 s,although the curve of collision force represents the strong mode of vibration,it increases along with the thickness of shell plate generally.But as the time is greater than 0.6 s,the impact force will present to be low level earlier as the thickness increases.Besides,we find that the maximum collision forces are also influenced by the thickness of theshell plate.As the thickness increases,the maximum value keeps an increase trend.The damage moment will be delayed accordingly.The maximum collision force and the cracking moment are presented in Tab.2.

    Fig.10 The time history of the collision force of Plan I

    Tab.2 The comparison of the maximum collision force and plate cracking moment under different cases

    3.3.3 Energy absorption

    Fig.11 and Tab.3 show the energy absorption of all the components under different cases. It can be found that the shell plate is the main energy absorption component and it takes up to 50%of the total energy absorption.At the beginning stage,the energy absorption increases quickly due to the influence of the elastic deformation of the structure.As collision continues, the shoulder component deforms and absorbs energy gradually.As a result,the energy absorption will keep increase.As the time becomes greater than 0.8 s,the energy absorption tends to be stable,since the direction of the movement of the ice starts to deviate and the ice breaks away from the hull structures.It can also be found that during the same period of collision time,for shell plate,the percentage of the total energy absorption decreases as the thickness increases.The discrepancies between different components become small.

    Fig.11 Curve of energy absorption under three cases in Plan I

    Tab.3 Energy absorption capacity of structural components

    3.4 Effects of diaphragm plate thickness on ship-ice collision

    From the previous results,we can see the percentage of the energy absorbing by the diaphragm plate is smaller than that of the shell plate.But it is still an important energy absorption component.

    3.4.1 Collision damage deformation

    Through the comparisons in Fig.12,it can be found that in terms of damage level of the shell plate and the internal side frame,the most serious damage occurs in Case A,followed by Cases B and C.Increasing the thickness of diaphragm plate can keep the shell plate from cracking,as well as stop the ice from embedding the inner structures of the ship hull.

    Fig.12 Damage deformation of hull structure under three cases in Plan II

    3.4.2 Collision force

    Fig.13 shows the time history of the collision force at different diaphragm plate thickness. Tab.4 presents the maximum collision force in different cases.It can be found that the maximum collision forces are influenced by the thickness of the diaphragm plate.As the thickness increases,the maximum value keeps an increase trend.This suggests that increasing the thickness of diaphragm plate can improve shell plate’s capability of resistance to deformation.Besides,at the same impact height,the structure possesses a better anti-collision capacity if it can endure a higher collision force.

    Tab.4 The comparison of the maximum collision force under different cases

    Fig.13 The time history of the collision force of Plan II

    3.4.3 Energy absorption

    Fig.14 and Tab.5 show the energy absorption of all the components under different cases. From the curve it can be found that as the thickness of diaphragm plate increases,the percentage of the energy absorption of the shell plate decreases,while the percentage of the diaphragmplate keeps a increase trend.The discrepancies between different components become smaller.Besides,the total energy absorption decreases as the thickness increases.As the diaphragm plate thickness increases from 12 mm to 22 mm,the total percentage of energy absorption decreases by 8%.In addition,based on the same increase of the mass of the diaphragm plate, the energy absorption of Case B decreases by 3.8%compared with Case A and for Case C,it drops by 4.2%compared with Case B.It can be concluded that in the present case study,the thickness of 17 mm is the optical design due to the collision damage deformation and energy absorption percentage.

    Fig.14 Curve of energy absorption under different cases in Plan II

    Tab.5 Energy absorption capacity of structural components

    3.5 Effects of frame spacing on ship-ice collision

    The frame can strengthen the shell plate with respect to resist the collision and squeeze from the ice.It can be found from energy absorption results that the frame is also a key energy absorption component,which takes about 15%of the total energy absorption.

    3.5.1 Collision damage deformation

    From Fig.15,we can see the damage of hull structure varies with the spacing of the frame. The crevasse area in Case B is bigger than that in Case A.For Case A and Case B,the spacing of the frame is 0.8 m and 2.4 m.Since the spacing of the frame of Case A is much smaller,the capacity of support from the frame could be increased.At the beginning of the collision, the surface of the ship is squeezed.Most area of side frames and stringers start to deform and fail.Since the spacing of the frame in Case A is small,the number of frames in the same impact area is larger and more collision energy will be absorbed.As a result,the crack of theshell plate and the failure of the other components can be restrained.As the collision continues, the shell plates start to crack and the failure mode transforms into cutting.The frames,diaphragm plates,platform plates,side longitudinals and other components will resist ice intrusion together.

    Fig.15 Damage deformation of hull structure under two cases in Plan III

    3.5.2 Collision force

    The time history of the collision force of Plan III is shown in Fig.16.It can be found that the peak value of collision force in Case A is higher than that in Case B.As the frame spacing increases,the supportability of the shell plate is reduced.Therefore,the shell plate could crack more easily.The maximum collision force and the cracking moment are presented in Tab.6.

    Fig.16 Time history of collision force of Plan III

    Tab.6 The comparison of the maximum collision force and plate cracking moment under different cases

    3.5.3 Energy absorption

    Fig.17 and Tab.7 show the energy absorption of all the components at different cases.It can be found that as the frame spacing decreases,the energy absorption percentage of the side frames increases and the discrepancies between different components become small.Besides, the total energy absorption increases as the frame spacing decreases.Compared Case B with Case A,the frame space increases from 0.8 m to 2.4 m,but the total energy absorption decreases by 18.5%.Therefore,when conducting the structural anti-collision design in the iceregion,the increase of the spacing of the frames can improve the total energy absorption.

    Fig.17 Curve of energy absorption under two cases in Plan III

    Tab.7 Energy absorption capacity of structural components

    In summary,the thickness of the shell plate was taken 19mm,22 mm and 25 mm,the thickness of the shell plate increases,the structure crashworthiness performance is better,the damage deformation is smaller;shell plate thickness variations have a significant effect on the impact force,in general,maximum impact force increases with the increasing of the thickness; shell plate is the main energy absorbing component,the structure of the total energy absorption is inversely proportional to the increase of plating thickness,It can be concluded that in the present case study,the thickness of 25 mm is a better choice.The thickness of the diaphragm plate was taken 12 mm,17 mm and 22 mm,with the increase of the thickness of diaphragm plate,the range of damage of the shell plate and the internal skeleton will decreases;The biggest impact value is proportional to the thickness of diaphragm plate size;The total energy absorption of the ship structure is inversely proportional to the increase of the diaphragm plate thickness,it can be concluded that in the present case study,the thickness of 17 mm is a better choice.Frame spacing was taken 0.8 m and 2.4 m,the smaller the frame spacing,the more frames,the damage of other components will be reduced;when the frame spacing is 0.8 m,the peak value of pounding force will be larger;as the frame spacing decreases,the overall structure absorption can be increased,it can be concluded that in the present case study,0.8 m is a better choice.

    4 Conclusions

    From the results presented in this paper,the following conclusions can be drawn:

    (1)The fluid-structure interaction method is an effective method to simulate the ship-ice collision in the water media.It can reflect the effect of the water on the ship motions.

    (2)Shell plate is the main energy absorbing component.It takes up around 50%of the total energy absorption,followed by the side supporting components,such as diaphragm plate, platform plate,side shell frame and side longitudinal.

    (3)The discrepancies of the energy absorption capacity among the diaphragms plate,platform plate,side longitudinal stiffness,side shell frame and shell plate can be reduced by increasing the shell plate and diaphragm plate thickness,as well as reducing the spacing of frames.

    (4)In the design of the anti-collision structures,we need to consider the relationship between the damage deformation and energy absorption,so as to find the optimal thickness of the shell plate in ice region.In the premise of the limitation of the structural weight,we can increase the thickness of the diaphragm plate and reduce the spacing of the frames to achieve a smaller deformation and better energy absorption.

    [1]Gagnon R.First results of numerical simulations of bergy-bit collisions with the CCGS terry fox icebreaker[C]//In:Proc. of the 18th Int.Ass.for the Hydo-Enviroment Eng.and Research(IAHR)on ice,Aug.28-Sep.1 2006,Sapporo,Japan, 2006.

    [2]Wang B,Yu H C.Ship and ice collison modeling and strength evaluation of LNG ship structure[C]//In:OMAE.shanghai, China,2008.

    [3]Liu Z H,Amdahl J,L?set S.A parametric study on the external mechanics of ship/iceberg collision[C]//In:Proceedings of the ASME 2010 29th International Conference on Ocean Offshore and Arctic Engineering.Shanghai,China,2010(4): 741-749.

    [4]Liu Z H,Amdahl J,L?set S.Integrated numerical analysis of an iceberg collision with a foreship structure[J].Marine Structures,2011(24):377-395.

    [5]Lee S G,Lee J S.Structural safety assessment in membrane-type CCS in LNG under iceberg collision[C].In:Royal Institution of Naval Architects International Conference-International Conference on Ship and Offshore Technology:Ice Class Vessels,2009:69-81.

    [6]Zhang J,Wan Z H,Chen C.Research on structure dynamic response of bulbous bow in ship-ice collision load[J].Journal of Ship Mechanics,2014,18(1):106-114.

    [7]Liu Z H,Amdahl J,L?set S.Plasticity based material modelling of ice and its application to ship-iceberg impacts[J]. Cold Region Science and Technology,2011,65:326-334.

    [8]Schulson E M.Brittle failure of ice[J].Engineering Fracture Mechanics,2011(68):1839-1887.

    [9]Gagnon R E,Gammon P H.Triaxial experiments on iceberg and glaciers ice[J].Glaciology,1995,41(1):39-54.

    [10]Zhang Jian,Wan Zhengquan,Huang Jinhao,Yin Qun.Research on numerical simulation and model test of collision between side grillage and icebergs[J].Journal of Ship Mechanics,2014,18(4):424-433.

    [11]Myhre S A.Analysis of accidental iceberg impacts with membrane tank LNG carriers[D].The Master Dissertation of Norwegian University of Science and Technology,2010:121-135.

    [12]DNV Technical Report 2006-0672.Ice Collision Scenario[R].2006.

    船冰碰撞載荷下船舶結(jié)構(gòu)加強(qiáng)方案研究

    張 健1,何文心1,元志明2,陳 聰1

    (1.江蘇科技大學(xué) 船舶與海洋工程學(xué)院,江蘇 鎮(zhèn)江212003;2.英國斯特拉斯克萊德大學(xué),英國,格拉斯哥)

    為了探究水介質(zhì)對(duì)船—冰碰撞結(jié)構(gòu)響應(yīng)的影響,文章首次對(duì)考慮水介質(zhì)中的船—冰碰撞問題進(jìn)行了研究。模擬了船舶與冰體在水介質(zhì)中的碰撞場(chǎng)景,研究船—水—冰三者共同耦合作用對(duì)船—冰碰撞的影響。將該計(jì)算方法應(yīng)用于船—冰碰撞的多種工況計(jì)算,對(duì)比分析了增加外板厚度、橫隔板厚度以及肋骨間距等多種加強(qiáng)方案對(duì)船體結(jié)構(gòu)響應(yīng)的差異。揭示了碰撞區(qū)域的損傷變形、碰撞力、結(jié)構(gòu)吸能隨外板厚度、橫隔板厚度和肋骨間距的變化而變化的規(guī)律,分析了船舶肩部各主要構(gòu)件對(duì)于抵抗冰載荷作用的能力及貢獻(xiàn)。所得結(jié)論對(duì)于進(jìn)行冰區(qū)船舶結(jié)構(gòu)設(shè)計(jì)具有參考作用。

    水介質(zhì)下船冰碰撞;冰模型;結(jié)構(gòu)響應(yīng);加固設(shè)計(jì)

    U661.4

    A

    張 ?。?977-),男,博士,江蘇科技大學(xué)副教授,研究生導(dǎo)師;

    U661.4

    A

    10.3969/j.issn.1007-7294.2016.06.008

    1007-7294(2016)06-0722-14

    何文心(1992-)女,江蘇科技大學(xué)碩士研究生;

    元志明(1985-),男,博士,英國斯特拉斯克萊德大學(xué)講師,研究生導(dǎo)師;

    陳 聰(1990-),男,江蘇科技大學(xué)碩士研究生。

    Received date:2016-03-30

    Founded item:Supported by National Natural Science Foundation of China(Grant No.51579121);Jiangsu Provincial Department of Education Natural Science Foundation of Major Project(Grant No.15KJA580002); Ship Industry Pre-research Foundation(13J1.5.1)

    Biography:ZHANG Jian(1977-),male,Ph.D.Associate professor of Jiangsu University of Science and Technology,

    E-mail:justzj@126.com;YUAN Zhi-ming(1985-),male,Ph.D.lecturer/tutor of University of Strathclyde,UK.

    猜你喜歡
    外板文心介質(zhì)
    信息交流介質(zhì)的演化與選擇偏好
    淬火冷卻介質(zhì)在航空工業(yè)的應(yīng)用
    側(cè)圍外板轉(zhuǎn)角深拉伸起皺缺陷研究
    模具制造(2019年10期)2020-01-06 09:13:02
    側(cè)圍外板尾燈處排料困難的解決方案
    模具制造(2019年10期)2020-01-06 09:13:00
    汽車側(cè)圍外板A柱起皺缺陷分析及處理方法
    模具制造(2019年4期)2019-06-24 03:36:42
    冬天來啦
    文心雜記
    請(qǐng)不斷修煉你的“文心”
    中國篆刻(2017年5期)2017-07-18 11:09:30
    某型快艇船體外板腐蝕點(diǎn)焊補(bǔ)工藝研究
    欧美乱妇无乱码| 色哟哟哟哟哟哟| 国产免费av片在线观看野外av| 久久99一区二区三区| 91麻豆精品激情在线观看国产 | 精品无人区乱码1区二区| 免费一级毛片在线播放高清视频 | 91在线观看av| ponron亚洲| 日韩av在线大香蕉| 每晚都被弄得嗷嗷叫到高潮| 亚洲一区中文字幕在线| 69精品国产乱码久久久| 在线永久观看黄色视频| 亚洲午夜理论影院| 性欧美人与动物交配| 亚洲熟女毛片儿| 国产一区二区三区在线臀色熟女 | 9191精品国产免费久久| av在线播放免费不卡| 18禁黄网站禁片午夜丰满| 欧美国产精品va在线观看不卡| 99在线人妻在线中文字幕| 中文亚洲av片在线观看爽| 欧美日韩精品网址| 国内久久婷婷六月综合欲色啪| 麻豆av在线久日| 日本黄色日本黄色录像| 中文亚洲av片在线观看爽| 久久精品成人免费网站| 91九色精品人成在线观看| 亚洲美女黄片视频| 在线观看日韩欧美| 欧美精品亚洲一区二区| 午夜福利免费观看在线| xxx96com| 一a级毛片在线观看| 99精品久久久久人妻精品| 日韩精品中文字幕看吧| 成人免费观看视频高清| 大码成人一级视频| 99国产极品粉嫩在线观看| 丝袜人妻中文字幕| 亚洲欧美一区二区三区久久| 日韩三级视频一区二区三区| 女人被狂操c到高潮| 亚洲全国av大片| 桃红色精品国产亚洲av| 国产成+人综合+亚洲专区| 在线免费观看的www视频| 久久久久国产一级毛片高清牌| 精品国产亚洲在线| av网站免费在线观看视频| 香蕉久久夜色| 女生性感内裤真人,穿戴方法视频| 91麻豆av在线| 久久久久久人人人人人| 一区二区三区精品91| 国产激情欧美一区二区| 亚洲成人国产一区在线观看| 国产单亲对白刺激| 91麻豆av在线| 成人特级黄色片久久久久久久| 国产精品乱码一区二三区的特点 | 精品国产乱码久久久久久男人| ponron亚洲| 国产精品九九99| 在线av久久热| 日韩精品青青久久久久久| 桃红色精品国产亚洲av| 久久国产精品影院| 高清黄色对白视频在线免费看| 性少妇av在线| 淫妇啪啪啪对白视频| 久久精品国产亚洲av香蕉五月| 精品欧美一区二区三区在线| 国产亚洲精品第一综合不卡| 亚洲国产欧美一区二区综合| 不卡av一区二区三区| 99精国产麻豆久久婷婷| 亚洲欧洲精品一区二区精品久久久| 久久精品aⅴ一区二区三区四区| 90打野战视频偷拍视频| 搡老乐熟女国产| 在线观看一区二区三区| 久久中文字幕一级| 欧美黄色淫秽网站| 亚洲三区欧美一区| 日韩欧美一区二区三区在线观看| 他把我摸到了高潮在线观看| 色婷婷久久久亚洲欧美| 深夜精品福利| 国产av一区在线观看免费| 久久精品91蜜桃| √禁漫天堂资源中文www| 少妇的丰满在线观看| 日韩av在线大香蕉| 午夜福利一区二区在线看| 在线观看66精品国产| 91精品国产国语对白视频| 欧美激情高清一区二区三区| 精品久久久精品久久久| 欧美日韩乱码在线| 精品福利观看| 90打野战视频偷拍视频| 国产成人啪精品午夜网站| 少妇被粗大的猛进出69影院| 高清在线国产一区| 天堂√8在线中文| 欧美黑人欧美精品刺激| 久久婷婷成人综合色麻豆| 成人精品一区二区免费| 最新在线观看一区二区三区| 免费观看人在逋| 美女大奶头视频| 精品欧美一区二区三区在线| 国产精品国产av在线观看| 多毛熟女@视频| 欧美激情 高清一区二区三区| 亚洲情色 制服丝袜| 亚洲精品美女久久av网站| 久久久国产一区二区| 亚洲精品一区av在线观看| aaaaa片日本免费| 欧美在线黄色| www.熟女人妻精品国产| 91麻豆精品激情在线观看国产 | 国产激情欧美一区二区| 免费高清在线观看日韩| 女人爽到高潮嗷嗷叫在线视频| 极品教师在线免费播放| 日韩免费高清中文字幕av| 国产有黄有色有爽视频| 亚洲午夜精品一区,二区,三区| 丝袜人妻中文字幕| 好男人电影高清在线观看| 三级毛片av免费| 国产精品自产拍在线观看55亚洲| 久久中文看片网| 国产一区二区在线av高清观看| 久久久久国产精品人妻aⅴ院| 国产成人精品久久二区二区91| 国产av在哪里看| 99re在线观看精品视频| 一区二区三区激情视频| 国产视频一区二区在线看| 夫妻午夜视频| 中文字幕色久视频| 亚洲avbb在线观看| 天堂√8在线中文| 又黄又爽又免费观看的视频| 日韩高清综合在线| 18禁观看日本| 久久久国产一区二区| 国产成年人精品一区二区 | 欧美国产精品va在线观看不卡| 1024视频免费在线观看| 欧美 亚洲 国产 日韩一| 亚洲精华国产精华精| 身体一侧抽搐| 欧美黄色淫秽网站| av欧美777| 国产激情欧美一区二区| 国产男靠女视频免费网站| 黑人猛操日本美女一级片| 久久香蕉国产精品| 欧美激情久久久久久爽电影 | 国产亚洲欧美98| 国产91精品成人一区二区三区| 亚洲第一欧美日韩一区二区三区| 精品一品国产午夜福利视频| 脱女人内裤的视频| 国产成人精品无人区| 一进一出抽搐gif免费好疼 | 久久这里只有精品19| 岛国视频午夜一区免费看| 亚洲久久久国产精品| 国产午夜精品久久久久久| cao死你这个sao货| 18禁裸乳无遮挡免费网站照片 | 19禁男女啪啪无遮挡网站| 欧美 亚洲 国产 日韩一| 嫁个100分男人电影在线观看| 村上凉子中文字幕在线| 国产av又大| 午夜影院日韩av| 9色porny在线观看| 51午夜福利影视在线观看| 天堂俺去俺来也www色官网| 亚洲成国产人片在线观看| 一a级毛片在线观看| 人妻丰满熟妇av一区二区三区| 欧美日韩乱码在线| 欧美日韩黄片免| 欧美中文综合在线视频| 欧美黑人欧美精品刺激| 久久久久久久久中文| 1024视频免费在线观看| 国产av在哪里看| 久久人妻av系列| 男女高潮啪啪啪动态图| 性欧美人与动物交配| 亚洲一卡2卡3卡4卡5卡精品中文| 深夜精品福利| 欧美在线黄色| 韩国精品一区二区三区| 精品一区二区三区av网在线观看| а√天堂www在线а√下载| 最好的美女福利视频网| 久久青草综合色| 日韩人妻精品一区2区三区| 精品人妻1区二区| 国产成人av激情在线播放| 亚洲成av片中文字幕在线观看| 亚洲国产欧美网| 老汉色av国产亚洲站长工具| 黑人巨大精品欧美一区二区蜜桃| 欧美日韩黄片免| 欧美日韩精品网址| 在线观看免费高清a一片| 天堂影院成人在线观看| 国产午夜精品久久久久久| 国产成人欧美| 黑人猛操日本美女一级片| 久久天堂一区二区三区四区| 高清av免费在线| av在线天堂中文字幕 | 黄色片一级片一级黄色片| 曰老女人黄片| 老司机午夜福利在线观看视频| 国产精品爽爽va在线观看网站 | 日日摸夜夜添夜夜添小说| 韩国av一区二区三区四区| 99久久99久久久精品蜜桃| 在线观看舔阴道视频| 精品久久久久久电影网| 首页视频小说图片口味搜索| 婷婷精品国产亚洲av在线| 嫁个100分男人电影在线观看| 精品国产一区二区三区四区第35| 国产亚洲欧美精品永久| 丝袜美腿诱惑在线| 久久人人精品亚洲av| 色综合婷婷激情| 波多野结衣一区麻豆| 精品久久久久久,| 亚洲欧美一区二区三区久久| 老鸭窝网址在线观看| 欧美日韩黄片免| 怎么达到女性高潮| 天堂中文最新版在线下载| 黄色丝袜av网址大全| 在线观看午夜福利视频| 美女午夜性视频免费| 久久国产亚洲av麻豆专区| 如日韩欧美国产精品一区二区三区| 日韩欧美一区视频在线观看| 中文字幕精品免费在线观看视频| 亚洲欧美激情综合另类| 国产精品自产拍在线观看55亚洲| 午夜免费成人在线视频| 欧美大码av| 桃红色精品国产亚洲av| 岛国视频午夜一区免费看| 国产单亲对白刺激| 国产有黄有色有爽视频| 亚洲熟女毛片儿| 国产精品一区二区三区四区久久 | 中文字幕另类日韩欧美亚洲嫩草| 激情视频va一区二区三区| 多毛熟女@视频| 视频在线观看一区二区三区| 午夜福利一区二区在线看| 亚洲熟女毛片儿| 50天的宝宝边吃奶边哭怎么回事| 亚洲五月婷婷丁香| 免费在线观看完整版高清| 亚洲精品在线美女| 法律面前人人平等表现在哪些方面| 久久久久亚洲av毛片大全| 多毛熟女@视频| 女生性感内裤真人,穿戴方法视频| 国产成人欧美| 丝袜美足系列| 后天国语完整版免费观看| 亚洲av成人不卡在线观看播放网| 精品卡一卡二卡四卡免费| 午夜免费激情av| 国产精品九九99| 村上凉子中文字幕在线| 曰老女人黄片| 成人18禁在线播放| 亚洲国产精品一区二区三区在线| 国产乱人伦免费视频| av网站在线播放免费| 又大又爽又粗| 久热这里只有精品99| 亚洲欧美激情在线| 亚洲色图综合在线观看| aaaaa片日本免费| av超薄肉色丝袜交足视频| 他把我摸到了高潮在线观看| 真人做人爱边吃奶动态| 欧美成人午夜精品| 丝袜在线中文字幕| 每晚都被弄得嗷嗷叫到高潮| ponron亚洲| 91字幕亚洲| 黑人巨大精品欧美一区二区蜜桃| 亚洲在线自拍视频| 一进一出抽搐动态| 人人妻人人爽人人添夜夜欢视频| 中文字幕最新亚洲高清| 女警被强在线播放| 国产欧美日韩一区二区三| 亚洲专区国产一区二区| 日日爽夜夜爽网站| av网站在线播放免费| 麻豆国产av国片精品| 午夜福利,免费看| 国产黄a三级三级三级人| 亚洲av日韩精品久久久久久密| av天堂在线播放| 99久久人妻综合| 少妇粗大呻吟视频| 精品国产亚洲在线| 亚洲男人天堂网一区| 五月开心婷婷网| 久久精品aⅴ一区二区三区四区| 午夜视频精品福利| 三上悠亚av全集在线观看| 伊人久久大香线蕉亚洲五| 久热这里只有精品99| 亚洲熟妇中文字幕五十中出 | 亚洲中文字幕日韩| 人人妻,人人澡人人爽秒播| 午夜激情av网站| 99国产精品99久久久久| 精品熟女少妇八av免费久了| 日韩欧美免费精品| 美女高潮到喷水免费观看| 国产av一区二区精品久久| 91字幕亚洲| 亚洲中文字幕日韩| 亚洲国产欧美网| 视频区图区小说| 90打野战视频偷拍视频| 日韩欧美国产一区二区入口| 国产精品1区2区在线观看.| 国产亚洲精品综合一区在线观看 | 久久香蕉激情| 国产精品av久久久久免费| 大型黄色视频在线免费观看| 亚洲人成网站在线播放欧美日韩| 91九色精品人成在线观看| 怎么达到女性高潮| 免费av中文字幕在线| 国产单亲对白刺激| cao死你这个sao货| 亚洲av日韩精品久久久久久密| 悠悠久久av| 欧美激情久久久久久爽电影 | 99国产精品99久久久久| 中文欧美无线码| 日本免费a在线| 成年人免费黄色播放视频| avwww免费| 日本欧美视频一区| 亚洲精品在线观看二区| 动漫黄色视频在线观看| 精品一区二区三区av网在线观看| 午夜久久久在线观看| 真人一进一出gif抽搐免费| 亚洲精品中文字幕一二三四区| 一区二区三区激情视频| 韩国av一区二区三区四区| 亚洲人成电影免费在线| 一夜夜www| 亚洲人成77777在线视频| 在线观看免费视频日本深夜| 亚洲免费av在线视频| 人人妻人人澡人人看| 黄色视频不卡| 婷婷六月久久综合丁香| 91字幕亚洲| 不卡一级毛片| 日韩中文字幕欧美一区二区| 三级毛片av免费| 午夜精品久久久久久毛片777| 两个人免费观看高清视频| 黑人巨大精品欧美一区二区蜜桃| 窝窝影院91人妻| 国产极品粉嫩免费观看在线| 国产在线观看jvid| 女同久久另类99精品国产91| 美女国产高潮福利片在线看| 黑人巨大精品欧美一区二区蜜桃| 19禁男女啪啪无遮挡网站| 精品熟女少妇八av免费久了| 无人区码免费观看不卡| 99在线视频只有这里精品首页| 午夜福利在线观看吧| 国产高清激情床上av| 欧美人与性动交α欧美软件| 成人精品一区二区免费| 欧美最黄视频在线播放免费 | 超碰成人久久| 91精品国产国语对白视频| 久久热在线av| 中文字幕色久视频| 可以免费在线观看a视频的电影网站| 在线天堂中文资源库| 色在线成人网| 18禁观看日本| 国产av在哪里看| 免费看a级黄色片| 在线视频色国产色| 国产99久久九九免费精品| 99在线视频只有这里精品首页| 国产成人影院久久av| 一进一出好大好爽视频| 在线观看一区二区三区| 国产有黄有色有爽视频| 国产成人精品无人区| 18禁黄网站禁片午夜丰满| 黑人猛操日本美女一级片| 黄色毛片三级朝国网站| 夫妻午夜视频| 啪啪无遮挡十八禁网站| 精品人妻在线不人妻| 精品国产一区二区三区四区第35| 999久久久精品免费观看国产| 日本免费一区二区三区高清不卡 | 热99re8久久精品国产| 久久久久国内视频| 久久久国产精品麻豆| 亚洲一码二码三码区别大吗| 狂野欧美激情性xxxx| 激情在线观看视频在线高清| av国产精品久久久久影院| 人妻丰满熟妇av一区二区三区| 欧美中文日本在线观看视频| 欧美不卡视频在线免费观看 | 亚洲专区国产一区二区| 精品久久久精品久久久| 一级a爱视频在线免费观看| 精品国产亚洲在线| 又紧又爽又黄一区二区| 日韩免费av在线播放| 亚洲国产精品sss在线观看 | 男女高潮啪啪啪动态图| 欧美最黄视频在线播放免费 | 久久久久久久久免费视频了| 99精品在免费线老司机午夜| 亚洲va日本ⅴa欧美va伊人久久| 精品少妇一区二区三区视频日本电影| 成人av一区二区三区在线看| 婷婷精品国产亚洲av在线| 黄色视频不卡| 大码成人一级视频| 国产精品一区二区在线不卡| 另类亚洲欧美激情| 国产精品九九99| 久久人妻福利社区极品人妻图片| 麻豆av在线久日| 午夜老司机福利片| 成熟少妇高潮喷水视频| 一级片'在线观看视频| 国产日韩一区二区三区精品不卡| xxxhd国产人妻xxx| 不卡av一区二区三区| 不卡一级毛片| 天天躁狠狠躁夜夜躁狠狠躁| 99精品欧美一区二区三区四区| 极品教师在线免费播放| 午夜免费激情av| 亚洲成国产人片在线观看| 久久人人爽av亚洲精品天堂| 国产精品乱码一区二三区的特点 | 如日韩欧美国产精品一区二区三区| 亚洲成a人片在线一区二区| xxx96com| 女人爽到高潮嗷嗷叫在线视频| 女同久久另类99精品国产91| 国产野战对白在线观看| 亚洲自拍偷在线| 80岁老熟妇乱子伦牲交| 亚洲精品美女久久av网站| 男女下面插进去视频免费观看| tocl精华| 久久伊人香网站| 一级毛片高清免费大全| 在线观看免费高清a一片| 欧美人与性动交α欧美软件| 丰满饥渴人妻一区二区三| 欧美日本亚洲视频在线播放| 久久天堂一区二区三区四区| 大陆偷拍与自拍| 久久精品成人免费网站| 久久久久九九精品影院| 少妇的丰满在线观看| 可以在线观看毛片的网站| 老司机在亚洲福利影院| 老熟妇仑乱视频hdxx| 国产欧美日韩一区二区三| 两性夫妻黄色片| 免费一级毛片在线播放高清视频 | 亚洲av成人一区二区三| 亚洲国产毛片av蜜桃av| 99国产精品一区二区蜜桃av| 国产熟女午夜一区二区三区| 高清av免费在线| 欧美日韩国产mv在线观看视频| 高清欧美精品videossex| 亚洲中文字幕日韩| 成人三级黄色视频| 国产一区二区在线av高清观看| 午夜福利在线观看吧| 在线av久久热| 精品高清国产在线一区| 亚洲熟女毛片儿| 丁香欧美五月| 午夜福利一区二区在线看| 51午夜福利影视在线观看| av福利片在线| 成人亚洲精品av一区二区 | 97超级碰碰碰精品色视频在线观看| 9色porny在线观看| 久久人人97超碰香蕉20202| 国产精品 国内视频| 啦啦啦在线免费观看视频4| 国产深夜福利视频在线观看| 淫秽高清视频在线观看| 久久这里只有精品19| 一本综合久久免费| avwww免费| 69精品国产乱码久久久| 日韩欧美在线二视频| 色婷婷久久久亚洲欧美| 在线观看www视频免费| 亚洲,欧美精品.| 亚洲一区二区三区欧美精品| 婷婷丁香在线五月| 成人永久免费在线观看视频| 国产深夜福利视频在线观看| 人人妻人人爽人人添夜夜欢视频| www.www免费av| 搡老乐熟女国产| 国产免费男女视频| 精品久久久久久,| 叶爱在线成人免费视频播放| 老司机午夜十八禁免费视频| 久久 成人 亚洲| 亚洲成人久久性| 99国产精品一区二区蜜桃av| 亚洲av成人一区二区三| 91精品三级在线观看| 欧美中文日本在线观看视频| 悠悠久久av| 中亚洲国语对白在线视频| 欧美精品啪啪一区二区三区| 欧美黄色片欧美黄色片| 欧美一级毛片孕妇| 伊人久久大香线蕉亚洲五| 国产精品成人在线| 亚洲 欧美一区二区三区| 最好的美女福利视频网| 一区在线观看完整版| 黄色成人免费大全| 免费在线观看完整版高清| 在线观看午夜福利视频| 日韩大尺度精品在线看网址 | 乱人伦中国视频| 一区在线观看完整版| 手机成人av网站| 88av欧美| 高清av免费在线| 日韩中文字幕欧美一区二区| 中文字幕人妻丝袜制服| 人妻丰满熟妇av一区二区三区| 日本 av在线| 日本免费a在线| 成人精品一区二区免费| 国产91精品成人一区二区三区| 一级作爱视频免费观看| 不卡av一区二区三区| 国产精品自产拍在线观看55亚洲| 最近最新中文字幕大全电影3 | 久久草成人影院| 国产亚洲精品久久久久久毛片| a级毛片在线看网站| 一本大道久久a久久精品| 美女高潮喷水抽搐中文字幕| 久久草成人影院| 久久国产乱子伦精品免费另类| 午夜精品在线福利| 九色亚洲精品在线播放| 亚洲性夜色夜夜综合| 国产亚洲欧美98| 亚洲成人久久性| 最近最新免费中文字幕在线| 久久精品亚洲精品国产色婷小说| 国产一区二区三区综合在线观看| 国产精品久久电影中文字幕| 99re在线观看精品视频| 欧美人与性动交α欧美精品济南到| 黄网站色视频无遮挡免费观看| 日韩av在线大香蕉| 色综合婷婷激情| 午夜福利,免费看| 精品无人区乱码1区二区| 母亲3免费完整高清在线观看| 又黄又粗又硬又大视频| 午夜福利在线免费观看网站| 黄片播放在线免费|