Table of Content

18 April 2006, Volume 27 Issue 4

Previous Issue    Next Issue

Articles

STUDY OF MECHANISM OF BREAKDOWN IN LAMINAR-TURBULENT TRANSITION
OF SUPERSONIC BOUNDARY LAYER ON FLAT PLATE

CAO Wei;HUANG Zhang-feng;ZHOU Heng

2006, 27(4):  425-434 .  doi:10.1007/s10483-006-0401-1

Abstract ( 1515 )   PDF (407KB) ( 416 )  

Related Articles | Metrics

Spatial mode direct numerical simulation has been applied to study the mechanism of breakdown in laminar-turbulent transition of a supersonic boundary layer on a flat plate with Mach number 4.5. Analysis of the result showed that, during the breakdown process in laminar-turbulent transition, the mechanism causing the mean flow profile to evolve swiftly from
laminar to turbulent was that the modification of mean flow profile by the disturbance, when they became larger, leads to remarkable change of its stability characteristics. Though the most unstable T-S wave was of second mode for laminar flow, the first mode waves played the key role in the breakdown process in laminar-turbulent transition.

EFFECTIVE SOLUTION METHOD OF CHEMICAL REACTION KINETICS WITH DIFFUSE

LÜHe-xiang;QIU Kun-yu;CHEN Jian-feng

2006, 27(4):  435-442 .  doi:10.1007/s10483-006-0402-z

Abstract ( 1514 )   PDF (244KB) ( 521 )  

Related Articles | Metrics

The time integration method with four-order accuracy, self-starting and implicit for the diffuse chemical reaction kinetics equation or the transient instantaneous temperature filed equation was presented. The examples show that both accuracy and stability are better than Runge-Kutta method with four-order. The coefficients of the equation are stored with sparse
matrix pattern, so an algorithm is presented which combines a compact storage scheme with reduced computation cost. The computation of the competitive and consecutive reaction in the rotating packed bed, taken as examples, shows that the method is effective.

DYNAMICAL FORMATION OF CAVITY FOR COMPOSED THERMAL HYPERELASTIC SPHERES IN NON-UNIFORM TEMPERATURE FIELDS

CHENG Chang-jun;MEI Bo

2006, 27(4):  443-452 .  doi:10.1007/s10483-006-0403-z

Abstract ( 1471 )   PDF (232KB) ( 507 )  

Related Articles | Metrics

Dynamical formation and growth of cavity in a sphere composed of two incompressible thermal-hyperelastic Gent-Thomas materials were discussed under the case of a non-uniform temperature field and the surface dead
loading. The mathematical model was first presented based on the dynamical theory of finite deformations. An exact differential relation between the void radius and surface load was obtained by using the variable transformation method. By numerical computation, critical loads and cavitation growth curves were obtained for different temperatures. The influence of the temperature and material parameters of the composed sphere on the void formation and
growth was considered and compared with those for static analysis. The results show that the cavity occurs suddenly with a finite radius and its evolvement with time displays a non-linear periodic vibration and that the critical load decreases with the increase of temperature and also the dynamical critical load is lower than the static critical load under the same conditions.

BIPARAMETRIC PERTURBATION SOLUTIONS OF LARGE DEFLECTION PROBLEM OF CANTILEVER BEAMS

HE Xiao-ting;CHEN Shan-lin

2006, 27(4):  453-460 .  doi:10.1007/s10483-006-0404-z

Abstract ( 1568 )   PDF (142KB) ( 486 )  

Related Articles | Metrics

The large deflection problem of cantilever beams was studied by means of the biparametric perturbation method and the first order derivative substitution from pseudolinear analysis approach. This kind of substitution can transform the basic equation, an integral differential equation into nonlinear algebraic ones, thus simplify computational process. Compared with present results, it indicates that the large deflection problem solved by using pseudolinear analysis can lead to simple and precise results.

FREE VIBRATION OF ANISOTROPIC RECTANGULAR PLATES BY GENERAL ANALYTICAL METHOD

HUANG Yan;LEI Yong-jun;SHEN Hui-jun

