A DQ based approach to simulate the vibrations of buckled beams

In this paper, the nonlinear planar response of a hinged–hinged buckled beam to a primary-resonance excitation of its first vibration mode is computed by a new numerical scheme. The beam is subjected to an axial force beyond the critical load of the first buckling mode and to a transverse harmonic excitation. The nonlinear dynamical problem is solved by deducing directly the discretized equations governing the problem thanks to a new approach, here called DQ based approach, since it is based on the application of the quadrature rules of the DQM. As it will be shown, for the problem here considered, the minimum number of degrees of freedom to be retained to limit the numerical errors is four. Computer simulations of the dynamic behaviour of the discretized system are conducted by means of the IDQ method, a method proposed and recently generalized by the author. A sequence of supercritical period-doubling bifurcations leading to chaos, snapthrough motions and quasi-periodic motions can be observed, similarly to some cases existing in literature.     

Efficient quadrature for a boundary element method to compute three-dimensional Stokes flow

A collocation-type boundary element method based on bilinear B-splines is used for the numerical solution of the Stokes Dirichlet problem in bounded domains D ? R³. The computation of the influence matrix requires the numerical evaluation of weakly singular integrals on the domain boundary if the usual double-layer potential ansatz is chosen. Here mostly standard methods with disjoint grids for collocation and integration are used. We develop a special integration scheme based on triangular co-ordinates near the singularity and show its efficiency compared with the method mentioned above.     

Quadrature method of moments for modeling multi-component spray vaporization

This article puts forward the quadrature method of moments (QMoM) for modeling droplet composition during the spray vaporization process. This method is implemented for solving the Continuous Thermodynamic Model (CTM) of multi-component droplet vaporization, an advantageous alternative to the classical Discrete Component Model (DCM) when the droplet is formed of a great number of components. The CTM approach consists in modeling the droplet’s composition using a probability density function (PDF). This method was first tried out for vaporizing droplets by Hallett, who assumed a Gamma-function for the PDF. However, Harstadt et al. underlined some problems in the case of vapor condensation on the droplet surface, since the Gamma-PDF model presumes the PDF’s mathematical form. The QMoM which does not require this hypothesis is studied in this article, according to Lage’s research dealing with QMoM application to phase equilibria. The numerical features of QMoM are investigated in detail, and then the method is implemented for the difficult test case of vapor condensation. The results are analyzed to illustrate the application of QMoM to multi-component droplet vaporization modeling and to provide a better understanding of the QMoM main advantages and limitations.     

Consistent vs. reduced integration penalty methods for incompressible media using several old and new elements

The frequently used reduced integration method for solving incompressible flow problems ‘a la penalty’ is critically examined vis-a-vis the consistent penalty method. For the limited number of quadrilateral and hexahedral elements studied, it is shown that the former method is only equivalent to the latter in certain special cases. In the general case, the consistent penalty method is shown to be more accurate. Finally, we demonstrate significant advantages of a new element, employing biquadratic (2-D) or triquadratic (3-D) velocity and linear pressure over that using the same velocity but employing bilinear (2-D) or trilinear (3-D) pressure approximation.     

Numerical breakdown at high Weissenberg number in non-Newtonian contraction flows

 Planar contraction flows of non-Newtonian fluids with integral constitutive models are studied to investigate the problem of numerical breakdown at high Weissenberg or Debrorah numbers. Spurious shear stress extrema are found on the wall downstream of the re-entrant corner for both sharp and rounded corners. Moreover, a non-monotonic relation between shear stress and strain rate is found when the Deborah number limit is approached, which correlates with these shear extrema. This strongly suggests that non-monotonicity between shear stress and strain rate may be responsible for the Deborah number limit problem in contraction flow simulations. This non-monotonicity is caused by the inaccuracy of the quadrature, using constitutive equations that do not have shear stress maxima when exactly evaluated. This conclusion agrees with recent analytical findings by others that inaccuracy of the integration along the streamlines – either by numerical integration or asymptotic approximation – makes the problem ill-conditioned, with spurious growth occurring on the wall downstream of the re-entrant corner. Received: 5 March 1999/Accepted: 1 September 1999     

A Simplified Quadrature Element Method to compute the natural frequencies of multispan beams and frame structures

In this paper a new version of the Modified Quadrature Element Method (MQEM) is proposed. Like MQEM, the proposed method overcomes the drawback of the distance δ of the Quadrature Element Method (QEM) without introducing further degrees of freedom at the ends of the element as in the Differential Quadrature Element Method (DQEM), but it makes the computational cost of the stiffness matrix (and the mass matrix) lighter and uses a general procedure to generate the sampling points distribution. The method here presented has been applied to compute the fundamental frequencies of some structures.     

