A nonlinear planar beam formulation with stretch and shear deformations under end forces and moments

A new nonlinear planar beam formulation with stretch and shear deformations is developed in this work to study equilibria of a beam under arbitrary end forces and moments. The slope angle and stretch strain of the centroid line, and shear strain of cross-sections, are chosen as dependent variables in this formulation, and end forces and moments can be either prescribed or resultant forces and moments due to constraints. Static equations of equilibria are derived from the principle of virtual work, which consist of one second-order ordinary differential equation and two algebraic equations. These equations are discretized using the finite difference method, and equilibria of the beam can be accurately calculated. For practical, geometrically nonlinear beam problems, stretch and shear strains are usually small, and a good approximate solution of the equations can be derived from the solution of the corresponding Euler–Bernoulli beam problem. The bending deformation of the beam is the only important one in a slender beam, and stretch and shear strains can be derived from it, which give a theoretical validation of the accuracy and applicability of the nonlinear Euler–Bernoulli beam formulation. Relations between end forces and moments and relative displacements of two ends of the beam can be easily calculated. This formulation is powerful in the study of buckling of beams with various boundary conditions under compression, and can be used to calculate post-buckling equilibria of beams. Higher-order buckling modes of a long slender beam that have complex configurations are also studied using this formulation.     

On electromagnetic forming processes in finitely strained solids: Theory and examples

The process of electromagnetic forming (EMF) is a high velocity manufacturing technique that uses electromagnetic (Lorentz) body forces to shape sheet metal parts. EMF holds several advantages over conventional forming techniques: speed, repeatability, one-sided tooling, and most importantly considerable ductility increase in several metals. Current modeling techniques for EMF processes are not based on coupled variational principles to simultaneously account for electromagnetic and mechanical effects. Typically, separate solutions to the electromagnetic (Maxwell) and motion (Newton) equations are combined in staggered or lock-step methods, sequentially solving the mechanical and electromagnetic problems.The present work addresses these issues by introducing a fully coupled Lagrangian (reference configuration) least-action variational principle, involving magnetic flux and electric potentials and the displacement field as independent variables. The corresponding Euler-Lagrange equations are Maxwell's and Newton's equations in the reference configuration, which are shown to coincide with their current configuration counterparts obtained independently by a direct approach. The general theory is subsequently simplified for EMF processes by considering the eddy current approximation.Next, an application is presented for axisymmetric EMF problems. It is shown that the proposed variational principle forms the basis of a variational integration numerical scheme that provides an efficient staggered solution algorithm. As an illustration a number of such processes are simulated, inspired by recent experiments of freely expanding uncoated and polyurea-coated aluminum tubes.     

A principle of virtual work for combined electrostatic and mechanical loading of materials

The equations governing mechanics and electrostatics are formulated for a system in which the material deformations and electrostatic polarizations are arbitrary. A mechanical/electrostatic energy balance is formulated for this situation in terms of the electric enthalpy, in which the electric potential and the electric field are the independent variables, and charge and electric displacement, respectively, are the conjugate thermodynamic forces. This energy statement is presented in the form of a principle of virtual work (PVW), in which external virtual work is equated to internal virtual work. The resulting expression involves an internal material virtual work in which (1) material polarization is work-conjugate to increments of electric field, and (2) a combination of Cauchy stress, Maxwell stress and a product of polarization and electric field is work-conjugate to increments of strain. This PVW is valid for all material types, including those that are conservative and those that are dissipative. Such a virtual work expression is the basis for a rigorous formulation of a finite element method for problems involving the deformation and electrostatic charging of materials, including electroactive polymers and switchable ferroelectrics. The internal virtual work expression is used to develop the structure of conservative constitutive laws governing, for example, electroactive elastomers and piezoelectric materials, thereby determining the form of the Maxwell or electrostatic stress. It is shown that the Maxwell or electrostatic stress has a form fully constrained by the constitutive law and cannot be chosen independently of it. The structure of constitutive laws for dissipative materials, such as viscoelastic electroactive polymers and switchable ferroelectrics, is similarly determined, and it is shown that the Maxwell or electrostatic stress for these materials is identical to that for a material having the same conservative response when the dissipative processes in the material are shut off. The form of the internal virtual work is used further to develop the structure of dissipative constitutive laws controlled by rearrangement of material internal variables.     

轴向运动导电薄板磁弹性耦合动力学理论模型

针对磁场环境中轴向运动导电薄板的动力学理论建模问题进行研究,得到较为完备的磁弹性耦合振动基本方程及相应的补充关系式。在考虑几何非线性效应下,给出薄板运动的动能、应变能以及外力虚功的表达式。应用哈密顿变分原理,推得磁场中轴向运动薄板的非线性磁弹性耦合振动方程,并得到力和位移满足的边界条件。基于麦克斯威尔电磁场方程,并考虑相应的电磁本构关系和电磁边界条件,推得任意磁场环境中轴向运动导电薄板满足的电动力学方程和所受电磁力表达式。分别针对纵向磁场环境、横向磁场环境、条形板等具体情形,给出了振动方程、电动力学方程和电磁力的简化形式。所得结果,可为此类问题的进一步求解和分析提供理论参考。     