2006, 27(4):  461-467 .  doi:10.1007/s10483-006-0405-y

Abstract ( 1547 )   PDF (152KB) ( 564 )  

Related Articles | Metrics

According to the differential equation for transverse displacement function of anisotropic rectangular thin plates in free vibration, a general analytical solution is established. This general solution, composed of the composite solutions of trigonometric function and hyperbolic function, can satisfy the problem of arbitrary boundary conditions along four edges. The algebraic polynomial with double sine series solutions can also satisfy the problem of boundary conditions at four corners. Consequently, this general solution can be used to solve the vibration problem of anisotropic rectangular plates with arbitrary boundaries accurately. The integral constants can be determined by boundary conditions of four edges and four corners. Each natural frequency and vibration mode can be solved by the determinate of coefficient matrix from the homogeneous linear algebraic equations equal to zero. For example, a composite symmetric angle ply laminated plate with four edges clamped has been calculated and discussed.

FINITE ELEMENT METHOD ON NUMERICAL SIMULATION OF STRATUM CORNEUM'S PENETRATION PROPERTY

LIU Yu-hong;QIAO Ai-ke;Dirk Feuchter;Gabriel Wittum;ZENG Yan-jun

2006, 27(4):  469-475 .  doi:10.1007/s10483-006-0406-y

Abstract ( 1521 )   PDF (274KB) ( 489 )  

Related Articles | Metrics

How the outer substance could penetrate through the skin lies in the stratum corneum, because it is the main barrier in the multi-layers of the skin. Supposing the keratin cell with a special geometry as tetrakaidecahedron, the penetration property of stratum corneum was the key problem which was numerically simulated with finite element method. At first the discretization of the stratum corneum region was given in two steps: first, the discretization of the keratin cell; second, the discretization of fattiness that surrounds the keratin. Then there was the work of numerical simulation. In this procedure, the finite element method and the multi-grid method were used. The former was to obtain the discretization of basic elements; the latter was to decrease the high frequency error. At last the visualization of the numerical simulation was shown.

NEW METHOD FOR LOW ORDER SPECTRAL MODEL AND ITS APPLICATION

CAO Jie;YOU Ya-lei

2006, 27(4):  477-484 .  doi:10.1007/s10483-006-0407-z

Abstract ( 1153 )   PDF (162KB) ( 451 )  

Related Articles | Metrics

In order to overcome the deficiency in classical method of low order spectral model, a new method for low order spectral model was advanced. Through
calculating the multiple correlation coefficients between combinations of different functions and the recorded data under the least square criterion, the truncated functions which can mostly reflect the studied physical phenomenon were objectively distilled from these data. The new method overcomes the deficiency of artificially selecting the truncated functions in the classical low order spectral model. The new method being applied to study the inter-annual variation of summer atmospheric circulation over
Northern Hemisphere, the truncated functions were obtained with the atmospheric circulation data of June 1994 and June 1998. The mechanisms for the two-summer atmospheric circulation variations over Northern Hemisphere were obtained with two-layer quasi-geostrophic baroclinic equation.

NORMAL EXPANSION THEORY FOR PENETRATION OF PROJECTILE AGAINST CONCRETE TARGET

GAO Shi-qiao;LIU Hai-peng;LI Ke-jie;HUANG Feng-lei;JIN Lei

2006, 27(4):  485-492 .  doi:10.1007/s10483-006-0408-y

Abstract ( 1396 )   PDF (241KB) ( 535 )  

Related Articles | Metrics

Based on the equations which describe the dynamic behavior of material under high-velocity and high-pressure shock, corresponding equations at shock front whose surface is general space curve surface were established. For concrete material, a normal expansion theory was proposed by which
some deceleration about time history of the projectile can be analytically given. This normal expansion theory is not only suitable for spherical and cylindrical-nose projectile, but also suitable for other general nose projectile, for example conical nose or ogive-nose. And it is not only suitable for perpendicular shock but also for oblique shock.

CONTINUOUS SELECTION THEOREMS FOR FAN-BROWDER MAPPINGS IN TOPOLOGICAL SPACES AND THEIR APPLICATIONS

