Eddy diffusivity derived from drifter data for dispersion model applications

Ocean transport and dispersion processes are at the present time simulated using Lagrangian stochastic models coupled with Eulerian circulation models that are supplying analyses and forecasts of the ocean currents at unprecedented time and space resolution. Using the Lagrangian approach, each particle displacement is described by an average motion and a fluctuating part. The first one represents the advection associated with the Eulerian current field of the circulation models while the second one describes the sub-grid scale diffusion. The focus of this study is to quantify the sub-grid scale diffusion of the Lagrangian models written in terms of a horizontal eddy diffusivity. Using a large database of drifters released in different regions of the Mediterranean Sea, the Lagrangian sub-grid scale diffusion has been computed, by considering different regimes when averaging statistical quantities. In addition, the real drifters have been simulated using a trajectory model forced by OGCM currents, focusing on how the Lagrangian properties are reproduced by the simulated trajectories.     

Non-linear finite element dynamic analysis of axisymmetric solids

The subject of this paper is the finite element analysis of axisymmetric solids subjected to axisymmetric static and dynamic loading, and taking into account material as well as geometric non-linearities. A general Lagrangian formulation forms the basis for the incremental equations of motion which are solved using direct integration methods. Solution accuracy is improved by applying equilibrium correction loads at each step. Finite element discretization is achieved through the use of quadrilateral plane stress and axisymmetric elements with incompatible modes added for improvement of the element flexural characteristics. Several numerical examples are presented to demonstrate the effectiveness of the developed computer program.     

A macro‐element model for inelastic building analysis

A three‐dimensional model for approximate inelastic analysis of buildings is presented herein. The model is based on a single macro‐element per building storey. The inelastic properties of the model are characterized by the so‐called ultimate storey shear and torque (USST) surfaces. Different algorithms for the construction of these surfaces, as well as their applications in building modelling, are presented and discussed. Two alternative procedures are developed to integrate the force‐deformation constitutive relationship of the macro‐elements. The first one follows the exact trajectory of the load path of the structure on the USST, and the second uses linear programming without ever forming the USST surface. The accuracy of the model and integration procedure is evaluated by means of the earthquake response of single‐storey systems. The model and integration procedure developed is finally used to compute the inelastic response of a seven‐storey R/C building. The results of this investigation show that the model proposed, although approximate, can be effective in estimating the inelastic deformation demand of a building. It also enables the engineer to capture and interpret important features of the three‐dimensional inelastic response of a structure even before performing any inelastic dynamic analysis. Copyright © 2000 John Wiley & Sons, Ltd.     

复杂加载路径下堆石料动力本构模型及数值模拟

采用所发展的亚塑性边界面模型，对堆石料在不同固结应力状态、排水条件和加载方式下的典型试验进行了数值模拟分析。在本构模型的数值实施中，联合应用自适应多步回退欧拉积分规律和局部迭代加快了迭代收敛的速度，提高了积分精度，进而编制了计算程序TESTROCK。数值模拟结果与已有试验成果的对比分析表明，该模型能够合理、有效地反映堆石料在复杂加载路径下的变形和强度特性，为深入揭示堆石料变形和孔隙水压力的发展特性和进行高混凝土面板堆石坝的非线性耦合静动力分析提供了理论基础。     

An ELLAM approximation for advective–dispersive transport with nonlinear sorption

We consider an Eulerian–Lagrangian localized adjoint method (ELLAM) applied to nonlinear model equations governing solute transport and sorption in porous media. Solute transport in the aqueous phase is modeled by standard advection and hydrodynamic dispersion processes, while sorption is modeled with a nonlinear local-equilibrium model. We present our implementation of finite volume ELLAM (FV-ELLAM) and finite element (FE-ELLAM) discretizations to the reactive transport model and evaluate their performance for several test problems containing self-sharpening fronts.     

Finite element analysis of the seismic response of anchored and unanchored liquid storage tanks

