Stokes＇ first problem for a viscoelastic fluid with the generalized Oldroyd-B model

The flow near a wall suddenly set in motion for a viscoelastic fluid with the generalized Oldroyd-B model is studied. The fractional calculus approach is used in the constitutive relationship of fluid model. Exact analytical solutions of velocity and stress are obtained by using the discrete Laplace transform of the sequential fractional derivative and the Fox H-function. The obtained results indicate that some well known solutions for the Newtonian fluid, the generalized second grade fluid as well as the ordinary Oldroyd-B fluid, as limiting cases, are included in our solutions. The project supported by the National Natural Science Foundation of China (10272067), the Doctoral Program Foundation of the Education Ministry of China (20030422046), the Natural Science Foundation of Shandong Province, China (Y2006A14) and the Research Foundation of Shandong University at Weihai. The English text was polished by Keren Wang.     

Plane surface suddenly set in motion in a viscoelastic fluid with fractional Maxwell model

The fractional calculus approach in the constitutive relationship model of viscoelastic fluid is introduced. The flow near a wall suddenly set in motion is studied for a non-Newtonian viscoelastic fluid with the fractional Maxwell model. Exact solutions of velocity and stress are obtained by using the discrete inverse Laplace transform of the sequential fractional derivatives. It is found that the effect of the fractional orders in the constitutive relationship on the flow field is significant. The results show that for small times there are appreciable viscoelastic effects on the shear stress at the plate, for large times the viscoelastic effects become weak. The project supported by the National Natural Science Foundation of China (10002003), Foundation for University Key Teacher by the Ministry of Education, Research Fund for the Doctoral Program of Higher Education     

Unsteady flows of a generalized second grade fluid with the fractional derivative model between two parallel plates

The fractional calculus approach in the constitutive relationship model of a generalized second grade fluid is introduced. Exact analytical solutions are obtained for a class of unsteady flows for the generalized second grade fluid with the fractional derivative model between two parallel plates by using the Laplace transform and Fourier transform for fractional calculus. The unsteady flows are generated by the impulsive motion or periodic oscillation of one of the plates. In addition, the solutions of the shear stresses at the plates are also determined. The project supported by the National Natural Science Foundation of China (10372007, 10002003) and CNPC Innovation Fund     

An exact solution for flow past an accelerated horizontal plate in a rotating fluid with the generalized Oldroyd-B model

In this paper, the generalized Oldroyd-B with fractional calculus approach is used. An exact solution in terms of Fox-H function for flow past an accelerated horizontal plate in a rotating fluid is obtained by using discrete Laplace transform method. A comparison among the influence of various parameters in the Oldroyd-B model and the angular velocity of the fluid on the velocity profiles is made through numerical method in graphic form.     

A mixed finite element scheme for viscoelastic flows with XPP model

A mixed finite element formulation for viscoelastic flows is derived in this paper, in which the FIC （finite incremental calculus） pressure stabilization process and the DEVSS （discrete elastic viscous stress splitting） method using the Crank-Nicolson-based split are introduced within a general framework of the iterative version of the fractional step algorithm. The SU （streamline-upwind） method is particularly chosen to tackle the convective terms in constitutive equations of viscoelastic flows. Thanks to the proposed scheme the finite elements with equal low-order interpolation approximations for stress-velocity-pressure variables can be successfully used even for viscoelastic flows with high Weissenberg numbers. The XPP （extended Pom-Pom） constitutive model for describing viscoelastic behaviors is particularly integrated into the proposed scheme. The numerical results for the 4：1 sudden contraction flow problem demonstrate prominent stability, accuracy and convergence rate of the proposed scheme in both pressure and stress distributions over the flow domain within a wide range of the Weissenberg number, particularly the capability in reproducing the results, which can be used to explain the ＂die swell＂ phenomenon observed in the polymer injection molding process.     

The transient fracture behavior for a functionally graded layered structure subjected to an in-plane impact load

The transient fracture behavior of a functionally graded layered structure subjected to an in-plane impact load is investigated. The studied structure is composed of two homogeneous layers and a functionally graded interlayer with a crack perpendicular to the boundaries. The impact load is applied on the face of the crack. Fourier transform and Laplace transform methods are used to formulate the present problem in terms of a singular integral equation in Laplace transform domain. Considering variations of parameters such as the nonhomogeneity constant, the thickness ratio and the crack length, the dynamic stress intensity factors (DSIFs) in time domain are studied and some meaningful conclusions are obtained.The project supported by the National Science Foundation for Excellent Young Investigators (10325208), the National Natural Science Foundation of China (10432030) and the China Postdoctoral Science Foundation (2004036018)The English text was polished by Ron Marshall.     

