Dynamic elastoplastic analysis by a hybrid BEM–FEM time-domain formulation

This study presents a hybrid BEM–FEM procedure for the dynamic analysis of elastoplastic models. In this hybrid approach, boundary node and internal point displacements are evaluated considering the time-domain BEM formulation (initial stress approach), and stresses are computed taking into account FEM techniques (domain discretization is only necessary where non-linear behaviour is expected to occur). This hybrid methodology is very appropriate to model infinite or semi-infinite elastoplastic models and, at the end of the paper, three numerical applications are presented, illustrating the potentialities of the proposed formulation.     

Stability boundaries of two-parameter non-linear elastic structures

The load-bearing capacity of structures can be influenced by variations in parameters, such as initial geometric defects, multi-parameter loadings, material specifications and temperature. This paper aims to introduce a new formulation to trace the stability boundaries of two-parameter elastic structures. The proposed procedure can find a set of critical points, both limit and bifurcation ones, via a modified Newton’s method. In the authors’ formulation, the residual force is set to zero, and a critically constraint is satisfied simultaneously. Numerical examples presented in this paper demonstrate the efficiency of the suggested method.     

Einige Aspekte der dynamischen Simulation von Stofftrennapparaten im Hinblick auf ihre Übertragung auf reale Anlagen

Because of the usual formulation of mathematical models for mass transfer apparatus in mol-units, serious differences can arise in their dynamic behaviour, when some independent operating variables are kept constant during the simulation-for example, product streams-in kmol/s, kg/s or in m³/s. The reason for this fact is that the average product molmass usually changes when a disturbance spreads. This article presents the dynamic transient behaviour of distillation columns in view of this problem. Here, two examples will be analyzed and discussed. The results show the importance of taking in consideration the measurement technology in the dynamic simulation of mass transfer apparatus, if the simulation results are to be transfered later to the production unit.     

A General Formulation and Solution Scheme for Fractional Optimal Control Problems

Accurate modeling of many dynamic systems leads to a set of Fractional Differential Equations (FDEs). This paper presents a general formulation and a solution scheme for a class of Fractional Optimal Control Problems (FOCPs) for those systems. The fractional derivative is described in the Riemann–Liouville sense. The performance index of a FOCP is considered as a function of both the state and the control variables, and the dynamic constraints are expressed by a set of FDEs. The Calculus of Variations, the Lagrange multiplier, and the formula for fractional integration by parts are used to obtain Euler–Lagrange equations for the FOCP. The formulation presented and the resulting equations are very similar to those that appear in the classical optimal control theory. Thus, the present formulation essentially extends the classical control theory to fractional dynamic system. The formulation is used to derive the control equations for a quadratic linear fractional control problem. An approach similar to a variational virtual work coupled with the Lagrange multiplier technique is presented to find the approximate numerical solution of the resulting equations. Numerical solutions for two fractional systems, a time-invariant and a time-varying, are presented to demonstrate the feasibility of the method. It is shown that (1) the solutions converge as the number of approximating terms increase, and (2) the solutions approach to classical solutions as the order of the fractional derivatives approach to 1. The formulation presented is simple and can be extended to other FOCPs. It is hoped that the simplicity of this formulation will initiate a new interest in the area of optimal control of fractional systems.     

延性动态损伤的细观描述

对动载荷下材料的延性动态损伤的细观力学研究现状进行了评述，详细讨论了现有的几类典型的动态延性损伤的微孔生长模型。针对目前该领域工作中存在的问题，指出了今后应开展的工作方向。     

Numerical study of PPE source term errors in the incompressible SPH models

In a recent work by Gui et al. 13 , an incompressible SPH model was presented that employs a mixed pressure Poisson equation (PPE) source term combining both the density‐invariant and velocity divergence‐free formulations. The present work intends to apply the model to a wider range of fluid impact situations in order to quantify the numerical errors associated with different formulations of the PPE source term in incompressible SPH (ISPH) models. The good agreement achieved between the model predictions and the documented data is taken as a further demonstration that the mixed source term formulation can accurately predict the fluid impact pressures and forces, both in the magnitude and in the spatial and temporal patterns. Furthermore, an in‐depth numerical analysis using either the pure density‐invariant or velocity divergence‐free formulation has revealed that the pure density‐invariant formulation can lead to relatively large divergence errors while the velocity divergence‐free formulation may cause relatively large density errors. As compared with these two approaches, the mixed source term formulation performs much better having the minimum total errors in all test cases. Although some recent studies found that the weakly compressible SPH models perform somewhat better than the incompressible SPH models in certain fluid impact problems, we have shown that this could be largely caused by the particular formulation of PPE source term in the previous ISPH models and a better formulation of the source term can significantly improve the accuracy of ISPH models. Copyright © 2014 John Wiley & Sons, Ltd.     