The seismic response of liquid-filled cylindrical storage tanks has been investigated using finite element techniques implemented in the general purpose structural analysis computer code ANSYS. Both added mass concepts and displacement-based fluid finite elements were employed to allow for the effects of the liquid. Simplified response spectrum modal analyses of a tank making use of the axisymmetric harmonic displacement patterns of the principal modes of deformation were found to give accurate predictions of the tank behaviour with a rigidly anchored base. Time history analyses of three-dimensional finite element models of unanchored and flexibly anchored tanks, with gap conditions between the tank base and the supporting floor to allow lift-off of the base, indicated that stresses in the tank and resultant loads on the floor can be much greater than for a rigidly restrained tank. These results demonstrate the importance of carefully considering the restraint conditions when performing seismic design calculations on storage tanks.     

Arbitrary Lagrangian Eulerian formulation for soil structure interaction problems

Simulation of soil structure interaction problems becomes more and more the focus of computational engineering in the last decades in the civil engineering community. In many applications in civil engineering, the Lagrangian formulation is the main approach for soil structure interaction problems. For the last decades, various approaches have been investigated in fluid and structure mechanics, to solve fluid structure interaction problems, including Arbitrary Eulerian Lagrangian (ALE) formulation and meshless methods as Smoothing Particle Hydrodynamics (SPH). These two approaches are more common in fluid mechanics than structure and civil engineering, where the classical Lagrangian formulation is still dominant for numerical simulations. The Lagrangian formulation, where the mesh moves with material is mainly used to solve problems in solid mechanics, and civil engineering for soil structure interaction problems. For small deformation, Lagrangian formulation can solve soil structure interface and material boundary accurately; the main limitation of the formulation is high mesh distortion for large deformation and moving structure. One of the commonly used approaches to solve these problems is the ALE formulation which has been used with success in the simulation of fluid structure interaction with large structure motion such as sloshing fuel tank in automotive industry and bird impact in aeronautic industry. In this paper, devoted to ALE formulation for civil engineering problems, the mathematical and numerical implementation of the ALE formulation and the contact algorithm are described. In order to validate the ALE algorithms, and prove its ability for civil engineering applications, we consider two different applications. First we perform the simulation of a shock wave propagation generated by a detonation of an explosive. In the second application, we consider a soil structure coupling problem, using an ALE formulation for the soil which undergoes large mesh deformation, and a Lagrangian formulation for the structure. A contact algorithm is performed at the soil structure interface.     

透射边界高频耦合失稳机理及稳定实现——SH波动

就大型近场波动的高效数值模拟而言,稳定实现高阶人工边界是一个尚未圆满解决的问题.本文针对使用多次透射公式的SH波动集中质量有限元模拟,依据GKS定理的群速度解释,进一步阐明了人工边界与内域离散格式耦合所导致高频失稳的机理,即两者支持群速度指向内域的外行高频平面谐波,波动能量自发地从人工边界进入內域,从而导致失稳,而这类谐波是由集中质量有限元离散引入的.本文提出了消除此种耦合失稳的一种方法:通过修改有限元刚度阵来改变内域离散格式,并保证修改格式的精度不低于原有格式的精度.理论分析和数值实验表明此法能稳定实现透射边界.本文研究结果具有推广应用前景.     

近地空中爆炸作用下钢框架结构冲击响应

采用多物质欧拉与拉格朗日耦合算法,对框架结构和混凝土地面采用Lagrange单元、空气和炸药采用多物质ALE单元,建立流固耦合有限元模型,对一个五层的钢框架结构在近地爆炸作用下冲击压力波的传播过程、结构冲击响应和变形以及破坏过程进行了数值模拟.数值模拟分析结果表明,拉格朗日-欧拉耦合算法较好地模拟了爆炸冲击波在介质中的传播和作用在结构上的爆炸荷载效应.在爆炸近区,空气爆炸冲击波的衰减速率快;爆炸产生强烈的冲击波首先使钢框架结构正面的柱和梁柱节点区产生很大的塑性变形,钢框架柱在发生爆炸的瞬间失去承载能力,最终导致建筑物整体倒塌.考虑结构变形和流体荷载间的相互影响可以较真实地模拟结构在爆炸荷载作用下的连续倒塌过程,为研究结构在爆炸荷载作用下的连续倒塌提供了有效手段.     

An Hermitian finite element solution of the two-dimensional saturated-unsaturated flow equation

