Free vibrations of simply supported nonhomogeneous isotropic rectangular plates of bilinearly varying thickness and elastically restrained edges against rotation using Rayleigh-Ritz method

This paper addresses the free transverse vibrations of thin simply supported nonhomogeneous isotropic rectangular plates of bilinearly varying thickness with elastically restrained edges against rotati...     

Symmetric modes of vibration for A circular plate using finite dynamic elements

The dynamic element method has been shown previously to provide a computational advantage over the ordinary finite element method for various beam elements. The Taylor expansions are computed here for the dynamic shape functions (two terms) and dynamic stiffness matrix (four terms) for the axisymmetric vibrations of an annular plate element. The complicated matrices which result are made more tractable by expressing them as power series in powers of the aspect ratio. The percentage error in the natural frequencies is then calculated using both the two- and the three-term dynamic stiffness matrix, demonstrating the increased accuracy for a given number of elements.     

New approach to analyzing soil-building systems

A new method of analyzing seismic response of soil-building systems is introduced. The method is based on the discrete-time formulation of wave propagation in layered media for vertically propagating plane shear waves. Buildings are modeled as an extension of the layered soil media by assuming that each story in the building is another layer. The seismic response is expressed in terms of wave travel times between the layers, and the wave reflection and transmission coefficients at layer interfaces. The calculation of the response is reduced to a pair of simple finite-difference equations for each layer, which are solved recursively starting from the bedrock. Compared with commonly used vibration formulation, the wave propagation formulation provides several advantages, including the ability to incorporate soil layers, simplicity of the calculations, improved accuracy in modeling the mass and damping, and better tools for system identification and damage detection.     

Vibration of continuous skew plates

This note presents vibration analysis of continuous skew plates with intermediate supports by the Rayleigh-Ritz method with B-spline functions. The accuracy of the present method is shown in a few typical cases. The results are compared with existing results based on other numerical methods and found to be in good agreement. Natural frequencies of continuous skew plates with several elastic point supports are also analysed for different skew angles and for different stiffnesses of point supports.     

G-S变换的快速算法

在电磁场瞬变响应的数值计算中 ,常采用G S变换法作逆拉氏变换 .它是纯实数运算 ,而且只需对较少的拉氏变换变量s值作计算 （通常对每一采样时间选用 1 2个s值 ） ,因而是一种计算速度较快的算法 .但是 ,要对大量采样时间作计算 ,其计算量仍太大 .本文基于拉氏变换的延迟定理 ,建立了一种新的G S变换算法 .数值检验结果表明 ,新算法可成级次地减少对大量采样时间作G S变换的计算量 ,显著提高电磁场瞬变响应的计算速度 .     

采用复合样条有限元方法求解板的动力初值问题

结构动力分析是工程设计中的重要组成部分，传统结构动力分析不能全面反映结构动力的初值特征，而Gurtin变分原理被认为是目前唯一能全面反映结构动力初值特征的变分原理。本文应用以位移为参变量的Gurtin变分原理，采用复合样条有限元的方法，即在时间域及空间域的y方向采用三次B样条函数，而在空间域的x方向采用多项式逼近广义位移，从而建立了精度较高的计算板动力初值问题的样条有限元模型。数值计算结果表明，本文所建立的方法能有效地求解板的动力初值问题，且计算精度高。     

Grid Refinement in Cartesian Coordinates for Groundwater Flow Models Using the Divergence Theorem and Taylor's Series

Grid refinement is introduced in a numerical groundwater model to increase the accuracy of the solution over local areas without compromising the run time of the model. Numerical methods developed for grid refinement suffered certain drawbacks, for example, deficiencies in the implemented interpolation technique; the non‐reciprocity in head calculations or flow calculations; lack of accuracy resulting from high truncation errors, and numerical problems resulting from the construction of elongated meshes. A refinement scheme based on the divergence theorem and Taylor's expansions is presented in this article. This scheme is based on the work of De Marsily (1986) but includes more terms of the Taylor's series to improve the numerical solution. In this scheme, flow reciprocity is maintained and high order of refinement was achievable. The new numerical method is applied to simulate groundwater flows in homogeneous and heterogeneous confined aquifers. It produced results with acceptable degrees of accuracy. This method shows the potential for its application to solving groundwater heads over nested meshes with irregular shapes.     

