Development of a finite-element analysis methodology for the propagation of delaminations in composite structures

Analysing the collapse of skin-stiffened structures requires capturing the critical phenomenon of skin-stiffener separation, which can be considered analogous to interlaminar cracking. This paper presents the development of a numerical approach for simulating the propagation of interlaminar cracks in composite structures. A degradation methodology was introduced in MSC.Marc, which involved the modelling of a structure with shell layers connected by user-defined multiple-point constraints (MPCs). User subroutines were written that employ the virtual crack closure technique (VCCT) to determine the onset of crack growth and modify the properties of the user-defined MPCs to simulate crack propagation. Methodologies for the release of failing MPCs are presented and are discussed with reference to the VCCT assumption of self-similar crack growth. The numerical results obtained by using the release methodologies are then compared with experimental data for a double-cantilever beam specimen. Based on this comparison, recommendations for the future development of the degradation model are made, especially with reference to developing an approach for the collapse analysis of fuselage-representative structures. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 43, No. 1, pp. 15–42, January–February, 2007.     

An improved linearization strategy for zero-one quadratic programming problems

We present a new linearized model for the zero-one quadratic programming problem, whose size is linear in terms of the number of variables in the original nonlinear problem. Our derivation yields three alternative reformulations, each varying in model size and tightness. We show that our models are at least as tight as the one recently proposed in [7], and examine the theoretical relationship of our models to a standard linearization of the zero-one quadratic programming problem. Finally, we demonstrate the efficacy of solving each of these models on a set of randomly generated test instances.     

A parameter-uniform numerical method for a Sobolev problem with initial layer

The present study is concerned with the numerical solution, using finite difference method of a one-dimensional initial-boundary value problem for a linear Sobolev or pseudo-parabolic equation with initial jump. In order to obtain an efficient method, to provide good approximations with independence of the perturbation parameter, we have developed a numerical method which combines a finite difference spatial discretization on uniform mesh and the implicit rule on Shishkin mesh(S-mesh) for the time variable. The fully discrete scheme is shown to be convergent of order two in space and of order one expect for a logarithmic factor in time, uniformly in the singular perturbation parameter. Some numerical results confirming the expected behavior of the method are shown.      

High order generalized upwind schemes and numerical solution of singular perturbation problems

High even order generalizations of the traditional upwind method are introduced to solve second order ODE-BVPs without recasting the problem as a first order system. Both theoretical analysis and numerical comparison with central difference schemes of the same order show that these new methods may avoid typical oscillations and achieve high accuracy. Singular perturbation problems are taken into account to emphasize the main features of the proposed methods. AMS subject classification (2000)  65L10, 65L12, 65L50     

Maximum-norm error analysis of a non-monotone FEM for a singularly perturbed reaction-diffusion problem

A non-monotone FEM discretization of a singularly perturbed one-dimensional reaction-diffusion problem whose solution exhibits strong layers is considered. The method is shown to be maximum-norm stable although it is not inverse monotone. Both a priori and a posteriori error bounds in the maximum norm are derived. The a priori result can be used to deduce uniform convergence of various layer-adapted meshes proposed in the literature. Numerical experiments complement the theoretical results. AMS subject classification (2000)  65L10, 65L50, 65L60     

Integral representations of nonnegative solutions for parabolic equations and elliptic Martin boundaries

We consider nonnegative solutions of a parabolic equation in a cylinder D×(0,T), where D is a noncompact domain of a Riemannian manifold. Under the assumption [IU] (i.e., the associated heat kernel is intrinsically ultracontractive), we establish an integral representation theorem: any nonnegative solution is represented uniquely by an integral on (D×{0})∪(M_∂D×[0,T)), where M_∂D is the Martin boundary of D for the associated elliptic operator. We apply it in a unified way to several concrete examples to explicitly represent nonnegative solutions. We also show that [IU] implies the condition [SP] (i.e., the constant function 1 is a small perturbation of the elliptic operator on D).     

Generalized 3G theorem and application to relativistic stable process on non-smooth open sets

Let G(x,y) and GD(x,y) be the Green functions of rotationally invariant symmetric α-stable process in Rd and in an open set D, respectively, where 0<α<2. The inequality GD(x,y)GD(y,z)/GD(x,z)?c(G(x,y)+G(y,z)) is a very useful tool in studying (local) Schrödinger operators. When the above inequality is true with c=c(D)∈(0,∞), then we say that the 3G theorem holds in D. In this paper, we establish a generalized version of 3G theorem when D is a bounded κ-fat open set, which includes a bounded John domain. The 3G we consider is of the form GD(x,y)GD(z,w)/GD(x,w), where y may be different from z. When y=z, we recover the usual 3G. The 3G form GD(x,y)GD(z,w)/GD(x,w) appears in non-local Schrödinger operator theory. Using our generalized 3G theorem, we give a concrete class of functions belonging to the non-local Kato class, introduced by Chen and Song, on κ-fat open sets. As an application, we discuss relativistic α-stable processes (relativistic Hamiltonian when α=1) in κ-fat open sets. We identify the Martin boundary and the minimal Martin boundary with the Euclidean boundary for relativistic α-stable processes in κ-fat open sets. Furthermore, we show that relative Fatou type theorem is true for relativistic stable processes in κ-fat open sets. The main results of this paper hold for a large class of symmetric Markov processes, as are illustrated in the last section of this paper. We also discuss the generalized 3G theorem for a large class of symmetric stable Lévy processes.     

On a new system of nonlinear A-monotone multivalued variational inclusions

In this paper, we introduce and study a new system of nonlinear A-monotone multivalued variational inclusions in Hilbert spaces. By using the concept and properties of A-monotone mappings, and the resolvent operator technique associated with A-monotone mappings due to Verma, we construct a new iterative algorithm for solving this system of nonlinear multivalued variational inclusions associated with A-monotone mappings in Hilbert spaces. We also prove the existence of solutions for the nonlinear multivalued variational inclusions and the convergence of iterative sequences generated by the algorithm. Our results improve and generalize many known corresponding results.     

q-positive definiteness and related operators

We are in progress of extending the family of ‘q-deformed operators’ considered in the previous papers by joining to them q-subnormal as well as q-formally subnormal ones. It turns out that q-positive definiteness, a notion generalizing Halmos' standard positive definiteness of bounded subnormal operators, is likewise central for our new scheme.     

Banach空间中含有无穷多个跳跃点的一阶脉冲积分-微分方程的无穷边值问题

通过建立一个新的比较引理,应用上下解方法和单调迭代技术,研究了Banach空间中含有无穷多个跳跃点的一阶脉冲积分-微分方程无穷边值问题在任意闭区间上最小解和最大解的存在性.