物质点法的理论和应用

物质点法采用质点离散材料区域, 用背景网格计算空间导数和求解动量方程,避免了网格畸变和对流项处理, 兼具拉格朗日和欧拉算法的优势, 非常适合模拟涉及材料特大变形和断裂破碎的问题. 本文详细论述了物质点法在基本理论、算法和软件开发方面的进展, 包括广义插值物质点法、接触算法、自适应算法、并行算法、与其他算法的杂交和耦合等. 系统地总结了物质点法在超高速碰撞、冲击侵彻、爆炸、动态断裂、流固耦合、多尺度分析、颗粒材料流动和岩土失效等一系列涉及材料特大变形问题中的应用,展示了其相对于传统数值计算方法的优势.     

RTM充模过程数值模拟的隐式有限元算法

建立了基于欧拉方法描述树脂传递模塑(RTM)工艺充模过程的基本数学方程,并采用有限元隐式时间积分方法对基本方程进行了数值求解．编制了基于隐式有限元算法及传统有限元控制体算法的程序,通过具体算例比较了这两种算法的优缺点．与传统的有限元控制体法相比,该文提出的隐式有限元算法能节省计算时间,特别适合于单元、节点数目多的情况．隐式有限元算法是一种纯有限元方法,不需要使用控制体积技术,采用该算法计算出的流动前沿与时间步长无关。     

多体系统Lagrange方程数值算法的研究进展

