用样条有限条法计算高层剪力墙结构

针对高层剪力墙结构的特殊性,运用样条有限条法对实体墙和开孔墙进行了计算,并通过计算实例予以验证。

三次梁样条函数及其应用

文中提出了一种求二阶导数连续的三次样条函数的新方法，构造了样条有限条方法。

Modeling Technology in Finite Element Analysis of Electrostatic Proximity Fuze Problem

In order to analyze the electrostatic field concerned with electrostatic proximity fuze problem using the available finite analysis software package, the technology to model the problem with a scale reduction object and boundary was presented. The boundary is determined by the maximum distance the sensor can detect. The object model is obtained by multiplying the terms in Poisson's equation with a scale reduction factor and the real value can be reconstructed with the same reverse process after software calculation. Using the finite element analysis program, the simulation value is close to the theoretical value with a little error. The boundary determination and scale reduction method is suitable to modeling the irregular electrostatic field around air targets, such as airplane, missile and so on, which is based on commonly used personal computer (PC). The technology reduces the calculation and storage cost greatly.

Improved absorbing boundary condition based on linear interpolation for ADI-FDTD method

With the linear interpolation method, an improved absorbing boundary condition（ABC）is introduced and derived, which is suitable for the alternating-direction-implicit finite- difference time-domain （ADI-FDTD） method. The reflection of the ABC caused by both the truncated error and the phase velocity error is analyzed. Based on the phase velocity estimation and the nonuniform cell, two methods are studied and then adopted to improve the performance of the ABC. A calculation case of a rectangular waveguide which is a typical dispersive transmission line is carried out using the ADI-FDTD method with the improved ABC for evaluation. According to the calculated case, the comparison is given between the reflection coefficients of the ABC with and without the velocity estimation and also the comparison between the reflection coefficients of the ABC with and without the nonuniform processing. The reflection variation of the ABC under different time steps is also analyzed and the acceptable worsening will not obscure the improvement on the absorption. Numerical results obviously show that efficient improvement on the absorbing performance of the ABC is achieved based on these methods for the ADI-FDTD.

Dynamic analysis of dam-reservoir-foundation interaction using finite difference technique

Time domain dynamic analysis of inclined dam-reservoir-foundation interaction was conducted using finite difference method （FDM）. The Timoshenko beam theory and the Euler-Bemoulli beam theory were implemented to draw out governing equation of beam. The interactions between the dam and the soil were modeled by using a translational spring and a rotational spring. A Sommerfeld＇s radiation condition at the infinity boundary of the fluid domain was adopted. The effects of the reservoir bottom absorption and surface waves on the dam-reservoir-foundation interaction due to the earthquake were studied. To avoid the instability of solution, a semi-implicit scheme was used for the discretization of the governing equation of dam and an explicit scheme was used for the discretization of the governing equation of fluid. The results show that as the slope of upstream dam increases, the hydrodynamic pressure on the dam is reduced. Moreover, when the Timoshenko beam theory is used, the system response increases.

Decomposition method for solving parabolic equations in finite domains

This paper presents a comparison among Adomian decomposition method (ADM), Wavelet-Galerkin method (WGM), the fully explicit (1,7) finite difference technique (FTCS), the fully implicit (7,1) finite difference method (BTCS), (7,7) Crank-Nicholson type finite difference formula (C-N), the fully explicit method (1,13) and 9-point finite difference method, for solving parabolic differential equations with arbitrary boundary conditions and based on weak form functionals in finite domains. The problem is solved rapidly, easily and elegantly by ADM. The numerical results on a 2D transient heat conducting problem and 3D diffusion problem are used to validate the proposed ADM as an effective numerical method for solving finite domain parabolic equations. The numerical results showed that our present method is less time consuming and is easier to use than other methods. In addition, we prove the convergence of this method when it is applied to the nonlinear parabolic equation.

3D prestack reverse time migration of ground penetrating radar data based on the normalized correlation imaging condition

