动力应力函数张量

本文将弹性静力学中的应力函数张量概念直接推广到连续介质动力学问题中,给出了动力应力函数张量的一般表达式.其形式在一般曲线坐标下,可写成而在Descartes坐标系下为.     

弹、粘性体动力学变分原理的Laplace变换形式、有限元构式及数值方法

作者在[1]文中提出了弹、粘动力学变分原理的谱分解形式,本文将其推广到Laplace变换形式,具体写出了薄壳动力学的混合变分原理以及弹-粘-孔隙介质力学的变分原理,并对后者作出了有限元构式. Laplace变换形式的变分原理具有简洁形式,为便于有限元法计算,当已知Laplace变换式的有限个值时,需求原时间函数的有限个值,对此当前尚无成熟方法,本文提供了求原函数的数值方法.从例题可见,这种数值方法是有效的. 结合以上两种理论:从变分原理进行有限元构式以及求Laplace反变换的数值方法,可以使相当广的一类固体动力学问题能够用电子计算机进行求解.     

折线强化弹塑性应力分析的有限元法

本文对材料的应力应变曲线用三段直线的折线拟合,按照弹塑性的简单加载理论,对以增量理论得出的完整应力应变关系进行简化,导出按位移求解的有限元的增量方程.其中弹塑性刚度矩阵可以从弹性刚度矩阵补充后得出,从而节省计算时间.根据von Mises屈服准则确定各次荷载的增量,引入迭代法进行求解,省去对弹塑性刚度矩阵的重复地三角分解,进一步减少计算时间.本文对于应用高次单元、偏离简单加载的荷载、卸载计算、曲线拟合以及荷载的估算问题,均作了说明.     

用张量导数法求横观各向同性弹性材料的弹性张量的一般形式

用对张量函数求导的方法导出了横观各向同性材料和各向同性材料的弹性张量的一般形式与应力-应变关系式.从推导过程可更清楚地看出为什么横观各向同性材料和各向同性材料分别有五个和两个独立的弹性常数,即材料有几个独立的弹性常数是由其应变能函数的形式所决定的.     

论连续体有限变形的位移协调条件

有限变形的协凋条件在文献中常以Riemann-Christoffel张量等于零表达.水文应用Cesaro方法和作者的非线性应变-转动张量分解定理证明上述条件仅是必要的,尚不充分保证位移场的单值性与连续;文中导出新的一般有限变形的位移协调条件.当应变与转动微小时,它化为Saint-Venant方程.     

热弹性材料的自由能表达式及其与物性系数的关系*

本文中,自由能的表达式展开成幂级数,其中,温度增量θ*取到三阶,应变张量γ_ij只取到二阶.由这个表达式可以导出物性系数随温度增量的变化规律.这些规律与参考文献中的实验图线是相符的.不过,自由能表达式中的常数须由实验数据确定.文中指出,变化的弹性模量E和剪切弹性模量G是彼此独立的,而其它的物性系数则与它们相关.     

钱氏定理在有限变形极矩弹性力学广义变分原理的应用

应用Lagrange乘子法和钱伟长证明的两类广义变分原理的等价定理,在本文中导出有限变形极矩弹性力学的广义变分原理.文中采用了在拖带坐标系描述法建立的有限变形应变张量(称为Biot有限变形应变定义的准确形式)和应变速率定义与拖带系应力张量构成完整的数学描述.     

多元Logistic回归模型的结构变点估计

多元Logistic回归模型是广义线性模型中的一种常见形式,在社会科学和生物医学等领域有着广泛的应用.本文首先基于Pearson卡方统计量对Logistic回归模型的结构变点进行估计,再结合二元分割方法将其推广到多变点的情形.数值模拟结果表明基于Pearson卡方统计量的二元分割方法能有效估计出变点,且当变点之间间隔的样本较多时,估计效果较好.最后将此方法应用于一组DNA数据上,说明方法的有效性.     

考虑损伤的修正梯度弹性理论

梯度弹性理论在描述材料微结构起主导作用的力学行为时具有显著优势,将其与损伤理论相结合,可在材料破坏研究中考虑微结构的影响.基于修正梯度弹性理论,将应变张量、应变梯度张量和损伤变量作为Helmholtz自由能函数的状态变量,并在自然状态附近对自由能函数作Taylor展开,进而由热力学基本定律,推导出修正梯度弹性损伤理论本构方程的一般形式.编制有限元程序,模拟土样损伤局部化带的发展演化过程.结果表明,修正梯度弹性损伤理论消除了网格依赖性;损伤局部化带不是与损伤同时发生,而是在损伤发展到一定程度后再逐渐显现出来.     

