The problem of the non-linear unsymmetrical bending for cylindrically orthotropic circular thin plate with variable thickness

I-Intr0ducti0nThebendingproblemofcircularthinplateisimriartantbothintheoryandinengineeringpractice,butbecausethebasicequationsoftheproblemsofnon-linearunsymmetricalbendingf0rcircularthinplate'arenonlineardifferentialequations,itisquitedifficulttoobtainana…     

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-)     

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.     

Computation of super-convergent nodal stresses of timoshenko beam elements by EEP method

The newly proposed element energy projection (EEP) method has been applied to the computation of super-convergent nodal stresses of Timoshenko beam elements. Generalformul as based on element projection theorem were derived and illustrative numerical examples using two typical elements were given. Both the analysis and examples show that EEP method also works very well for the problems with vector function solutions. The EEP method gives super-convergent nodal stresses, which are well comparable to the nodal displacements in terms of both convergence rate and error magnitude. And in addition, it can overcome the “ shear locking“ difficulty for stresses even when the displacements are badly affected. This research paves the way for application of the EEP method to general onedimensional systems of ordinary differential equations.     

Lower bound limit analysis of three-dimensional elastoplastic structures by boundary element method

IntroductionThelimitanalysisofstructuresisoneofthemostpracticalandusefulbranchesinplasticity .Ithasimportantapplicationbackgroundforproblemssuchasthedeterminationofloadcarryingcapacityandplasticformingofmetal.Thepurposeofthelimitanalysisofstructuresistoprovidereliabletheoreticalbasesforengineeringdesignandsafetyassessment.Asasimplifiedmethodforelastoplasticproblems,limitanalysisneednotrequirethehistoryofloadandcancomputethelimitloadsdirectlyinsteadofelastoplasticincrementalcomputationwhichisus…     

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

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

Periodic solutions of a class of higher order neutral type equations

In this paper, by ssing the theory of Fourier series, some necessary and sufficient conditions of existence and uniqueness of periodic solutions of a class of higher order neutral type equations are obtained. The main results by Shi Jianguo in “Discussion on the periodic solutions for linear equation of neutral type with constant coefficients” are improved, i. e., the condition |b₀|≠1 instead of the condtion |b₀|<1/2 of Theorem 1 by Shi Jianguo is given. Other theorems by Shi are rebuilt and improved according to the new assumption. Foundation item: the Natural Science Foundation of Yunnan Province, China (97A012G)     

STREAMLINE-BASED MATHEMATICAL MODEL FOR CO2 MISCIBLE FLOODING

According to the research theory of improved black oil simulator, a practical mathematical model for C02 miscible flooding was presented. In the model, the miscible process simulation was realized by adjusting oil/gas relative permeability and effective viscosity under the condition of miscible flow. In order to predict the production performance fast, streamline method is employed to solve this model as an alternative to traditional finite difference methods. Based on streamline distribution of steady-state flow through porous media with complex boundary confirmed with the boundary element method (BEM), an explicit total variation diminishing (TVD) method is used to solve the one-dimensional flow problem. At the same time, influences of development scheme, solvent slug size, and injection periods on CO2 drive recovery are discussed. The model has the advantages of less information need, fast calculation, and adaptation to calculate CO2 drive performance of all kinds of patterns in a random shaped porous media with assembly boundary. It can be an effective tool for early stage screening andmiscible oil field.reservoir dynamic management of the CO2 miscible oil field.     

Stress analysis for an infinite strip weakned by periodic cracks

Stress analysis for an infinite stripcracks were assumed in a horizontal position,weakened by periodic cracks is studied. The and the strip was applied by tension “p“ in y-direction. The boundary value problem can be reduced into a complex mixed one. It is found that the EEVM ( eigenfunction expansion variational method) is efficient to solve the problem. The stress intensity factor at the crack tip and the T-stress were evaluated. From the deformation response under tension the cracked strip can be equivalent to an orthotropic strip without cracks. The elastic properties in the equivalent orthotropic strip were also investigated. Finally, numerical examples and results were given.     