On virtual work and stress in granular media

A granular medium can be treated as an equivalent continuum. Appropriate representative stresses can be derived from the virtual work principle. However, the expression of virtual work is not unique and therefore may lead to different results of stress expressions in terms of discrete quantities—contact forces, contact moments, and branch vectors. In this paper, we introduced a generalized expression of virtual work that includes the restriction of boundary conditions. To show the advantages of the current expression, the virtual work expression is applied to derive expressions for stress, couple stress, a higher-order stress, and the stress moment. A distinction is made between the average stress within a granular volume and the representative stress that is conjugate with the representative strain of the volume. The current work is compared with that of [International Journal of Solids and Structures 38 (2) (2001) 353–367], and the current stress expressions are shown to satisfy three essential conditions of a stress measure.     

耗散系统的虚功变分不等式及其应用

对于保守系统,能量变分原理为推导力学系统控制方程提供了简洁的途径.对于耗散系统,控制方程的建立往往需要引入经验的或理性的假定,增大了建模的难度.针对耗散系统,引入系统局部稳定的概念,并在此基础上,提出一类虚功变分不等式.这一不等式事实上揭示了耗散系统的一类虚功不等原理.该原理的物理含义为:使系统状态稳定的必要条件是,在该状态附近所有可能的虚拟路径上系统释放的势能不大于系统耗散的能量.研究表明:仅需结合虚功不等原理和能量守恒原理,即可导出准静态系统力学状态量的全部控制方程.作为应用,文章重新讨论了塑性力学,结合虚功不等原理与能量守恒原理,导出经典塑性力学的全部控制方程,并证明了经典的最大塑性耗散原理可以作为虚功不等原理的推论导出;同时,以Mohr-Coulomb强度准则为例,讨论了虚功不等原理在强度理论中的应用,说明基于应力的强度准则可以是基于能量的稳定性准则的推论.上述例子说明了虚功不等原理的广泛适用性和在建立耗散系统控制方程中的有效性.     

Modeling of fiber-reinforced cement composites: Discrete representation of fiber pullout

The discrete modeling of individual fibers in cement-based materials provides several advantages, including the ability to simulate the effects of fiber dispersion on pre- and post-cracking composite performance. Recent efforts in this direction have sought a balance between accurate representation of fiber behavior and computational expense. This paper describes a computationally efficient approach to representing individual fibers, and their composite behavior, within lattice models of cement-based materials. Distinguishing features of this semi-discrete approach include: (1) fibers can be positioned freely in the computational domain, irrespective of the background lattice representing the matrix phase; (2) the pre- and post-cracking actions of the fibers are simulated with little computational expense, since the number of system degrees of freedom is independent of fiber count. Simulated pullouts of single fibers are compared with theory and test results for the cases of perfectly-plastic and slip-hardening behavior of the fiber–matrix interface. To achieve objective results with respect to discretization of the matrix, pullout forces are distributed along the embedded lengths of fibers that bridge a developing crack. This is in contrast to models that lump the pullout force at the crack surfaces, which can lead to spurious break-off of matrix particles as the discretization of the matrix is refined. With respect to fracture in multi-fiber composites, the proposed model matches theoretical predictions of post-cracking strength and pullout displacement corresponding to the load-free condition. The work presented herein is a significant step toward the modeling of strain-hardening composites that exhibit multiple cracking.     

Introducing unsteady non‐uniform source terms into the lattice Boltzmann model

Taking body forces into account is not new for the lattice Boltzmann method, yet most of the existing approaches can only treat steady and uniform body forces. To manage situations with time‐ and space‐dependent body forces or source terms, this paper proposes a new approach through theoretical derivation and numerical verification. The method by attaching an extra term to the lattice Boltzmann equation is still used, but the expression of the extra term is modified. It is the modified extra term that achieves the particularity of the new approach. This approach can not only introduce unsteady and non‐uniform body forces into momentum equations, but is also able to add an arbitrary source term to the continuity equation. Both the macroscopic equations from multi‐scale analysis and the simulated results of typical examples show that the accuracy with second‐order convergence can be guaranteed within incompressible limit. Copyright © 2007 John Wiley & Sons, Ltd.     

Inverse Dynamics of Servo-Constraints Based on the Generalized Inverse

