Three-dimensional discontinuous deformation analysis (3-D DDA) using a new contact resolution algorithm     

S. Amir Reza Beyabanaki  Roozbeh Geraili Mikola  Kianoosh Hatami   《Computers and Geotechnics》2008,35(3):346-356

This paper presents a new contact calculating algorithm for contacts between two polyhedra with planar boundaries in the three-dimensional discontinuous deformation analysis (3-D DDA). In this algorithm, all six type contacts in 3-D (vertex-to-face, vertex-to-edge, vertex-to-vertex, face-to-face, edge-to-edge, and edge-to-face) are simply transformed into the form of point-to-face contacts. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. In this paper, formulations of normal contact, shear contact and frictional force submatrices based on the new method are derived and the algorithm has been programmed in VC⁺⁺. Examples are provided to demonstrate the new contact rule between two blocks.  相似文献   

胶州湾三维变动边界的潮流数值模拟   总被引：35，自引：3，他引：35  

孙英兰  张越美 《海洋与湖沼》2001,32(4):355-362

基于Blumberg等（1986）的河口、陆架和海洋模式，引入变边界处理技术，建立了胶州湾三维变动边界的潮流模型，模型以σ-正交曲线坐标下三维非线性潮波方程为基本方程，引入2.5阶瑞封闭方程组，采用分裂算子法数值求解方程组，利用湍封闭方程求解垂直紊动粘滞系数，采用干湿网格方法模拟潮流漫滩过程，三维变动边界潮流模型计算结果与实测值吻合良好。  相似文献   

Comparison of two three-dimensional hydrodynamic modeling systems for coastal tidal motion     

Q.Y. Zhang   《Ocean Engineering》2006,33(2):137-151

A comparison of two three-dimensional numerical modeling systems for tidal elevations and velocities in the coastal waters is presented. The two modeling systems are: (1) the Princeton Ocean Model (POM) and (2) the MIKE 3 flow model. The model performance results for Singapore's coastal waters show that the predicted tidal elevations from the two hydrodynamic modeling systems are almost identical and are in very good agreement with field measurement data. The simulated tidal current velocities match well with field measurement data at the selected stations, but it seems that the POM provides the slightly better simulation, compared to the MIKE 3 flow model. The depth profiles of the velocities obtained from the two modeling systems may be greatly different at some time, due to the vertical diffusion coefficient calculated from different turbulent sub-models in the two modeling systems. The POM generally predicts larger peak tidal velocities. The maximum speed differences for the model results from the two modeling systems occur in the top and differ from time to time and from location to location, reaching up to 20%.  相似文献   

渤海、黄海强温跃层数值预报模式     

王宗山  徐伯昌  金梅兵  邹娥梅  李繁华 《海洋科学进展》1996,(3)

文中对前人关于水温垂直结构计算、模拟和预报研究成果进行了评价，并在此基础上建立了考虑海面吸收辐射和透射辐射、地形、风、海流、界面摩擦及其引起的混合作用的强温跃层三维数值预报模式。该模式运行的驱动量仅为风和气温场资料，因而具有较好的实用性；此外，从试报的结果看，效果是令人满意的。  相似文献   

Numerical simulation of migration–accumulation of oil resources     

Yirang Yuan  Yuji Han 《Computational Geosciences》2008,12(2):153-162

Numerical simulation of oil migration and accumulation is to describe the history of oil migration and accumulation in basin evolution. It is of great value to the evaluation of oil resources and to the determination of the location and amount of oil deposits. This thesis discusses the characteristics of petroleum geology and permeation fluid mechanics. For the three-dimensional problems of Dongying hollow of Shengli Petroleum Oil Field, it puts forward a new model and a kind of modified method of upwind finite difference fractional steps implicit interactive scheme. For the famous hydraulic experiment of secondary migration–accumulation, the numerical simulation test has been done, and both the computational and experimental results are basically identical. For the actual problem of Dongying hollow, the numerical simulation test and the actual conditions are basically coincident. Thus, the well-known problem has been solved.  相似文献   

弹性波的三维有限元模拟   总被引：5，自引：0，他引：5       下载免费PDF全文

邓玉琼  张之立 《地球物理学报》1990,33(1):44-53

本文讨论了弹性波瞬态传播问题的三维有限元计算方法及当前存在的实际困难.针对要求计算机内存大和计算时间长的问题,采取了改进措施.由于采用了集中质量矩阵和修正的中心差分时间积分显格式相结合的方法,可以使计算机内存和计算时间大为减少;由于采用结点定位法,最适合用于目前发展的并行计算机系统,可使计算速度大大增快;还采用了有效激发震源法,有效激发区是随时间步进的增加而逐步增大,这不仅能节省计算时间,而且使波场的传播过程一目了然,本文计算了由两种介质组成的三维楔形问题,得到若干典型剖面的瞬时波场图及随时间变化的合成地震图.  相似文献   

