Reconstructing the transient,dissipative dynamics of a bistable Duffing oscillator with an enhanced averaging method and Jacobian elliptic functions

Systems characterized by the governing equation of the bistable, double-well Duffing oscillator are ever-present throughout the fields of science and engineering. While the prediction of the transient dynamics of these strongly nonlinear oscillators has been a particular research interest, the sufficiently accurate reconstruction of the dissipative behaviors continues to be an unrealized goal. In this study, an enhanced averaging method using Jacobian elliptic functions is presented to faithfully predict the transient, dissipative dynamics of a bistable Duffing oscillator. The analytical approach is uniquely applied to reconstruct the intrawell and interwell dynamic regimes. By relaxing the requirement for small variation of the transient, averaged parameters in the proposed solution formulation, the resulting analytical predictions are in excellent agreement with exact trajectories of displacement and velocity determined via numerical integration of the governing equation. A wide range of system parameters and initial conditions are utilized to assess the accuracy and computational efficiency of the analytical method, and the consistent agreement between numerical and analytical results verifies the robustness of the proposed method. Although the analytical formulations are distinct for the two dynamic regimes, it is found that directly splicing the inter- and intrawell predictions facilitates good agreement with the exact dynamics of the full reconstructed, transient trajectory.     

内聚滑移爆轰加载下金属圆管压缩过程的工程分析

对内聚滑移爆轰加载下圆管压缩运动过程进行工程分析,分析解与数值模拟结果吻合。利用工程模型得到了不同管道参数（包括内径、厚度）下管道封闭所需要的炸药量和封闭时间表达式,可为相关研究工作提供参考。     

Numerical evaluation of generalized stress-intensity factors in multi-layered composites

The problem of the evaluation of the generalized stress-intensity factors for re-entrant corners in multi-layered structural components is addressed. An approximate analytical model based on the theory of multi-layered beams is presented. This approach provides a simple closed-form solution for the direct computation of the Mode I stress-intensity factor for the general problem of a re-entrant corner symmetrically meeting a bi-material interface. For the self-consistency of the theory, re-entrant corners in homogeneous materials and cracks perpendicular to bi-material interfaces can also be gained as limit cases of this formulation. According to this approach, the effects of the elastic mismatch parameters, the value of the notch angle and the thicknesses of the layers on the stress-intensity factor are carefully quantified and the results are compared with FE solutions. FE results are obtained by applying a combination of analytical and numerical techniques based on the knowledge a priori of the asymptotic stress field for re-entrant corners perpendicular to a bi-material interface and on the use of generalized isoparametric singular finite elements at the notch tip. A good agreement between approximate and analytical/numerical predictions is achieved, showing the effectiveness of this approach.     

Analytical and numerical study of buoyancy-driven convection in a vertical enclosure filled with nanofluids

This paper presents an analytical and numerical study of natural convection of nanofluids contained in a rectangular enclosure subject to uniform heat flux along the vertical sides. Governing parameters of the problem under study are the thermal Rayleigh number Ra, the Prandtl number Pr, the aspect ratio of the cavity A and the solid volume fraction of nanoparticles, Φ. Three types of nanoparticles are taken into consideration: Cu, Al₂O₃ and TiO₂. Various models are used for calculating the effective viscosity and thermal conductivity of nanofluids. In the first part of the analytical study, a scale analysis is made for the boundary layer regime situation. In the second part, an analytical solution based on the parallel flow approximation is reported for tall enclosures (A ≫ 1). In the boundary layer regime a good agreement is obtained between the predictions of the scale analysis and those of the analytical solution. Solutions for the flow fields, temperature distributions and Nusselt numbers are obtained explicitly in terms of the governing parameters of the problem. A numerical study of the same phenomenon, obtained by solving the complete system of the governing equations, is also conducted. A good agreement is found between the analytical predictions and the numerical simulations.     

分数阶达芬振子的超谐与亚谐联合共振