Free axisymmetric vibration of transversely isotropic laminated circular plates

Based on the three dimensional elastic theory, the state equation of the axisymmetric free vibration of transversely isotropic circular plates is established. Taking the advantage of finite Hankel transform, two exact solutions are derived for two boundary conditions, i. e. the rigid-slipping boundary and elastic simply supported boundary. Finally, numerical results are presented and compared with those of FEM. The project is supported by the National Natural Science Foundation of China and the Zhejiang Provincial Natural Science Foundation.     

Effect of microstructure on the hardening and softening behaviors of polycrystalline shape memory alloys Part II: Numerical simulation under axisymmetrical loading

Based on the microstructure-based constitutive model established in Part I, a detailed numerical investigation on the role of each microstructure parameter in the kinematical and kinetic evolution of polycrystalline SMA under axisymmetrical tension loading is performed. Some macroscopic constitutive features of stress-induced martensite transformation are discussed. The subject supported by the Research Grant Committee (RGC) of Hong Kong SAR, the National Natural Science Foundation of China and the Provincial Natural Science Foundation of Jiangxi Province of China     

Numerical study of an oscillatory turbulent flow over a flat plate

Oscillatory turbulent flow over a flat plate is studied using large eddy simulation (LES) and Reynolds-average Navier-Stokes (RANS) methods. A dynamic subgrid-scale model is employed in LES and Saffman's turbulence model is used in RANS. The flow behaviors are discussed for the accelerating and decelerating phases during the oscillating cycle. The friction force on the wall and its phase shift from laminar to turbulent regime are also investigated for different Reynolds numbers. The project supported by the Youngster Funding of Academia Sinica and by the National Natural Science Foundation of China     

广义Maxwell黏弹性流体在两平板间的非定常流动

将分数阶微积分运算引入Maxwell黏弹性流体的本构方程，研究了黏弹性流体在两平板问的非定常流动．对于广义Maxwell黏弹性流体的分数阶导数模型，导出了对时间具有分数阶导数的特殊运动方程，利用分数阶微积分的Laplace变换理论，得到了流动的解析解．     

Investigation of the dynamic behavior of two collinear anti-plane shear cracks in a piezoelectric layer bonded to two half spaces by a new method

The dynamic behavior of two collinear anti-plane shear cracks in a piezoelectric layer bonded to two half spaces subjected to the harmonic waves is investigated by a new method. The cracks are parallel to the interfaces in the mid-plane of the piezoelectric layer. By using the Fourier transform, the problem can be solved with two pairs of triple integral equations. These equations are solved by using Schmidt’s method. This process is quite different from that adopted previously. Numerical examples are provided to show the effect of the geometry of cracks, the frequency of the incident wave, the thickness of the piezoelectric layer and the constants of the materials upon the dynamic stress intensity factor of cracks.     

The behavior of a crack in functionally graded piezoelectric/piezomagnetic materials under anti-plane shear loading

Summary In this paper, the behavior of a crack in functionally graded piezoelectric/piezomagnetic materials subjected to an anti-plane shear loading is investigated. To make the analysis tractable, it is assumed that the material properties vary exponentially with the coordinate parallel to the crack. By using a Fourier transform, the problem can be solved with the help of a pair of dual integral equations in which the unknown variable is the jump of the displacements across the crack surfaces. These equations are solved using the Schmidt method. The relations among the electric displacement, the magnetic flux and the stress field near the crack tips are obtained. Numerical examples are provided to show the effect of the functionally graded parameter on the stress intensity factors of the crack.The authors are grateful for financial support from the Natural Science Foundation of Hei Long Jiang Province (A0301), the National Natural Science Foundation of China (50232030, 10172030), the Natural Science Foundation with Excellent Young Investigators of Hei Long Jiang Province(JC04-08) and the National Science Foundation with Excellent Young Investigators (10325208).     

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     

A finite element analysis of the flange earrings of strong anisotropic sheet metals in deep-drawing processes

Flange earrings of strong anisotropic sheet metals in deep-drawing process are numerically analyzed by the elastic-plastic large deformation finite element formulation based on a discrete Kirchhoff triangle plate shell element model. A Barlat-Lian anisotropic yield function and a quasi-flow corner theory are used in the present formulation. The numerical results are compared with the experimental ones of cylindrical cup drawing process. The focus of the present researches is on the numerical analysis and the constraining scheme of the flange earring of circular sheets with strong anisotropy in square cup drawing process. The project supported by the National Natural Science Foundation of China (19832020) and Provincial Natural Science Foundation of Jilin, China (200000519)     

