常微分方程技术及其在固体力学计算中的应用

常微分方程(ODE——Ordinary Differential Equation)边值问题的最新计算求解技术的迅速发展推出了一批高质高效的通用软件,而工程中大量的ODE问题并非呈现为这些求解器(Solver)所接受的标准形式。然而,运用一些简单的ODE变换技巧可以将大量的不同类型的特殊问题转化为标准形式。本文列举了若干常用的变换技巧,并广泛地应用于各种固体力学问题的计算中,使大量的ODE问题在形式上得到统一,得以用标准的ODE Solver方便有效地求解。     

Explicit solutions of the problems of elastostatics for an elastic circle with double porosity

In this paper the Aifantis linear theory of consolidation with double porosity is considered. The boundary value problems (BVPs) of elastostatics for an elastic circle are formulated and the uniqueness theorems for regular (classical) solutions are proved. The explicit solutions of these BVPs are constructed by means of absolutely and uniformly convergent series.     

Plane Waves and Uniqueness Theorems in the Coupled Linear Theory of Elasticity for Solids with Double Porosity

This paper concerns with the coupled linear dynamical theory of elasticity for solids with double porosity. Basic properties of plane harmonic waves are established. Radiation conditions of regular vectors are given. Basic internal and external boundary value problems (BVPs) of steady vibrations are formulated, and finally, uniqueness theorems for regular (classical) solutions of these BVPs are proved.     

Dimension Reduction in the Context of Structured Deformations

In the present paper the basic boundary value problems (BVPs) of the full coupled linear theory of elasticity for triple porosity materials are investigated by means of the potential method (boundary integral equation method) and some basic results of the classical theory of elasticity are generalized. In particular, the Green’s identities and the formula of Somigliana type integral representation of regular vector and regular (classical) solutions are presented. The representation of Galerkin type solution is obtained and the completeness of this solution is established. The uniqueness theorems for classical solutions of the internal and external BVPs are proved. The surface (single-layer and double-layer) and volume potentials are constructed and their basic properties are established. Finally, the existence theorems for classical solutions of the BVPs are proved by means of the potential method and the theory of singular integral equations.     

Uniqueness theorems in the theory of thermoelasticity for solids with double porosity

In this paper the coupled linear theory of thermoelasticity for solids with double porosity is considered. The governing system of field equations of this theory is based on motion equations, conservation of fluid mass, constitutive equations, extended Darcy’s law for materials with double porosity and Fourier’s law for heat conduction. A wide class of the basic internal and external boundary value problems (BVPs) of steady vibrations is formulated and uniqueness theorems for regular (classical) solutions of these BVPs are proved.     

Potential Method in the Linear Theory of Viscoelastic Materials with Voids

In the present paper the linear theory of viscoelasticity for Kelvin–Voigt materials with voids is considered and some basic results of the classical theory of elasticity are generalized. Indeed, the basic properties of plane harmonic waves are established. The explicit expression of fundamental solution of the system of equations of steady vibrations is constructed by means of elementary functions. The Green’s formulas in the considered theory are obtained. The uniqueness theorems of the internal and external basic boundary value problems (BVPs) are proved. The representation of Galerkin type solution is obtained and the completeness of this solution is established. The formulas of integral representations of Somigliana type of regular vector and regular (classical) solution are obtained. The Sommerfeld-Kupradze type radiation conditions are established. The basic properties of elastopotentials and singular integral operators are given. Finally, the existence theorems for classical solutions of the internal and external basic BVPs of steady vibrations are proved by using of the potential method (boundary integral method) and the theory of singular integral equations.     

Thermoelastostatic equilibrium of a spatial quadrant: Green’s function and solution in integrals

In this study, a new Green??s function and a new Green-type integral formula for a 3D boundary value problem (BVP) in thermoelastostatics for a quarter-space are derived in closed form. On the boundary half-planes, twice mixed homogeneous mechanical boundary conditions are given. One boundary half-plane is free of loadings and the normal displacements and the tangential stresses are zero on the other one. The thermoelastic displacements are subjected by a heat source applied in the inner points of the quarter-space and by mixed non-homogeneous boundary heat conditions. On one of the boundary half-plane, the temperature is prescribed and the heat flux is given on the other one. When the thermoelastic Green??s function is derived, the thermoelastic displacements are generated by an inner unit point heat source, described by ??-Dirac??s function. All results are obtained in elementary functions that are formulated in a special theorem. As a particular case, when one of the boundary half-plane of the quarter-space is placed at infinity, we obtain the respective results for half-space. Exact solutions in elementary functions for two particular BVPs for a thermoelastic quarter-space and their graphical presentations are included. They demonstrate how to apply the obtained Green-type integral formula as well as the derived influence functions of an inner unit point body force on volume dilatation to solve particular BVPs of thermoelasticity. In addition, advantages of the obtained results and possibilities of the proposed method to derive new Green??s functions and new Green-type integral formulae not for quarter-space only, but also for any canonical Cartesian domain are also discussed.     