研究了含分数阶微分项的达芬(Duffing)振子的超谐与亚谐联合共振.采用平均法得到了系统的一阶近似解析解,提出了超、亚谐联合共振时等效线性阻尼和等效线性刚度的概念.建立了联合共振定常解幅频曲线的解析表达式,并对联合共振幅频响应的近似解析解和数值解进行了比较,二者吻合良好,证明了求解过程及近似解析解的正确性.然后,将等效线性阻尼和等效线性刚度的概念与传统整数阶系统进行比较,证明分数阶微分项不仅起着阻尼的作用同时还起着刚度的作用.最后,通过数值仿真研究了不同的分数阶微分项系数和阶次对联合共振幅频曲线多值性和跳跃现象的影响,并与单一频率下超谐共振或亚谐共振进行了对比.研究发现,分数阶微分项系数与阶次不仅影响着系统的响应幅值、共振频率,同时还对系统的周期解个数、发生区域面积、发生先后等有重要影响.并且,在不同的基本参数下该系统分别表现出单独超谐共振、单独亚谐共振以及超谐共振和亚谐共振同时存在的现象.这些结果对系统动力学特性的研究具有重要意义.     

Development of an approach to constitutive modelling of concrete: Isotropic damage coupled with plasticity

The paper presents an approach to constitutive modelling of concrete using damage mechanics and plasticity theory. The thermodynamic formulation, and parameter identification of a non-local coupled damage-plasticity model are presented in this study. The particular focus is the calibration of model parameters. It is shown that both the local parameters and the parameters governing the non-local interaction can be determined from experimental data reliably and consistently. A novel procedure is developed for parameter identification, using the separation of total dissipation energy into additive parts corresponding to different dissipation mechanisms. The relationship between the local and non-local parameters is also addressed, helping to obtain model responses consistent with the fracture energy of the material. The application of the model and the calibration procedure proposed in this study to the numerical failure analysis of concrete structures is illustrated through a series of real structural tests, showing both the performance of the model and the consistency of the proposed calibration procedure.     

计及表面能的应变梯度理论本构参数识别

本文提出利用静态位移信息对一种计及表面能的应变梯度理论本构参数进行识别的求解策略.基于Vardoulakis和Sulem的计及表面能的简单线性应变梯度理论,文献[13]给出了伯努力-欧拉梁弯曲问题的正演解析模型,本文将其反演归结为两个带有不等式约束的非线性规划问题.在此基础上,采用黄金分割一维搜索方法进行求解,给出了数值验证,讨论了信息误差对反演结果的影响.结果显示,这种方法可以用来对应变梯度理论本构参数进行识别,即使在体积和表面能常数非常小的情况下,仍然能够得到满意的结果.     

Analytical solutions for resonant response of suspended cables subjected to external excitation

In this study, two analytical methods are applied to study the primary resonances response of suspended cables subjected to external excitation. We choose four different sag-to-span ratios and the first two modes to investigate the differences in nonlinear responses obtained with analytical methods. First, we summarize the equations of motion by applying the Hamilton’s principle and quasi-static assumption, and then these equations are discretized by the Galerkin procedure. Second, the multiple-scale method and homotopy analysis method are adopted to obtain the approximate solutions. Moreover, numerical integrations are introduced in order to verify the obtained approximate results. The numerical results show that frequency response curves obtained by different analytical methods show different quantitative predictions in some cases of motion, modes, and particular sag-to-span ratios. Finally, the differences in displacement fields and axial tension forces are compared and analyzed.     

Analytical and numerical study of double diffusive convection in a vertical enclosure

This paper presents an analytical and numerical study of natural convection of a double-diffusive fluid contained in a rectangular slot subject to uniform heat and mass fluxes along the vertical sides. Governing parameters of the problem under study are the thermal Rayleigh number, Ra _T ; buoyancy ratio, N; Lewis number, Le; Prandtl number, Pr and aspect ratio of the cavity, A. In the first part of the analytical study a scale analysis is applied to the two extreme cases of heat-transfer and mass-transfer-driven flows. In the second part, an analytical solution, based on the parallel flow approximation, is reported for tall enclosures (A?1). Solutions for the flow fields, temperature and concentration distributions and Nusselt and Sherwood numbers are obtained in terms of the governing parameters of the problem. In the limits of heat-driven and solute-driven flows a good agreement is obtained between the prediction of the scale analysis and those of the analytical solution. The numerical solutions are based on the complete governing equations for two-dimensional flows, and cover the range 1≤Ra _T ≤10⁷, 0≤N≤10⁵, 10^-3≤Le≤10³, 1≤A≤20 and Pr=7. A good agreement is found between the analytical predictions and the numerical simulation.     

