共查询到20条相似文献，搜索用时 31 毫秒
1.
B. Nath 《地震工程与结构动力学》1981,9(2):117-131
A novel Finite Element Method (FEM) is proposed for the analysis of the uncoupled hydrodynamic pressures generated on arch dams due to a steady-state ground acceleration. In this method the equation governing hydrodynamic pressures and also the prescribed boundary conditions are all transformed from the Cartesian space to a logarithmically condensed cylindrical polar space; in this process the physical configuration of the reservoir-dam is also mapped into an ‘image’ domain. The transformed governing equation is then solved in the image domain, subject to the transformed boundary conditions, using standard finite elements. Because physical dimensions are logarithmically condensed in the image space, the proposed method is particularly suitable for dealing with large or very large aspect-ratio reservoir-dam systems, economically and efficiently. The high degree of accuracy which the proposed method is capable of, and also the simple way in which it can be applied to complex reservoir-dam shapes, have been demonstrated by means of examples. The method has also been applied to study the uncoupled hydrodynamic pressures on the upstream face of a cylindrical arch dam, generated by a steady-state vertical ground acceleration. 相似文献
2.
B. Nath 《地震工程与结构动力学》1981,9(1):33-51
A novel method has been proposed for the finite element solution of the steady-state scalar wave equation in three-dimensions. In this the governing equation and the prescribed boundary conditions in the physical space are transformed into a spherical polar space in which the radial direction is logarithmically condensed; the physical problem domain is also mapped into the new space. The transformed equation is then solved in the mapped domain using conventional finite elements. Because physical dimensions of the problem are logarithmically condensed in the proposed spherical polar space, the method is particularly suitable for solving truly three-dimensional problems in which the aspect ratio(s) is large or very large. A number of illustrative examples considered show that the proposed method is capable of a high degree of accuracy, achieved efficiently and economically. A hybrid scheme has also been proposed for dealing with awkward-shaped domains. 相似文献
3.
Traditional Ensemble Kalman Filter (EnKF) data assimilation requires computationally intensive Monte Carlo (MC) sampling, which suffers from filter inbreeding unless the number of simulations is large. Recently we proposed an alternative EnKF groundwater-data assimilation method that obviates the need for sampling and is free of inbreeding issues. In our new approach, theoretical ensemble moments are approximated directly by solving a system of corresponding stochastic groundwater flow equations. Like MC-based EnKF, our moment equations (ME) approach allows Bayesian updating of system states and parameters in real-time as new data become available. Here we compare the performances and accuracies of the two approaches on two-dimensional transient groundwater flow toward a well pumping water in a synthetic, randomly heterogeneous confined aquifer subject to prescribed head and flux boundary conditions. 相似文献
4.
A detailed comparison between fully dynamic and kinematic plate formulations has been made in models of mantle convection. Plate velocity is computed self-consistently from fully dynamic plate models with temperature- and stress-dependent viscosity and preexisting mobile faults. In fully dynamic models, the flow is driven solely by internal buoyancy, while in kinematic models the flow is driven by a combination of the prescribed surface velocity and internal buoyancy. Only a temperature-dependent viscosity, close to the effective viscosity determined from the fully dynamic models, is used in the kinematic models. The two types of models give very similar temperature structures and slab evolutionary histories when the effective viscosity and surface velocity are nearly identical. In kinematic plate models, the additional work introduced by the prescribed velocity boundary condition is apparently dissipated within the lithosphere and has little influence on the convection under the lithosphere. In models with periodic lateral boundary conditions, slabs sink into the lower mantle at an oblique angle and this contrasts with the vertical sinking which occurs with reflecting boundary conditions. Models show that we can simulate fully dynamic models with kinematic models under either periodic boundary conditions or reflecting boundary conditions. 相似文献
5.
A novel finite element method has been described in this paper for the solution of the steady-state scalar wave equation in two dimensions. In this method the physical domain of the problem is mapped into an image domain; the governing equation and the prescribed boundary conditions are also appropriately transformed. Because logarithmic mapping functions are used, the physical domain is logarithmically condensed in the image plane. The method is therefore particularly suitable for the efficient and economical solution of large or very large aspect ratio problems. The high degree of accuracy which the method is capable of is demonstrated by means of two typical examples. Possible extension of the method to complex problems including non-linearity, multiplicity, etc. is also discussed. 相似文献
6.
This paper presents a time-dependent semi-analytical artificial boundary for numerically simulating elastic wave propagation problems in a two-dimensional homogeneous half space. A polygonal boundary is considered in the half space to truncate the semi-infinite domain, with an appropriate boundary condition imposed. Using the concept of the scaled boundary finite element method, the wave equation of the truncated semi-infinite domain is represented by the partial differential equation of non-constant coefficients. The resulting partial differential equation has only one spatial coordinate variable and time variable. Through introducing a few auxiliary functions at the truncated boundary, the resulting partial differential equations are further transformed into linear time-dependent equations. This allows an artificial boundary to be derived from the time-dependent equations. The proposed artificial boundary is local in time, global at the truncated boundary and semi-analytical in the finite element sense. Compared with the scaled boundary finite element method, the main advantage in using the proposed artificial boundary is that the requirement for solving a matrix form of Lyapunov equation to obtain the unit-impulse response matrix is avoided, so that computer efforts are significantly reduced. The related numerical results from some typical examples have demonstrated that the proposed artificial boundary is of high accuracy in dealing with time-dependent elastic wave propagation in two-dimensional homogeneous semi-infinite domains. 相似文献
7.
The distinct element method (DEM) has been used successfully for the dynamic analysis of rigid block systems. One of many
difficulties associated with DEM is modeling of damping. In this paper, new procedures are proposed for the damping modeling
and its numerical implementation in distinct element analysis of rigid multi-block systems. The stiffness proportional damping
is constructed for the prescribed damping ratio, based on the non-zero fundamental frequency effective during the time interval
while the boundary conditions remain essentially constant. At this time interval, the fundamental frequency can be estimated
without complete eigenvalue analysis. The damping coefficients will vary while the damping ratio remains the same throughout
the entire analysis. A new numerical procedure is developed to prevent unnecessary energy loss that can occur during the separation
phases. These procedures were implemented in the development of the distinct element method for the dynamic analyses of piled
multi-block systems. The analysis results for the single-block and two-block systems were in a good agreement with the analytic
predictions. Applications to the seismic analyses of piled fourblock systems revealed that the new procedures can make a significant
difference and may lead to much-improved results. 相似文献
8.
A modified method for simulating non-stationary multi-point earthquake ground motion 总被引：1，自引：1，他引：0
A spectral-representation-based algorithm is proposed to simulate non-stationary and stochastic processes with evolutionary power,according to a prescribed non-stationary cross-spectral density matrix. Non-stationary multi-point seismic ground motions at different locations on the ground surface are generated for use in engineering applications. First,a modified iterative procedure is used to generate uniformly modulated non-stationary ground motion time histories which are compatible with the prescribed power spectrum. Then,ground motion time histories are modeled as a non-stationary stochastic process with amplitude and frequency modulation. The characteristic frequency and damping ratio of the Clough-Penzien acceleration spectrum are considered as a function of time in order to study the frequency time variation. Finally,two numerical examples are presented to validate the efficiency of the proposed method,and the results show that this method can be effectively applied to the dynamic seismic analysis of long and large scale structures. 相似文献
9.
10.
David P Thambiratnam 《地震工程与结构动力学》1986,14(3):475-485
Herein the propagation of transient waves in a rod due to impulsive end loading is investigated using the method of wavefront expansion. The Mindlin–Herrmann equations which include the Poisson effect are used to model the rod. Since this rod theory is restricted to lower values of Poisson's ratio (< 0.283), the present work also has the same restriction. The impulsive end loading on the rod can be prescribed in the form of stress, strain, velocity or acceleration boundary conditions. The analysis is based on the concept of a wave as a carrier of discontinuities in the field variables and their derivatives. These discontinuities are determined from a set of recurrence relations which are in turn generated by the use of time harmonic asymptotic series solutions to the equations of motion. Transients due to both velocity and strain (or stress) boundary conditions are treated. Numerical examples are presented to illustrate the method of solution. The results indicate the distortion of the transient waves due to the dispersive nature of the governing equations. Transients due to other boundary conditions at the rod end can be treated by using the results presented in the paper. 相似文献
11.
Online numerical simulation: A hybrid simulation method for incomplete boundary conditions 
下载免费PDF全文
下载免费PDF全文 《地震工程与结构动力学》2018,47(4):889-905
Hybrid simulation is a powerful and cost‐effective simulation technique to evaluate structural dynamic performance. However, it is sometimes rather difficult to guarantee all the boundaries on the physical substructures, especially when the boundary conditions are very complex, due to limited laboratory resources. Lacking of boundary conditions is bound to change the stress state of the structure and eventually result in an inaccurate evaluation of structural performance. A model updating‐based online numerical simulation method is proposed in this paper to tackle the problem of incomplete boundary conditions. In the proposed method, 2 sets of finite element models with the same constitutive model are set up for the overall analysis of the whole structure and the constitutive model parameter estimation of the physical substructure, respectively. The boundary conditions are naturally satisfied because the response is calculated from the overall structural model, and the accuracy is improved as the material constitutive parameters are updated. The effectiveness of the proposed method is validated via numerical simulations and actual hybrid tests on a RC frame structure, and the results show that the negative effect of incomplete boundary conditions is almost eliminated and the accuracy of hybrid simulation is very much improved. 相似文献
12.
在多源直流电阻率法有限元三维数值模拟中,传统混合边界条件由于刚度矩阵与源点位置相关,无法实现线性方程组的快速回代求解,目前常用齐次边界条件或无限元边界进行替代,虽然实现了快速回代求解,但同时也降低了数值模拟的精度.为了实现快速回代求解,并确保数值模拟的计算精度,本文提出了一种近似边界条件方法.其核心思想是将与场源位置相关的边界系数矩阵从刚度矩阵中分离出来,使得分离后的刚度矩阵与场源位置无关.并将边界系数矩阵与边界处一次电场向量的乘积移到线性方程组右端源项中,当场源位置发生改变时,只需要重新计算右端源项就可以实现快速回代求解.理论模型数值计算表明,在水平地形条件下,本文边界条件数值精度优于混合边界条件;在起伏地形条件下,与齐次边界条件相比,本文边界条件数值结果与混合边界条件吻合度更高. 相似文献
13.
N. Bhanja 《Pure and Applied Geophysics》1969,75(1):42-46
Stresses and displacements in a transversely isotropic cylindrical shell when it vibrates radially under prescribed boundary loadings have been obtained with the help of finite Hankel transforms [1]. 相似文献
14.
The energy transmitting boundary used in programs such as FLUSH and ALUSH is a very accurate and useful technique for the earthquake response analysis of soil–structure interaction systems. However, it is applicable only to linear analyses or equivalent linear analyses, because it can be calculated only in the frequency domain. The author has proposed methods for transforming frequency-dependent impedance into the time domain. In this paper, an earthquake response analysis method for a soil–structure interaction system, using the energy transmitting boundary in the time domain, is proposed. First, the transform of the transmitting boundary matrices to the time domain using the methods proposed by the author is studied. Then, linear and nonlinear time history earthquake response analyses using the boundary are performed. Through these studies, the validity and efficiency of the proposed methods are confirmed. 相似文献
15.
本文讨论了在波动方程偏移问题中构造吸收边界条件方程的原则,并提出了求解具有吸收边界条件的偏移问题的一种新方法--矩阵拟合分解方法。这一方法的特点在于所构造的差分方程的拟合分解格式可以通过三点褶积和正、反递归被快速求解。 通过合成记录和野外实际资料的偏移对这一方法进行了验证。理论分析和实际计算都表明,该方法可以有效地消除来自网格边界的人工反射,减少了因增加边道或其他方法所需要的辅助计算而带来的计算工作量,提高了计算效率,因而是一种有效方法。 相似文献
16.
A numerical scheme is developed in the paper for calculating torsional, vertical, horizontal, coupling and rocking impedances in frequency domain for axial-symmetric foundations embedded in layered media. In the scheme, the whole soil domain is divided into interior and exterior domains. For the exterior domain, the analytic solutions with unknown coefficients are obtained by solving three-dimensional (3D) wave equations in cylindrical coordinates satisfying homogeneous boundary conditions. For the interior domain, the analytical solutions are also obtained by solving the same 3D wave equations satisfying the homogeneous boundary conditions and the prescribed boundary conditions. The prescribed conditions are the interaction tractions at the interfaces between embedded foundation and surrounding soil. The interaction tractions are assumed to be piecewise linear. The piecewise linear tractions at the bottom surface of foundation will be decomposed into a series of Bessel functions which can be easily fitted into the general solutions of wave equations in cylindrical coordinates. After all the analytic solutions with unknown coefficients for both interior and exterior domains are found, the variational principle is employed using the continuity conditions (both displacements and stresses) at the interfaces between interior and exterior domains, interior domain and foundation, and exterior domain and foundation to find impedance functions. 相似文献
17.
在地震波场数值模拟过程中,边界反射是影响其模拟结果的一个重要因素。实际地下介质具有各向异性特征,传统的完全匹配层边界（PML）对于小入射角地震波具有良好效果,但该方法并不能有效地吸收低频波和大角度入射波。针对VTI介质边界反射的问题,本文提出在频率-空间域有限差分法数值模拟中采用卷积完全匹配层（CPML）和特征分析法的组合边界条件,并对该组合边界条件进行数值模拟实验和边界反射吸收效果分析,验证所提方法是一种可靠的人工吸收边界条件,能够有效地压制波场模拟过程中产生的边界反射。 相似文献
18.
在地震波场的波动方程数值模拟中,由于计算量的限制,必须加入人为的边界,使模拟计算可以在一定的空间范围内进行. 由于边界节点上的波场值不能像模拟区域内部的节点一样使用中心差分来计算,使其计算精度大大降低,从而产生边界反射. 为了消除边界反射,本文提出了延迟边界方法,根据弹性波在传播方向上等距离质点的等相位延迟性质和振幅衰减特性,由内部波场的时空分布,推算出边界波场的相位延迟的大小和振幅衰减系数,从而提高边界节点上的波场值计算精度,消除边界反射的产生. 相似文献
19.
Review on localized boundary integral equation:Discrete wavenumber method for 2D irregular layers 下载免费PDF全文
下载免费PDF全文 The pioneer study of simulating the wave field in media with irregular interface belongs to Aki and Larner.Since that many numerical methods on the subject have been developed,such as pure numerical techniques,ray method and boundary method.The boundary method based on boundary integral equation is a semi-analytical method which is suitable to modeling wave field induced by irregular border.According to the property of the applied Green's function the boundary methods can be sorted into space domain boundary method and wavenumber domain boundary method.For both of them it is necessary to solve a large equation,which means much computation is needed.Thus,it is difficult for the boundary methods to be applied in simulating wave field with high frequency or in large range.To develop a new method with less computation is meaningful.For this purpose,localized boundary integral equation,i.e.,discrete wavenumber method is proposed.It is rooted in the Bouchon-Campillo method,an important wavenumber domain boundary method.Firstly the force on interface is separated into two parts:one is on flat part and the other on irregular part of the interface.Then Fourier transform is applied to identify their relation,the unknown distributes only on irregular part.Consequently computation efficiency is dramatically improved.Importantly its accuracy is the same as that of Bouchon-Campillo. 相似文献
20.
Modified Mixed Lagrangian–Eulerian Method Based on Numerical Framework of MT3DMS on Cauchy Boundary 下载免费PDF全文
下载免费PDF全文 Heejun Suk 《Ground water》2016,54(4):508-520
MT3DMS, a modular three‐dimensional multispecies transport model, has long been a popular model in the groundwater field for simulating solute transport in the saturated zone. However, the method of characteristics (MOC), modified MOC (MMOC), and hybrid MOC (HMOC) included in MT3DMS did not treat Cauchy boundary conditions in a straightforward or rigorous manner, from a mathematical point of view. The MOC, MMOC, and HMOC regard the Cauchy boundary as a source condition. For the source, MOC, MMOC, and HMOC calculate the Lagrangian concentration by setting it equal to the cell concentration at an old time level. However, the above calculation is an approximate method because it does not involve backward tracking in MMOC and HMOC or allow performing forward tracking at the source cell in MOC. To circumvent this problem, a new scheme is proposed that avoids direct calculation of the Lagrangian concentration on the Cauchy boundary. The proposed method combines the numerical formulations of two different schemes, the finite element method (FEM) and the Eulerian–Lagrangian method (ELM), into one global matrix equation. This study demonstrates the limitation of all MT3DMS schemes, including MOC, MMOC, HMOC, and a third‐order total‐variation‐diminishing (TVD) scheme under Cauchy boundary conditions. By contrast, the proposed method always shows good agreement with the exact solution, regardless of the flow conditions. Finally, the successful application of the proposed method sheds light on the possible flexibility and capability of the MT3DMS to deal with the mass transport problems of all flow regimes. 相似文献