Axisymmetric dynamic response of the multi-layered transversely isotropic medium

By virtue of the precise integration method (PIM) and the technique of mixed variable formulations, solutions for the dynamic response of the multi-layered transversely isotropic medium subjected to the axisymmetric time-harmonic forces are presented. The planes of cross anisotropy are assumed to be parallel to the horizontal surface of the stratified media. Four kinds of vertically acting axisymmetric loads are prescribed either at the external surface or in the interior of the soil system. Thicknesses and number of the medium strata are not limited. Employing the Hankel integral transform in cylindrical coordinate, the axisymmetric governing equations in terms of displacements of the multi-layered media are uncoupled. Applying mixed variable formulations, more concise first-order ordinary differential matrix equations from the uncoupled motion equations can be obtained. Solutions of the ordinary differential matrix equations in the transformed domain are acquired by utilizing the approach of PIM. Since PIM is highly accurate to solve the sets of first-order ordinary differential equations, any desired accuracy of the solutions can be achieved. All calculations are based on the corresponding algebraic operations and computational efforts can be reduced to a great extent. Comparisons with the existing numerical solutions are made to confirm the accuracy of the present solutions proposed by this procedure. Several examples are illustrated to explore the influences of the type and degree of material anisotropy, the frequency of excitation and loading positions on the dynamic response of the stratified medium.     

An efficient,adaptive algorithm for large-scale random vibration analysis

An efficient, adaptive and robust algorithm is proposed to reduce the cost of large-scale stationary and transient random vibration analysis of structures excited by multiple partially correlated nodal and or base excitations. The cost saving is accomplished by computing integrals selectively, and yet attempting to maintain a level of accuracy desired by the analyst. Recently, proposed closed-form solutions for fully coherent propagating band-limited white-noise excitation are used to rank approximately the terms in the modal covariance matrix. Terms are then evaluated starting from the most important one, and the computations are terminated in such a way that the accuracy level requested by the user is satisfied in an approximate sense. Two variations of the algorithm are proposed: the first one is more robust and is preferred, and is recommended when computing a relatively small number (hundreds) of response quantities; the second one is more efficient when computing a very large number (thousands) of response quantities. Both variations are adaptive, and consider the closeness, damping and participation of all modes explicitly, while the first method also considers the mode shapes. The efficiency and accuracy of the algorithm is investigated by using it to compute the stationary and transient seismic response of the Golden Gate suspension bridge.     

A practical method for estimating dynamic soil stiffness on surface of multi‐layered soil

It is important to estimate the influence of layered soil in soil–structure interaction analyses. Although a great number of investigations have been carried out on this subject, there are very few practical methods that do not require complex calculations. In this paper, a simple and practical method for estimating the horizontal dynamic stiffness of a rigid foundation on the surface of multi‐layered soil is proposed. In this method, waves propagating in the soil are traced using the conception of the cone model, and the impulse response function can be calculated directly and easily in the time domain with a good degree of accuracy. The characteristics of the impedance, that is the transformed value to the frequency domain of the obtained impulse response, are studied using two‐ to four‐layered soil models. The cause of the fluctuation of impedance is expressed clearly from its relation to reflected waves from the lower layer boundary in the model. Copyright © 2005 John Wiley & Sons, Ltd.     

Restoring forces on vertical circular cylinders forced by earthquakes

An eigenfunction solution is presented for the dynamic response of vertical circular cylinders to earthquake excitation in a compressible fluid of finite depth. This single eigenseries expansion eliminates the need for a double summation over both the eigenfunctions and the trial functions as required by Rayleigh-Ritz methods. Revised definitions for the added mass and hydrodynamic radiation damping coefficients per unit length are derived from the hydrodynamic fluid pressures. Based on comparisons between these newly defined coefficients, the compressibility of the fluid is found to be relatively more important at dimensionless frequencies greater than unity (ω > 1.0) when analysing both rigid and flexible cylinders having relatively large diameter to water depth ratios,r₀/h > 0.25 (squatty type). This conclusion regarding the relative importance of the fluid compressibility is derived from a comparison between the relative magnitudes and the vertical distributions over depth of both the added mass and radiation damping coefficients per unit length for both rigid and flexible squatty cylinders. From additional comparisons with Rayleigh-Ritz solutions that require trial functions, the results for totally immersed flexible slender cylinders (r₀/h< 0.10) are shown to be equivalent; but the results for totally immersed flexible squatty cylinders (r₀/h > 0.25) are not. The reason for this difference appears to be in the truncation of the trial function series in the Rayleigh-Ritz methods, which excludes the higher mode shapes, and in the definitions of the added mass coefficients. Comparisons with laboratory data for both rigid and flexible cylinders confirm the accuracy of the solutions obtained by the eigenseries in the limited frequency interval above the highest frequency for surface gravity waves (f > 1.0 Hz) and below the first dimensionless cut-off frequency for acoustic waves (ω< 1.0).     

