A two‐step Taylor–Galerkin algorithm applied to shock wave propagation in soils

This paper presents a two‐step Taylor–Galerkin algorithm formulated in terms of velocities and stresses which can be applied to solid dynamics problems requiring good resolution of small wavelengths, such as propagation of shocks. The proposed model is both fast, as it uses simple linear elements (triangles in 2D and tetrahedra in 3D), and accurate. It avoids locking and mesh alignment problems, and therefore can be applied to localized failure computations without the limitations exhibited by the classical displacement formulations. Copyright © 2003 John Wiley & Sons, Ltd.     

倾角函数展开及其在分析法轨道预报中的应用

田谐项摄动是分析法轨道预报中的重要部分,其中包含大量倾角函数及其偏导数的计算.由于具有精度更高、速度更快的优点,倾角函数一般通过递推方法计算.以文献中提出的改进Gooding方法为基础,将其给出的程序稍加改进,在计算2–50阶倾角函数时缩短了约24%的计算时间.考虑到分析法预报过程中轨道平倾角变化很小,以泰勒展开式计算倾角函数,可极大提高计算速度,较大程度地减小分析法预报耗时,且引力场阶次越高,减小幅度越大,取50阶时预报耗时缩短了48%.另一方面,以2阶展开式计算倾角函数时,与改进Gooding法相比,分析法预报星历偏差很小.对于500 km高度的低轨卫星,分别以改进Gooding法和2阶泰勒展开式计算倾角函数,预报3天,当地球引力场阶次不高于50时,二者预报星历偏差RMS (Root Mean Square)低于1 mm,且随着轨道高度的增加,预报星历偏差RMS逐渐减小.     

FINITE ELEMENT MODELLING FOR LAND DISPLACEMENTS DUE TO PUMPING

Equations of equilibrium (force balance) and flow in multidimensions were coupled in this paper to describe land displacements due to pressure decline in aquifers. A Galerkin finite element model based on these equations was developed. The saturated/unsaturated behaviour and the isotropic/anisotropic properties of permeability and elasticity were considered when the model was formulated. This model was verified by comparing its simulation results with those of known analytical solutions for simplified cases. The simulation of displacements due to pressure decline in unsaturated media was also performed. Those results demonstrated that the choice of boundary ranges for an aquifer with infinite domain may significantly affect the estimated horizontal and vertical displacements. To obtain a good estimation of land displacements, the boundary ranges should be carefully chosen. The displacements occurring in unconfined aquifers are caused by the drop of the water table and the change in body force in the dewatering zone. Simulation results also indicated that the change in body force should be considered once an unconfined aquifer has been pumped. Otherwise, the horizontal and vertical displacements in unconfined aquifers would be overestimated and underestimated, respectively. The behaviour of land displacements due to pumping was shown to be affected by changes in the total stresses in aquifers.     

随机波浪下泰勒离散系数的时域解

利用Wolk提出的粒子追踪方程,通过等分频率法划分不规则波谱,利用MATLAB做粒子运动模拟计算,得到无因次化泰勒离散系数K/D随时间t变化的曲线;通过与Huang等得到的P-M谱的泰勒离散系数K/D计算结果比较证明了本计算方法的可靠性。采用该方法研究了不规则波条件下,波序列(同一谱型不同波面序列)和谱型(谱峰周期、有效波高、谱峰升高因子)对波浪离散系数的影响;计算结果表明:同一谱型不同波序列对泰勒纵向离散系数稳定值和稳定时间无影响;不规则波谱峰周期越大,纵向离散系数K/D越小,稳定时间越短;有效波高越大,纵向离散系数K/D越大,稳定时间越长;谱峰升高因子越大,泰勒离散系数K/D越大,稳定时间越长;与规则波相比,不规则波的泰勒离散系数K/D的值略小10%~30%。     

没水倾斜板式防波堤消波性能分析

提出了一种将波浪中倾斜板问题等效化简为波浪中水平板单元组的方法,该方法建立在使用分离变量法求解水工结构边值问题的基础上,并使用伽辽金法精确求解连续边界条件,确定考虑衰减波态的速度势函数,从而求解没水倾斜板结构的消波性能。等效化简法计算精度于边界元法相当,且计算单元数量少、开销低。基于二维线性势波理论,对没水倾斜板式防波堤消波性能分析显示,没水板的倾斜角度、没水深度与板长是结构消波性能的控制因素:没水倾斜板防波堤的消波性能优于没水水平板防波堤的消波性能,随着没水板结构的倾斜角度增大,没水倾斜板结构的波浪透射系数显著减小,且长板优于短板,浅板优于深板;与前人的水槽实验对比显示,相对没水深度与波陡影响结构的消波性能,且波陡造成的波浪破碎贡献了显著的波能消耗。该结论对板式防波堤的结构配置、优化设计有重要意义。     

Estimation of Land Subsidence Based on Groundwater Flow Model

This article presents an approach for estimating land subsidence due to withdrawal of groundwater. The proposed method calculates the groundwater seepage in 3-D-condition and calculates the land subsidence one-dimensionally. The governing equation on groundwater seepage is based on the three-dimensional mass conservation law and the principle of effective stress. The land subsidence calculation method is derived based on the following assumptions: (1) displacements occur only in the vertical direction, and (2) in vertical direction the total stresses do not change. The governing equation is solved by numerical method, i.e., finite element method (FEM) in spatial discretization and finite difference method (FDM) in time series discretization. In FEM Galerkin method is adopted and in FDM, lumped matrix method is employed. The proposed method is calibrated via analyzing 1-D consolidation problem and the results are compared with those from Terzaghi's one-dimensional consolidation theory and oedometer test. The proposed method is employed to analyze the consolidation of a soft layer due to withdrawal of groundwater from an aquifer under it. Moreover, this method is also applied to a field case of land subsidence due to groundwater pumping in a gas production field in Japan. The analytical results are compared with the field observed data. The results show that this approach simulates the field case well.     

Accelerating the discontinuous Galerkin method for seismic wave propagation simulations using multiple GPUs with CUDA and MPI

We have successfully ported an arbitrary high-order discontinuous Galerkin method for solving the three-dimensional isotropic elastic wave equation on unstructured tetrahedral meshes to multiple Graphic Processing Units (GPUs) using the Compute Unified Device Architecture (CUDA) of NVIDIA and Message Passing Interface (MPI) and obtained a speedup factor of about 28.3 for the single-precision version of our codes and a speedup factor of about 14.9 for the double-precision version. The GPU used in the comparisons is NVIDIA Tesla C2070 Fermi, and the CPU used is Intel Xeon W5660. To effectively overlap inter-process communication with computation, we separate the elements on each subdomain into inner and outer elements and complete the computation on outer elements and fill the MPI buffer first. While the MPI messages travel across the network, the GPU performs computation on inner elements, and all other calculations that do not use information of outer elements from neighboring subdomains. A significant portion of the speedup also comes from a customized matrix–matrix multiplication kernel, which is used extensively throughout our program. Preliminary performance analysis on our parallel GPU codes shows favorable strong and weak scalabilities.     

���ڵڶ�����Բ���ֵ������߻��������·���

????????????????????????????????????????????????????·??????÷?????黯γ??????????????????????????????????????????????????????????????????????????????????????t??????