基于多弹簧模型的空间梁柱单元(Ⅱ):参数分析

为模拟地震作用时杆件在变动轴力和双向弯曲的耦合特点,基于简化纤维模型的多弹簧单元模型近年被广泛应用.本文作者在相关论文中提出了考虑剪切变形弹塑性刚度影响的多弹簧模型空间梁柱单元,杆件的计算模拟精度除了取决于单元模型之外,还决定于两个重要参数:弹簧位置及个数布置、塑性区长度,本文将针对其进行参数分析并提出合理建议取值.本文最后还提出了模拟杆件残余应力的合理弹簧布置方式.     

用于正交各向异性织物仿真的弹簧质点模型

针对正交各向异性织物,本文提出了一种含附加节点的三角形弹簧质点模型．通过冗余自由度的引入和缩聚,能够实现有限元模型和弹簧质点模型三角形单元刚度矩阵的精确相等;据此,针对一般正交各向异性材料常数和单元形状,推导出弹簧质点模型中弹簧刚度系数及附加节点位置的解析表达式,有效解决了传统弹簧质点模型在正交各向异性织物仿真中存在的计算精度不高、弹簧参数选取随意等缺陷;并进一步完善了基于弹簧质点模型的柔性织物空间变形算法．数值模拟表明本文提出的弹簧质点模型和变形算法在进行织物仿真计算时,具有精度好、通用性强等优点．     

模拟火炮后坐运动的非线性弹簧－阻尼单元

根据非线性有限元法的基本理论，建立了火炮发射动力学的有限元运动方程，推导出单元耗散为向量和单元阻尼矩阵的计算公式。提出采用非线性弹簧－阻尼单元来模拟反后坐装置对炮身后坐的阻尼作用，并给出了该单元的刚度矩阵和阻尼矩阵。根据炮身的后坐过程的三个阶段，提出采用分段拟合的方法确定非线性单元的参数以模拟火炮后坐阻力－位移曲线。采用非线性有限元单元模拟后坐装置，能够较准确地描述炮身部分的后坐运动，同时在大简化了全炮模型的建立，为全炮的动力学仿真提出了一条新的思路。     

三角形截面梁柱构件的三维等效弹簧计算模型

传统有限元在分析梁柱构件时一般采用常应变单元和双线性单元,但此类单元在梁柱构件受弯分析中计算精度不是很高.本文根据梁柱构件的力学性能,在三维连续介质体受到轴向变形和弯曲变形的状态下,利用其轴向变形和弯曲挠度相同,得到具有相同宽度和高度且刚度等效的超静定桁架力学模型.然后,通过桁架杆的截面参数求得弹簧的刚度系数,从而得等效弹簧元模型.本文提出的等效弹簧模型计算方法简单,便于扩展到更为复杂的构件分析中.     

完备化ANCF薄板单元及在钢板弹簧动力学建模中的应用

绝对节点坐标方法已在多体系统动力学研究中广泛应用, 但常用来描述板壳类结构的薄板单元, 由于梯度不完备而无法直接用于带有初始弯曲参考构型的柔性体变形描述. 为避免全参数板单元建立车辆钢板弹簧模型时存在的严重截面闭锁问题, 拟采用薄板单元用于板簧建模. 为此, 探索了将现有绝对节点坐标薄板单元纳入一般连续介质力学弹性力表达的方法, 采用中面上单位法向量作为单元厚度方向的梯度向量, 从而得到了完备化的薄板单元及其描述初始弯曲构型时消除初应变的方法. 在此基础上通过定义簧片的未变形构型, 在钢板弹簧中引入可控的预应力, 实现对钢板弹簧装配过程的准确模拟. 通过数值算例验证了本方法的正确性. 建立了车辆钢板弹簧模型, 通过建立在簧片上的局部坐标系实现接触点的跨单元搜索, 并采用惩罚函数法和平滑化的库伦摩擦模型施加簧片间的接触力. 引入参考节点的概念建立了整合车身与吊耳及其机构运动关系的刚柔耦合模型.}}      

基于单弹簧联结单元法的钢筋混凝土梁损伤模拟

