Nonlinear Schrdinger equation with combined power-type nonlinearities and harmonic potential

This paper discusses a class of nonlinear Schrdinger equations with combined power-type nonlinearities and harmonic potential. By constructing a variational problem the potential well method is applied. The structure of the potential well and the properties of the depth function are given. The invariance of some sets for the problem is shown. It is proven that, if the initial data are in the potential well or out of it, the solutions will lie in the potential well or lie out of it, respectively. By the convexity method, the sharp condition of the global well-posedness is given.     

General existence theorems for unilateral problems in continuum mechanics

The problem of minimizing a possibly non-convex and non-coercive functional is studied. Either necessary or sufficient conditions for the existence of solutions are given, involving a generalized recession functional, whose properties are discussed thoroughly. The abstract results are applied to find existence of equilibrium configurations of a deformable body subject to a system of applied forces and partially constrained to lie inside a possibly unbounded region.     

Stability,instability and interaction of solitary pulses in a composite medium

The questions of a dynamical stability and instability of soliton-like solutions (solitary pulses) of the Hamiltonian equations, describing planar waves in nonlinear elastic composites are considered, both in the presence as well as in the absence of the anisotropy. In the anisotropic case one has the slow and the fast two-parametric soliton families on the background of the quiescent state. In the absence of the anisotropy these two families coalesce into the unique three parametric family. It was shown recently that solitary pulses of the slow family in the anisotropic composite and pulses in the isotropic composite are stable when their speeds lie inside a certain range, the so-called range of stability. In the present paper, on the basis of numerical solving of the Cauchy problem for the basic governing equations, the classification is given of the types of instability of solitary pulses from the fast family for all range of speeds as well as in the case of the slow family and in the isotropic case, when the speeds of the pulses lie without the range of stability. The first type of instability is the blow-up instability for the slow anisotropic and isotropic pulses, living without the range of stability and also for high amplitude fast anisotropic pulses. The second type of instability is the instability resulting in energy exchange between the components of strain tensor for low amplitude fast anisotropic solitary pulses. The reasons of the both types of instability are discussed in detail.The interaction between the pairs of solitary pulses of different nature is investigated both analytically as well as numerically. It is found out that solitary pulses having the different polarization, i.e. different sign of amplitudes, can form bound states, oscillating about the common center, subjected to a motion with a constant speed, approximately equal to the average of speeds of two pulses when they are far apart.     

Simplified Gurtin-type generalized variational principles for fully dynamic magneto-electro-elasticity with geometrical nonlinearity

The fundamental equations, governing all the variables of the initial boundary value problem in fully dynamic magneto-electro-elasticity with geometrical nonlinearity, are expressed in covariant differential form. The generalized principle of virtual work is given in terms of convolutions for the present problem. Two simplified Gurtin-type generalized variational principles, directly leading to all the fundamental equations, are deduced by using He’s semi-inverse method instead of Laplace transforms. By enforcing some fundamental equations as constraint conditions, one of various constrained variational principles is given as an example. By simply dropping out selected field functions, several reduced variational principles are obtained as special forms for piezoelectricity, elastodynamics, and electromagnetics, respectively. This paper aims at providing a more complete theoretical foundation for the finite element applications for the discussed problem.     

关于可变区域问题的变分法以及广义数学规划问题

本文具体研究了定义在可变区域上或可变边界表面上的泛函的一阶、二阶变分问题，得到了与经典变分法相对应的关于可变区域问题的变分法。并用该变分法讨论了具有可变区域的弹性系统的势能原理。另一方面，与传统的数学规划相对应，研究了可变区域上泛函的约束极值问题——广义数学规划问题，给出了相应的广义Ｋｕｈｎ－Ｔｕｃｋｅｒ条件。     

Stability of a nonorthogonal stagnation flow to three-dimensional disturbances

A similarity solution for a low Mach number nonorthogonal flow impinging on a hot or cold plate is presented. For the constant-density case, it is known that the stagnation point shifts in the direction of the incoming flow and that this shift increases as the angle of attack decreases. When the effects of density variations are included, a critical plate temperature exists; above this temperature the stagnation point shifts away from the incoming stream as the angle is decreased. This flow field is believed to have applications to the reattachment zone of certain separated flows or to a lifting body at a high angle of attack. Finally, we examine the stability of this nonorthogonal flow to self-similar, three-dimensional disturbances. Stability characteristics of the flow are given as a function of the parameters of this study: ratio of the plate temperature to that of the outer potential flow and angle of attack. In particular, it is shown that the angle of attack can be scaled out by a suitable definition of an equivalent wave number and temporal growth rate, and the stability problem for the nonorthogonal case is identical to the stability problem for the orthogonal case. By use of this scaling, it can be shown that decreasing the angle of attack decreases the wave number and the magnitude of the temporal decay rate, thus making nonlinear effects important. For small wave numbers, it is shown that cooling the plate decreases the temporal decay of the least-stable mode, while heating the plate has the opposite effect. For moderate to large wave numbers, density variations have little effect except that there exists a range of cool plate temperatures for which these disturbances are extremely stable.This work was supported by the National Aeronautics and Space Administration under NASA Contract NAS1-18605 while the authors were in residence at the Institute for Compute Applications in Science and Engineering, NASA Langley Research Center, Hampton, VA 23665, U.S.A.     