A General Formulation and Solution Scheme for Fractional Optimal Control Problems

Accurate modeling of many dynamic systems leads to a set of Fractional Differential Equations (FDEs). This paper presents a general formulation and a solution scheme for a class of Fractional Optimal Control Problems (FOCPs) for those systems. The fractional derivative is described in the Riemann–Liouville sense. The performance index of a FOCP is considered as a function of both the state and the control variables, and the dynamic constraints are expressed by a set of FDEs. The Calculus of Variations, the Lagrange multiplier, and the formula for fractional integration by parts are used to obtain Euler–Lagrange equations for the FOCP. The formulation presented and the resulting equations are very similar to those that appear in the classical optimal control theory. Thus, the present formulation essentially extends the classical control theory to fractional dynamic system. The formulation is used to derive the control equations for a quadratic linear fractional control problem. An approach similar to a variational virtual work coupled with the Lagrange multiplier technique is presented to find the approximate numerical solution of the resulting equations. Numerical solutions for two fractional systems, a time-invariant and a time-varying, are presented to demonstrate the feasibility of the method. It is shown that (1) the solutions converge as the number of approximating terms increase, and (2) the solutions approach to classical solutions as the order of the fractional derivatives approach to 1. The formulation presented is simple and can be extended to other FOCPs. It is hoped that the simplicity of this formulation will initiate a new interest in the area of optimal control of fractional systems.     

A model of an elastic-plastic medium with delayed yield

Stress-strain relationships for metals at high strain rates have long been studied, but no really reliable and generally accepted theory has emerged. It is sometimes assumed that the dynamic stress-strain diagram is largely insensitive to the rate over a certain range. Another approach is to insert derivatives of the stress and strata with respect to time. One difficulty in establishing the actual reIationships is that experiment provides only indirect evidence (direct tests are usually impossible). Any real dynamic experiment tends to produce complicated effects, which can be interpreted only if the basic equations are taken as known. The best that experiment can then do is to confirm or reject some prior assumptions.Many experimental studies deal with mechanical characteristics such as breaking strength and yield point as functions of strain rate; however, strain rate characterizes a range of conditions rather than defines a parameter. We therefore have to use simple models that allow formulation and solution of definite mechanical problems in relation to the dynamics of elastic-plastic media.     

Static and dynamic studies for coupling discrete and continuum media; Application to a simple railway track model

In this paper, we present a formulation for coupling discrete and continuum models for both dynamic and static analyses. This kind of formulation offers the possibility of carrying out better simulations of material properties than the discrete calculations, and with both larger length scales and longer times. Using only a discrete approach to simulate a large medium composed of many degrees of freedom seems very difficult in terms of calculation and implementation. Moreover, using only a continuum approach does not give an accurate solution in a zone where particular and localized phenomena can occur. A direct application of our coupling approach to the case of railway track models subjected to an external load, is proposed for its validation.     

A constitutive model accounting for strain ageing effects on work-hardening. Application to a C–Mn steel

One of the most successful models for describing the Portevin–Le Chatelier effect in engineering applications is the Kubin–Estrin–McCormick model (KEMC). In the present work, the influence of dynamic strain ageing on dynamic recovery due to dislocation annihilation is introduced in order to improve the KEMC model. This modification accounts for additional strain hardening rate due to limited dislocation annihilation by the diffusion of solute atoms and dislocation pinning at low strain rate and/or high temperature. The parameters associated with this novel formulation are identified based on tensile tests for a C–Mn steel at seven temperatures ranging from 20 °C to 350 °C. The validity of the model and the improvement compared to existing models are tested using 2D and 3D finite element simulations of the Portevin–Le Chatelier effect in tension.     

Nonlinear dynamics of a cable–pulley system using the absolute nodal coordinate formulation