The acceleration form of constraint equations is utilized in this paper to solve for the inverse dynamics of servo-constraints. A condition for the existence of control forces that enforce servo-constraints is derived. For overactuated dynamical systems, the generalized Moore-Penrose inverse of the constraint matrix is used to parameterize the solutions for these control forces in terms of free parameters that can be chosen to satisfy certain requirements or optimize certain criterions. In particular, these free parameters can be chosen to minimize the Gibbsian (i.e., the acceleration energy of the dynamical system), resulting in ideal control forces (those satisfying the principle of virtual work when the virtual displacements satisfy the servo-constraint equations). To achieve this, the nonminimal nonholonomic form recently derived by the authors in the context of Kanes method is used to determine the accelerations of the system, and hence to determine the forces to be generated by the redundant manipulators. Finally, an extension to inverse dynamics of servo-constraints involving control variables is made. The procedures are illustrated by two examples.     

Configurational forces in discrete elastic systems

The concept of configurational forces is applied to a simple, one-dimensional problem that is solved by finite elements. Both the exact solution and its finite-element approximation are provided in closed form. The total energy according to the approximate solution depends on the choice of the nodes. Any virtual shift of a node results in a virtual change of energy, which can be interpreted as the virtual work done by a configurational force acting on that node. It is shown that, in the case of equidistant nodes, the configurational forces acting on the interior nodes vanish. Also, the relation between the nodal configurational forces and the Eshelby stress resultant along the rod is investigated.     

MAGNETO-THERMO-ELASTIC COUPLED DYNAMICS EQUATIONS OF AXIALLY MOVING CARRY CURRENT PLATE IN MAGNETIC FIELD

针对磁场环境中轴向运动载电流导电板磁热弹性耦合动力学建模问题进行研究. 考虑几何非线性和热效应条件下, 给出薄板运动的动能、应变能以及外力虚功的表达式.应用哈密顿变分原 理, 推得力、运动、电、磁和热效应相互作用下轴向运动导电板的非线性磁热弹性耦合振动方程.基于麦克斯韦电磁场方程, 考虑相应的电磁本构关系和电磁边界条件, 推得磁场环境中轴向运动载电流导电板满足的电动力学方程和所受电磁力表达式, 并给出焦耳热作用下耦合形式的热传导方程. 算例表明, 磁场等参量对动力学系统分岔特性有明显影响.所得结果可为此类问题的进一步求解和分 析提供理论参考.     

The inverse identification problem and its technical application

This paper presents an overview of a loading force identification technique. Load identification methods are based on the solution of the inverse identification problem. Many different approaches for linear systems have been developed in this area. For both linear and nonlinear systems, methods based on the minimization of assumed objective functions are formulated. The least square error between the simulated and measured system responses is mainly used as the objective function. The dynamic programming optimization method formulated by Bellman is commonly used for the minimization of the objective function to estimate the excitation forces. The inverse identification problem in most practical cases is ill-posed because not all the state variables or initial conditions are known. Ill-posed inverse identification problems can be solved using several techniques, the most useful of which are: the generalized cross-validation method, the dynamic programming technique and Tikhonov’s method. This article presents the theoretical background and main limits to the application of inverse identification methods. Numerical and experimental tests on a laboratory rig were made to verify the formulated procedures. The method is applied to the identification of wheel–rail contact forces during rail vehicle operation. The method can be applied for indirect measurements of contact forces in railway equipment testing.     

The role of inter-particle collisions on elbow erosion

Erosion due to particle impingement, which can occur in a variety of practical cases, is often the key factor in pipeline failure. Parts such as elbows, for instance, are particularly prone to erosion issues. In this work, the Unsteady Reynolds Averaged Navier–Stokes (URANS) equations are combined with a stochastic Lagrangian particle tracking scheme considering all relevant elementary processes (drag and lift forces, particle rotation, inter-particle collisions, particle-wall interactions, coupling between phases) to numerically predict the erosion phenomenon on a 90° elbow pipe. After a detailed validation of the erosion model based on the experimental data of Solnordal et al. (2015), several cases regarding the wall roughness and static and dynamic coefficients of friction are analysed to elucidate the nature of the erosive process. For such analysis, more fundamental variables related to particle-wall interactions (impact velocity, impact angle, impact frequency) were used to scrutinize the basic erosion mechanisms. Finally, to prove the importance of inter-particle collision on elbow erosion, different mass loadings are additionally simulated. Especially for the high mass loading cases, interesting results about the role of the inter-particle collisions on elbow erosion are enlightened.     

An investigation on the external inversion of thin-walled tubes using a die