大型陀螺特征值问题的广义Arnoldi减缩算法

基于Ａｒｎｏｌｄｉ法，建立陀螺特征值问题的广义Ａｒｎｏｌｄｉ格式，并利用系统矩阵的反对称特性，得到极其简洁的甚至比对称矩阵Ｌａｎｃｚｏｓ法更为简单的递推格式，可称为陀螺Ａｒｎｏｌｄｉ减缩算法。     

Phonon dispersion analysis of carbon nanotubes based on inter-belt model and symplectic solution method

The objective of the study is to develop a totally new theory and structural mechanics model for phonon dispersion analysis of carbon nanotubes. The fundamental theory and computational algorithm for phonon dispersion analysis of carbon nanotubes are developed based on the symplectic theory and algorithm established in applied mechanics in recent years. Carbon nanotubes are simulated by two kinds of structural mechanics models, i.e. the conventional sub-structure model and the inter-belt model. The variational principle for wave dissipation analysis of periodic structure is given on the basis of the symplectic-mathematical theory. Numerical examples are carried out to demonstrate the validity of the theory and algorithm developed. By the comparison of the results obtained by the two kinds of structural mechanics models, it can be found that the inter-belt model has more advantages than the conventional sub-structure model in the calculation of phonon spectra of nanotubes. As a basic research work, the present study illustrates well the potential of the symplectic-mathematical theory as well as the inter-belt model and is valuable for the further research in computational nanomechanics.     

ANALYSIS OF THE LARGE DEFLECTION DYNAMIC RESPONSE OF RIGID-PLASTIC CIRCULAR PLATE RESTING ON FLUID

A study is presented for the large deflection dynamic response of rigid-plastic circular plate resting on potential fluid under a rectangular pressure pulse load.By virtue of Hankel integral transform technique,this interaction problem is reduced toa problem of dynamic plastic response of the plate in vacuum.The closed-formsolutions are derived for both middle and high pressure loads by solving the equationsof motion with the large deflection in the range where both bending moments andmembrane forces are important.Some numerical results are given.     

虚拟测试的理论体系与系统技术研究

从工程背景的一般性测试问题出发,提出虚拟测试的概念,将虚拟现实与测试技术结合在一起,形成虚拟测试研究领域。指出虚拟测试研究的理论、技术和应用背景在于虚拟现实和可视化与一般测试问题的结合。界定其研究内容和技术特征,并阐明其研究方法和实验系统。对于虚拟现实环境下的测试系统、测试对象和测试信息智能模型的建立与合成方法的研究,以及虚拟现实的其它相关应用领域具有推广意义。     

置于液面上的刚塑性圆板大挠度动力响应分析

分析了置于无旋不可压理想流体流面上的简支刚塑性圆板受矩形脉冲载荷作用的大挠度动力响应，借助Ｈａｎｋｅｌ变换，将液－固耦合作用为在空气中的圆板塑性动力响应问题，进而求解弯矩和膜力联合作用的大挠度运动方程，得到了中载及高载下各相运动的完全解，并提供了数值算例。     

Hamiltonian structures of dynamics of a gyrostat in a gravitational field

The dynamics of a gyrostat in a gravitational field is a fundamental problem in celestial mechanics and space engineering. This paper investigates this problem in the framework of geometric mechanics. Based on the natural symplectic structure, non-canonical Hamiltonian structures of this problem are derived in different sets of coordinates of the phase space. These different coordinates are suitable for different applications. Corresponding Poisson tensors and Casimir functions, which govern the phase flow and phase space structures of the system, are obtained in a differential geometric method. Equations of motion, as well as expressions of the force and torque, are derived in terms of potential derivatives. We uncover the underlying Lie group framework of the problem, and we also provide a systemic approach for equations of motion. By assuming that the gravitational field is axis-symmetrical and central, SO(2) and SO(3) symmetries are introduced into the general problem respectively. Using these symmetries, we carry out two reduction processes and work out the Poisson tensors of the reduced systems. Our results in the central gravitational filed are in consistent with previous results. By these reductions, we show how the symmetry of the problem affects the phase space structures. The tools of geometric mechanics used here provide an access to several powerful techniques, such as the determination of relative equilibria on the reduced system, the energy-Casimir method for determining the stability of equilibria, the variational integrators for greater accuracy in the numerical simulation and the geometric control theory for control problems.     

伸缩虚拟边界元法解二维Helmholtz外问题

以位势理论为基础，提出了求解Helmholtz外问题的伸缩虚拟边界元法．给出了该方法在全波数域内获得唯一解的严格数学证明，其核心是通过伸缩虚拟边界使对偶内问题的特征频率(本征值)避开与波数重合，从而保证了解的唯一性，同以往前人提出的几种解法途径相比，该法简单得多；通过诸多边界曲线形状和不同边界量的声辐射算例，从计算精度、稳定性以及克服解的非唯一性等方面，对该方法进行了检验．计算结果表明：对远场或近场辐射声压，该方法都具有非常高的效率和精度．     

