考虑地基-结构-散粒体相互作用时贮仓结构的静、动力研究[Ⅱ]--有限元分析

对贮仓结构的静、动力问题进行了系统的分析计算：考虑到地基—结构—散粒体间的相互作用，引入新的计算模式，对不同地基上的贮仓结构模型进行了系统的有限元静、动力分析计算，并与作者所完成的试验结果进行了比较。结果表明，所提出的计算模式及有限元计算模型是正确的。     

近似不可压缩粘弹性结构动力响应的有限元分析

采用近似不可压缩有限元方法,分析了粘弹性结构的动力响应.针对粘弹性材料近似不可压缩性质和准静态分析的局限性,基于Hamilton变分原理,推导出了一种计算近似不可压缩粘弹性结构动力响应的增量有限元方法,考虑不同工况和松弛模量,计算了固体火箭发动机药柱结构的动力响应.理论分析和算例表明该方法简单实用、通用性强,有重要的工程价值.     

散粒体-贮仓结构-地基的动力特性分析

以筒仓模型结构为例，考虑到散粒体-结构-地基的相互作用，对弹性地基上筒仓内的散粒体的不同计算模型进行了多种工况、系统的有限元动力分析计算。通过与筒仓模型动力实验结果进行了比较，得出结论：在对地基施以水平激励时，弹性地基上筒仓的动力响应大于刚性地基上筒仓的动力响应，散粒体与仓壁的相对运动对筒仓结构有减振作用。     

复合结构振动分析的数值方法

本文提出了弹性-粘弹性复合结构动力特性和动力响应的数值计算方法。为了给动力特性分析提供初值,定义了复合结构的“对偶保守结构”概念。用“对偶保守结构”的实模态参数确定复特征值初值,逐个迭代计算复合结构的复模态参数,同时还给出了两种数值方法计算复合结构的时域响应,可适应于一般载荷的响应计算,最后给出了部分计算与实验结果,并进行了比较,结果令人满意。     

精细积分时域平均法和随机扩阶系统法

讨论含随机参数结构的动力响应的计算问题,发展了精细积分时域平均法（ＴＡＰＩＭ）,它可以用来计算确定性系统受到随机激励时的动力响应;结合随机扩阶系统方法与随机有限元法,将ＴＡＰＩＭ方法应用于计算随机参数结构的动力响应,取得了较好的结果。结出了数值算例,结果表明随机扩阶系统法,随机有限元法与精细积分时域平均法的结合是计算 随机参数结构动力响应的有效方法。     

基础振动动力响应的边界单元分析

本文应用边界单元法对基础振动的动力响应进行了数值求解。结构的弹性动力微分方程在通过Laplace积分变换后,可以得到弹性动力的基本边界积分方程。然后在变换空间内划分边界单元进行数值求解。最后通过Laplace的数值逆变换求得时间域内的动力响应值。文中对刚性的动力基础,在简谐荷载的作用下,对于不同频率、不同压缩层厚度和基础埋深等动力响应进行了计算与探讨。     

基于梁平截面假设的复杂细长结构动力模型简化方法

针对纵向尺度显著大于横向的复杂细长结构,提出了一个基于梁平截面假设的模型简化方法,建立了具有物理意义的基向量。同时考虑这类结构截面变形的翘曲模式,增加了翘曲基向量,得到精度较高的简化动力模型。这样的简化模型可以用于原结构在不同边界条件的动力分析。结构频率计算的两个具体算例表明了简化方法的有效性。     

周期结构动力响应的高效数值方法

提出一种计算周期结构动力响应的高效率算法. 以精细积分方法为基础, 利用周期结构的对称性和动力问题的物理特性, 分析了周期结构对应矩阵指数的特殊结构, 并基于此给出一种计算周期结构对应矩阵指数的高效率方法. 在高效和精确计算周期结构对应矩阵指数的基础上, 得到了周期结构动力响应的高效率和高精度算法. 数值算例表明, 该方法效率高且节省存储要求.      

随机结构动力分析的递归聚食缩算法

首先介绍一种新型的随机结构动力分析方法，扩阶系统方法，然后建议了用于扩阶系统动力分析的递归聚缩算法，在不降低计算精度的条件下，递归聚缩算法可以大幅度地提高随机结构动力反应分析的速度，使关于扩阶系统动力分析的计算速度接近于相应的确定性结构系统的计算速度，从而，为随机结构分析的扩阶系统方法进入到实用了阶段铺平了道路。     