It is reasonable to develop models and to investigate the dynamic behaviour of systems composed of cables since cable vibration can have an important effect on the motion of these mechanical systems. This paper deals with the application of the nonlinear formulation for flexible body dynamics called the absolute nodal coordinate formulation (ANCF). It is used for modelling the systems composed of cables, pulleys, other rigid bodies and a motor with prescribed motion. The ANCF was chosen as a suitable approach, which that can allow to consider a detailed interaction of the cable and the pulley with its nonlinear dynamical behaviour. The ANCF uses absolute positions of nodes (reference vectors) and slopes (reference vector derivations) as a set of nodal coordinates. An in-house modelling tool in the MATLAB system was created based on the proposed modelling methodology and two case studies were performed. A simple system containing a pulley and a cable with two attached bodies was used in order to test the simulation tool based on the proposed modelling methodology with respect to different parameters. A more complex mechanical system composed of a driven weight joined with a motor by a cable led over a pulley was numerically and also experimentally investigated. The comparison of obtained numerical and experimental results shows sufficient agreement and proves that the proposed modelling approach can be used for dynamic analyses of such systems.     

Impact of the interplanetary magnetic field on thewave flow pattern in the neighborhood of the Earth’s bow shock under sharp variations in the solar wind dynamic pressure

The impact of the interplanetary magnetic field on transformation and disintegration of the Earth’s bow shock into a system of magnetohydrodynamic (MHD) shock waves, rotational discontinuities and rarefaction waves under the action of abrupt variations in the solar wind dynamic pressure is simulated in the three-dimensional non-plane-polarized formulation within the framework of the ideal magnetohydrodynamic model using the solution of the MHD Riemann problem of breakdown of an arbitrary discontinuity. This discontinuity arises when a contact discontinuity, on which the solar wind density increases or decreases suddenly and which travels together with the solar wind, impinges on the Earth’s bow shock and propagates along its surface. The interaction pattern is constructed in the quasisteady- state formulation as a mosaic of exact solutions obtained on computer using an original MHD Riemann solver. The wave flow patterns are found for all elements of the surface of the bow shock as functions of their latitude and longitude for various jumps in the density on the contact discontinuity and characteristics parameters of the solar wind and interplanetary magnetic field at the Earth’s orbit. It is found that when the solar wind dynamic pressure increases, a fast MHD shock wave, that first penetrates into the magnetosheath, is always formed. When the solar wind dynamic pressure decreases, the influence of the interplanetary magnetic field can lead to the development of the leading fast MHD shock wave in certain zones on the surface of the Earth’s bow shock. The solution obtained can be used to interpret measurements on spacecraft in the solar wind at the libration point and in the neighborhood of the Earth’s magnetosphere.     

Two non-probabilistic set-theoretical models for dynamic response and buckling failure measures of bars with unknown-but-bounded initial imperfections

This paper is concerned with the problem of comparison of two non-probabilistic set-theoretical models for dynamic response and buckling failure measures of bars with unknown-but-bounded initial imperfections. Two kinds of non-probabilistic set-theoretical models are convex models and interval analysis models. In convex models and interval analysis models, the uncertain quantities are considered to be unknown except that they belong to given sets in an appropriate vector space. In this case, all information about the dynamic response and buckling failure measures of bars is provided by the set of dynamic responses and buckling failure measures consistent with the constraints on the uncertain quantities. The dynamic response estimate is actually a set in appropriate response space rather than a single vector. The set estimate is the smallest calculable set which contains the uncertain dynamic response, but it is usually impractical to calculate this set. Two set estimate methods are developed which can calculate the time varying box or hyperrectangle, i.e. interval vector in the response space that contains the true dynamic response. Comparison between convex models and interval analysis models in mathematical proofs and numerical calculations shows that, under the condition of the outer enclosed ellipsoid from a hyperrectangle or an interval vector, the set dynamic response predicted by interval analysis models is smaller than that yielded by convex models; under the condition of the outer enclosed hyperrectangle or an interval vector from an ellipsoid, the dynamic response set calculated by convex models is smaller than that obtained by interval analysis models.     

Regularization modeling for LES of separated boundary layer flow