Anti-plane crack moving along the interface of dissimilar piezoelectric materials

The problem of an anti-plane Griffith crack moving along the interface of dissimilar piezoelectric materials is solved by using the integral transform technique. It is shown from the result that the intensity factors of anti-plane stress and electric displacement around the crack tip are dependent on the speed of the Griffith crack as well as the material coefficients. When the two piezoelectric materials are identical, the present result will be reduced to the result for the problem of an anti-plane moving Griffith crack in homogeneous piezoelectric materials. Supported by the National Natural Science Foundation and the National Post-doctoral Science Foundation of China.     

The scattering of harmonic elastic anti-plane shear waves by two collinear cracks in the piezoelectric plate by using the non-local theory

In this paper, the dynamic interaction between two collinear cracks in a piezoelectric material plate under anti-plane shear waves is investigated by using the non-local theory for impermeable crack surface conditions. By using the Fourier transform, the problem can be solved with the help of two pairs of triple integral equations. These equations are solved using the Schmidt method. This method is more reasonable and more appropriate. Unlike the classical elasticity solution, it is found that no stress and electric displacement singularity is present at the crack tip. The non-local dynamic elastic solutions yield a finite hoop stress at the crack tip, thus allowing for a fracture criterion based on the maximum dynamic stress hypothesis. The project supported by the Natural Science Foundation of Heilongjiang Province and the National Natural Science Foundation of China(10172030, 50232030)     

Asymmetric vortices flow over slender body and its active control at high angle of attack

The studies of asymmetric vortices flow over slender body and its active control at high angles of attack have significant importance for both academic field and engineering area. This paper attempts to provide an update state of art to the investigations on the fields of forebody asymmetric vortices. This review emphasizes the correlation between micro-perturbation on the model nose and its response and evolution behaviors of the asymmetric vortices. The critical issues are discussed, which include the formation and evolution mechanism of asymmetric multi-vortices; main behaviors of asymmetric vortices flow including its deterministic feature and vortices flow structure; the evolution and development of asymmetric vortices under the perturbation on the model nose; forebody vortex active control especially discussed micro-perturbation active control concept and technique in more detail. However present understanding in this area is still very limited and this paper tries to identify the key unknown problems in the concluding remarks. The project supported by the National Natural Science Foundation of China (10172017), Aeronautical Science Foundation of China (02A51048) and Foundation of National Key Laboratory of Aerodynamic Design and Research (51462020504HK0101)     

Integrability for peaffian constrained systems: A geometrical theory

There exists an Ehresmann connection on the fibred constrained sub-manifold defined by Pfaffian differential constraints. It is proved that curvature of the connection is closely related to the d-σ commutation relation in the classical nonholonomic mechanics. It is also proved that conditions of complete integrability for Pfaffian systems in Frobenius sense are equivalent to the three requirements upon the conditional variations in the classical calculus of variations: (1) the variations belong to the constrained manifold, (2) variational operators commute with differential operators, (3) variations satisfy the Chetaev's conditions. Thus this theory verifies the conjecture or experience of researchers of mechanics on the integrability conditions in terms of variation calculus. The project supported by the National Natural Science Foundation of China     

INVESTIGATION OF BEHAVIOR OF MODE-I INTERFACE CRACK IN PIEZOELECTRIC MATERIALS BY USING SCHMIDT METHOD

The behavior of a Mode-Ⅰinterface crack in piezoelectric materials was investigated under the assumptions that the effect of the crack surface overlapping very near the crack tips was negligible. By use of the Fourier transform, the problem can be solved with the help of two pairs of dual integral equations. To solve the dual integral equations, the jumps of the displacements across the crack surfaces were expanded in a series of Jacobi polynomials. It is found that the stress and the electric displacement singularities of the present interface crack solution are the same as ones of the ordinary crack in homogenous materials. The solution of the present paper can be returned to the exact solution when the upper half plane material is the same as the lower half plane material.     

Poincaré-cartan integral invariants of Birkhoffian systems

Based on modern differential geometry, the symplectic structure of a Birkhoffian system which is an extension of conservative and nonconservative systems is analyzed. An integral invariant of Poincaré-Cartan’s type is constructed for Birkhoffian systems. Finally, one-dimensional damped vibration is taken as an illustrative example and an integral invariant of Poincaré’s type is found.