钢筋混凝土梁的损伤模拟的重难点在于描述两种材料非线性的本构关系以及钢筋混凝土粘结滑移关系非线性的描述,采用单弹簧联结单元法模拟钢筋和混凝土相互作用,结合混凝土四参数损伤模型,模拟钢筋混凝土结构多轴应力状态下从开始加载到结构破坏的全过程。数值算例不仅体现该方法的有效性,也具有重要的工程实践意义。     

基于多弹簧模型的空间梁柱单元(I):理论模型

引入多轴应力状态下的塑性应力-应变关系理论,在单元模型中考虑了弹塑性区域剪切变形对单元的弹塑性刚度影响,推导了考虑剪切变形弹塑性刚度影响的多弹簧模型的空间梁柱单元刚度矩阵,为整体钢框架结构的弹塑性动力分析奠定了单元基础。     

基于多弹簧模型的空间梁柱单元(Ⅰ):理论模型

引入多轴应力状态下的塑性应力-应变关系理论,在单元模型中考虑了弹塑性区域剪切变形对单元的弹塑性刚度影响,推导了考虑剪切变形弹塑性刚度影响的多弹簧模型的空间粱柱单元刚度矩阵,为整体钢框架结构的弹塑性动力分析奠定了单元基础.     

金属橡胶材料基于微弹簧组合变形的细观本构模型

研究了金属橡胶材料细观本构模型。依据金属橡胶材料变形的主要特征,分析了微弹簧沿轴向和径向变形的规律,分别推导了应力应变关系的解析表达式,由两种微弹簧的组和变形构造了代表性体积单元-微元体。为了反映因铺层引起的工艺各向异性,本构方程中引入了一个新的材料参数铺层比例系数iβ,该参数可以较好地表征细观结构中微弹簧取向的分布状况。在大量试验的基础上,提出了确定铺层系数的实验方法,进而建立了包含金属橡胶材料细观结构信息的本构方程。理论与实验结果比较表明,本文建立的本构模型能较好的反映材料的力学行为。     

基于离散元的杆件断裂全过程数值模拟

杆件的断裂会涉及到大变形、非线性以及不连续等问题,通常的数值计算方法模拟这种复杂力学行为具有局限性。本文基于颗粒离散元法DEM,将接触粘结处的分布式弹簧用梁纤维进行等效,提出了一种适于结构弹塑性分析的DEM纤维梁模型,然后在此基础上构建了构件断裂模拟算法以及纤维破环准则。将该模型应用于悬臂梁结构,模拟了悬臂梁从弹性到弹塑性阶段,再到断裂破坏的全过程,数值模拟得到的结构响应和截面开裂破坏形态均较合理。最后将该方法应用于单层网壳倒塌破坏模拟,并与网壳振动台倒塌试验进行对比,结果表明,数值模拟得到的杆件断裂过程及结构倒塌模式与试验现象一致,验证了该模型的正确性和适用性。     

等径通道挤压过程三维有限元模拟

利用三维有限元模型对等径通道挤压过程中变形分布的不均匀性进行了模拟,通过对不同摩擦系数下截面等效塑性应变分布比较发现:沿三个方向截面上的变形分布都不是均匀的,这说明二维有限元模型不能真实模拟等径通道挤压过程中试样中的变形分布。此外,摩擦对等效塑性应变分布及挤压力的影响较大,截面变形不均匀参数和最大等效塑性应变出现的位置随摩擦系数的增大而变化。压力-位移曲线的变化可以结合试样的变形过程来解释。     

Semi-active hydropneumatic spring and damper system

This paper is concerned with the development of a semi-active hydropneumatic spring and damper system, comprising of a two state hydropneumatic spring and a two state hydraulic damper. The system was specifically developed to improve the ride comfort and handling of large off-road vehicles. The suspension requirements for good ride comfort and handling for heavy off-road vehicles are discussed with special reference to the advantages of semi-active hydropneumatic springs and semi-active dampers. The layout and functioning of an experimental spring and damper unit used for laboratory tests are discussed. Spring and damper characteristics, as well as valve response times for both the semi-active spring and semi-active damper were determined. A single degree of freedom test rig with a sprung mass of 3 tons was used to perform first order ride comfort tests. Tests include step response and random input response. The test rig was also used to evaluate semi-active control strategies for both spring and damper as well as a control strategy for implementing ride height adjustment without using an external hydraulic pump.     

Three dimensional combined fracture–plastic material model for concrete