Regularization models for the turbulent stress tensor are applied to mixing and separated boundary layers. The Leray and the NS-α models in large-eddy simulation (LES) are compared to direct numerical simulation (DNS) and (dynamic) eddy-viscosity models. These regularization models are at least as accurate as the dynamic eddy-viscosity model, and can be derived from an underlying dynamic principle. This allows one to maintain central transport properties of the Navier-Stokes equations in the model and to extend systematically toward complex applications. The NS-α model accurately represents the small-scale variability, albeit at considerable resolution. The Leray model was found to be much more robust, allowing simulations at high Reynolds number. Leray simulations of a separated boundary layer are shown for the first time. The strongly localized transition to turbulence that arises under a blowing and suction region over a flat plate was captured accurately, quite comparable to the dynamic model. In contrast, results obtained with the Smagorinsky model, either with or without Van Driest damping, yield considerable errors, due to its excessive dissipation.     

多柔体系统动力学建模与优化研究进展

多柔体系统是由柔性部件和运动副组成的力学系统,在航空、航天、车辆、机械与兵器等众多工程领域具有广泛的应用前景, 其典型的代表包括柔性机械臂、直升机旋翼、卫星的可展开天线、太阳帆航天器等. 近年来,随着工程技术的发展,多柔体系统动力学问题日益突出,尤其是含变长度柔性部件的多柔体系统,不仅涉及其动力学 建模与计算,还涉及其动力学优化设计. 事实上,部件柔性对多柔体系统的动力学行为影响很大,直接影响到优化结果,因此需要发展基于多柔体系统动力学的优化设计方法. 本文首先阐述了多柔体系统动力学优化的研究背景及意义,简要回顾了多柔体系统动力学建模的3类方法:浮动坐标方法、几何 精确方法和绝对节点坐标方法,并介绍了含变长度柔性部件的多柔体系统动力学建模方法. 系统概述了多柔体系统动力学响应优化、动力学特性优化和动力学灵敏度分析3个方面的研究进展,并从尺寸优化、形状优化和 拓扑优化 3 个方面综述了多柔体系统部件优化的研究进展. 本文最后提出了在多柔体系统动力学优化研究中值得关注的若干问题.      

A generalized variational formulation for convective heat transfer

The Scope of this paper is to develop the basic equations for a variational formulation which can be used to solve problems related to convection and/or diffusion dominated flows. The formulation is based on the introduction of a generalized quantity defined as the hear displacement. The governing equation is expressed in terms of this quantity and a variational formulation is developed which leads to a system of equations similar in form to Lagrange's equations of mechanics. These equations can be used for obtaining approximate solutions, though they are of particular interest for application of the finite element method. As an example of the formulation two finite element models are derived for solving convectiondiffusion boundary value problems. The performance of the two models is investigated and numerical results are given for different cases of convection and diffusion with two types of boundary conditions. The applications of the developed formulations are not limited to convection-diffusion problems but can also be applied to other types of problems such as mass transfer, hydrodynamics and wave propagation.     

考虑人载的汽车多自由度3维动力学模型

在目前的汽车柔顺性分析中,通常使用的2/4/7自由度动力学模型一般都未考虑车中人载质量,而在有人载与无人载时,汽车的动力学模型是不同的,它们代表了不同的动力学系统。因此未考虑人载质量的汽车动力学模型不能全面、准确地给出汽车的动力学特性,不能准确地反映柔顺性所包含的车身垂直方向、俯仰与侧倾振动,以及车轮跳动的加速度等物理量的变化。本文提出了考虑车中人载质量-座椅附加系统的汽车振动分析的3维动力学模型,给出了系统相应的动力学方程,并以轿车为例,给出了不同模型的动力学特性与响应的仿真分析结果,说明了汽车在考虑了人载-座椅附加系统后的多自由度3维动力学模型的正确性与必要性。     

A Pedagogical Approach to the Thermodynamically Constrained Averaging Theory

The thermodynamically constrained averaging theory (TCAT) is an evolving approach for formulating macroscale models that are consistent with both microscale physics and thermodynamics. This consistency requires some mathematical complexity, which can be an impediment to understanding and efficient application of this model-building approach for the non-specialist. To aid understanding of the TCAT approach, a simplified model formulation approach is developed and used to show a more compact, but less general, formulation compared to the standard TCAT approach. This new simplified model formulation approach is applied to the case of binary species diffusion in a single-fluid-phase porous medium system, clearly showing a TCAT approach that is applicable to many other systems as well. Recent extensions to the TCAT approach that enable a priori parameter estimation, and approaches to leverage available TCAT modeling building results are also discussed.