拉萨台地磁资料的深部电性反演

本文分析了拉萨台长周期成分的感应尺度Ｃ值。周期范围从０．３－６７ｄ，超过两个量级。用Ｍｏｎｔｅ－Ｃａｒｌｏ法对所得的Ｃｒｅ（Ｔ）和Ｃｉｍ（Ｔ）结果进行反演计算，反演结果给出了深部电导率结构可能的多解范围和最佳拟合的电性模型。所谓最佳拟合是指目标函数值最小的一种拟合。     

The complete SRSS modal combination rule

The complete Square‐Root‐of‐Sum‐of‐Squares (c‐SRSS) modal combination rule is presented. It expresses the structural response in terms of uncoupled SDOF modal responses, yet accounting fully for modal response variances and cross‐covariances. Thus, it is an improvement over the classical SRSS rule which neglects contributions from modal cross‐covariances. In the c‐SRSS rule the spectral moments of the structural response are expressed rigorously in terms of the spectral moments of uncoupled modal responses and of some coefficients that can be computed straightforwardly as a function of modal frequencies and damping, without involving the computation of cross‐correlation coefficients between modal responses. An example shows an application of the c‐SRSS rule for structural systems with well separated and closely spaced modal frequencies, subjected to wide‐band and narrow‐band excitations. Comparisons with response calculations using the SRSS and the Complete Quadratic Combination rules are given and discussed in detail. Based on the c‐SRSS rule a response spectrum formulation is introduced to estimate the maximum structural response. An example considering a narrow‐band excitation from the great Mexico earthquake of September 19, 1985, is given and the accuracy of the response spectrum formulation is examined. Copyright © 2010 John Wiley & Sons, Ltd.     

Errors in response calculations for non-classically damped structures

The classical normal mode method of determining response is extremely useful for practical calculations, but depends upon the damping matrix being orthogonal with respect to the modal vectors. Approximations that allow the method to be used when this condition is not satisfied have been suggested; the simplest approach is to neglect off-diagonal terms in the triple matrix product formed from the damping and modal matrices. In this paper the errors in response caused by this approximation are determined for several simple structures for a wide range of damping parameters and different types of excitation. Based on these results a criterion, relating modal damping and natural frequencies, is formulated; if this is satisfied, the errors in response caused by this diagonalization procedure are within acceptable limits.     

Stochastic characterization of earthquakes through their response spectrum

If earthquakes are modelled by a stochastic process, it is possible to interpret the associated response spectrum in terms of the statistics of extreme values of oscillator response to the process. For a stationary earthquake model this interpretation leads to a relationship between the power spectral density function of the process and the response spectrum. This relationship is examined in this paper and forms the basis for two methods presented to obtain the power spectrum of the earthquake process from its response spectrum. One of these methods is approximate but leads to an explicit representation of the power spectral density function in terms of the response spectrum. The other method is exact wherein an iterative scheme for the solution of the problem is established. An example problem is solved to illustrate the use of the two methods and it is shown that for small values of damping, the approximate derivation yields a fairly accurate solution.     

Axisymmetric seismic response of a thick circular plate supporting many rods by modal synthesis