This paper describes a Galerkin-type finite element solution of the two-dimensional saturated-unsaturated flow equation. The numerical solution uses an incomplete (reduced) set of Hermitian cubic basis functions and is formulated in terms of normal and tangential coordinates. The formulation leads to continuous pressure gradients across interelement boundaries for a number of well-defined element configurations, such as for rectangular and circular elements. Other elements generally lead to discontinuous gradients; however, the gradients remain uniquely defined at the nodes. The method avoids calculation of second-order derivatives, yet retains many of the advantages associated with Hermitian elements. A nine-point Lobatto-type integration scheme is used to evaluate all local element integrals. This alternative scheme produces about the same accuracy as the usual 9- or 16-point Gaussian quadrature schemes, but is computationally more efficient.     

Newtonian nudging for a Richards equation-based distributed hydrological model

The objective of data assimilation is to provide physically consistent estimates of spatially distributed environmental variables. In this study a relatively simple data assimilation method has been implemented in a relatively complex hydrological model. The data assimilation technique is Newtonian relaxation or nudging, in which model variables are driven towards observations by a forcing term added to the model equations. The forcing term is proportional to the difference between simulation and observation (relaxation component) and contains four-dimensional weighting functions that can incorporate prior knowledge about the spatial and temporal variability and characteristic scales of the state variable(s) being assimilated. The numerical model couples a three-dimensional finite element Richards equation solver for variably saturated porous media and a finite difference diffusion wave approximation based on digital elevation data for surface water dynamics. We describe the implementation of the data assimilation algorithm for the coupled model and report on the numerical and hydrological performance of the resulting assimilation scheme. Nudging is shown to be successful in improving the hydrological simulation results, and it introduces little computational cost, in terms of CPU and other numerical aspects of the model’s behavior, in some cases even improving numerical performance compared to model runs without nudging. We also examine the sensitivity of the model to nudging term parameters including the spatio-temporal influence coefficients in the weighting functions. Overall the nudging algorithm is quite flexible, for instance in dealing with concurrent observation datasets, gridded or scattered data, and different state variables, and the implementation presented here can be readily extended to any of these features not already incorporated. Moreover the nudging code and tests can serve as a basis for implementation of more sophisticated data assimilation techniques in a Richards equation-based hydrological model.     

面向目标自适应海洋可控源电磁三维矢量有限元正演

本文实现了一种面向目标自适应海洋可控源电磁三维矢量有限元方法.为满足三维复杂电性结构模拟的需求,网格剖分采用非结构化六面体.在组装刚度矩阵之后,形成的大型复数线性方程组分解为等价的实数形式,利用带预条件的广义最小残差法进行求解.在获得微分方程的解之后,为提高解的准确性,通过面向目标的自适应误差估计来指示网格细化,重点加密能使观测点数值模拟精度提高的网格.对于大规模三维数据,为了使模型空间的并行计算达到均衡负载的效果,我们使用METIS函数库来进行网格计算任务量的划分.最后,通过对比一维解析解与三维自适应矢量有限元计算结果,验证了程序的正确性;通过自适应过程中误差指示子的分布,验证了面向目标自适应的有效性;通过对三维复杂模型进行均衡负载下的并行计算,测试了程序的可扩展性.     

A new efficient Eulerian–Lagrangian localized adjoint method for solving the advection–dispersion equation on unstructured meshes

A new and high efficient scheme is developed for the Eulerian–Lagrangian Localized Adjoint Method (ELLAM) to solve the Advection–Dispersion transport Equation (ADE) on unstructured triangular meshes. To obtain accurate results, the new method requires a very limited number of integration points (usually 1 per element).     

Nonlinear analysis of masonry structures using fiber-section line elements

For seismic analysis of unreinforced masonry (URM) buildings characterized by a box-like behavior, a widely accepted model is based on the equivalent frame idealization of walls. The equivalent frame model uses 1D elements to represent the vertical piers and horizontal spandrels which are connected by rigid nodes. The mechanical characterization of the elements is one of the crucial aspects to predict reasonably the building seismic behavior. Through the comparison with pseudo-static and dynamic experimental tests performed on two-story full-scale buildings, this paper validates the frame modeling in the OpenSees framework, which includes a fiber-section force-based beam element for the axial-flexural behavior, coupled with a cyclic shear-deformation phenomenological law.     

Mass lumping models of the linear diffusion equation