Winkler-Pasternak弹性地基梁自由振动的二维弹性分析Two-dimensional elastic analysis for free vibration of beams set on winkler-pasternak elastic foundations

基于二维线弹性理论,应用Hamilton原理,获得Winkler-Pasternak弹性地基梁自由振动的控制微分方程,应用微分求积法(DQM)数值研究了梁自由振动的无量纲频率特性。计算结果与已有的结果(Bernoulli-Euler梁和Timoshenko梁)比较表明,本文的分析方法对弹性地基长梁和短梁自由振动的研究都有效。最后考虑了几何参数对梁频率的影响,以及不同边界条件下地基系数对频率的影响和收敛性。     

Winkler-Pasternak弹性地基梁自由振动的二维弹性分析

基于二维线弹性理论,应用Hamilton原理,获得Winkler-Pasternak弹性地基梁自由振动的控制微分方程,应用微分求积法(DQM)数值研究了梁自由振动的无量纲频率特性。计算结果与已有的结果(Bernoulli-Euler梁和Timoshenko梁)比较表明,本文的分析方法对弹性地基长梁和短梁自由振动的研究都有效。最后考虑了几何参数对梁频率的影响,以及不同边界条件下地基系数对频率的影响和收敛性。     

基于改进降维法的可靠度分析

单变量维数缩减法可以高效、准确地进行结构响应矩的分析.与传统的一阶可靠度算法FORM(First Order Reliability Method),二阶可靠度算法SORM(Second Order Reliability Method)相比,该方法不需要响应的导数,也不需要迭代搜索最可能失效点.然而近期的研究发现,该...     

Non-axisymmetric bending of thin annular plates due to circumferentially distributed moments

The non-linear deformation of a thin annular plate subjected to circumferentially distributed bending moments is studied. A von Kármán plate model is adopted to formulate the equations of motion. Free–free boundary conditions have been applied at both inner and outer edges. The matrix formulation of the Differential Quadrature Method is used to discretise and solve the governing equations. Linear analysis predicts that the annular disk deforms axisymmetrically into a spherical dome. However, the proposed non-linear analysis shows that a symmetry breaking bifurcation may occur after which the linear solution becomes unstable and the plate transitions into a non-axisymmetric cylindrical deformation. This is the case when at least one of two parameters reaches a critical value. These parameters are the non-dimensionalised ratio between applied moment and bending stiffness and the ratio between inner and outer radius. Furthermore, it is noted that free-free boundary conditions and circumferentially distributed bending moments do not break the circular symmetry of the annular disk. Hence, the principal axes of curvature in the deformed configuration do not have a preferred orientation. Therefore, the present work describes a shell possessing infinite identical equilibria, having different yet no favoured direction and, hence, links to previous researches on neutrally stable structures.     

A moment method for splashing and evaporation processes of polydisperse sprays

A new moment method for the modelling of polydisperse sprays is proposed that simultaneously takes into account the dispersion in droplet size and droplet velocity. For the derivation of this Eulerian method the kinetic spray equation is used which constitutes a partial differential equation for the probability density function of droplets. To reduce the complex kinetic spray equation to a form that can be managed with the available numerical procedures, moment transforms with respect to the droplet velocity and the droplet size are conducted. The resulting moment equations are closed by choosing an approximate probability density function which applies to polydisperse sprays. The method is successfully tested for configurations in which a polydisperse spray is either splashed, evaporated or effected by a Stokes drag force. The tests are organised in such a way that crossing of two spray distributions is always included. The new method is able to capture the polydisperse nature of sprays as well as the bi-(or multi-) modal character of the droplet velocity distribution function, for example, when droplets cross each other.     

Flexural vibration analysis of an eccentric annular Mindlin plate

A weak-form quadrature element method is presented to study the flexural vibrations of an eccentric annular Mindlin plate. Typical combinations of boundary conditions are considered and the natural frequencies are obtained for both thin and moderately thick plates. All results are verified using the commercial computer code ANSYS. Excellent agreement is reached in all cases. Comparison of the present predictions with other available results for thin plates is also made.     

三阶非线性发展方程的线性化有理逼近方法

研究了三阶非线性发展方程的初边值问题的解。采用基于Sinc函数的微分求积法发展了线性化有理逼近方法。通常的配点法不适用于上述三阶问题的求解。本文把提出的方法用于求解KdV方程,取得了良好的效果。     

微分求积区域分裂法在裂缝问题上的应用

