Analytic solutions of a class of nonlinear partial differential equations

An approach is presented for computing the adjoint operator vector of a class of nonlinear （that is, partial-nonlinear） operator matrices by using the properties of conjugate operators to generalize a previous method proposed by the author. A unified theory is then given to solve a class of nonlinear （partial-nonlinear and including all linear） and non-homogeneous differential equations with a mathematical mechanization method. In other words, a transformation is constructed by homogenization and triangulation, which reduces the original system to a simpler diagonal system. Applications are given to solve some elasticity equations.     

Initial-boundary value problem of one class of nonlinear Schrödinger equations described in molecular crystals

In this paper, we study the initiial-boundary value problem of one class of nonlinear Schrodinger equations described in molecular crystals. Furthermore, the existence of the global solution is obtained by means of interpolation inequality and a priori estimation.     

Superconvergence analysis of lower order anisotropic finite element

The convergence analysis of the lower order nonconforming element pro- posed by Park and Sheen is applied to the second-order elliptic problem under anisotropic meshes.The corresponding error estimation is obtained.Moreover,by using the interpo- lation postprocessing technique,a global superconvergence property for the discretization error of the postprocessed discrete solution to the solution itself is derived.Numerical results are also given to verify the theoretical analysis.     

Hamilton原理四种表达形式的异同

Ｈａｍｉｌｔｏｎ原理四种表达形式的异同陈滨（北京大学力学系北京１００８７１）１．引言Ｈａｍｉｌｔｏｎ原理是经典力学中的重要原理。它是以泛函（Ｈａｍｉｌｔｏｎ作用量）驻值的变分形式给出的力学系统动力学原理,并可以作为整个力学系统动力学研究的基础,Ｈａｍ...     

CAVITATION BIFURCATION FOR COMPRESSIBLE ANISOTROPIC HYPERELASTIC MATERIALS

I. INTRODUCTION The cavitation bifurcation problem, sudden formation and growth of voids in solid materials, haslong attracted much attention because of the fundamental role it plays in the local failure and fractureof materials. For hyperelastic materi…     

损伤体有效弹性性质的细观分析和不变性描述—— 一个考虑微缺陷相互作用的一般理论模式

通过引进微缺陷相互作用张量，建立了一个二维情况下考虑微缺陷（微裂纹或微孔洞）间相互作用的损伤固体有效弹性性质的一般理论模式模型中考虑了微缺陷的几何形状、取向分布和空间分布所造成的有效柔度张量的各向异性和材料中微缺陷之间的相互作用所引起的损伤柔度张量的高阶效应针对微椭圆孔、微圆孔和微裂纹问题，求得了相互作用张量的解析形式     

各向异性对平板裂纹扩展路径的影响研究

近场动力学理论克服了经典连续介质力学在模拟裂纹扩展方面的不足。本文基于微梁键构建了一般各向异性近场动力学平面模型,并探讨了各向异性对于裂纹扩展路径的影响。首先,建立一般各向异性的微梁键,并基于插值法建立了一般各向异性近场动力学平面模型的应变能密度表达式;然后,利用应变能密度互等原理求解出近场动力学平面本构模型的一般各向异性参数;最后,通过带初始切口的平板断裂实验研究了不同各向异性参数对于平板裂纹扩展路径的影响。     

Micromechanics of breakage in sharp-edge particles using combined DEM and FEM

By combining DEM （Discrete Element Method） and FEM （Finite Element Method）, a model is established to simulate the breakage of twodimensional sharp-edge particles, in which the simulated particles are assumed to have no cracks. Particles can, however, crush during different stages of the numerical analysis, if stress-based breakage criteria are fulfilled inside the particles. With this model, it is possible to study the influence of particle breakage on macro- and micro-mechanical behavior of simulated angular materials. Two series of tests, with and without breakable particles, are simulated under different confining pressures based on conditions of biaxial tests. The results, presented in terms of micromechanical behavior for different confining pressures, are compared with macroparameters. The influence of particle breakage on microstructure of sharp-edge materials is discussed and the related confining pressure effects are investigated. Breakage of particles in rockfill materials are shown to reduce the anisotropy coefficients of the samples and therefore their strength and dilation behaviors.     

A streamline diffusion nonconforming finite element method for the time-dependent linearized Navier-Stokes equations

A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming （P1）2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.     

Least-squares mixed finite element method for a class of stokes equation

ＩｎｔｒｏｄｕｃｔｉｏｎＡｇｅｎｅｒａｌｔｈｅｏｒｙｏｆｔｈｅｌｅａｓｔ＿ｓｑｕａｒｅｓｍｅｔｈｏｄｈａｓｂｅｅｎｄｅｖｅｌｏｐｅｄｂｙＡＫＡｚｉｚ,ＲＢＫｅｌｌｏｇｇａｎｄＡＢＳｔｅｐｈｅｎｓｉｎ[1].Ｔｈｅｍｏｓｔｉｍｐｏｒｔａｎｔａｄｖａｎｔａｇｅｌｅａｄｓｔｏａｓｙｍｍｅｔｒｉｃｐｏｓｉｔｉｖｅｄｅｆｉｎｉｔｅｐｒｏｂｌｅｍ.ＪＨＢｒａｍｂｌｅａｎｄＪＡＮｉｔｓｈｅｐｒｅｓｅｎｔｅｄａｌｅａｓｔ＿ｓｑｕａｒｅｓｍｅｔｈｏｄｆｏｒＤｉｒｉｃｈｌｅｔｐｒｏｂｌｅｍｓｉｎ[2].Ｔｈｅｍｅｔｈｏｄｇｅ…     