一类分两段收取保费的相关风险模型

本文研究了一类索赔时间间隔与索赔量相关且分两段收取保费的风险模型.利用微分的方法得到了罚金折现期望满足的积分-微分方程,并给出了此方程解的一个表达式.     

Basic equations of the theory of reinforced media

The author derives the basic equations of the theory of composite elastic media obtained by reinforcing some elastic medium with a large number of linear or planar elastic elements with high strength and deformation resistance. The argument is based on macrostructural considerations. The stress-strain state of each of the reinforcing elements is considered with allowance for interaction with the matrix material. In addition, the "smoothing" principle introduced in [1–3] is applied. This corresponds to approximating the reinforced medium with some equivalent quasi-homogeneous anisotropic medium.The case of a fibrous medium in which the reinforcing elements are rods or filaments [4] is discussed in detail. Allowance for moment effects leads to equations analogous to the equations of the Voight-Cosserat moment theory and its later generalizations. Similar equations are obtained for the case of laminated media, where the reinforcing elements are membranes or plates. On the basis of the viscoelastic analogy [7], the equations of the theory of reinforced media are extended to include the case in which the matrix and/or reinforcing materials are linear viscoelastic.Mekhanika Polimerov, Vol 1, No. 2, pp. 27–37, 1965     

On the integral representation formula for a two‐component elastic composite

The aim of this paper is to derive, in the Hilbert space setting, an integral representation formula for the effective elasticity tensor for a two‐component composite of elastic materials, not necessarily well‐ordered. This integral representation formula implies a relation which links the effective elastic moduli to the N‐point correlation functions of the microstructure. Such relation not only facilitates a powerful scheme for systematic incorporation of microstructural information into bounds on the effective elastic moduli but also provides a theoretical foundation for inverse‐homogenization. The analysis presented in this paper can be generalized to an n‐component composite of elastic materials. The relations developed here can be applied to the inverse‐homogenization for a special class of linear viscoelastic composites. The results will be presented in another paper. Copyright © 2005 John Wiley & Sons, Ltd.     

非线性弹性体的弹性力学变分原理

作者自1978年以后,曾发表了一系列有关弹性力学的变分原理和广义变分原理的文章如[1](1978),[6](1980),[2]、[3](1983),[4]、[5](1984),都是指线性应力应变关系的线性弹性体的.在1985年出版的广义变分原理中,初步推广至非线性弹性体,但并未进行较全面的探讨.本文特别讨论非线性应力应变关系的弹性体的变分原理和广义变分原理,这里有不少问题是值得注意的,有时,它对线性弹性体的变分原理,有指导意义.当应变很小,其高次项可以略去时,本文所得结论,都能近似地化简为通常线性理论的结果.     

On the propagation of elastic waves in two-phase media

At the present time a number of papers has been already devoted to the dynamics of two-phase media. One may mention the papers by Frenkel' [1], Rakhmatulin [2], Biot [3,4], Zwikker and Kosten [5], and others. However, the basic problem of the setting up of the equations of motion in two-phase media still cannot be considered solved and requires additional study and experimental verification.

This paper is concerned with the study of the simplest case of motion, which is the propagation of elastic waves in a homogeneous isotropic medium consisting of a solid and a fluid phase. The problems of the reflection of plane waves and surface waves at the free boundary of the half-space are solved. It is shown that the stress-strain relations established by Frenkel' are equivalent to the analogous relations proposed by Biot and that the equations of motion of the latter are more general.     

Formulating an <Emphasis Type="Italic">n</Emphasis>-person noncooperative game as a tensor complementarity problem

In this paper, we consider a class of n-person noncooperative games, where the utility function of every player is given by a homogeneous polynomial defined by the payoff tensor of that player, which is a natural extension of the bimatrix game where the utility function of every player is given by a quadratic form defined by the payoff matrix of that player. We will call such a problem the multilinear game. We reformulate the multilinear game as a tensor complementarity problem, a generalization of the linear complementarity problem; and show that finding a Nash equilibrium point of the multilinear game is equivalent to finding a solution of the resulted tensor complementarity problem. Especially, we present an explicit relationship between the solutions of the multilinear game and the tensor complementarity problem, which builds a bridge between these two classes of problems. We also apply a smoothing-type algorithm to solve the resulted tensor complementarity problem and give some preliminary numerical results for solving the multilinear games.     

