岩土材料应变局部化的有限元分析方法

采用有限单元法分析岩土材料的应变局部化时经常会遇到单元尺寸敏感性问题和网格锁定问题。自适应网格技术能够有效地解决网格锁定问题,但仍然无法克服计算结果对单元尺寸的依赖性,尽管在一维情况下被证明是可行的。复合体理论(均匀化理论)和弱非连续有限元方法可以成功地解决岩土材料的单元尺寸敏感性问题,在一维情况下两类方法实际上是一致的。本文针对岩土材料应变局部化的有限元新技术所存在的若干问题进行了详细的讨论,并给出了有关算例。     

GRCA(1)模型中误差方差自加权估计的渐近分布

考虑随机系数自回归模型Y_t =Φ_ty_t-1+ u_t,其中 Φ_t为随机系数,u_t为随机误差。在允许Φ_t与u_t相依以及Eu_t无穷的条件下,构造了误差方差的自加权估计,并证明了该估计的渐近正态性。最后通过数值模拟,说明 自加权估计的稳健和有效性。     

Global behavior of a discrete haematopoiesis model with several delays

In this paper, we study the difference equation:

     

1-forms over the moduli space of irreducible connections defined by the spectrum of Dirac operators

Let (P) be the moduli space of irreducible connections of a G-principal bundle P over a closed Riemannian spin manifold M. Let D_A be the Dirac operator of M coupled to a connection A of P and f a smooth function on M. We consider a smooth variation A(u) of A with tangent vector ω and denote T_ω:= (D_A(u)−f) (u=0. The coefficients of the asymptotic expansion of trace (T_ω · e^-t(D_A−f)2) near t=0 define 1-forms a^(k)_f, K=0, 1, 2, … on (P). In this paper we calculate aa⁽⁰⁾_f, a⁽¹⁾_f, a⁽²⁾_f and study some of their properties. For instance using the 1-form a⁽²⁾_f for suitable functions f we obtain a foliation of codimension 5 of the space of G-instantons of S⁴.     

Identification of non-minimum phase transfer function using higher-order spectrum

The present paper treats the identification of parametric nonminimum phase transfer function. We propose a method of identification based on the inner outer factorization of stable transfer function. It consists of identifying the outer and inner parts of a transfer function separately. The outer part is identified by the use of the second-order spectral estimate from the observed linear process, while the inner part is identified by the use of a higher-order cumulant spectral estimate from the observed process. Respective parameter estimators are determined in the light of asymptotic efficiency. In order to estimate the order of the inner part of a transfer function, a criterion is proposed. It is introduced based on the same principle as in the case of Akaike's AIC.     

Asymptotic theorems for estimating the distribution function under random truncation

Representation theorem and local asymptotic minimax theorem are derived for nonparametric estimators of the distribution function on the basis of randomly truncated data. The convolution-type representation theorem asserts that the limiting process of any regular estimator of the distribution function is at least as dispersed as the limiting process of the product-limit estimator. The theorems are similar to those results for the complete data case due to Beran (1977, Ann. Statist., 5, 400–404) and for the censored data case due to Wellner (1982, Ann. Statist., 10, 595–602). Both likelihood and functional approaches are considered and the proofs rely on the method of Begun et al. (1983, Ann. Statist., 11, 432–452) with slight modifications.Division of Biostatistics, School of Public Health, Columbia Univ.     

Asymptotic expansions of the distributions of some test statistics for Gaussian ARMA processes

Let {X_t} be a Gaussian ARMA process with spectral density f_θ(λ), where θ is an unknown parameter. The problem considered is that of testing a simple hypothesis H:θ = θ₀ against the alternative A:θ ≠ θ₀. For this problem we propose a class of tests , which contains the likelihood ratio (LR), Wald (W), modified Wald (MW) and Rao (R) tests as special cases. Then we derive the χ² type asymptotic expansion of the distribution of T up to order n⁻¹, where n is the sample size. Also we derive the χ² type asymptotic expansion of the distribution of T under the sequence of alternatives A_n: θ = θ₀ + /√n, ε > 0. Then we compare the local powers of the LR, W, MW, and R tests on the basis of their asymptotic expansions.     

球面介质作用下的测度值分枝过程

我们考虑空间上一粒子系统，当其受到分布于求面上的介质作进行粒子分枝和衍生，产生新了体，而新粒子仍按原粒子的运动规则继续空间运动。通过合理的假设和极限过程，粒子在空间的散布一测度值分枝过程来刻划。     

Rado Partition Theorem for Random Subsets of Integers

For an l x k matrix A = (a_ij) of integers, denote by L(A) thesystem of homogenous linear equations a_i1x₁ + ... + a_ikx_k =0, 1 i l. We say that A is density regular if every subsetof N with positive density, contains a solution to L(A). Fora density regular l x k matrix A, an integer r and a set ofintegers F, we write if for any partition F = F₁ ... F_r there exists i {1, 2,..., r} and a column vector x such that Ax = 0 and all entriesof x belong to F_i. Let [n]_N be a random N-element subset of{1, 2, ..., n} chosen uniformly from among all such subsets.In this paper we determine for every density regular matrixA a parameter = (A) such that lim_n P([n]_N (A)_r)=0 if N =O(n) and 1 if N = (n). 1991 Mathematics Subject Classification:05D10, 11B25, 60C05     

Numerical algorithms for uniform Airy-type asymptotic expansions

Airy-type asymptotic representations of a class of special functions are considered from a numerical point of view. It is well known that the evaluation of the coefficients of the asymptotic series near the transition point is a difficult problem. We discuss two methods for computing the asymptotic series. One method is based on expanding the coefficients of the asymptotic series in Maclaurin series. In the second method we consider auxiliary functions that can be computed more efficiently than the coefficients in the first method, and we do not need the tabulation of many coefficients. The methods are quite general, but the paper concentrates on Bessel functions, in particular on the differential equation of the Bessel functions, which has a turning point character when order and argument of the Bessel functions are equal.