YANG Ming-ge;DENG Lei

2006, 27(4):  493-500 .  doi:10.1007/s10483-006-0409-z

Abstract ( 1443 )   PDF (160KB) ( 419 )  

Related Articles | Metrics

The concept of Fan-Browder mappings was first introduced in topological spaces without any convex structure. Then a new continuous selection theorem was obtained for the Fan-Browder mapping with range in a topological space without any convex structure and noncompact domain. As applications, some fixed point theorems, coincidence theorems and a nonempty intersection theorem were given. Both the new concepts and results unify and extend many known results in recent literature.

ADAPTIVE REGULATION OF HIGH ORDER NONHOLONOMIC SYSTEMS

MU Xiao-wu;YU Ji-min;CHENG Gui-fang

2006, 27(4):  501-507 .  doi:10.1007/s10483-006-0410-z

Abstract ( 1454 )   PDF (158KB) ( 594 )  

Related Articles | Metrics

The problem of adaptive regulation of a class of high-order parametric nonholonomic systems in chained-form was discussed. Using adding a power integrator technique and state scaling with discontinuous projection technique, a discontinuous adaptive dynamic controller was constructed. The controller guarantees the estimated value of unknown parameter is in the prescribed extent.

CRACK PROPAGATION IN POLYCRYSTALLINE ELASTIC-VISCOPLASTIC MATERIALS USING COHESIVE ZONE MODELS

WU Yan-qing;ZHANG Ke-shi

2006, 27(4):  509-518 .  doi:10.1007/s10483-006-0411-y

Abstract ( 1371 )   PDF (388KB) ( 425 )  

Related Articles | Metrics

Cohesive zone model was used to simulate two-dimensional plane strain crack propagation at the grain level model including grain boundary zones. Simulated results show that the original crack-tip may not be separated firstly in an elastic-viscoplastic polycrystals. The grain interior's material properties (e.g. strain rate sensitivity) characterize the competitions between plastic and cohesive energy dissipation mechanisms. The higher the strain rate sensitivity is, the larger amount of the external work is transformed into plastic dissipation energy than into cohesive energy, which delays the cohesive zone rupturing. With the strain rate sensitivity decreased, the material property tends to approach the elastic-plastic responses. In this case, the plastic dissipation energy decreases and the cohesive dissipation energy increases which accelerates the cohesive zones debonding. Increasing the cohesive strength or the critical separation displacement will reduce the stress triaxiality at grain interiors and grain boundaries. Enhancing the cohesive zones ductility can improve the matrix materials resistance to void damage.

A DOMAIN DECOMPOSITION ALGORITHM WITH FINITE ELEMENT-BOUNDARY ELEMENT COUPLING

YAN Bo;DU Juan;HU Ning;SEKINE Hideki

2006, 27(4):  519-525 .  doi:10.1007/s10483-006-0412-y

Abstract ( 1511 )   PDF (212KB) ( 666 )  

Related Articles | Metrics

A domain decomposition algorithm coupling the finite element and the boundary element was presented. It essentially involves subdivision of the analyzed domain into sub-regions being independently modeled by two methods, i.e., the finite element method (FEM) and the boundary
element method (BEM). The original problem was restored with continuity and equilibrium conditions being satisfied on the interface of the two sub-regions using an iterative algorithm. To speed up the convergence rate of the iterative algorithm, a dynamically changing relaxation parameter during iteration was introduced. An advantage of the proposed algorithm is that the
locations of the nodes on the interface of the two sub-domains can be inconsistent. The validity of the algorithm is demonstrated by the consistence of the results of a numerical example obtained by the proposed method and those by the FEM, the BEM and a present finite element-boundary element (FE-BE) coupling method.

VISCOPLASTIC SOLUTION TO FIELD AT STEADILY PROPAGATING CRACK TIP IN LINEAR-HARDENING MATERIALS

JIA Bin;WANG Zhen-qing;LI Yong-dong;LIANG Wen-yan

2006, 27(4):  527-533 .  doi:10.1007/s10483-006-0413-1