This paper describes a combined fracture–plastic model for concrete. Tension is handled by a fracture model, based on the classical orthotropic smeared crack formulation and the crack band approach. It employs the Rankine failure criterion, exponential softening, and it can be used as a rotated or a fixed crack model. The plasticity model for concrete in compression is based on the Menétrey–Willam failure surface, the plastic volumetric strain as a hardening/softening parameter and a non-associated flow rule based on a nonlinear plastic potential function. Both models use a return-mapping algorithm for the integration of constitutive equations. Special attention is given to the development of an algorithm for the combination of the two models. The suggested combination algorithm is based on a recursive substitution, and it allows for the two models to be developed and formulated separately. The algorithm can handle cases when failure surfaces of both models are active, but also when physical changes such as crack closure occur. The model can be used to simulate concrete cracking, crushing under high confinement and crack closure due to crushing in other material directions. The model is integrated in a general finite element package ATENA and its performance is evaluated by comparisons with various experimental results from the literature.     

A theoretical model of collision between soft-spheres with Hertz elastic loading and nonlinear plastic unloading

This paper presents a theoretical model on the normal (head-on) collision between soft-spheres on the basis of elastic loading of the Hertz contact for compression process and a nonlinear plastic unloading for restitution one, in which the parameters all are determined in terms of the material and geometric ones of the spheres, and the behaviors of perfect elastic, inelastic, and perfect plastic collisions appeared in the classical mechanics are fully described once a value of coefficient of restitution is specified in the region of 0 ≤ ε ≤ 1. After an empirical formula of the coefficient of restitution dependent on the impact velocity is suggested to fit the existing experimental measurements by means of the least square method, the predictions of the dependency and the collision duration are in well quantitative agreement with their experimental measurements. It is found that the measurable quantities are dependent on both the impact velocity and the parameters of spheres. Following this model, finally, an approach to determine the spring coefficient in the linear viscoelastic model of the collision is also displayed. These results obtained here will be significantly beneficial for the applications where a collision model is requested in the simulations of relevant grain flows and impact dynamics etc..     

含裂纹梁刚塑性动态断裂的非线性线弹簧模型分析

用非线性线弹簧模型分析了带裂纹梁的刚塑性动态断裂问题.在塑性势理论基础上,建立了全塑性状态下的弹簧本构关系,并用此关系导出带裂纹梁刚塑性动态断裂分析的基本方程,计算了在冲击载荷作用下,裂纹梁的动态断裂响应.     

塑性各向异性板料本构关系及冲压成形的数值模拟

推导了具有一般屈服函数形式的弹塑性速率型本构关系;给出了用于板料成形的Hill塑性各向异性屈服模型下本构关系的具体形式;用有限元动力显式计算程序MSC/DYTRAN模拟了金属板料的冲压成形;通过算例分析,考察了塑性各向异性对凸耳形成和大小以及对成形模拟结果准确性的影响;数值结果和实验结果表明：各向（厚向）异性本构模型比各向同性本构模型更真实地反映了板料的成形性。     

板材多点成形过程的有限元分析

多点成形过程采用静力隐式格式进行数值模拟是比较合适的。本文建立了用于多点成形过程分析的静力隐式弹塑性大变形有限元方法 ,给出了对稳定迭代收敛过程效果较好的板壳有限单元模型、处理多点不连续接触边界的接触单元方法以及增量变形过程中应力及塑性应变计算的多步回映计算方法。基于这些方法编制了计算软件 ,应用该软件进行了矩形板的液压胀形过程及球形模具拉伸成形过程的有限元分析 ,数值计算结果与典型的实验结果及计算结果吻合很好。最后给出了球形、圆柱形目标形状的实际多点成形过程的数值模拟结果。     

On plastic anisotropy of constitutive model for rate-dependent single crystal

An algorithm for single crystals was developed and implemented to simulate plastic anisotropy using a rate-dependent slip model. The proposed procedure was a slightly modified form of single crystal constitutive model of Sarma and Zacharia. Modified Euler method, together with Newton-Raphson method was used to integrate this equation which was stable and efficient. The model together with the developed algorithm was used to study three problems. First, plastic anisotropy was examined by simulating the crystal deformation in tension and plane strain compression, respectively. Secondly, the orientation effect of some material parameters in the model and applied strain rate on plastic anisotropy for single crystal also is investigated. Thirdly, the influence of loading direction on the active slip system was discussed.