The role of plastic deformation on the impact behaviour of high aspect ratio aluminium foam-filled sections

The main objective of this work is to investigate the role of the plastic deformation of metal foams on the dynamic behaviour of aluminium foam-filled columns with respect to their energy absorbing capabilities. The influence of the cross-section shape as well as other parameters is thoroughly studied. A comparison with correspondent hollow-sections is performed concerning the dissipation of kinetic energy and the obtained deformed profiles. For this particular purpose, three-dimensional finite element modelling dynamic analyses are carried out using ABAQUS/Explicit in order to achieve an in-depth study of the structural crash behaviour, during which energy needs to be absorbed in a controlled manner. A comprehensive numerical study of the crush behaviour of aluminium foam-filled sections undergoing axial compressive loading is performed. The results obtained are also analysed with respect to the reduction in the length of the structural element and impact time, the effect of friction between the foam and the outer skin, the energy decomposition, the role of plastic deformation, the influence of the skin material and impact velocity, and the influence of the shape of the cross-section on the impact behaviour. A comparison with existing analytical expressions is made in order to corroborate the numerical results.     

Numerical Calculation of Mass Transfer from Elliptical Pools in Uniform Flow Using the Boundary Element Method

The three-dimensional problem of advection-dispersion associated with an elliptical non-aqueous-phase liquid (NAPL) pool is addressed using the boundary element method. The boundary condition on the plane of the pool is such that over the pool the concentration is equal to the saturation concentration while a no flux boundary condition is imposed in the region not covered by the pool. The numerical results are verified by asymptotic analytical solutions obtained in the limits of diffusion-dominated and convection-dominated mass transport. For cases of practical interest an empirical expression is obtained for the Sherwood number that matches the numerical results over a wide range of the relevant parameters. Comparison with experimental results suggests that the corresponding numerical results predict a higher overall mass transfer coefficient.     

Multi-frequency oscillations in self-excited systems

A self-excited three-mass chain system is considered here. For a self-excitation of van der Pol type, the possibility of multi-frequency oscillations is investigated. Both analytical approximate solutions and numerical simulation are used. The averaging method is used to establish existence and stability of the normal modes, the two-frequency modes as well as the three-frequency oscillations solutions. We found at first that the single mode seems to prevail. However a three-frequency solution can be stabilised by adapting the system slightly. A generic bifurcation diagram is given where all the possible phase portraits are sketched. The flow turns out to be quite predictable. There is no “room” for chaos or strange attractors. This behaviour is not typical for systems of coupled oscillators but turns out to be partly related to the involved symmetries as well as the particular choice of the system parameters.     

四节点二十四自由度平板壳单元几何刚度矩阵显式解析式的推演算法研究

几何刚度矩阵的推演是结构几何非线性有限元分析的重点和难点之一。推导几何刚度矩阵显式解析表达式成为简化非线性有限元列式,提高分析效率的关键。本文在协同转动法框架下,基于刚体运动法则对四节点二十四自由度的平板壳单元几何刚度矩阵显式解析式进行了推导和讨论;分析了悬臂梁大转动、不同壁厚条件下简支圆柱形屋顶空间大变位两个经典算例。研究结果表明：（1）几何刚度矩阵的显式计算公式不仅为板壳结构几何非线性列式提供了方便而且具有良好的精度;（2）推导的几何刚度矩阵适用于各类型四边形二十四自由度平板壳单元模型;（3）与数值积分相比,采用解析形式的几何刚度矩阵可以显著提高非线性响应计算效率。     

Complex variable approach in studying modified polarization saturation model in two-dimensional semipermeable piezoelectric media