Dynamic response of a thick, horizontal, circular plate supporting a large number of slender rods subjected to uniform boundary motion in the vertical direction has been studied by synthesizing component modes of continuous substructures. The excitation considered corresponds to the vertical component of boundary movement produced by earthquake disturbances and the axisymmetric response problem was solved. Mindlin theory was used to formulate the component equations of the plate which is treated as the main component in a modal synthesis technique. The slender rods, which are attached vertically to the plate, are handled as branch components. Vibration modes of a classical thin plate were used as the initial displacement functions for the Mindlin plate. These functions were subsequently modified by a component mode improvement process to obtain plate modes. System modes were generated by combining the improved plate modes with component modes of rods. Numerical results for the natural frequencies and time-history response of the coupled system are compared with those given by a three-dimensional finite element method.     

基于Bootstrap方法确定岩土参数的标准值

提出一种基于Bootstrap重抽样技术的岩土参数标准值确定方法,将置信概率为95%对应的重抽样均值作为标准值,并通过一组实测土体的抗剪强度参数阐述标准值的确定方法和过程。对两个实例进行研究,分别探讨Bootstrap重抽样次数(即Bootstrap样本集大小)和原始数据的变异性对标准值精度和收敛性的影响。结果表明,基于Bootstrap方法确定的岩土参数标准值具有较好的收敛性和精度,当重抽样次数不小于100时,Bootstrap重抽样次数和原始样本的变异性对标准值的收敛性和精度不会产生显著影响。该方法能够有效地处理岩土参数的不确定性问题。     

A parametric study of the finite element eigenfunction method for the scattering of SH-waves

The performance of a coupled Finite Element—Analytic method in solving elastic wave scattering in infinite spaces is studied for plane SH-waves. The influence of several parameters on the accuracy of the results is investigated. These parameters are the number of terms in the series solution, the number of compatibility points across an imaginary circular interface (i.e. the points for which the continuity of both stresses and displacements are enforced), and the ratio of element dimension to wavelength in the direction of propagation. It was found that the method converges rapidly with the number of terms. Furthermore, both element dimension to wavelength ratio and the number of compatibility points significantly influence the accuracy of the approximate solution.     

A method for stochastic seismic response analysis of non-classically damped structures

A method to calculate the stationary random response of a non-classically damped structure is proposed that features clearly-defined physical meaning and simple expression. The method is developed in the frequency domain, The expression of the proposed method consists of three terms, i.e., modal velocity response, modal displacement response, and coupled (between modal velocity and modal displacement response), Numerical results from the parametric study and three example structures reveal that the modal velocity response term and the coupled term are important to structural response estimates only for a dynamic system with a tuned mass damper. In typical cases, the modal displacement term can provide response estimates with satisfactory accuracy by itself, so that the modal velocity term and coupled term may be ignored without loss of accuracy, This is used to simplify the response computation of non-classically damped structures. For the white noise excitation, three modal correlation coefficients in closed form are derived. To consider the modal velocity response term and the coupled term, a simplified approximation based on white noise excitation is developed for the case when the modal velocity response is important to the structural responses. Numerical results show that the approximate expression based on white noise excitation can provide structural responses with satisfactory accuracy~     

Fast and accurate three‐dimensional controlled source electromagnetic modelling†

We present a fast approximate method for three‐dimensional low frequency controlled source electro‐magnetic modeling. We apply the method to a synthetic model in a typical marine controlled source electromagnetic scenario, where conductivity and permittivity are different from the known background medium. For 3D configurations, fast computational methods are relevant for both forward and inverse modelling studies. Since this problem involves a large number of unknowns, it has to be solved efficiently to obtain results in a timely manner, without compromising accuracy. For this reason, the Born approximation, extended Born approximation and iterative extended Born approximation are implemented and compared with the full solution of the conjugate gradient fast Fourier transformation method. These methods are based on an electric field domain integral equation formulation. It is shown here how well the iterative extended Born approximation method performs in terms of both accuracy and speed with different configurations and different source positions. The improved accuracy comes at virtually no additional computational cost. With the help of this method, it is now possible to perform sensitivity analysis using 3D modelling in a timely manner, which is vital for controlled source electromagnetic applications. For forward modeling the solution at the sea‐bottom is of interest, because that is where the receivers are usually located. For inverse modeling, the accuracy of the solution in the target zone is important to obtain reasonably accurate conductivity values from the inversion using this approximate solution method. Our modelling studies show that the iterative extended Born approximation method is fast and accurate for both forward and inverse modelling. Sensitivity analysis as a function of the source position and different reservoir sizes validate the accuracy of the iterative extended Born approximation.