Lagrange方法是建立多体系统动力学方程的普遍方法之一,其方程的形式为常微分方程组或微分 - 代数方程组,数值计算与数值分析是研究多体系统动力学特性的重要方法.本文简要介绍了多体系统动力学方程的第一、二类Lagrange方程和修正的Lagrange方程的基本形式及这些方程的正则形式,着重介绍了正则方程在数值计算中的特点,就多体系统Lagrange方程的隐式算法、辛算法和多体系统动力学特性的数值分析方法(包括数值仿真、Poincar'e映射和Lyapunov指数的计算方法)的研究现状进行了综述.     

含初缺陷裂纹损伤梁的冲击动力屈曲

由Hamilton原理导出考虑初始缺陷及横向剪切变形时裂纹梁的动力屈曲控制方程;应用断裂力学中常用的线弹簧模型将裂纹引入到屈曲控制方程中;基于B-R动力屈曲判断准则,采用数值方法求解了受轴向冲击载荷作用时裂纹梁的动力屈曲;对比讨论了不同冲击速度、初始几何缺陷大小以及分布形式等因素对梁冲击动力屈曲的影响。     

基于SPH法的爆炸螺栓冲击特性研究

由于爆炸螺栓解锁过程具有瞬时性、变形大等特点,同时接触爆炸过程中交界面处介质阻抗相差 较大,因此基于计及材料强度的SPH 方法,在控制方程的粒子近似过程中将其转化为粒子密度及质量对物理 量近似过程影响较小的形式,建立了爆炸螺栓三维SPH 数值模型,分析了爆炸螺栓解锁过程中的冲击特性。 同时,通过与LS-DYNA软件的计算结果进行对比,验证了本文所建立数值模型的准确性和可靠性。     

有限变形问题过程最优控制算法

在非增量算法的基础上,提出了用最优控制变分原理形成过程最优控制迭代求解的基本思路,并给出求解的基本控制方程。这一工作为有限变形力学问题的数值求解提供了一个新的处理方法。     

多体系统Lagrange方程数值算法的研究进展

Lagrange方法是建立多体系统动力学方程的普遍方法之一, 其方程的形式为常微分方程组或微分-代数方程组,数值计算与数 值分析是研究多体系统动力学特性的重要方法。本文简要介绍了多 体系统动 力学方程的第一、二类Lagrange方程和修正的Lagrange方 程的基本形式及这些方程的正则形式,着重介绍了正则方程在数值 计算中的特点,就多体系统Lagrange方程的隐式算法、辛算法和多 体系统动力学特性的数值分析方法(包括数值仿真、 Poincarè映射 和Lyapunov指数的计算方法)的研究现状进行了综述。     

弹塑性球形薄壳在冲击载荷作用下的动力分析

通过曲面弯曲的等度量变换，给出了受冲击球壳的变形模态；接着，分别假定材料力弹性或刚塑性，基地能量守恒，得到了壳本和撞击体在运动过程中控制方程；最后，对所得到的控制方程进行了数值求工与实验数据作了比较，发展二者具有较好的一致性。     

单晶有限变形的热力耦合计算模型

以弹性变形梯度作为基本变量,结合热力学理论构造了单晶有限变形的热、力耦合计算模型。该模型考虑了温度、变温速率以及塑性耗散等条件对单晶有限变形的影响,相对于传统的以弹性变形梯度为基本变量的晶体塑性模型,算法能够体现温度效应的影响。采用隐式的积分方法对建立的控制方程进行计算以保证求解过程的稳定。以1100Al单晶为例计算了不同升温、降温速率,以及不同应变率影响下的材料应力-应变的响应。结果表明,模型能较好地反映变温过程中,单晶各向异性性质的演化以及应力、应变之间关系的变化。     

岩体渗透结构类型的划分及其渗透特性研究

提出了一种计算二维有限变形弹塑性摩擦接触问题形状设计灵敏度的算法. 采用主动 集策略和mortar方法处理接触边线上的约束条件. 在mortar接触边线的切线和法线方向上 采用相同的名义罚函数,提出基于名义罚函数的移动摩擦锥算法来正则化接触约束条件,发 展了一种新的二维多体有限变形摩擦接触算法. 在此基础上, 通过将离散形式的摩擦接触问题 控制方程对形状设计变量微分,得到了该路径相关问题的直接微分法解析设计灵敏度 计算格式, 其节点位移灵敏度方程在每个增量步不用迭代、直接求解. 与国际上现有 的二维多体有限变形摩擦接触问题的解析设计灵敏度算法相比,本算法不需分 解为法向和切向推导,表达式较简洁,便于编程实现. 数值算例验证了算法的精度 和有效性.     

泡沫铝夹层结构抗冲击性能的近场动力学模拟分析

针对某光学舱所采用的泡沫铝夹层防护结构在破片冲击下的抗冲击性能问题,采用Monte-Carlo方法创建了泡沫铝结构的二维细观模型,在常规态型近场动力学理论中引入了Mises屈服准则和线性各向同性强化模型,建立了近场动力学塑性本构的数值计算框架。基于近场动力学计算程序模拟了低速冲击作用下泡沫铝夹层结构的塑性变形以及有机玻璃背板的裂纹扩展形态,分析了泡沫铝芯材孔隙率对该夹层结构抗冲击性能和损伤模式的影响规律。结果表明：泡沫铝夹层结构良好的塑性变形能力是其发挥缓冲与防护作用的主要因素,并且在一定范围内,泡沫铝芯材孔隙率越高,则夹层结构具有更好的抗冲击性能;当泡沫铝孔隙率从0.4提升到0.7时,泡沫铝对冲击物的动能吸收率从90%提高到99%;模拟结果与实验结果具有较好的一致性,验证了模拟结果的准确性和分析结论的有效性。通过数值模拟,预测了有机玻璃背板的裂纹扩展形态,发现提高泡沫铝的孔隙率能获得更好的防护效果。     

A variational differential quadrature solution to finite deformation problems of hyperelastic shell-type structures: a two-point formulation in Cartesian coordinates

A new numerical approach is presented to compute the large deformations of shell-type structures made of the Saint Venant-Kirchhoff and Neo-Hookean materials based on the seven-parameter shell theory. A work conjugate pair of the first Piola Kirchhoff stress tensor and deformation gradient tensor is considered for the stress and strain measures in the paper. Through introducing the displacement vector, the deformation gradient, and the stress tensor in the Cartesian coordinate system and by means of the chain rule for taking derivative of tensors, the difficulties in using the curvilinear coordinate system are bypassed. The variational differential quadrature (VDQ) method as a pointwise numerical method is also used to discretize the weak form of the governing equations. Being locking-free, the simple implementation, computational efficiency, and fast convergence rate are the main features of the proposed numerical approach. Some well-known benchmark problems are solved to assess the approach. The results indicate that it is capable of addressing the large deformation problems of elastic and hyperelastic shell-type structures efficiently.     

NONLINEAR DYNAMIC ANALYSIS OF FLEXIBLE MULTIBODY SYSTEM

The nonlinear dynamic equations of a multibody system composed of ?exible beams are derived by using the Lagrange multiplier method. The nonlinear Euler beam theory with inclusion of axial deformation e?ect is employed and its deformation ?eld is …     

Theory and model implementation for analyzing line structures subject to dynamic motions of large deformation and elongation using the absolute nodal coordinate formulation (ANCF) approach

This paper presents the theory development and numerical implementation of a new gradient-deficient-based ANCF (Absolute Nodal Coordinate Formulation) model applied to perform the nonlinear dynamic analysis of elastic line structures subject to large stretching and deformation. The derivations of model equations, introduced numerical approaches, and result validations are the focuses of this study. Different from the traditional rod theory for small stretching consideration, the present model implements the line structures’ large elongation concepts into both the control mechanisms of constitutive formulations and equations of motion. The effect of external hydrodynamic forces on structures is also included in the model formulations. Based on the conservation of energy, the line model developed in this study covers the variation in strain and takes a full account of the bending effect with large stretching. A finite-element-based implicit scheme according to a modified Newmark-beta method is employed to solve the assembled system equations with unknown variables of nodal position vectors, their tangential derivatives, and strains. Selected cases with dynamic motions, such as nonlinear oscillation of a compound pendulum, free falling of a horizontal elastic beam in air with two different settings of gravity, free falling of a submerged horizontal tether with and without an attached concentrated mass, and a submerged vertical tether with a prescribed translational motion, are simulated to verify the developed model by comparing the results with analytical solutions and published experimental data and numerical results. It is found the present ANCF model, as noticed with good matched results with analytical solutions, measurements and other published solutions, is demonstrated to be able to provide converged and reasonably accurate predictions on the responses of line structures subject to large dynamic motions.

     

柔性梁刚柔耦合碰撞动力学及变形传播研究

运用柔性多体系统刚柔耦合动力学理论,研究了作大范围回转运动柔性梁的碰撞动力学问题.考虑柔性梁的横向变形,以及横向变形引起的纵向缩短项即非线性耦合变形项.采用基于Hertz接触理论及非线性阻尼理论的非线性弹簧阻尼模型来求解碰撞过程中产生的碰撞力,运用第二类拉格朗日方程建立了系统的刚柔耦合碰撞动力学方程.编制仿真软件进行动力学仿真计算,得到了碰撞力和系统动力学响应,对比分析了不同动力学模型对系统动力学响应的影响.同时研究了碰撞导致的柔性梁横向变形传播的波动特性.     

Mathematical Models and Methods of the Mechanics of Stochastic Composites

The basic approaches used in mathematical models and general methods for solution of the equations of the mechanics of stochastic composites are generalized. They can be reduced to the stochastic equations of the theory of elasticity of a structurally inhomogeneous medium, to the equations of the theory of effective elastic moduli, to the equations of the theory of elastic mixtures, or to more general equations of the fourth order. The solution of the stochastic equations of the elastic theory for an arbitrary domain involves substantial mathematical difficulties and may be implemented only rather approximately. The construction of the equations of the theory of effective moduli is associated with the problem on the effective moduli of a stochastically inhomogeneous medium, which can be solved by the perturbation method, by the method of moments, or by the method of conditional moments. The latter method is most appropriate. It permits one to determine the effective moduli in a two-point approximation and nonlinear deformation properties. In the structure of equations, the theory of elastic mixtures is more general than the theory of effective moduli; however, since the state equations have not been strictly substantiated and the constants have not been correctly determined, theoretically or experimentally, this theory cannot be used for systematic designing composite structures. A new model of the nonuniform deformation of composites is more promising. It is constructed by performing strict mathematical transformations and averaging the output stochastic equations, all the constants being determined. In the zero approximation, the equations of the theory of effective moduli follow from this model, and, in the first approximation, fourth-order equations, which are more general than those of the theory of mixtures, follow from it     

无网格MPM法三维爆炸焊接数值模拟

物质点法MPM(Material Point Method)是无网格方法之一.它是在质点网格法(PIC)基础上发展而来的一种新数值方法,它利用了欧拉法和拉格朗日法两者的优点,计算物质点在冲击载荷下的应力和应交;通过物质点来跟踪材料体的变形和破损,而在整个计算过程中背景网格始终固定不变,避免了重新划分网格.本文应用MPM法计算三维爆炸焊接问题,在爆轰载荷作用下的飞板和基板的金属动态变形过程进行了三维数值模拟,并且对飞板的碰撞点速度和爆轰压力变化进行了计算分析.     

固体介质中SPH方法的拉伸不稳定性问题研究进展

光滑粒子流体动力学法(smoothed particle hydrodynamics, SPH)是一种基于核估计的无网格Lagrange数值方法.它用粒子方程离散流体动力学的连续方程, 既可以处理有限元难于处理的大变形和严重扭曲问题, 又可以处理有限差分法不易处理的自由边界和材料界面的问题, 在固体力学中的冲击、爆炸和裂纹模拟中具有广阔的发展前景.但是, 该算法的拉伸不稳定性(tensile instability)问题是它在固体力学领域中应用的最大障碍.对SPH稳定性分析表明, 算法不稳定性的条件仅与应力状态和核函数的2阶导数有关.目前, 应力点法(stress points)、Lagrange核函数法、人工应力法(artificialstress)、修正光滑粒子法(corrective smoothed particle method, CSPM)和守恒光滑法(conservativesmoothing)以及其他一些方法成功地改善了SPH的拉伸不稳定性, 但是每一种方法都不能彻底解决SPH的拉伸不稳定性问题.本文介绍了SPH法的方程和Von Neumann稳定性分析的思想, 以及国内外在这几个方面的研究成果及其最新进展, 同时指出目前研究中存在的问题和研究的方向.      

AN ANALYTICAL SOLUTION OF G. I. TAYLOR’S THEORY OF PLASTIC DEFORMATION IN IMPACT OF CYLINDRICAL PROJECTILES AND ITS IMPROVEMENT

The theory of plastic deformation in the impact of cylindrical projectiles on rigid targets was first introduced by G. I. Taylor(1948)^[1]. The importance of this theory lies in the fact that the dynamic yield strength of the materials can be determined from the measurement of the plastic deformation of flat-ended cylindrical projectiles. From the experimental results^[2] we find that the dynamic yield strength is independent of impact velocity, and that it is higher than the static yield strength in general, and several times higher than the static yield strength in certain cases. This gives an important foundation for the study of elastoplastic impact problems in general. However, it is well known that the complexity of differential equations in Taylor’s theory compelled us to use the troublesome numerical solution. In this paper, the analytical solution of all the equations in Taylor’s theory is given in parametrical form and the results are discussed in detail.In the latter part of this paper, the method of calculation of impulse of impact is improved by considering the processes of radial’ movement of materials. The analytical solution of the improved theory shows that it gives better agreement with the experimental results than that of original Taylor’s theory.