A modified polarization saturation model is proposed and addressed mathematically using a complex variable approach in two-dimensional(2 D) semipermeable piezoelectric media. In this model, an existing polarization saturation(PS) model in 2D piezoelectric media is modified by considering a linearly varying saturated normal electric displacement load in place of a constant normal electric displacement load, applied on a saturated electric zone. A centre cracked infinite 2D piezoelectric domain subject to an arbitrary poling direction and in-plane electromechanical loadings is considered for the analytical and numerical studies. Here, the problem is mathematically modeled as a non-homogeneous Riemann-Hilbert problem in terms of unknown complex potential functions representing electric displacement and stress components. Having solved the Hilbert problem, the solutions to the saturated zone length, the crack opening displacement(COD), the crack opening potential(COP), and the local stress intensity factors(SIFs) are obtained in explicit forms. A numerical study is also presented for the proposed modified model, showing the effects of the saturation condition on the applied electrical loading, the saturation zone length, and the COP. The results of fracture parameters obtained from the proposed model are compared with the existing PS model subject to electrical loading, crack face conditions, and polarization angles.     

Fractal porous media

The transport properties of continuous deterministic fractals are reviewed. The method of construction, the fractal dimension, and the major features of transport are summarized. Then the major single-phase transports are addressed; attention is focused on the numerical results and on the analytical arguments which may be used to derive these results in a simple way, whenever it is possible.     

结构损伤识别的非线性解析法

本文在直接解析法的基础上,进一步提出结构损伤识别的非线性解析法.它与直接解析法不同之处是:考虑了频率和振型对损伤参数的二阶偏导数项,以此构造损伤参数为未知量的超定非线性方程组.求解该方程组可以一次性得到全部损伤参数值.本方法不仅具有直接解析法的优点,而且由于构建方程组时增加了损伤参数的二次项,极大地提高了损伤识别的精度.最后,通过对一个五层框架结构进行数值模拟分析验证了该方法的可行性,同时经过比较计算结果表明:在损伤识别中,非线性解析法在精度及节省计算工作量上明显优于直接解析法.     

Interaction between a submerged evacuated cylindrical shell and a shock wave—Part II: Numerical aspects of the solution

The interaction between a submerged elastic circular cylindrical shell and an external shock wave is addressed. A linear, two-dimensional formulation of the problem is considered. A semi-analytical solution is obtained using a combination of the classical analytical approach based on the use of the Laplace transform and separation of variables, and finite difference methodology. The study consists of two parts. Part I focuses on the simulation and analysis of the acoustic fields induced during the interaction. Both the diffraction (absolutely rigid cylinder) and complete diffraction–radiation (elastic shell) are considered. Special attention is paid to the lower-magnitude shell-induced waves representing radiation by the elastic waves circumnavigating the shell. The focus of Part II is on the numerical analysis of the solution. The convergence of the series solution and finite-difference scheme is analysed. The computation of the response functions of the problem is discussed as well, as is the effect of the bending stiffness on the acoustic field. The membrane model of the shell is considered to analyse such effect, which, in combination with the models addressed in Part I, allows for the analysis of the evolution of the acoustic field around the structure as its elastic properties change from an absolutely rigid cylinder to a membrane. The results of the numerical simulations are compared to available experimental data, and a good agreement is observed.     

Interaction between a submerged evacuated cylindrical shell and a shock wave—Part I: Diffraction–radiation analysis

The interaction between a submerged elastic circular cylindrical shell and an external shock wave is addressed. A linear, two-dimensional formulation of the problem is considered. A semi-analytical solution is obtained using a combination of the classical analytical approach based on the use of the Laplace transform and separation of variables, and finite difference methodology. The study consists of two parts. Part I focuses on the simulation and analysis of the acoustic fields induced during the interaction. Both the diffraction (absolutely rigid cylinder) and complete diffraction–radiation (elastic shell) are considered. Special attention is paid to the lower-magnitude shell-induced waves representing radiation by the elastic waves circumnavigating the shell. The focus of Part II is on the numerical analysis of the solution. The convergence of the series solution and finite-difference scheme is analysed. The computation of the response functions of the problem is discussed as well, as is the effect of the bending stiffness on the acoustic field. The membrane model of the shell is considered to analyse such an effect, which, in combination with the models addressed in Part I, allows for the analysis of the evolution of the acoustic field around the structure as its elastic properties change from an absolutely rigid cylinder to a membrane. The results of the numerical simulations are compared to available experimental data, and a good agreement is observed.