Solution for dissimilar elastic inclusions in a finite plate using boundary integral equation method

This paper studies the boundary value problem for a finite plate containing two dissimilar inclusions. The matrix and the two inclusions have different elastic properties. The loadings applied along the outer boundary are in equilibrium. The mentioned problem is decomposed into three boundary value problems (BVPs). Two of them are interior BVP for the elastic inclusions, while the other is a BVP for the triply-connected region. Three problems are connected together through the common displacements and tractions along the interface boundaries. Explicit form for the complex variable boundary integral equation (CVBIE) is derived. After discretization of relevant BIEs, the solutions are evaluated numerically. Three numerical examples for different elastic constant combinations are provided.     

分段常微分方程边值问题的改进COLSYS算法

对分段定义的常微分方程边值问题的 COLSYS算法进行了改进 ,提高了其计算效率及数值稳定性 ,并且可以节省内存 ,这些优势随着分段数的增加而迅速加大。     

Peristaltic motion of third grade fluid in curved channel

＂ Analysis is performed to study the slip effects on the peristaltic flow of non-Newtonian fluid in a curved channel with wall properties. The resulting nonlinear partial differential equations are transformed to a single ordinary differential equation in a stream function by using the assumptions of long wavelength and low Reynolds number. This differential equation is solved numerically by employing the built-in routine for solving nonlinear boundary value problems （BVPs） through the software Mathematica. In addition, the analytic solutions for small Deborah number are computed with a regular perturbation technique. It is noticed that the symmetry of bolus is destroyed in a curved channel. An intensification in the slip effect results in a larger magnitude of axial velocity. Further, the size and circulation of the trapped boluses increase with an increase in the slip parameter. Different from the case of planar channel, the axial velocity profiles are tilted towards the lower part of the channel. A comparative study between analytic and numerical solutions shows excellent agreement.     

具有待定形函数的有限元线法

有限元线法是一种优良的半解析法，但其解函数存在解析方向和离散方向上的精度不相称的问题。本文提出了一种基于有限元线法的双向解析法——具有待定形函数的有限元线法。该方法设形函数为未知待定并借助能量原理对形函数进行解析性求解，使解函数在两个方向上的精度达到相称，从而大幅度提升了解答的整体质量和精度。文中，以二维Poisson方程问题为例，具体给出了待定形函数的有限元线法的算法和算例，用以表明本法的可行性和有效性。     

Determination of the deformed position of a thin shell from surface strains and height function

This paper focuses on the following problem: given the strain tensor of a deformed reference surface of a thin shell and the distances of the points on this surface from some arbitrarily fixed reference plane (the so-called height function) find the position of this reference surface. Two alternative procedures supplying the solution are developed. The first one follows from the ideas developed by Darboux (Leçons sur la Théorie Général des Surfaces, Troisiéme Partie, Gauthier-Villars, Paris, 1894), whereas the second one is based on the polar decomposition theorem and techniques developed in continuum mechanics. These procedures are purely kinematic, valid for arbitrary surface geometry and for unrestricted surface strains. Szwabowicz (Deformable surfaces and almost inextensional deflections of thin shells, Habilitation Thesis, 1999) proposed a relatively simple non-linear boundary-value problem (BVP) for thin elastic shells, which was expressed in three surface strains and the height function as basic independent field variables. The results of this paper suggest that this approach to the non-linear problems of thin shells may be an attractive alternative to other BVPs developed in the literature.     

Plane Elastic Boundary Value Problem Posed by Displacement and Force Orientations on a Closed Contour

A boundary value problem (BVP) of the plane elasticity posed in terms of the orientations of forces and displacements is considered. The main aim of the present paper is to investigate the solvability of BVPs of this kind. Firstly, analysis of two cases is performed: the case of a circle with special orientations of force and displacement vectors on the circumference and the case of an arbitrary contour with coaxial orientations of these vectors. The solutions obtained indicate that the problem can have a certain number of solutions or be unsolvable. Then the BVP is reduced to a boundary integral equation and its solvability is investigated for the general case of a smooth simple-connected closed contour. As a result, the number of linearly independent solutions is determined. This number only depends upon the angle between the force and displacement vectors. This revised version was published online in July 2006 with corrections to the Cover Date.     

Second-gradient homogenized model for wave propagation in heterogeneous periodic media