Initial layer phenomena for a class of singular perturbed nonlinear system with slow variables

The initial layer phenomena for a class of singular perturbed nonlinear system with slow variables are studied. By introducing stretchy variables with different quantity levels and constructing the correction term of initial layer with different “ thickness“, the Norder approximate expansion of perturbed solution concerning small parameter is obtained, and the “ multiple layer“ phenomena of perturbed solutions are revealed. Using the fixed point theorem, the existence of perturbed solution is proved, and the uniformly valid asymptotic expansion of the solutions is given as well.     

K-ε-T model of dense liquid-solid two-phase turbulent flow and its application to the pipe flow

To predict the characteristics of dense liquid-solid two-phase flow, K-ε-T model is established, in which the turbulent flow of fluid phase is described with fluid turbulent kinetic energy K_f and its dissipation rate ε_f, and the particles random motion is described with particle turbulent energy K_p and its dissipation rate ε_p and pseudothermal temperature T_p. The governing equations are also derived. With K-ε-T model, numerical study of dense liquid-solid two-phase turbulent up-flow in a pipe is performed. The calculated results are in good agreement with experimental data of Alajbegovic et al. (1994), and some flow features are captured.     

The nonlinear nonlocal singularly perturbed problems for reaction diffusion equations

A class of nonlinear nonlocal for singularly perturbed Robin initial boundary value problems for reaction diffusion equations is considered. Under suitable conditions, firstly, the outer solution of the original problem is obtained, secondly, using the stretched variable, the composing expansion method and the expanding theory of power series the initial layer is constructed, finally, using the theory of differential inequalities the asymptotic behavior of solution for the initial boundary value problems are studied and educing some relational inequalities the existence and uniqueness of solution for the original problem and the uniformly valid asymptotic estimation is discussed. Foundation items: the National Natural Science Foundation of China (10071048); the “Hunfred Talents Project” by Chinese Academy of Sciences Biography: Mo Jia-qi (1937−)     

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)     

THE PROBLEMS OF THE NONLINEAR UNSYMMETRICAL BENDING FOR CYLINDRICALLY ORTHOTROPIC CIRCULAR PLATE（II）

ＴＨＥＰＲＯＢＬＥＭＳＯＦＴＨＥＮＯＮＬＩＮＥＡＲＵＮＳＹＭＭＥＴＲＩＣＡＬ．ＢＥＮＤＩＮＧＦＯＲＣＹＬＩＮＤＲＩＣＡＬＬＹＯＲＴＨＯＴＲＯＰＩＣＣＩＲＣＵＬＡＲＰＬＡＴＥ（ＩＩ）ＨｕａｎｇＪｉａｙｉｎ（黄家寅）;ＱｉｎＳｈｅｎｇｌｉ（秦圣立）;Ｘｉ...     

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.     

Dynamic flight stability of hovering insects

The equations of motion of an insect with flapping wings are derived and then simplified to that of a flying body using the “rigid body” assumption. On the basis of the simplified equations of motion, the longitudinal dynamic flight stability of four insects (hoverfly, cranefly, dronefly and hawkmoth) in hovering flight is studied (the mass of the insects ranging from 11 to 1,648 mg and wingbeat frequency from 26 to 157 Hz). The method of computational fluid dynamics is used to compute the aerodynamic derivatives and the techniques of eigenvalue and eigenvector analysis are used to solve the equations of motion. The validity of the “rigid body” assumption is tested and how differences in size and wing kinematics influence the applicability of the “rigid body” assumption is investigated. The primary findings are: (1) For insects considered in the present study and those with relatively high wingbeat frequency (hoverfly, drone fly and bumblebee), the “rigid body” assumption is reasonable, and for those with relatively low wingbeat frequency (cranefly and howkmoth), the applicability of the “rigid body” assumption is questionable. (2) The same three natural modes of motion as those reported recently for a bumblebee are identified, i.e., one unstable oscillatory mode, one stable fast subsidence mode and one stable slow subsidence mode. (3) Approximate analytical expressions of the eigenvalues, which give physical insight into the genesis of the natural modes of motion, are derived. The expressions identify the speed derivative M _u (pitching moment produced by unit horizontal speed) as the primary source of the unstable oscillatory mode and the stable fast subsidence mode and Z _w (vertical force produced by unit vertical speed) as the primary source of the stable slow subsidence mode. The project supported by the National Natural Science Foundation of China (10232010 and 10472008).     

