A semi-inverse algorithm in application of computer algebra

ProjectsupportcdbytheNationalNaturalScicnceFoundationofChinaandShanghaiMunicipalNaturaIScicnceF0undationI.Introducti0nCQmputeralgebra(i.e.,symbolic'computation)hasfounditswideapp1icationinscientificcomputationandengineeringana1ysis(cf.Refs.[1~31).Aswaspoi…     

Microstructure observation and mechanical behavior modeling for limnetic nacre

In the present research, microstructure of akind of limnetic shell （Hyriopsis cumingii） is observed and measured by using the scanning electron microscopy, and mechanical behavior experiments of the shell nacre are carried out by using bending and tensile tests. The dependence of mechanical properties of the shell nacre on its microstructure is analyzed by using a modified shear-lag model, and the overall stress-strain relation is obtained. The experimental results reveal that the mechanical properties of shell nacre strongly depend on the water contents of the limnetic shell. Dry nacre shows a brittle behavior, whereas wetting nacre displays a strong ductility. Compared to the tensile test, the bending test overestimates the strength and underestimates the Young＇s modulus. The modified shear-lag model can characterize the deformation features of nacre effectively.     

一类生物材料界面的结构及其裂纹阻力

 通过对天然珍珠母材料有机基质界面的微结构及其力学性能 的研究，简要地分析了珍珠母有机界面的弹性模量以及裂纹阻力与 其微结构的联系，由此说明在珍珠母所表现出来的优异力学性能中 微结构所起的重要作用.     

MIXED FINITE ELEMENT METHODS FOR THE SHALLOWWATER EQUATIONS INCLUDING CURRENT AND SILTSEDIMENTATION （Ⅱ）——THE DISCRETE-TIMECASE ALONG CHARACTERISTICS

The mixed finite element(MFE) methods for a shallow water equation systemconsisting of water dynamics equations, silt transport equation, and the equation of bottomtopography change were derived. A fully discrete MFE scheme for the discrete-time alongcharacteristics is presented and error estimates are established. The existence andconvergence of MFE solution of the discrete current velocity, elevation of the bottomtopography, thickness of fluid column, and mass rate of sediment is demonstrated.     

Uniformly valid asymptotic solutions of the nonlinear unsymmetrical bending for orthotropic rectangular thin plate of four clamped edges with variable thickness

By using “the method of modified two-variable”,“the method of mixing perturbation ”and introducing four small parameters, the problem of the nonlinear unsymmetrical bending for orthotropic rectangular thin plate with linear variable thickness is studied. And the uniformly valid asymptotic solution of Nth- order for ε1 and Mth- order for ε2 of the deflection functions and stress function are obtained.     

Numerical study of vortex reconnection for two anti-parallel vortex tubes

By direct numerical simulation of the Navier-Stokes equations we investigate the reconnection of two antiparallel vortex tubes. A new type of perturbation of the initial vorticity field is given which is different from that presented in Refs. [8] and [9]. The formation and the evolution of the “curved vortex belts”, their mutual action with the “bridges” are found. These are important phenomena not studied by others. The project supported by the LNM of Institute of Mechanics, Academia Sinica and The National Natural Science Foundation of China     

The coupling dynamical modeling theory of flexible multibody system

Based on the deformation theory of elastic beams, the coupling effect between the coupling displacements of a point on the middle line of beam and large overall motion is presented. The “coupling matrix library” and Jourdain's variation principle and single direction recursive formulation method are used to establish the general coupling dynamical equations of flexible multibody system. Two typical examples show the coupling effect between coupling displacements and large overall motion on the dynamics of flexible multibody system consisting of beams. The project supported by the National Natural Science Foundation of China (No. 19832040).     

Formation Mechanism of Cracks in Saturated Sand

The formation mechanism of “water film” (or crack) in saturated sand is analyzed theoretically and numerically. The theoretical analysis shows that there will be no stable “water film” in the saturated sand if the strength of the skeleton is zero and no positions are choked. It is shown by numerical simulation that stable water films initiate and grow if the choking state keeps unchanged once the fluid velocities decrease to zero in the liquefied sand column. The developments of “water film” based on the model presented in this paper are compared with experimental results.The project supported by the National Natural Science Foundation of China (40025103 and 10202024) and Key Laboratory of Mountain Hazards and Surface Process, Chinese Academy of Sciences. The English text was polished by Keren Wang.     

Particulate size effects in the particle-reinforced metal-matrix composites

The influences of particle size on the mechanical properties of the particulate metal matrix composite are obviously displayed in the experimental observations. However, the phenomenon can not be predicted directly using the conventional elastic-plastic theory. It is because that no length scale parameters are involved in the conventional theory. In the present research, using the strain gradient plasticity theory, a systematic research of the particle size effect in the particulate metal matrix composite is carried out. The roles of many composite factors, such as: the particle size, the Young's modulus of the particle, the particle aspect ratio and volume fraction, as well as the plastic strain hardening exponent of the matrix material, are studied in detail. In order to obtain a general understanding for the composite behavior, two kinds of particle shapes, ellipsoid and cylinder, are considered to check the strength dependence of the smooth or non-smooth particle surface. Finally, the prediction results will be applied to the several experiments about the ceramic particle-reinforced metal-matrix composites. The material length scale parameter is predicted. The project supported by the National Natural Science Foundation of China (19891180, 19925211) and by the Chinese Academy of Sciences (KJ951-1-201) and “Bai Ren” plan     

Comprehensive mathematical model of microcirculatory dynamics (II)—Calculation and the results

The mathematical model described in Part I was solved using “influence line method” combining analytical method and finite element method. Many important aspects of microcirculatory dynamics were analyzed and discussed. It show that interstitial fluid pressure changes its sign twice within one arteriolar vasomotion period and it is therefore not important that interstitial fluid pressure is a little higher or lower than atmospheric pressure; arteriolar vasomotion can periodically result in lymph formation and interstitial total pressure plays an important role in this procedure; local regulation of microcirculation can meet metabolic need some extent in the form of dynamic equilibrium. The property of arteriole as a “resistant vessel” and the efficiency of microvascular network as heat exchanger are also shown. These results show that the comprehensive mathematical model developed in Part I is physiologically resonable. Foundation item: the Natural Science Foundation of Sichuan Province, P R China Biography: Guo Zhongsan (1947-)     

Secondary buckled states of rectangular sandwich plates

In this paper, for a rectangular sandwich plate with edges simply supported and subjected to a constant compressive thrust λ along two opposite edges; the secondary bifurcation points and the secondary buckled states that bifurcate from the primary buckled states are determined by a perturbation method. These results are useful for their numerical calculation and can be used to explain the phenomenon of “mode-jumping”. Project Supported by National Natural Science Foundation of China.     

Block pivot methods for solving frictional contact problems

Based on elementary group theory, the block pivot methods for solving two-dimensional elastic frictional contact problems are presented in this paper. It is proved that the algorithms converge within a finite number of steps when the friction coefficient is “relative small”. Unlike most mathematical programming methods for contact problems, the block pivot methods permit multiple exchanges of basic and nonbasic variables. The project supported by the National Natural Science Foundation of China     

MIXED FINITE ELEMENT METHODS FOR THE SHALLOW WATER EQUATIONS INCLUDING CURRENT AND SILT SEDIMENTATION （Ⅰ）——THE CONTINUOUS-TIME CASE

An initial-boundary value problem for shallow equation system consisting of water dynamics equations, silt transport equation, the equation of bottom topography change, and of some boundary and initial conditions is studied, the existence of its generalized solution and semidiscrete mixed finite element (MFE) solution was discussed, and the error estimates of the semidiscrete MFE solution was derived. The error estimates are optimal.     

THE MULTI-SYMPLECTIC ALGORITHM FOR “GOOD” BOUSSINESQ EQUATION

The multi-symplectic formulations of the “Good” Boussinesq equation were considered. For the multi-symplectic formulation, a new fifteen-point difference scheme which is equivalent to the multi-symplectic Preissman integrator was derived. The numerical experiments show that the multi-symplectic scheme have excellent long-time numerical behavior. Foundation items: the Foundation for Key Laboratory of Scientific/Engineering Computing Institute of Computational Mathematics and Scientific/Engineering Computing, Chinese Academy of Sciences; the Natural Science Foundation of Huaqiao University. Biography: ZENG Wen-ping (1940-), Professor (E-mail: qmz@1sec.cc.ac.cn)     

A Simple Fast Method in Finding Particular Solutions of Some Nonlinear PDE

1　ＡＴｒｉａｌＦｕｎｃｔｉｏｎａｎｄａＲｏｕｔｉｎｅｔｏＦｉｎｄＡｎａｌｙｔｉｃａｌＳｏｌｕｔｉｏｎｏｆＴｗｏＴｙｐｅｓｏｆＮｏｎｌｉｎｅａｒＰＤＥ　　Ｗｅｔｒｅａｔｔｈｅｎｏｎｌｉｎｅａｒｅｖｏｌｕｔｉｏｎｅｑｕａｔｉｏｎ ,ｗｈｉｃｈｉｓｆｏｒｍｅｄｂｙａｄｄｉｎｇｈｉｇｈｏｒｄｅｒｄｅｒｉｖａｔｉｖｅｔｅｒｍｓａｎｄｎｏｎｌｉｎｅａｒｔｅｒｍｓｔｏｔｈｅＢｕｒｇｅｒｓｅｑｕａｔｉｏｎ ｕ ｔ ｕ ｕ ｘ … ｕｐ ｕ ｘｑ α1 ｕ ｘ … αｎ ｎｕ ｘｎ =0 ,( 1)ｗｈｉｃｈｐ ,ｑ ,ｎａｎｄαｉ(ｉ =1,2…     

Numerical investigation on the flowfield of ＂swallowtail＂ cavity for supersonic mixing enhancement

A ＂swallowtail＂ cavity for the supersonic combustor was proposed to serve as an efficient flame holder for scramjets by enhancing the mass exchange between the cavity and the main flow. A numerical study on the ＂swallow- tail＂ cavity was conducted by solving the three-dimensional Reynolds-averaged Navier-Stokes equations implemented with a k-e turbulence model in a multi-block mesh. Turbu- lence model and numerical algorithms were validated first, and then test cases were calculated to investigate into the mechanism of cavity flows. Numerical results demonstrated that the certain mass in the supersonic main flow was sucked into the cavity and moved spirally toward the combustor walls. After that, the flow went out of the cavity at its lateral end, and finally was efficiently mixed with the main flow. The comparison between the ＂swallowtail＂ cavity and the conventional one showed that the mass exchanged between the cavity and the main flow was enhanced by the lateral flow that was induced due to the pressure gradient inside the cavity and was driven by the three-dimensional vortex ring generated from the ＂swallowtail＂ cavity structure.     

Mechanical behavior of amorphous polymers in shear

Based on the non-equilibrium thermodynamic theory, a new thermo-viscoelastic constitutive model for an incompressible material is proposed. This model can be considered as a kind of generalization of the non-Gaussian network theory in rubber elasticity to include the viscous and the thermal effects. A set of second rank tensorial internal variables was introduced, and in order to adequately describe the evolution of these internal variables, a new expression of the Helmholtz free energy was suggested. The mechanical behavior of the thermo-viscoelastic material under simple shear deformation was studied, and the “ viscous dissipation induced“ anisotropy due to the change of orientation distribution of molecular chains was examined. Influences of strain rate and thermal softening produced by the viscous dissipation on the shear stress were also discussed. Finally, the model predictions were compared with the experimental results performed by G‘ Sell et al. , thus the validity of the proposed model is verified.     

IMPACT MODEL RESOLUTION ON PAINLEVE’S PARADOX

Painlevé‘s paradox is one of the basic difficulties for solving LCP of dynamic systems subjected to unilateral constraints. A bi-nonlinear parameterized impact model, consistent with dynamic principles and experimental results, is established on the localized and quasi-static impact model theory. Numerical simulations are carried out on the dynamic motion of Painlevé‘s example. The results confirm ““““impact without collision““““ in the inconsistent states of the system. A ““““critical normal force““““ which brings an important effect on the future movement of the system in the indeterminate states is found. After the motion pattern for the impact process is obtained from numerical results, a rule of the velocity‘s jump that incorporates the tangential impact process is deduced by using an approximate impulse theory and the coefficient of restitution defined by Stronge. The results of the jump rule are quite precise if the system rigidity is big enough.     

Microgravity experiments of two-phase flow patterns aboard MIR space station

A first experimental study on two-phase flow patterns at a long-term, steady microgravity condition was conducted on board the Russian Space Station “MIR” in August 1999. Carbogal and air are used as the liquid and the gas phase, respectively. Bubble, slug, slug-annular transitional, and annular flows are observed. A new region of annular flow with lower liquid superficial velocity is discovered, and the region of the slug-annular transitional flow is wider than that observed by experiments on board the parabolic aircraft. The main patterns are bubble, slug-annular transitional and annular flows based on the experiments on board MIR space station. Some influences on the two-phase flow patterns in the present experiments are discussed. The project supported by the National Natural Science Foundation of China (19789201), the Ministry of Science and Technology of China (95-Yu-34), and the Post-doctoral Science Foundation of China     

Experimental investigation of nanomechanics of the mineral-protein interface in nacre

Nacre, the inner layer of molluscan seashells is a model biomimetic nanocomposite system that has been an inspiration for design of novel composites. The organic phase present in nacre is believed to play an important role in enhancing the toughness of nacre. The understanding of the adhesion forces of the organic matrix on to the mineral phase is essential for the fundamental assessment of the toughening properties in nacre. Our prior work using molecular dynamics simulations revealed that a very large force is needed to pull the protein molecules when in close proximity with the aragonite phase. In the current work, we have experimentally described the mechanical response of the organic phase in proximity of aragonite using force mode atomic force microscopy. Our results indicate that a very large force is required (>5–6 nN) to pull the proteins away from the aragonite. Our experiments show that the molecular interactions at the organic–inorganic interface in nacre are substantial and may play a significant role on the overall toughness of nacre. Thus, molecular interactions albeit weak and non-bonded play a significant role on the mechanics of hybrid nanocomposite systems.