Abstract ( 1371 )   PDF (197KB) ( 573 )  

Related Articles | Metrics

An elastic-viscoplastic constitutive model was adopted to analyze asymptotically the tip-field of moving crack in linear-hardening materials under plane strain condition. Under the assumption that the artificial viscosity coefficient was in inverse proportion to power law of the rate of effective plastic strain, it is obtained that stress and strain both possess power law singularity and the singularity exponent is uniquely determined by the power law exponent of the rate of effective plastic strain. Variations of zoning structure according to each material parameter were discussed by means of numerical computation for the tip-field of mode II dynamic propagating crack, which show that the structure of crack tip field is dominated by hardening coefficient rather than viscosity coefficient. The secondary plastic zone can be ignored for weak hardening materials while the secondary plastic zone and the secondary elastic zone both have important influence on crack tip field for strong hardening materials. The dynamic solution approaches to the corresponding quasi-static solution when the crack moving speed goes to zero, and further approaches to the HR (Hui-Riedel) solution when the hardening coefficient is equal to zero.

NUMERICAL STUDY OF PARTICLE DISTRIBUTION IN WAKE OF LIQUID-PARTICLE FLOWS PAST A CIRCULAR CYLINDER USING DISCRETE VORTEX METHOD

HUANG Yuan-dong;WU Wen-quan

2006, 27(4):  535-542 .  doi:10.1007/s10483-006-0414-1

Abstract ( 1483 )   PDF (298KB) ( 506 )  

Related Articles | Metrics

Particle-laden water flows past a circular cylinder were numerically investigated. The discrete vortex method (DVM) was employed to evaluate the unsteady water flow fields and a Lagrangian approach was applied for tracking individual solid particles. A dispersion function was defined to represent the dispersion scale of the particle. The wake vortex patterns, the distributions and the time series of dispersion functions of particles with different Stokes numbers were obtained. Numerical results show that the particle distribution in the wake of the circular cylinder is closely related to the particle's Stokes number and the structure of wake vortices: (1) the intermediate sized particles with Stokes numbers, St, of 0.25, 1.0 and 4.0 can not enter the vortex cores and concentrate near the peripheries of the vortex structures, (2) in the circular cylinder wake, the dispersion intensity of particles decreases as St is increased from 0.25 to 4.0.

RECURRENT NEURAL NETWORK MODEL BASED ON PROJECTIVE OPERATOR AND ITS APPLICATION TO OPTIMIZATION PROBLEMS

MA Ru-ning;CHEN Tian-ping

2006, 27(4):  543-554 .  doi:10.1007/s10483-006-0415-z

Abstract ( 1586 )   PDF (281KB) ( 511 )  

Related Articles | Metrics

The recurrent neural network (RNN) model based on projective operator was studied. Different from the former study, the value region of projective operator in the neural network in this paper is a general closed convex subset of n-dimensional Euclidean space and it is not a compact convex set
in general, that is, the value region of projective operator is probably unbounded. It was proved that the network has a global solution and its solution trajectory converges to some equilibrium set whenever objective function satisfies some conditions. After that, the model was applied to continuously differentiable optimization and nonlinear or implicit complementarity problems. In addition, simulation experiments confirm the efficiency of the RNN.

DYNAMIC BEHAVIOR OF THIN RECTANGULAR PLATE ATTACHED TO MOVING RIGID

XIAO Shi-fu;CHEN Bin

2006, 27(4):  555-566 .  doi:10.1007/s10483-006-0416-1

Abstract ( 1480 )   PDF (187KB) ( 503 )  

Related Articles | Metrics

A nonlinear dynamic model of a thin rectangular plate attached to a moving rigid was established by employing the general Hamilton's variational principle. Based on the new model, it is proved theoretically that both phenomena of dynamic stiffening and dynamic softening can occur in the plate when the rigid undergoes different large overall motions including overall translational and rotary motions. It was also proved that dynamic
softening effect even can make the trivial equilibrium of the plate lose its stability through bifurcation. Assumed modes method was employed to validate the theoretical result and analyze the approximately critical bifurcation value and the post-buckling equilibria.