The production of sound thin-walled tubular parts by external inversion using a die is normally limited to components having geometrical features within a compact range. Recent published works in the field by Sun and Yang [Int. J. Mach. Tools Manuf. 42 (2002) 15] and Sekhon et al. [J. Mater. Proc. Technol. 133 (2003) 24] present new results relevant to process and die design. However, gaps of knowledge can still be found in understanding the influence of interface friction on the material flow and the effect of strain path and material damage on the occurrence of fracture. The aim of this paper is to contribute to an extension of the actual knowledge of the process by means of a comprehensive theoretical and experimental investigation. Theoretical investigation is accomplished by using virtual prototyping modeling techniques based on the finite element method. Experimental work is performed on Aluminum Al6060 industrial tubes and it is mainly utilized for supporting and validating the theoretical investigation.     

基于互补理论和扩展有限元法的摩擦接触裂纹应力强度因子计算

从虚功方程出发,结合扩展有限元离散技术与接触条件的非线性互补表述,建立了摩擦接触裂纹问题的扩展有限元非线性互补模型,将不等式接触条件转化为非线性互补类的非光滑方程组,并采用基于广义导数的非光滑阻尼牛顿法求解方程组,无需引入任何额外人工变量以及迭代求解。以含中心倾斜裂纹平板和边裂纹平板为例,运用相互作用积分法计算摩擦接触裂纹的应力强度因子,将其结果与理论解进行对比分析,该方法都能给出精确的计算结果;基于扩展有限元方法对单轴压缩作用下倾斜裂纹扩展过程进行了数值模拟,计算结果表明,受压裂纹数值结果与实验结果比较吻合,从而验证了本文方法的有效性与正确性。     

滑轮绳索系统中动态节点绳索单元

解除了传统有限元方法中单元节点与物质点固结的假设, 建立了空间点的速度和加速度与相应物质点的速度和加速度之间的数学关系, 强调了虚功率原理中出现的速度和加速度皆为物质速度和物质加速度. 在此基础上构造了单元节点既不与空间坐标固定也不与物质坐标固定的动态节点绳索单元. 根据滑轮绳索系统的特点, 选取绳索出入绳点的弧长坐标、出入绳角、面外摆角以及拉伸应变等空间参数描述了滑轮绳索系统的运动. 将绳索与滑轮以及绳索与卷筒之间的相互作用合理简化为物质速度间的约束条件, 避免了传统方法中接触力计算不收敛、效率低等缺点. 所提方法可精细求解绳索与滑轮间接触边界点位置和卷筒入绳点在卷筒上的运动、滑轮的中心和其连体基的运动、绳索出入滑轮和卷筒时空间方位的变化以及绳索上任意点的拉力变化等细节. 可为含绳索机械系统的力学分析提供新的理论基础. 所用的解除单元节点与物质点绑定的理论具有一定的普适性, 可为有限元方法的理论和应用提供参考.      

MULTIPLE TIME SCALES OF AEOLIAN AND FLUVIAL PROCESSES AND DEPTH-AVERAGED INTEGRAL MODELLING

针对风沙与水沙动力学研究的共性,提出风沙与水沙运动统一深度积分模式;在此基础上,将水沙运动过程多重时间尺度理论扩展至风沙运动,比较研究风沙与水沙运动向平衡状态调整的时间尺度特征. 数值算例与风洞实验结果对比表明,积分平均模型能够基本准确地捕捉风沙运动的输沙率变化等基本特征;风沙与水沙运动的积分时间尺度比较分析表明,风沙运动与水沙运动类似,推移质运动能够很快地调整到平衡状态,而悬移质运动调整到平衡状态则需要相对较长的过程,在湍流悬移质实验和数值模拟研究中应考虑恢复平衡过程的影响.     

多目标优化算法在弹头侵彻明胶运动模拟中的应用

为使弹道模型更好地模拟弹头侵彻明胶的运动规律,从而有助于揭示弹头对人体组织的致伤机理,以弹头水平位移、竖直位移、侧向位移、俯仰角以及偏航角理论值的均方根误差作为目标函数组,以弹头的待定初始运动参数和弹道模型的待定力学系数作为优化变量,通过Gamultiobj多目标遗传算法获得了的优化变量的Pareto解集,并通过TOPSIS综合分析方法获得了Pareto解集的最优解。研究结果表明该方法能够快速获得可信性较强的最优解,使最优解能较好地模拟7.62 mm步枪弹侵彻明胶的运动过程。

     

Lubricated revolute joints in rigid multibody systems

The main purpose of this work is to present a general methodology for modeling lubricated revolute joints in constrained rigid multibody systems. In the dynamic analysis of journal-bearings, the hydrodynamic forces, which include both squeeze and wedge effects, generated by the lubricant fluid, oppose the journal motion. The hydrodynamic forces are obtained by integrating the pressure distribution evaluated with the aid of Reynolds’ equation, written for the dynamic regime. The hydrodynamic forces built up by the lubricant fluid are evaluated from the system state variables and included into the equations of motion of the multibody system. Numerical examples are presented in order to demonstrate the use of the methodologies and procedures described in this work.