随机结构动力分析的递归聚缩算法

首先介绍了一种新型的随机结构动力分析方法：扩阶系统方法，然后建议了用于扩阶系统动力分析的速归聚缩算法．在不降低计算精度的条件下，递归聚缩算法可以大幅度地提高随机结构动力反应分析的速度，使关于扩阶系统动力分析的计算速度接近于相应的确定性结构系统的计算速度．从而，为随机结构分析的扩阶系统方法进入到实用阶段铺平了道路．     

Effective elastic properties of matrix composites with transversely-isotropic phases

The present work addresses the problem of calculation of the macroscopic effective elastic properties of composites containing transversely isotropic phases. As a first step, the contribution of a single inhomogeneity to the effective elastic properties is quantified. Relevant stiffness and compliance contribution tensors are derived for spheroidal inhomogeneities. The limiting cases of spherical, penny-shaped and cylindrical shapes are discussed in detail. The property contribution tensors are used to derive the effective elastic moduli of composite materials formed by transversely isotropic phases in two approximations: non-interaction approximation and effective field method. The results are compared with elastic moduli of quasi-random composites.     

The boundaries of the effective elastic moduli for inhomogeneous solids

The calculation of the effective elastic moduli of inhomogeneous solids, which connect the stresses and strains averaged for the material, is accompanied by certain mathematical difficulties owing to correlation relationships of arbitrary orders. Neglect of correlation relationships leads to average elastic moduli, where averaging according to Voigt and Reuss establishes boundaries containing the effective elastic moduli [1]. Approximate values of the latter can be found by taking into account the correlation relationships of the second order in both calculation schemes [2, 3]. Another method of evaluating the true moduli consists of narrowing the boundaries of Voigt and Reuss on the basis of model representations [4-6]. The approximate effective elastic moduli for a series of polycrystals with various common-angle values are presented in [7]. An analysis of the effect of the correlation relationships between the grains of a mechanical mixture of isotropic components on the effective elastic moduli is carried out in [8], although in all the papers just mentioned the use of correlative corrections to narrow the range of elastic moduli is not investigated. Below it is shown that the calculation of the correlation corrections in the second approximation allows the range for the effective moduli to be narrowed.     

拉压不同弹性模量薄板上圆孔的应力分析

采用弹性理论研究了拉压不同弹性模量薄板上圆孔的孔边应力集中问题．采用广义虎克定律推导出了拉压不同弹性模量薄板上圆孔边的应力平衡方程,并联合利用应力函数及边界条件得到了拉压不同弹性模量薄板上圆孔边的应力表达式．算例分析表明,当薄板材料的拉压弹性模量相差较大时,采用经典弹性理论研究薄板上圆孔的孔边应力是不合适的,当经典弹性理论与拉压不同弹性模量弹性理论的计算结果间的差别超过工程允许误差5%时,应该采用拉压不同弹性模量弹性理论进行计算．     

An Energy Method Applied to Large Elastic Deflection of a Thin Plate of Bimodulus Material

Abstract

Some composite materials and high-polymers are known to behave differently in simple tension and compression under static loads. The present paper is concerned with a method of analysis of the bending of bimodulus elastic plates employing Ambartsumyan-Khachatryan's model for isotropic bimodulus materials. This problem may be reduced to the conventional problem of minimizing the potential energy of the plate as a whole. A simply supported thin square plate subjected to lateral load is analyzed numerically by a simplex method. Results of the calculation show that the effect of the difference between the tensile and compressive elastic moduli on the deformation of the plate may be substantial     

Analytical procedure for determining the linear and nonlinear effective properties of the elastic composite cylinder

In this work we consider a cylindrical structure composed of a nonlinear core (inhomogeneity) surrounded by a different nonlinear shell (matrix). We elaborate a technique for determining its linear elastic moduli (second order elastic constants) and the nonlinear elastic moduli, which are called Landau coefficients (third order elastic constants). Firstly, we develop a nonlinear perturbation method which is able to turn the initial nonlinear elastic problem into a couple of linear problems. Then, we prove that only the solution of the first linear problem is necessary to calculate the linear and nonlinear effective properties of the heterogeneous structure. The following step consists in the exact solution of such a linear problem by means of the complex elastic potentials. As result we obtain the exact closed forms for the linear and nonlinear effective elastic moduli, which are valid for any volume fraction of the core embedded in the external shell.     

Simple algebraic approximations for the effective elastic moduli of cubic arrays of spheres