The nodal domain integration method is used to develop a numerical model of the linear diffusion equation. The nodal domain integration approach is shown to represent an infinity of finite element mass matrix lumping schemes including the Galerkin and subdomain integration versions of the weighted residual method and an integrated finite difference method. Neumann, Dirichlet and mixed boundary conditions are accommodated analogous to the Galerkin finite element method. In order to reduce the overall integrated approximation relative error, a mass matrix lumping formulation is developed which is based on the Crank-Nicolson time advancement approximation. The optimum mass lumping factors are found to be strongly related to the model timestep size.     

一致激励与多点激励对悬索桥地震响应影响分析

利用大型通用有限元软件ANSYS的APDL参数化语言,建立某悬索桥结构三维有限元模型,通过加速度时程积分曲线,获得相应的位移时程曲线,在桥台及桥墩处施加位移时程荷载,分析桥梁结构多点激励地震响应,并分析一致激励与多点激励下桥梁关键位置位移、内力差异,研究不同激励方式的影响。     

Analysis of GPS Measurements in Eastern Canada Using Principal Component Analysis

Continuous Global Positioning System (CGPS) position time series from eastern North America constrain the pattern and magnitude of regional crustal deformation. Initial analysis delineates consistent uplift patterns, as expected from glacial isostatic adjustment (GIA) predictions, but the associated horizontal deformation is not definitive, in part due to the short time periods for a significant number of the available stations. We employ an eigenpattern decomposition in order to define a unique, finite set of deformation patterns for this continuous GPS data. Similar in nature to the empirical orthogonal functions historically employed in the analysis of atmospheric and oceanographic phenomena, the method derives the eigenvalues and eigenstates from the diagonalization of the correlation matrix using a Karhunen–Loeve expansion (KLE). The KLE technique is used to identify the important modes in both time and space for the CGPS data, modes that potentially include signals such as horizontal and vertical GIA, tectonic strain, and seasonal effects. Here we filter both the vertical and horizontal velocity patterns on different spatiotemporal scales in order to study the potential geophysical sources, after the removal of correlated and random noise. The method is successful in disaggregating the linear vertical signal from more variable and less spatially correlated signals. The vertical and horizontal results are compared to the predictions of the ICE-3G GIA loading model with a number of plausible mantle viscosity profiles. The horizontal velocity analysis allows for qualitative differentiation between several potential GIA models and suggests that, with longer time series, this technique can be employed to remove correlated noise and improve estimates of crustal strain patterns and their sources.     

Nodal domain integration model of one-dimensional advection—diffusion

The nodal domain integration method is applied to a one-dimensional advection—diffusion mathematical model without a source term. Comparison of the resulting numerical model to the well known Galerkin finite element, subdomain, and finite difference domain models indicates that a single numerical statement can be developed which includes the Galerkin finite element, subdomain, and finite difference models as special cases.     

粘弹性人工边界单元及地震动输入方法比较研究

本文引用一致粘弹性边界的概念,推导了与一致粘弹性边界等效的粘弹性边界单元,得到了粘弹性边界单元的刚度矩阵和阻尼矩阵。通过对比粘弹性边界单元阻尼系数不同取值的计算精度,给出了阻尼系数推荐取值。在输入地震动时,采用一种基于局部人工边界的子结构输入方法,实现了和波动法相同的输入效果,并且避免了复杂的应力计算和荷载作用方向判断等问题,应用更加方便。采用自由场的算例,比较了人工边界和不同输入方式的组合在SV波垂直入射时的模拟精度。为了说明边界单元的变形对计算结果的影响,提出了使用底边界输入位移幅值u_max和边界单元宽度b的比值ε=u_max/b作为变形指标,调节边界单元边长尺寸。通过算例给出了ε的推荐取值。最后,选取合适的ε后建立一个圆形孔洞散射场地模型,对散射场地模型的模拟结果验证了粘弹性边界单元和子结构输入方式应用于散射场地时仍具有较高精度。     

徐淮地区煤矿井筒破裂原因初探

本用震源力学理论和方法研究了徐淮地区从１９７０年以来构造应力场的方向和强度的时空变化过程。结果表明：以唐山地震为分界线，本区的应力场Ｐ轴取向由震前平均６１．８°变为震后平均７７．７°。如果将本区以宿北断裂为界分为南区和北区两个部分，则北区的Ｐ轴取向从６８．１°变７１．２°，而南区的Ｐ轴取向由６２．５°变到８３．６°，南区的变化明显于北区。