微分求积法DQM在处理裂缝问题时,会产生很大的误差。因此,本文用微分求积法结合不带重叠的区域分裂法DQDDM来求解。通过本文的讨论,可以看到DQDDM在处理裂缝问题时,在节点数目不多的条件下获得比较精确的解,同时计算量又不大。     

FGM环扇形板的面内自由振动分析

假定功能梯度材料(FGM)的物性参数沿环扇形板径向按照幂律梯度变化,基于平面线弹性理论,建立了FGM环扇形板面内自由振动的运动控制微分方程。采用二维微分求积法(DQM)对FGM环扇形板面内自由振动的无量纲运动控制微分方程进行离散,数值求解了不同边界条件下FGM环扇形板面内自由振动的无量纲固有频率,同时也给出了FGM环扇形板扇形角为!/4时有限元商用软件ANSYS的部分计算结果,验证了本文方法的正确性。结果表明,在相应边界条件下,FGM环扇形板的梯度指标、内外半径比以及扇形角对无量纲固有频率均有影响,其计算结果和分析方法可供设计和研究参考。     

Continuous thermodynamics for droplet vaporization: Comparison between Gamma-PDF model and QMoM

The Continuous Thermodynamics Model (CTM) (Cotterman et al., 1985) is a suitable method to reduce computational cost of multi-component vaporization models. The droplet composition is described by a probability density function (PDF) rather than tens of components in the classical Discrete Component Model (DCM). In the first CTM method developed for this application, the PDF was assumed to be a Γ-function (Hallett, 2000), but some problems had appeared in the case of vapor condensation at the droplet surface (Harstadt et al., 2003). The method put forward in this article, the Quadrature Method of Moments (QMoM), enables one to avoid any assumption on the PDF mathematical form. Following Lage who has developed this method for phase equilibria (Lage, 2007), this article widens the scope of QMoM to the modelling of multi-component droplet vaporization. The different CTM approaches are presented in the first part and the results obtained for a vapor condensation test case are then compared and analysed to illustrate improvements made by QMoM. To cite this article: C. Laurent et al., C. R. Mecanique 337 (2009).     

A dual-porosity model for gas reservoir flow incorporating adsorption behaviour—Part II. Numerical algorithm and example applications

In the present study, an algorithm is presented for the dual-porosity model formulated in Part I of this series. The resultant flow equation with the dual-porosity formulation is of an integro-(partial) differential equation involving differential terms for the Darcy flow in large fractures and integrals in time for diffusion within matrix blocks. The algorithm developed here to solve this equation involves a step-by-step finite difference procedure combined with a quadrature scheme. The quadrature scheme, used for the integral terms, is based on the trapezoidal method which is of second-order precision. This order of accuracy is consistent with the first- and second-order finite difference approximations used here to solve the differential terms in the derived flow equation. In an approach consistent with many petroleum reservoir and groundwater numerical flow models, the example formulation presented uses a first-order implicit algorithm. A two-dimensional example is also demonstrated, with the proposed model and numerical scheme being directly incorporated into the commercial gas reservoir simulator SIMED II that is based on a fully implicit finite difference approach. The solution procedure is applied to several problems to demonstrate its performance. Results from the derived dual-porosity formulation are also compared to the classic Warren–Root model. Whilst some of this work confirmed previous findings regarding Warren–Root inaccuracies at early times, it was also found that inaccuracy can re-enter the Warren–Root results whenever there are changes in boundary conditions leading to transient variation within the domain.     

A dynamical strategy for approximation methodsUne stratégie dynamique pour les méthodes d'approximation

The numerical result provided by an approximation method is affected by a global error, which consists of both a truncation error and a round-off error. Let us consider the converging sequence generated by successively dividing by two the step size used. If computations are performed until, in the convergence zone, the difference between two successive approximations is only due to round-off errors, then the global error on the result obtained is minimal. Furthermore its significant bits which are not affected by round-off errors are in common with the exact result, up to one. To cite this article: F. Jézéquel, C. R. Mecanique 334 (2006).     

The 7th International Conference on Fluid Mechanics May 24-27,2015,Qingdao,China

正http://www.icfm7.org First Announcement and Call for PapersThe objective of International Conference on Fluid Mechanics(ICFM)is to provide a forum for researchers to exchange new ideas and recent advances in the fields of theoretical,experimental,computational Fluid Mechanics as well as interdisciplinary subjects.It was successfully convened by the Chinese Society of Theoretical and Applied Mechanics(CSTAM)in Beijing(1987,     

INFORMATION FOR CONTRIBUTORS

Contributions： The Journal, Acta Mechanica Solida Sinica, is pleased to receive papers from engineers and scientists working in various aspects of solid mechanics. All contributions are subject to critical review prior to acceptance and publication.