基于新修正偶应力理论的Mindlin层合板热稳定性分析

以新修正偶应力理论为基础,首次提出了机械载荷与热载荷共同作用下的微尺度Mindlin层合板热稳定性模型,该模型只引入一个材料尺度参数,通过虚功原理推导出了控制方程和边界条件,以四边简支方板为例,进行了热稳定性分析,应用纳维叶解法得到解析解。结果表明,所建模型可以捕捉到尺度效应。材料尺度参数值越大,屈曲临界温度越高;当跨厚比增大时,屈曲临界温度下降;随着板几何参数的增大,模型将退化为宏观模型;温度变化量越大,考虑热载荷作用下的屈曲临界载荷越大,尺度效应体现越显著。     

Extended Kantorovich method for local stresses in composite laminates upon polynomial stress functions

The extended Kantorovich method is employed to study the local stress concentrations at the vicinity of free edges in symmetrically layered composite laminates subjected to uniaxial tensile load upon polynomial stress functions. The stress fields are initially assumed by means of the Lekhnitskii stress functions under the plane strain state. Applying the principle of complementary virtual work, the coupled ordinary differential equations are obtained in which the solutions can be obtained by solving a generalized eigenvalue problem. Then an iterative procedure is estab-lished to achieve convergent stress distributions. It should be noted that the stress function based extended Kantorovich method can satisfy both the traction-free and free edge stress boundary conditions during the iterative processes. The stress components near the free edges and in the interior regions are calculated and compared with those obtained results by finite element method (FEM). The convergent stresses have good agreements with those results obtained by three dimensional (3D) FEM. For generality, various layup configurations are considered for the numerical analysis. The results show that the proposed polynomial stress function based extended Kan-torovich method is accurate and efficient in predicting the local stresses in composite laminates and computationally much more efficient than the 3D FEM.     

An eight noded composite plate element for the dynamic analysis of spatial mechanism systems

The development of a shear-deformable laminated plate element, based on the Mindlin plate theory, for use in large reference displacement analysis is presented. The element is sufficiently general to accept an arbitrary number of layers and an arbitrary number of orthotrophic material property sets. Coordinate mapping is utilized so that non-rectangular elements may be modeled. The Gauss quadrature method of numerical integration is utilized to evaluate volume integrals. A comparative study is done on the use of full Gauss quadrature, reduced Gauss quadrature, mixed Gauss quadrature, and closed form integration techniques for the element. Dynamic analysis is performed on the RSSR (Revolute-Spherical-Spherical-Revolute) mechanism, with the coupler modeled as a flexible plate. The results indicate the differences in the dynamic response of the transverse shear deformable eight-noded element as compared to a four-noded plate element. Dynamically induced stresses are examined, with the results indicating that the primary deformation mode of the eight-noded Mindlin plate model being bending.     

Fully computable a posteriori error bounds for stabilised FEM approximations of convection–reaction–diffusion problems in three dimensions

Fully computable upper bounds are developed for the discretisation error measured in the natural (energy) norm for convection–reaction–diffusion problems in three dimensions. The upper bounds are genuine upper bounds in the sense that the numerical value of the estimated error exceeds the actual numerical value of the true error regardless of the coarseness of the mesh or the nature of the data for the problem. All constants appearing in the bounds are fully specified. Examples show the estimator to be reliable and accurate even in the case of complicated three‐dimensional problems. Copyright © 2013 John Wiley & Sons, Ltd.     

Solution of the Navier–Stokes equations in velocity–vorticity form using a Eulerian–Lagrangian boundary element method

This paper describes the Eulerian–Lagrangian boundary element model for the solution of incompressible viscous flow problems using velocity–vorticity variables. A Eulerian–Lagrangian boundary element method (ELBEM) is proposed by the combination of the Eulerian–Lagrangian method and the boundary element method (BEM). ELBEM overcomes the limitation of the traditional BEM, which is incapable of dealing with the arbitrary velocity field in advection‐dominated flow problems. The present ELBEM model involves the solution of the vorticity transport equation for vorticity whose solenoidal vorticity components are obtained iteratively by solving velocity Poisson equations involving the velocity and vorticity components. The velocity Poisson equations are solved using a boundary integral scheme and the vorticity transport equation is solved using the ELBEM. Here the results of two‐dimensional Navier–Stokes problems with low–medium Reynolds numbers in a typical cavity flow are presented and compared with a series solution and other numerical models. The ELBEM model has been found to be feasible and satisfactory. Copyright © 2000 John Wiley & Sons, Ltd.     

Bubble‐based stabilized finite element methods for time‐dependent convection–diffusion–reaction problems