Recently, Cohen and Bergman (Phys. Rev. B 68 (2003a) 24104) applied the method of elastostatic resonances to the three-dimensional problem of nonoverlapping spherical isotropic inclusions arranged in a cubic array in order to calculate the effective elastic moduli. The leading order in this systematic perturbation expansion, which is related to the Clausius-Mossotti approximation of electrostatics, was obtained in the form of simple algebraic expressions for the elastic moduli. Explicit expressions were derived for the case of a simple cubic array of spheres, and comparison was made with some accurate results. Here, we present explicit expressions for the effective elastic moduli of base-centered and face-centered cubic arrays as well, and make a comparison with other estimates and with accurate numerical results. The simple algebraic expressions provide accurate results at low volume fractions of the inclusions and are good estimates at moderate volume fractions even when the contrast is high.     

Ultrasonic attenuation in polycrystals using a self-consistent approach

The self-consistent method of averaging elastic moduli to define the effective medium of a polycrystal is used to investigate the dynamic problem of wave propagation. An alternative covariance tensor describing the elastic moduli fluctuations of the polycrystal containing self-consistent elastic properties is derived and found to be significantly smaller than the covariance tensor formed through traditional Voigt averaging. Attenuation curves are generated using the self-consistent elastic moduli and covariance tensors and these results are compared with previous Voigt-averaged estimates. The second-order polycrystalline dispersion relation for the self-consistent scheme is found for cases of low and high crystallite anisotropy. The attenuation coefficients and dispersion relations derived through the self-consistent scheme are considerably different than previous estimates. Experimentally measured longitudinal attenuation coefficients support the use of the self-consistent scheme for estimation of attenuation.     

On a Displacement Based Solution to an Antiplane Crack Problem for Inhomogeneous Anisotropic Elastic Materials

This paper employs a displacement based method to examine an antiplane crack problem for an inhomogeneous elastic material in which the elastic moduli vary continuously with the spatial coordinates. Expressions for the crack tip stress intensity factors and the crack displacement are obtained in terms of Chebyshev polynomials. Numerical results are obtained for some particular inhomogeneous elastic materials.     

Elastic properties of inhomogeneous transversely isotropic rocks

The problem to determine the effective elastic moduli and velocities of elastic wave propagation in transversely isotropic solid containing aligned spheroidal inhomogeneities (solid grains, vugs and micro-cracks) has been solved using the self-consistent scheme known as effective medium approximation (EMA). Since a solution of so-called one-particle problem is a base for each self-consistent method, we solved this problem as a first step for spheroidal inhomogeneity in a transversely isotropic medium. In contrast to the known solution of this problem by Lin and Mura we obtained the expressions for the strain field inside inclusion in the explicit form (without quadratures). The obtained solution was used then in the symmetric variant of the EMA where each component of the system was considered as spheroid with its own aspect ratio. This approach was applied to simulate the properties of the rocks containing isolated pores and micro-cracks. For connected fluid-filled pores we used the anisotropic variant of the Gassmann theory. The results of the calculations, obtained for the effective elastic moduli, have been compared with the experimental data and theoretical simulations of the other authors. Unlike many other rock mechanics theories, EMA approximation gives correct elastic moduli values even in the nondilute concentration of inhomogeneities. The comparison of the experimental data for oriented crack system with the EMA predictions indicates their good correspondence.     

Determination of elastic moduli of composite medium containing bimaterial matrix and non-uniform inclusion concentrations

Reservoir porous rocks usually consist of more than two types of matrix materials,forming a randomly heterogeneous material.The determination of the bulk modulus of such a medium is critical to the elastic wave dispersion and attenuation.The elastic moduli for a simple matrix-inclusion model are theoretically analyzed.Most of the efforts assume a uniform inclusion concentration throughout the whole single-material matrix.However,the assumption is too strict in real-world rocks.A model is developed to estimate the moduli of a heterogeneous bimaterial skeleton,i.e.,the host matrix and the patchy matrix.The elastic moduli,density,and permeability of the patchy matrix differ from those of the surrounding host matrix material.Both the matrices contain dispersed particle inclusions with different concentrations.By setting the elastic constant and density of the particles to be zero,a double-porosity medium is obtained.The bulk moduli for the whole system are derived with a multi-level effective modulus method based on Hashin's work.The proposed model improves the elastic modulus calculation of reservoir rocks,and is used to predict the kerogen content based on the wave velocity measured in laboratory.The results show pretty good consistency between the inversed total organic carbon and the measured total organic carbon for two sets of rock samples.