On the representation formula for well‐ordered elastic composites: a convergence of measure approach

The aim of this paper is to derive an integral representation formula for the effective elasticity tensor for a two‐component well‐ordered composite of elastic materials without using a third reference medium and without assuming the completeness of the eigenspace of the operator ? defined in (2.16) in (J. Mech. Phys. Solids 1984; 32 (1):41–62). As shown in (J. Mech. Phys. Solids 1984; 32 (1):41–62) and (Math. Meth. Appl. Sci. 2006; 29 (6):655–664), this integral representation formula implies a relation which links the effective elastic moduli to the N‐point correlation functions of the microstructure. Such relation not only facilitates a powerful scheme for systematic incorporation of microstructural information into bounds on the effective elastic moduli but also provides a theoretical foundation for de‐homogenization. The analysis presented in this paper can be generalized to an n‐component composite of elastic materials. The relations developed here can be applied to the de‐homogenization for a special class of linear viscoelastic composites. The results will be presented in another paper. Copyright © 2006 John Wiley & Sons, Ltd.     

CP decomposition for tensors via alternating least squares with QR decomposition

The CP tensor decomposition is used in applications such as machine learning and signal processing to discover latent low-rank structure in multidimensional data. Computing a CP decomposition via an alternating least squares (ALS) method reduces the problem to several linear least squares problems. The standard way to solve these linear least squares subproblems is to use the normal equations, which inherit special tensor structure that can be exploited for computational efficiency. However, the normal equations are sensitive to numerical ill-conditioning, which can compromise the results of the decomposition. In this paper, we develop versions of the CP-ALS algorithm using the QR decomposition and the singular value decomposition, which are more numerically stable than the normal equations, to solve the linear least squares problems. Our algorithms utilize the tensor structure of the CP-ALS subproblems efficiently, have the same complexity as the standard CP-ALS algorithm when the input is dense and the rank is small, and are shown via examples to produce more stable results when ill-conditioning is present. Our MATLAB implementation achieves the same running time as the standard algorithm for small ranks, and we show that the new methods can obtain lower approximation error.     

Effect of random irregularities of the creep of reinforced layered plastics

The authors investigate the creep of inhomogeneous materials consisting of a large number of stiff orthotropic elastic layers alternating with layers of linear isotropic viscoelastic material. The elastic layers are assumed to be almost plane; the functions describing the irregularities (curvature) form a random field. The averaged characteristics of the medium are found together with the variation of the averaged displacements and strains in time. An analogous problem was previously considered in [1, 6] on the assumption that the binder layers are elastic. The present paper is based on the equations of [1] and the elastic-viscoelastic correspondence principle [4]. When the correlation scales of the irregularities are small as compared with the dimensions of the body and the characteristic distances over which the averaged parameters of the stress-strain state vary appreciably is considered in detail. A relation is established between the creep functions for simple cases of the state of stress and the parameters characterizing the properties of the components, the properties of the random field of initial irregularities, etc. The development of perturbations with different wave numbers is investigated. The theory is used to describe the creep of reinforced layered plastics.Mekhanika Polimerov, Vol. 2, No. 5, pp. 755–762, 1966     

Quasistatic anti-plane motion in the simplest theory of nonlinear viscoelasticity

We consider a very simple model in the framework of differential viscoelastic materials which are isotropic and incompressible. In this model the Cauchy stress tensor is split in an elastic part and a dissipative part. The elastic part is derived from a strain-energy density function only of the first invariant of the Cauchy–Green strain tensor. The dissipative part is like the Navier–Stokes equations: linear in the stretching tensor with a constant viscosity parameter. For this model we provide some time and spatial estimates in the quasistatic approximations for the equations governing anti-plane shear motions. Several explicit examples for specific form of the strain energy are produced. Our results impose analytical restrictions on the mathematical properties of the strain energy to ensure a physical behavior in the creep and recovery experiments. Moreover, we show polynomial decay for the spatial behavior in the class of stress-hardening (or strain-stiffening) materials. For stress-softening materials a Phragmen–Lindelof alternative is proved.     

Three-dimensional travelling waves for nonlinear viscoelastic materials with memory

In this paper, we show the existence of nontrivial travelling wave solutions, which propagating speeds are between ones determined by the equilibrium and instantaneous elastic tensor, respectively, to the nonlinear three-dimensional viscoelastic system with fading memory.