In this paper, we propose a numerical algorithm for time‐dependent convection–diffusion–reaction problems and compare its performance with the well‐known numerical methods in the literature. Time discretization is performed by using fractional‐step θ‐scheme, while an economical form of the residual‐free bubble method is used for the space discretization. We compare the proposed algorithm with the classical stabilized finite element methods over several benchmark problems for a wide range of problem configurations. The effect of the order in the sequence of discretization (in time and in space) to the quality of the approximation is also investigated. Numerical experiments show the improvement through the proposed algorithm over the classical methods in either cases. Copyright © 2016 John Wiley & Sons, Ltd.     

Finite element simulation of turbulent Couette–Poiseuille flows using a low Reynolds number k–ϵ model

Developing Couette–Poiseuille flows at Re=5000 are studied using a low Reynolds number k–ϵ two‐equation model and a finite element formulation. Mesh‐independent solutions are obtained using a standard Galerkin formulation and a Galerkin/least‐squares stabilized method. The predictions for the velocity and turbulent kinetic energy are compared with available experimental results and to the DNS data. Second moment closure's solutions are also compared with those of the k–ϵ model. The deficiency of eddy viscosity models to predict dissymmetric low Reynolds number channel flows has been demonstrated. Copyright © 1999 John Wiley & Sons, Ltd.     

Elastic–plastic near field solution of an eccentric crack under shear in a finite width plate

The near crack line analysis method is used to investigate an eccentric crack loaded by shear forces in a finite width plate, and the analytical solution is obtained in this paper. The solution includes: the unit normal vector of the elastic–plastic boundary near the crack line, the elastic–plastic stress fields near crack line, variations of the length of the plastic zone along the crack line with an external loads, and the bearing capacity of a finite plate with a centric crack loaded by shear stress in the far field. The results obtained in this paper are sufficiently precise near the crack line because the assumptions of small scale yielding theory have not been made and no other assumptions have been taken. Subsequently, the present results are compared with the traditional line elastic fracture mechanical solutions and elastoplastic near field solutions under small scale yielding condition. On the basis of the minimum strain energy density (SED) theory, the minimum values of SED in the vicinity of the crack tip are determined, the initial growth orientation of crack are determined. It is found that the normalized load under large scale yielding condition is higher than those under small scale yielding condition when the length of the plastic zone is the same.     

Numerical modelling of step–strain for stretched filaments

This article addresses the modelling of filament-stretching/step–strain deformation under viscoelastic capillary break-up configurations of the CaBER-type. Start-up, prior to step–strain, is conducted under constant stretch-rate synchronous plate retraction with impulsive sessation of plate motion. The study encompasses variation in material rheology, appealing to Oldroyd, Geisekus and Phan-Thien/Tanner-type models, which display differences in shear and extensional viscosity properties (shear thinning/extension hardening). Two different viscosity ratio settings are considered to reflect high- and low-solvent viscosity constituent components; the former representing typical Boger fluids, the latter high-polymer concentration fluids. We compare and contrast results for three alternative filament aspect ratios at the onset of step–strain. Throughout the step–strain period, we have been able to successfully capture such physical features as drainage to the filament feet, necking at the filament centre, and periods with travelling waves through the axial filament length. In addition, we have identified the suppressive influence that larger capillary forces have upon radial fluctuations, and the minor impact that gravitational forces have upon the ensuing deformation. From this study, estimates for rheometrical data have been derived in terms of characteristic material relaxation time and apparent extensional viscosity. The computational techniques employed include a compressed-mesh (CM) procedure, an Arbitrary Lagrangian–Eulerian scheme (ALE) and a free-surface particle tracking technique. Spatial discretisation of the problem is accomplished through a hybrid finite element/finite volume algorithm implemented in the form of a time-stepping incremental pressure-correction formulation.     

Investigation on a coupled Navier–Stokes–Direct Simulation Monte Carlo method for the simulation of plume flowfield of a conical nozzle

To assess the plume effects of space thrusters, the accurate plume flowfield is indispensable. The plume flow of thrusters involves both continuum and rarefied flow regimes. Coupled Navier–Stokes–Direct Simulation Monte Carlo (NS–DSMC) method is a major approach to the simulation of continuum‐rarefied flows. An axisymmetric coupled NS–DSMC solver, possessing adaptive‐interface and two‐way coupling features, is investigated in this paper for the simulation of the nozzle and plume flows of thrusters. The state‐based coupling scheme is adopted, and the gradient local Knudsen number is used to indicate the breakdown of continuum solver. The nitrogen flows in a conical nozzle and its plume are chosen as the reference case to test the coupled solver. The threshold value of the continuum breakdown parameter is studied based on both theoretical kinetic velocity sampling and coupled numerical tests. Succeeding comparisons between coupled and full DSMC results demonstrate their conformities, meanwhile, the former saves 58.8% computational time. The pitot pressure evaluated from the coupled simulation result is compared with the experimental data proposed in literatures, revealing that the coupled method makes precise predictions on the experimental pitot pressure. Copyright © 2014 John Wiley & Sons, Ltd.     