Wave equations and reaction-diffusion equations with several nonlinear source terms

The initial boundary value problem of wave equations and reaction-diffusion equations with several nonlinear source terms in a bounded domain is studied by potential well method.The invariance of some sets under the flow of these problems and the vac- uum isolation of solutions are obtained by introducing a family of potential wells.Then the threshold result of global existence and nonexistence of solutions are given.Finally, the problem with critical initial conditions are discussed.     

Periodic set of limited electrically permeable interface cracks with contact zones

Plane problem for an infinite space composed of two different piezoelectric or piezoelectric/dielectric semi-infinite spaces with a periodic set of limited electrically permeable interface cracks is considered. Uniformly distributed electromechanical loading is applied at infinity. The frictionless contact zones at the crack tips are taken into account. The problem is reduced to the combined Dirichlet–Riemann boundary value problem by means of the electromechanical factors presentation via sectionally analytic functions, assuming that the electric flux is uniformly distributed inside the cracks. An exact solution of the problem is proposed. It permits to find in a closed form all necessary electromechanical characteristics at the interface and to formulate the equation for the determination of the electric flux value. Analysis of this equation confirms the correctness of the assumption concerning the uniform distribution of the electric flux in the crack region.Formulae for stresses, electric displacement vector, elastic displacements and electric potential jump at the interface as well as the intensity factors at the crack tips are given. Equation for the contact zone length determination is presented. Calculations for certain material combinations are carried out. The influence of electric permeability of cracks on electromechanical fields and the fracture mechanical parameters is analyzed.     

Thermal stresses in an elastic space with a perfectly rigid flat inclusion under perpendicular heat flow

This paper examines the three-dimensional problem of finding thermal stresses due to an insulated rigid sheet-like inclusion (anticrack) in an elastic space under a uniform perpendicular heat flow. By using appropriate harmonic potentials, a general method of solving this problem is presented. The resulting boundary-value problems are reduced to classical mixed problems of potential theory. For the purpose of illustration, a complete solution in terms of elementary functions for a rigid circularly shaped inclusion is given and discussed from the point of view of material failure.     

Modelling the cracks produced by settlements in masonry structures

The present work is concerned with the prediction of the crack pattern produced by given kinematical data (settlements/distortions) in masonry constructions. By adopting the simplified model of Heyman, extending it to masonry structures treated as continuous bodies, we restrict the analysis to the search of displacement fields which are piecewise rigid. Restricting to small strains and displacements we look for the solution of the kinematical problem by minimizing the potential energy of the structure. A variational approximation of the minimum problem is obtained by considering a fixed finite element subdivision of the structure into rigid blocks. Two case studies are presented to illustrate the way in which a particular fracture pattern can be identified as the one associated to the minimum of the energy in this restricted class of piecewise rigid displacements.     

Two dimensional half-space problems for the Broadwell discrete velocity model

Stationary boundary value problems for the Broadwell model in a half-space and in a half-infinite channel are considered. By means of the analogy between the stationary boundary value problems for the Broadwell equations and the initial-boundary value problem of Carleman's system, solutions are found for various situations. Uniqueness and non-uniqueness of solutions is discussed as well. The non-uniqueness problem in the channel leads to the investigation of the initial value problem for Carleman's equation with partly negative initial densities. Some new results for this problem are given. Received January 20, 1996     

Local and parallel finite element algorithms based on two-grid discretization for steady Navier-Stokes equations

Local and parallel finite element algorithms based on two-grid discretization for Navier-Stokes equations in two dimension are presented. Its basis is a coarse finite element space on the global domain and a fine finite element space on the subdomain. The local algorithm consists of finding a solution for a given nonlinear problem in the coarse finite element space and a solution for a linear problem in the fine finite element space, then droping the coarse solution of the region near the boundary. By overlapping domain decomposition, the parallel algorithms are obtained. This paper analyzes the error of these algorithms and gets some error estimates which are better than those of the standard finite element method. The numerical experiments are given too. By analyzing and comparing these results, it is shown that these algorithms are correct and high efficient.     

Three-dimensional dynamic Green’s functions in transversely isotropic bi-materials

By virtue of a complete representation using two displacement potentials, an analytical derivation of the elastodynamic Green’s functions for a linear elastic transversely isotropic bi-material full-space is presented. Three-dimensional point-load Green’s functions for stresses and displacements are given in complex-plane line-integral representations. The formulation includes a complete set of transformed stress–potential and displacement–potential relations, within the framework of Fourier expansions and Hankel integral transforms, that is useful in a variety of elastodynamic as well as elastostatic problems. For numerical computation of the integrals, a robust and effective methodology is laid out which gives the necessary account of the presence of singularities including branch points and pole on the path of integration. As illustrations, the present Green’s functions are analytically degenerated to the special cases such as half-space, surface and full-space Green’s functions. Some typical numerical examples are also given to show the general features of the bi-material Green’s functions.