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法和异常曲线校正法对气田异常试井资料进行了重新处理,最终获得气井产能方程,为同类气藏产能试井可能出现的异常资料处理解释探索了新的途径.     

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衍射效率的有效途径在于提高声光频率比并给出计算声光器件最优长度的计算公式。     

连续梁的简便计算

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

用核磁共振测井确定渗透率的原理和方法分析

对常规测井和核磁共振测井的渗透率公式进行了比较 ,发现各渗透率公式存在相互联系。几种核磁共振测井渗透率公式的计算结果和实验室渗透率的对比表明核磁共振参数和渗透率有很好的相关性 ;改进的渗透率公式可以准确地求取岩心和地层的渗透率。     

Network topology and the theory of rubber elasticity

The earliest investigations on rubber elasticity, commencing in the 19th century, were necessarily limited to phenomenological interpretations. The realisation that polymers consist of very long molecular chains. commencing c. 1930, gave impetus to the molecular theory of rubber elasticity (1932-). according to which the high deformability of an elastomer, and the elastic force generated by deformation, stem from the configurations accessible to long molecular chains. Theories of rubber elasticity put forward from 1934-1946 relied on the assumption that the junctions of the rubber network undergo displacements that are affine in macroscopic strain. The theory of James and Guth (1947) dispensed with this premise, and demonstrated instead that the mean positions of the junctions of a ‘phantom’ network consisting of Gaussian chains devoid of material properties are affine in the strain. The vital significance of the distinction between the actual distribution of chain vectors in a network and their distribution if the junctions would be fixed at their mean positions went unnoticed for nearly 30 years. Experimental investigations, commencing with the incisive work of Gee in 1946. revealed large departures from the relationship of stress to strain predicted by the theories cited. This discrepancy prompted extensive studies, theoretical and experimental, during succeeding years. Inquiry into the fundamentals of polymer networks, formed for example by interlinking very long polymer molecules, exposed the need to take account of network imperfections, typically consisting of chains attached at only one end to a network junction. Various means were advocated to make corrections for these imperfections. The cycle rank ζ of the network has been shown (1976) to be the fundamental measure of its connectivity, regardless of the junction functionality and pattern of imperfections. Often overlooked is the copious interpenetration of the chains comprising typical elastomeric networks. Theories that attempt to represent such networks on a lattice are incompatible with this universal feature. Moreover, the dense interpenetration of chains may limit the ability of junctions in real networks to accommodate the fluctuations envisaged in the theory of phantom networks. It was suggested in 1975 that departures from the form predicted for the elastic equation of state are due to constraints on the fluctuations of junctions whose effect diminishes with deformation and with dilation. Formulation of a self-consistent theory based on this suggestion required recognition of the non-affine connection between the chain vector distribution function and the macroscopic strain in a real network, which may partake of characteristics of a phantom network in some degree. Implementation of the idea was achieved through postulation of domains of constraint affecting the equilibrium distribution of fluctuations of network junctions from their mean positions. This led in due course to a theory that accounts for the relationship of stress to strain virtually throughout the ranges of strain accessible to measurement. The theory establishes connections between structure and elastic properties. This is achieved with utmost frugality in arbitrary parameters.