Development of a convection–diffusion–reaction model for solving Maxwell's equations in frequency domain

We propose in this study a numerically accurate and computationally efficient convection–diffusion–reaction finite difference scheme to discretize the full‐vector and semi‐vector optical waveguide equations. The scheme formulated in a grid stencil of five nodal points for solving the three‐dimensional waveguide equations employs the locally analytic solution. In this three‐dimensional study, calculations were carried out for the investigation of wave propagation in diffused channel, rectangular and rib types of optical waveguide. Copyright © 2011 John Wiley & Sons, Ltd.     

A robust and adaptive recovery‐based discontinuous Galerkin method for the numerical solution of convection–diffusion equations

In this paper, we introduce and test the enhanced stability recovery (ESR) scheme. It is a robust and compact approach to the computation of diffusive fluxes in the framework of discontinuous Galerkin methods. The scheme is characterized by a new recovery basis and a new procedure for the weak imposition of Dirichlet boundary conditions. These features make the method flexible and robust, even in the presence of highly distorted meshes. The implementation is simplified with respect to the original recovery scheme (RDG1x). Furthermore, thanks to the proposed approach, a robust implementation of p‐adaptive algorithms is possible. Numerical tests on unstructured grids show a convergence rate equal to p + 1, where p is the reconstruction order. Comparisons are shown with the original recovery scheme RDG1x and the widely used BR2 method. Results show a significantly larger stability region for the proposed discretization when explicit Runge–Kutta time integration is employed. Interestingly, this advantage grows quickly when the reconstruction order is increased. The proposed procedure for the weak imposition of Dirichlet boundary conditions does not need the introduction of ghost cells, and it is truly local because it does not require data exchange with other elements. It can be easily used with curvilinear wall elements. Several test cases are considered. They include some benchmark tests with the heat equation and compressible Navier–Stokes equations, with test cases designed also to evaluate the behaviour of the scheme with very stretched elements and separated flows. Copyright © 2014 John Wiley & Sons, Ltd.     

Near crack line elastic–plastic analysis for a infinite plate loaded by two pairs of point tensile forces

The near crack line analysis method is used to investigate a center crack loaded by two pairs of point tensile forces in an infinite plate in an elastic–perfectly plastic solid, and the analytical solutions are obtained in this paper. These solutions include: the elastic–plastic stress field near the crack line, the law that the length of the plastic zone along the crack line is varied with an external loads and the bearing capacity of an infinite plate with a center crack. The results of this paper are sufficiently precise near the crack line because the assumptions of the small scale yielding theory have not been used and no other assumptions have been taken.     

Simulation of two‐phase flow–body interaction problems using direct forcing/fictitious domain–level set method

In the present paper, a direct forcing/fictitious domain (DF/FD)–level set method is proposed to simulate the twophase flow–body interaction. The DF/FD does not sacrifice accuracy and robustness by employing a discrete δ (Dirac delta) function to transfer quantities between the Eulerian nodes and Lagrangian points explicitly as the immersed boundary method. The advantages of this approach are the simple concept, the easy implementation and the utilization of original governing equation without modification. The main idea is to combine DF/FD method with the level set method in the Cartesian coordinates. We present the results of a number of test cases to illustrate the effectiveness of the proposed method for single‐phase flow–body interaction problem and the two‐phase flows with a stationary body. Eventually, the simulations of various water entry problems have been conducted to validate the capability and the accuracy of the present method on solving the twophase flow–body interaction. Consequently, the present results are found to be in good agreement with those of previous studies. Copyright © 2013 John Wiley & Sons, Ltd.     

A class of higher order compact schemes for the unsteady two‐dimensional convection–diffusion equation with variable convection coefficients

A class of higher order compact (HOC) schemes has been developed with weighted time discretization for the two‐dimensional unsteady convection–diffusion equation with variable convection coefficients. The schemes are second or lower order accurate in time depending on the choice of the weighted average parameter μ and fourth order accurate in space. For 0.5?μ?1, the schemes are unconditionally stable. Unlike usual HOC schemes, these schemes are capable of using a grid aspect ratio other than unity. They efficiently capture both transient and steady solutions of linear and nonlinear convection–diffusion equations with Dirichlet as well as Neumann boundary condition. They are applied to one linear convection–diffusion problem and three flows of varying complexities governed by the two‐dimensional incompressible Navier–Stokes equations. Results obtained are in excellent agreement with analytical and established numerical results. Overall the schemes are found to be robust, efficient and accurate. Copyright © 2002 John Wiley & Sons, Ltd.