Three-dimensional model and late stage warping of the Plattengneis Shear Zone in the Eastern Alps     

Martin Putz  Kurt Stüwe  Mark Jessell  Philippe Calcagno   《Tectonophysics》2006,412(1-2):87-103

The Plattengneis shear zone is a 250–600 m thick, flat lying, Cretaceous, eclogite facies, mylonitic shear zone, with north-over-south transport direction, that is exposed over almost 1000 km² in the Koralpe region along the eastern margin of the Alps. Although the shear zone is one of the largest in the Alps, its role in the Eoalpine metamorphic evolution and the subsequent exhumation of the region, remain enigmatic and its large-scale geometry is not well understood. The outcrop pattern suggests that the shear zone is made up of a single sheet that is folded into a series of open syn- and antiforms with wavelengths of about 10 km. Eclogite bodies occur above, within and below the shear zone and there is no metamorphic grade change across the shear zone. In the south, the fold axes strike east–west and plunge shallowly to the east. In the north, the fold axes are oriented in north–south direction and form a dome shaped structure of the shear zone. Total shortening during this late stage warping event was of the order of 5%. Indirect evidence constrains this folding event to have occurred between 80 and 50 Ma and the fold geometry implies that the final exhumation in the Koralpe occurred somewhat later than further north. Interestingly, the shear zone appears to strike out of the topography in the south and dip into the topography in the north, so that north of the shear zone only hanging-wall rocks are exposed and south of it only foot-wall rocks. Possibilities for the geometric relationship of the Plattengneis shear zone with the surrounding south dipping detachments are discussed.  相似文献   

笛卡儿坐标系的双多普勒天气雷达三维风场反演技术   总被引：27，自引：9，他引：27  

周海光  张沛源 《气象学报》2002,60(5):585-593

文中研究了笛卡尔坐标系下双多普勒天气雷达三维风场反演技术 ,提出了包括雷达原始资料的预处理 ,空间插值以及可靠性检验的新方法 ,提高了反演结果的可信度和精确度。使用模拟的双多普勒雷达体扫资料进行了反演试验 ,结果表明 :本文的方法能够比较真实地反映风场的三维结构 ,可以用于真实风场的反演。  相似文献   

VORTEX FLOW FIELD IN A SCOUR HOLE AROUND ABUTMENTS     

Abdul Karim BARBHUIYA  Subhasish DEY 《国际泥沙研究》2003,18(4)

The three-dimensional flow field in a scour hole around different abutments under a clear water regime was experimentally measured in a laboratory flume, using the Acoustic Doppler Velocimeter (ADV). Three types of abutments used in the experiments were vertical-wall (rectangular section), 45° wing-wall (45° polygonal section) and semicircular. The three-dimensional time-averaged velocity components were detected at different vertical planes for vertical-wall abutment and azimuthal planes for wing-wall and semicircular abutments. The velocity components were also measured at different horizontal planes. In the upstream, presentations of flow field through vector plots at vertical / azimuthal and horizontal planes show the existence of a primary vortex associated with the downflow inside the scour hole. On the other hand, in the downstream, the flow field shows a reversed flow near the abutments having a subsequent recovery with a passage of flow as a part of the main flow. The data presented in this p  相似文献   

Field vibration mitigation by honeycomb WIB for pile foundations of a high-speed train viaduct   总被引：5，自引：0，他引：5  

Hirokazu Takemiya   《Soil Dynamics and Earthquake Engineering》2004,24(1):69-87

Environmental vibrations from recent high-speed trains are becoming a special concern in the civil and environmental engineering field since they can give detrimental effects to residents, sensitive equipments and high-tech production facilities in the vicinity of train tacks. Herein, aiming at the vibration mitigation for a specific high-tech industrial area, the low-frequency vibrations from a train viaduct are targeted over an anticipated range. A theoretically designed innovative countermeasure, called honeycomb wave impeding barrier (WIB) for a wave impeding barrier, is introduced and its effects are investigated by computer simulation. The present WIB is based on the wave dispersion phenomenon that can modulate the incoming wavelengths into the shorter wavelengths, creating an apparent wave cut-off characteristic in the wave field across WIB installation area. The shorter wavelengths are further impeded due to the impedance ratio of the WIB walls and in-fill materials and absorbed by the in-fills more. The three-dimensional FEM simulation demonstrates a dramatic reduction effect that is difficult to achieve by conventional measures.  相似文献