The paper is focused on a homogenization procedure for the analysis of wave propagation in materials with periodic microstructure. By a reformulation of the variational-asymptotic homogenization technique recently proposed by Bacigalupo and Gambarotta (2012a), a second-gradient continuum model is derived, which provides a sufficiently accurate approximation of the lowest (acoustic) branch of the dispersion curves obtained by the Floquet–Bloch theory and may be a useful tool for the wave propagation analysis in bounded domains. The multi-scale kinematics is described through micro-fluctuation functions of the displacement field, which are derived by the solution of a recurrent sequence of cell BVPs and obtained as the superposition of a static and dynamic contribution. The latters are proportional to the even powers of the phase velocity and consequently the micro-fluctuation functions also depend on the direction of propagation. Therefore, both the higher order elastic moduli and the inertial terms result to depend by the dynamic correctors. This approach is applied to the study of wave propagation in layered bi-materials with orthotropic phases, having an axis of orthotropy parallel to the direction of layering, in which case, the overall elastic and inertial constants can be determined analytically. The reliability of the proposed procedure is analysed by comparing the obtained dispersion functions with those derived by the Floquet–Bloch theory.     

Free vibrations of a beam with elastic end restraints subject to a constant axial load

A general model for vibration of beams restrained with two transversal and two rotational elastic springs subject to a constant axially load is presented. The frequency equations and the shape functions are derived analytically. The proposed model can be employed for simulating the dynamic responses of elastically supported beams in tension or compression for most classical boundary conditions. Some simplifications in the degenerate cases are deduced to evaluate the effectiveness of the model. Numeric examples are given for engineering applications. This model unifies most of the previous vibration models and provides a convenient tool for the analyses of various beam vibrations in tension and compression conditions.     

On the treatment of soft soil parameter uncertainties in planetary rover mobility simulations

Nowadays soft soil–wheel contact models are widely used for predicting the mobility of rovers in off-road applications. However, most of the contact models used in computer simulations are based on semi-empirical laws for which soil parameters can be assessed only with large uncertainty. This lack of knowledge results in significant uncertainty on the rover position predictions. Applied to a planetary rover model, this paper illustrates probabilistic and non-probabilistic techniques for efficient treatment of soil parameter uncertainty for rover position predictions.     

Recent developments in axial tomography for heat transfer and fluid flow studies

The need for three-dimensional, nonintrusive field measurements in the area of heat transfer and fluid flows grows very rapidly, and at this time there are very few experimental techniques that can be used for such a purpose. Axial tomography is a promising technique for accurate quantitative measurements for a variety of heat transfer and fluid flow problems. It has already been tested and practiced in a number of applications, including medical scanning technology. The technique is based on Radon's original work and reconstructs a three-dimensional field from its two-dimensional projections (integrated measurements) obtained at different view angles. This review highlights some of the most important methods used in tomography and their applications to experimental heat transfer and fluid flow studies.     

DEVELOPMENT AND APPLICATIONS OF THE ELASTO-PLASTIC CELLULAR AUTOMATON

The paper presents the advancement and applications of the elasto-plastic cellular automaton (EPCA),a simulator for rock mechanics and rock engineering.The most significant feature of EPCA lies in its ...     

单纤维实验技术的现状和展望

本文综述了目前流行的测量纤维／树脂界面力学性能的单纤维实验技术。重点介绍了三种单纤维实验技术——单纤维拔出法、单纤维段埋入法和单纤维压出法的基本原理。讨论了三种单纤维实验技术在理论和应用上的发展重点     

A thermomechanical constitutive model for cemented granular materials with quantifiable internal variables. Part II – Validation and localization analysis

We study the mechanical failure of cemented granular materials (e.g., sandstones) using a constitutive model based on breakage mechanics for grain crushing and damage mechanics for cement fracture. The theoretical aspects of this model are presented in Part I: Tengattini et al. (2014), A thermomechanical constitutive model for cemented granular materials with quantifiable internal variables, Part I – Theory (Journal of the Mechanics and Physics of Solids, 10.1016/j.jmps.2014.05.021). In this Part II we investigate the constitutive and structural responses of cemented granular materials through analyses of Boundary Value Problems (BVPs).The multiple failure mechanisms captured by the proposed model enable the behavior of cemented granular rocks to be well reproduced for a wide range of confining pressures. Furthermore, through comparison of the model predictions and experimental data, the micromechanical basis of the model provides improved understanding of failure mechanisms of cemented granular materials. In particular, we show that grain crushing is the predominant inelastic deformation mechanism under high pressures while cement failure is the relevant mechanism at low pressures. Over an intermediate pressure regime a mixed mode of failure mechanisms is observed. Furthermore, the micromechanical roots of the model allow the effects on localized deformation modes of various initial microstructures to be studied. The results obtained from both the constitutive responses and BVP solutions indicate that the proposed approach and model provide a promising basis for future theoretical studies on cemented granular materials.