Questioning numerical integration methods for microsphere (and microplane) constitutive equations

In the last few years, more and more complex microsphere models have been proposed to predict the mechanical response of various polymers. Similarly than for microplane models, they consist in deriving a one-dimensional force vs. stretch equation and to integrate it over the unit sphere to obtain a three-dimensional constitutive equation. In this context, the focus of authors is laid on the physics of the one-dimensional relationship, but in most of the case the influence of the integration method on the prediction is not investigated.Here we compare three numerical integration schemes: a classical Gaussian scheme, a method based on a regular geometric meshing of the sphere, and an approach based on spherical harmonics. Depending on the method, the number of integration points may vary from 4 to 983,040! Considering simple quantities, i.e. principal (large) strain invariants, it is shown that the integration method must be carefully chosen. Depending on the quantities retained to described the one-dimensional equation and the required error, the performances of the three methods are discussed. Consequences on stress–strain prediction are illustrated with a directional version of the classical Mooney–Rivlin hyperelastic model. Finally, the paper closes with some advices for the development of new microsphere constitutive equations.     

A Constitutive Equation of Thermoelasticity for Solid Materials

A new constitutive equation of thermoelasticity for crystals is presented based on the interatomic potential and solid mechanics at finite temperature. Using the new constitutive equation, the calculations for crystal copper and graphene are carried out under different loading paths at different temperatures. The calculated results are in good agreement with those of the previous thermoelasticity constitutive equation based on quantum mechanics, which clearly indicates that our new constitutive equation of thermoelasticity is correct. A lot of comparisons also show that the present theory is more concise and efficient than the previous thermal stress theory in the practical application.     

大牛地气田异常修正等时试井资料求取产能方法

大牛地气田为低渗-致密气田,因存在启动压力梯度等因素,修正等时试井曲线按常规方法处理出现异常,无法准确获取气井的产能方程.采用LIT法和异常曲线校正法对气田异常试井资料进行了重新处理,最终获得气井产能方程,为同类气藏产能试井可能出现的异常资料处理解释探索了新的途径.     

Stress-strain relations in elastoplastic solids with Dugdale-type cracks

Mechanical behavior of a two-dimensional elastoplastic solid with rectilinear cracks is investigated. Plastic strip model is used to reduce plasticity problem to the equivalent linear elasticity formulation. Two realizations of the mixed mode plastic strip model are considered: in-line plastic strips as proposed by Becker and Gross [Int. J. Fract. 37 (1988) 163], and inclined plastic strips of Panasyuk and Savruk [Appl. Mech. Rev. 47 (1994) 151]. The effective mechanical response predictions are based on the procedure presented in Kachanov et al., [Appl. Mech. Rev. 47 (1994) 151]. Stress-strain relations are obtained for parallel and randomly oriented non-interacting cracks. Results are compared with known elastic solutions.     

A comparison of transformation methods for evaluating two‐dimensional weakly singular integrals

Accurate numerical evaluation of integrals arising in the boundary element method is fundamental to achieving useful results via this solution technique. In this paper, a number of techniques are considered to evaluate the weakly singular integrals which arise in the solution of Laplace's equation in three dimensions and Poisson's equation in two dimensions. Both are two‐dimensional weakly singular integrals and are evaluated using (in a product fashion) methods which have recently been used for evaluating one‐dimensional weakly singular integrals arising in the boundary element method. The methods used are based on various polynomial transformations of conventional Gaussian quadrature points where the transformation polynomial has zero Jacobian at the singular point. Methods which split the region of integration into sub‐regions are considered as well as non‐splitting methods. In particular, the newly introduced and highly accurate generalized composite subtraction of singularity and non‐linear transformation approach (GSSNT) is applied to various two‐dimensional weakly singular integrals. A study of the different methods reveals complex relationships between transformation orders, position of the singular point, integration kernel and basis function. It is concluded that the GSSNT method gives the best overall results for the two‐dimensional weakly singular integrals studied. Copyright © 2002 John Wiley & Sons, Ltd.     

波动方程数值模拟的隐式差分法

在有限差分波动方程数值模拟中,通常采用高阶差分方法来提高空间导数的数值逼近精度,以实现降低数值频散,提高数值模拟精度的目的。首先对差分频散进行了理论分析;然后讨论了估计一阶空间导数的隐式差分格式,并与通常采用的高阶精度显式差分格式进行了对比分析,结果表明,隐式差分格式能够在更宽的波数范围使差分频散控制在可接受的水平,如8阶精度的显式差分格式所适应的波数带宽约为O．55kmax,而隐式差分格式所适应的波数带宽约为0．7kmax;最后通过模型试算,对隐式差分格式的有效性进行了验证。模拟结果表明,用隐式差分格式在一定程度上降低了差分频散,提高了模拟精度。     

A multigrid accelerated hybrid unstructured mesh method for 3D compressible turbulent flow

 A cell vertex finite volume method for the solution of steady compressible turbulent flow problems on unstructured hybrid meshes of tetrahedra, prisms, pyramids and hexahedra is described. These hybrid meshes are constructed by firstly discretising the computational domain using tetrahedral elements and then by merging certain tetrahedra. A one equation turbulence model is employed and the solution of the steady flow equations is obtained by explicit relaxation. The solution process is accelerated by the addition of a multigrid method, in which the coarse meshes are generated by agglomeration, and by parallelisation. The approach is shown to be effective for the simulation of a number of 3D flows of current practical interest. Sponsored by The Research Council of Norway, project number 125676/410 Dedicated to the memory of Prof. Mike Crisfield, a respected colleague     

A note on the Gauss-Jacobi quadrature formulae for singular integral equations of the second kind

A fast and efficient numerical method based on the Gauss-Jacobi quadrature is described that is suitable for solving Fredholm singular integral equations of the second kind that are frequently encountered in fracture and contact mechanics. Here we concentrate on the case when the unknown function is singular at both ends of the interval. Quadrature formulae involve fixed nodal points and provide exact results for polynomials of degree 2n − 1, where n is the number of nodes. Finally, an application of the method to a plane problem involving complete contact is presented.     

声光效应拉曼-内斯方程矩阵级数解法及声光调制器优化设计

提出一种求解正常声光相互作用拉曼 内斯 (Raman Nath)方程的矩阵级数解法 ,该解法直观方便且具有普遍性。计算结果表明 ,对Q =4 1π ,Bragg衍射的效率只有 97 5 % ;对非对称入射 ,以往的Raman Nath近似解误差较大 ;指出提高Bragg衍射效率的有效途径在于提高声光频率比并给出计算声光器件最优长度的计算公式。     

连续梁的简便计算

利用三转角方程，对连续梁的内力进行精确的计算，该方法简捷明了，便于应用。