波浪作用下斜坡上护面块体内部应力分布数值模拟
|
摘要
基于3D FEMDEM方法建立三维原型尺度数值模型,模拟波浪荷载作用下斜坡上护面块体内部的应力分布。波浪作用下结构物的水动力荷载采用微幅波理论模拟,护面块体之间的运动、碰撞接触以及块体内部的应力变化采用3D FEMDEM方法模拟。块体之间的接触力采用基于势函数的罚函数法计算,有限元的变形采用中心差分的显式方法求解。应用该数值模型与ANSYS软件程序对自重作用下混凝土扭王字块的内部应力分布特性进行了比较分析,验证了数值模型应力计算的可行性和计算精度。通过数值模拟计算给出了波浪作用下斜坡上护面块体之间的相对运动和块体内部的应力分布及应力历时曲线,探讨了块体内部应力变化特性。
The paper,based on a 3D FEMDEM method,presents a three-dimensional prototype scale numerical model to simulate internal stress distribution of armor blocks on the slopes under wave action. The wave-induced hydrodynamic loads of the structures are treated using the small amplitude wave theory,the armor blocks movement and collision or contact between the armor blocks and the stress changes within the blocks are simulated using the 3D FEMDEM method. The contact force between blocks is calculated by the penalty function method based on potential function. The deformation of finite element is solved by the central difference explicit method. By the application of the numerical model and ANSYS software program,the distribution characteristics of the concrete accropode blocks internal stress under gravity are compared and analyzed,which verifies the feasibility and accuracy of stress calculation of the numerical model. Numerical calculation gives the stress distribution within the blocks and the relative movement between armor blocks on the slopes under wave action and a stress duration curve. The stress properties within the blocks are deeply analyzed.
The paper,based on a 3D FEMDEM method,presents a three-dimensional prototype scale numerical model to simulate internal stress distribution of armor blocks on the slopes under wave action. The wave-induced hydrodynamic loads of the structures are treated using the small amplitude wave theory,the armor blocks movement and collision or contact between the armor blocks and the stress changes within the blocks are simulated using the 3D FEMDEM method. The contact force between blocks is calculated by the penalty function method based on potential function. The deformation of finite element is solved by the central difference explicit method. By the application of the numerical model and ANSYS software program,the distribution characteristics of the concrete accropode blocks internal stress under gravity are compared and analyzed,which verifies the feasibility and accuracy of stress calculation of the numerical model. Numerical calculation gives the stress distribution within the blocks and the relative movement between armor blocks on the slopes under wave action and a stress duration curve. The stress properties within the blocks are deeply analyzed.
引文
[1]Harlow E H.Large rubble-mound breakwater failures[J].Journal of the Waterway Port Coastal and Ocean Division,1980,106(2):275-278.
[2]Oumeraci H.Review and analysis of vertical breakwater failures—lessons learned[J].Coastal Engineering,1994,22(1):3-29.
[3]Van der Meer J W.Deterministic and probabilistic design of breakwater armor layers[J].ASCE Journal of Waterway,Port,Coastal,and Ocean Engineering,1988.114(1):66-80.
[4]Munjiza A,Owen D R J,Bicanic N.A combined finite-discrete element method in transient dynamics of fracturing solids[J].Engineering Computations,1995,12(2):145-174.
[5]Munjiza A.The combined finite-discrete element method[M].England:John Wiley&Sons Ltd.,2004.
[6]Latham J P,Munjiza A,Mindel J,et al.Modelling of massive particulates for breakwater engineering using coupled FEMDEM and CFD[J].Particuology,2008,6(6):572-583.
[7]Gotoh H,Ikari H,Yasuoka T,et al.Particle method for simulating wave overtopping on stepped seawall with drainage[C]//Proceeding of the Coastal Engineering.2008:3071-3083.
[8]李增志.抛石防波堤稳定性的离散单元法分析[D].天津:天津大学,2003.(LI Zeng-zhi.Stablity analysis of rubble mound breakwater by discrete element method[D].Tianjin:Tianjin University,2003.(in Chinese))
[9]任冰,金钊,高睿,等.波浪与斜坡堤护面块体相互作用的SPH-DEM数值模拟[J].大连理工大学学报,2013,53(2):241-248.(REN Bing,JIN Zhao,GAO Rui,et al.Numerical simulation of wave interaction with armor block of mound breakwater by SPH-DEM model[J].Journal of Dalian University of Technology,2013,53(2):241-248.(in Chinese))
[10]REN B,JIN Z,GAO R,et al.SPH-DEM modeling of the hydraulic stability of 2D blocks on a slope[J].Journal of Waterway,Port,Coastal,and Ocean Engineering,2013:DOI:10.1061/(ASCE)WW.1943-5460.0000247
[11]Munjiza A,Andrews K R F.Penalty function method for combined finite-discrete element systems comprising large number of separate bodies[J].International Journal for Numerical Methods in Engineering,2000,49(11):1377-1396.
[12]JTS154-1-2011,防波堤设计与施工规范[S].2011.(JTS 154-1-2011,Design and construction of breakwaters[S].2011.(in Chinese))
[13]Xiang J,Munjiza A,Latham J P.Finite strain,finite rotation quadratic tetrahedral element for the combined finite-discrete element method[J].International Journal for Numerical Methods in Engineering,2009,79(8):946-978.
[14]Munjiza A,Andrews K R F.NBS contact detection algorithm for bodies of similar size[J].International Journal for Numerical Methods in Engineering,1998,43(1):131-149.
[2]Oumeraci H.Review and analysis of vertical breakwater failures—lessons learned[J].Coastal Engineering,1994,22(1):3-29.
[3]Van der Meer J W.Deterministic and probabilistic design of breakwater armor layers[J].ASCE Journal of Waterway,Port,Coastal,and Ocean Engineering,1988.114(1):66-80.
[4]Munjiza A,Owen D R J,Bicanic N.A combined finite-discrete element method in transient dynamics of fracturing solids[J].Engineering Computations,1995,12(2):145-174.
[5]Munjiza A.The combined finite-discrete element method[M].England:John Wiley&Sons Ltd.,2004.
[6]Latham J P,Munjiza A,Mindel J,et al.Modelling of massive particulates for breakwater engineering using coupled FEMDEM and CFD[J].Particuology,2008,6(6):572-583.
[7]Gotoh H,Ikari H,Yasuoka T,et al.Particle method for simulating wave overtopping on stepped seawall with drainage[C]//Proceeding of the Coastal Engineering.2008:3071-3083.
[8]李增志.抛石防波堤稳定性的离散单元法分析[D].天津:天津大学,2003.(LI Zeng-zhi.Stablity analysis of rubble mound breakwater by discrete element method[D].Tianjin:Tianjin University,2003.(in Chinese))
[9]任冰,金钊,高睿,等.波浪与斜坡堤护面块体相互作用的SPH-DEM数值模拟[J].大连理工大学学报,2013,53(2):241-248.(REN Bing,JIN Zhao,GAO Rui,et al.Numerical simulation of wave interaction with armor block of mound breakwater by SPH-DEM model[J].Journal of Dalian University of Technology,2013,53(2):241-248.(in Chinese))
[10]REN B,JIN Z,GAO R,et al.SPH-DEM modeling of the hydraulic stability of 2D blocks on a slope[J].Journal of Waterway,Port,Coastal,and Ocean Engineering,2013:DOI:10.1061/(ASCE)WW.1943-5460.0000247
[11]Munjiza A,Andrews K R F.Penalty function method for combined finite-discrete element systems comprising large number of separate bodies[J].International Journal for Numerical Methods in Engineering,2000,49(11):1377-1396.
[12]JTS154-1-2011,防波堤设计与施工规范[S].2011.(JTS 154-1-2011,Design and construction of breakwaters[S].2011.(in Chinese))
[13]Xiang J,Munjiza A,Latham J P.Finite strain,finite rotation quadratic tetrahedral element for the combined finite-discrete element method[J].International Journal for Numerical Methods in Engineering,2009,79(8):946-978.
[14]Munjiza A,Andrews K R F.NBS contact detection algorithm for bodies of similar size[J].International Journal for Numerical Methods in Engineering,1998,43(1):131-149.