Least-squares finite element methods for compressible Euler equations

A method based on backward finite differencing in time and a least-squares finite element scheme for first-order systems of partial differential equations in space is applied to the Euler equations for gas dynamics. The scheme minimizes the L²-norm of the residual within each time step. The method naturally generates numerical dissipation proportional to the time step size. An implicit method employing linear elements has been implemented and proves robust. For high-order elements, computed solutions based on the L²-method may have oscillations for calculations at similar time step sizes. To overcome this difficulty, a scheme which minimizes the weighted H¹-norm of the residual is proposed and leads to a successful scheme with high-degree elements. Finally, a conservative least-squares finite element method is also developed. Numerical results for two-dimensional problems are given to demonstrate the shock resolution of the methods and compare different approaches.     

Mathematical analysis of EEP method for one-dimensional finite element postprocessing

For a class of two-point boundary value problems, by virtue of one- dimensional projection interpolation, it is proved that the nodal recovery derivative obtained by Yuan's element energy projection (EEP) method has the accuracy O(hmin{2k,k 4}). The theoretical analysis coincides the reported numerical results.     

Perturbational solutions for fuzzy-stochastic finite element equilibrium equations (FSFEEE)

I.IntroductionInthefieldofengineeringstructureanalysis,suchparametersasmodulusofmaterial-Poissonratio,geometricalsizesofstructure,boundaryconditionsandexternalloadingsprobablyhavefuzzinessor/andrandomness.Thestochasticfiniteelementmethod(SFEM)1'1disposeso…     

Non-oscillatory shock-capturing finite element methods for the one-dimensional compressible Euler equations

A class of shock-capturing Petrov–Galerkin finite element methods that use high-order non-oscillatory interpolations is presented for the one-dimensional compressible Euler equations. Modified eigenvalues which employ total variation diminishing (TVD), total variation bounded (TVB) and essentially non-oscillatory (ENO) mechanisms are introduced into the weighting functions. A one-pass Euler explicit transient algorithm with lumped mass matrix is used to integrate the equations. Numerical experiments with Burgers' equation, the Riemann problem and the two-blast-wave interaction problem are presented. Results indicate that accurate solutions in smooth regions and sharp and non-oscillatory solutions at discontinuities are obtainable even for strong shocks.     

Three dimensional finite element analysis for Rayleight-Benard convection in rectangular enclosures

The pattern of Rayleigh-Benard convection of air in a rectangular box heated-from-below is studied by numerically solving the three-dimensional time-dependent Navier-Stokes equations under the Boussinesq approximation. Slightly supercritical Rayleigh number was adopted to track the evolutions of flow structure as a function of enclosure's aspect ratio (A=L/H). The flow will asymptotically evolve to different patterns, among which, two possible types of flow pattern are found. One consists of the pair of straight vortex rolls and the other appears as closed vortex rings. The transition between the flow patterns indicates that there exists a flow bifurcation with the variation of container's aspect ratio. In addition, both steady and oscillatory flows have been observed, corresponding to the pair of straight vortex rolls and the vortex ring, respectively. The complexity of flow structure tends to increase with the increasing aspect ratio of the rectangular enclosure. The project supported by the National Natural Science Foundation of China (19889210), the National Distinguished Young Fund (10125210), the Hundred Talents Program of CAS, and the Training Program for the Trans-Century Outstanding Young of MOE     

A parallel two-level finite element method for the Navier-Stokes equations

Based on domain decomposition, a parallel two-level finite element method for the stationary Navier-Stokes equations is proposed and analyzed. The basic idea of the method is first to solve the Navier-Stokes equations on a coarse grid, then to solve the resulted residual equations in parallel on a fine grid. This method has low communication complexity. It can be implemented easily. By local a priori error estimate for finite element discretizations, error bounds of the approximate solution are derived. Numerical results are also given to illustrate the high efficiency of the method.     

A finite element method for the one‐dimensional extended Boussinesq equations

A new finite element method for Nwogu's (O. Nwogu, ASCE J. Waterw., Port, Coast., Ocean Eng., 119 , 618–638 (1993)) one‐dimensional extended Boussinesq equations is presented using a linear element spatial discretisation method coupled with a sophisticated adaptive time integration package. The accuracy of the scheme is compared to that of an existing finite difference method (G. Wei and J.T. Kirby, ASCE J. Waterw., Port, Coast., Ocean Eng., 121 , 251–261 (1995)) by considering the truncation error at a node. Numerical tests with solitary and regular waves propagating in variable depth environments are compared with theoretical and experimental data. The accuracy of the results confirms the analytical prediction and shows that the new approach competes well with existing finite difference methods. The finite element formulation is shown to enable the method to be extended to irregular meshes in one dimension and has the potential to allow for extension to the important practical case of unstructured triangular meshes in two dimensions. This latter case is discussed. Copyright © 1999 John Wiley & Sons, Ltd.