The reverse time migration(RTM)of ground penetrating radar(GPR)is usually implemented in its two-dimensional(2D)form,due to huge computational cost.However,2D RTM algorithm is difficult to focus the scattering signal and produce a high precision subsurface image when the object is buried in a complicated subsurface environment.To better handle the multi-off set GPR data,we propose a three-dimensional(3D)prestack RTM algorithm.The high-order fi nite diff erence time domian(FDTD)method,with the accuracy of eighth-order in space and second-order in time,is applied to simulate the forward and backward extrapolation electromagnetic fi elds.In addition,we use the normalized correlation imaging condition to obtain pre-stack RTM result and the Laplace fi lter to suppress the low frequency noise generated during the correlation process.The numerical test of 3D simulated GPR data demonstrated that 3D RTM image shows excellent coincidence with the true model.Compared with 2D RTM image,the 3D RTM image can more clearly and accurately refl ect the 3D spatial distribution of the target,and the resolution of the imaging results is far better.Furthermore,the application of observed GPR data further validates the eff ectiveness of the proposed 3D GPR RTM algorithm,and its fi nal image can more reliably guide the subsequent interpretation.

Numerical Solution of Some Diffusion Problems in 3-Layered 3D Domain

In this paper we consider averaging and finite difference methods for solving the 3-D boundary-value problem in a multilayered domain. We consider the metal concentration in the 3 layered peat blocks. Using experimental data the mathematical model for calculating the concentration of metal at different points in peat layers is developed. A specific feature of these problems is that it is necessary to solve the 3-D boundary-value problems for the partial differential equations （PDEs） of the elliptic type of second order with piece-wise diffusion coefficients in the three layer domain. We develop here a finite-difference method for solving a problem of the above type with the periodical boundary condition in x direction. This procedure allows reducing the 3-D problem to a system of 2-D problems by using a circulant matrix.

On the Existence for Some Special Primitive Elements in Finite Fields

Let F_q be a finite field of characteristic p. In this paper, by using the index sum method the authors obtain a sufficient condition for the existence of a primitive elementα∈ F_(q^n) such that α + α^(-1)is also primitive or α + α^(-1)is primitive and α is a normal element of F_(q^n) over F_q.

Analysis of thin plates by the weak form quadrature element method

The recently proposed weak form quadrature element method （QEM） is applied to flexural and vibrational analysis of thin plates The integrals involved in the variational description of a thin plate are evaluated by an efficient numerical scheme and the par- tial derivatives at the integration sampling points are then approximated using differential quadrature analogs. Neither the grid pattern nor the number of nodes is fixed, being adjustable according to convergence need. The C~ continuity conditions char- acterizing the thin plate theory are discussed and the robustness of the weak form quadrature element for thin plates against shape distortion is examined. Examples are presented and comparisons with analytical solutions and the results of the finite element method are made to demonstrate the convergence and computational efficiency of the weak form quadrature element method. It is shown that the present formulation is applicable to thin plates with varying thickness as well as uniform plates.

Effects of first-order chemical reactions on the dispersion coefficient associated with laminar flow through the lungs

This paper aims to look into the determination of effective area-average concentration and dispersion coefficient associated with unsteady flow through a small-diameter tube where a solute undergoes first-order chemical reaction both within the fluid and at the boundary. The reaction consists of a reversible component due to phase exchange between the flowing fluid and the wall layer, and an irreversible component due to absorption into the wall. To understand the dispersion, the governing equations along with the reactive boundary conditions are solved numerically using the Finite Difference Method. The resultant equation shows how the dispersion coefficient is influenced by the first-order chemical reaction. The effects of various dimensionless parameters e.g. Da （the Damkohler number）, a （phase partitioning number） and F （dimensionless absorption number） on dispersion are discussed. One of the results exposes that the dispersion coefficient may approach its steady-state limit in a short time at a high value of Damkohler number （say Da 〉 10） and a small but nonzero value of absorption rate （say P 〈0.5）.

Orientation-preservation conditions on an iso-parametric FEM in cavitation computation

The orientation-preservation condition, i.e., the Jacobian determinant of the deformation gradient det Vu being required to be positive, is a natural physical constraint in elasticity as well as in many other fields. It is well known that the constraint can often cause serious difficulties in both theoretical analysis and numerical computation, especially when the material is subject to large deformations. We derive a set of necessary and sufficient conditions for the quadratic iso-parametric finite element interpolation functions of cavity solutions to be orientation preserving on a class of radially symmetric large expansion accommodating triangulations. The result provides a practical quantitative guide for meshing in the neighborhood of a cavity and shows that the orientation-preservation can be achieved with a reasonable number of total degrees of freedom by the quadratic iso-parametric finite element method.