Difference scheme for two-phase flow

A numerical method for two-phase flow with hydrodynamics behavior was considered. The nonconservative hyperbolic governing equations proposed by Saurel and Gallout were adopted. Dissipative effects were neglected but they could be included in the model without major difficulties. Based on the opinion proposed by Abgrall that “a two phase system, uniform in velocity and pressure at t = 0 will be uniform on the same variable during its temporal evolution“, a simple accurate and fully Eulerian numerical method was presented for the simulation of multiphase compressible flows in hydrodynamic regime. The numerical method relies on Godunov-typescheme, with HLLC and Lax-Friedrichs type approximate Riemann solvers for the resolution of conservation equations, and nonconservative equation. Speed relaxation and pressure relaxation processes were introduced to account for the interaction between the phases. Test problem was presented in one space dimension which illustrated that our scheme is accurate, stable and oscillation free.     

THE PROBLEMS OF THE NONLINEAR UNSYMMETRICAL BENDING FOR CYLINDRICALLY ORTHOTROPIC CIRCULA RPLATE（I）

ＴＨＥＰＲＯＢＬＥＭＳＯＦＴＨＥＮＯＮＬＩＮＥＡＲＵＮＳＹＭＭＥＴＲＩＣＡＬＢＥＮＤＩＮＧＦＯＲＣＹＬＩＮＤＲＩＣＡＬＬＹＯＲＴＨＯＴＲＯＰＩＣＣＩＲＣＵＬＡＲＰＬＡＴＥ（Ｉ）ＱｉｎＳｈｅｎｇ－Ｉｉ（秦圣立）ＨｕａｎｇＪｉａ－ｙｉｎ（黄家寅）（Ｑｕｆ...     

The problems of nonlinear bending for orthotropic rectangular plate with four clamped edges

（黄家寅）（秦圣立）ＴＨＥＰＲＯＢＬＥＭＳＯＦＮＯＮＬＩＮＥＡＲＢＥＮＤＩＮＧＦＯＲＯＲＴＨＯＴＲＯＰＩＣＲＥＣＴＡＮＧＵＬＡＲＰＬＡＴＥＷＩＴＨＦＯＵＲＣＬＡＭＰＥＤＥＤＧＥＳ￥ＨｕａｎｇＪｉａｙｉｎ;ＱｉｎＳｈｅｎｇｌｉ（ＱｕｆｕＮｏｒｍａｌＵｎ...     

Modeling the interaction of solitary waves and semi-circular breakwaters by using unsteady reynolds equations

A vertical 2 -D numerical wave model was developed based on unsteady Reynolds equations. In this model, the k-epsilon models were used to close the Reynolds equations, and volume of fluid( VOF) method was used to reconstruct the free surface. The model was verified by experimental data. Then the model was used to simulate solitary wave interaction with submerged, alternative submerged and emerged semi-circular breakwaters. The process of velocity field, pressure field and the wave surface near the breakwaters was obtained. It is found that when the semi-circular breakwater is submerged, a large vortex will be generated at the bottom of the lee side wall of the breakwater ; when the still water depth is equal to the radius of the semi-circular breakwater, a pair of large vortices will be generated near the shoreward wall of the semi-circular breakwater due to wave impacting, but the velocity near the bottom of the lee side wall of the breakwater is always relatively small. When the semi-circular breakwater is emerged, and solitary wave cannot overtop it, the solitary wave surface will run up and down secondarily during reflecting from the breakwater. It can be further used to estate the diffusing and transportation of the contamination and transportation of suspended sediment.