不同力学模型的胶泥阻尼器比较与分析.docxVIP

? ? 不同力学模型的胶泥阻尼器比较与分析 ? ? 刘昕运，马吉胜，吴大林 (陆军工程大学，河北 石家庄 050003) 1 引言 胶泥阻尼器出现以来，被广泛用于建筑、车辆、军工等行业。胶泥阻尼器内部含有独特的粘弹性材料，这种材料属于高分子硅氧有机聚合物，具有弹性和粘性双重特性。相比于常见的液压阻尼器和橡胶阻尼器，胶泥阻尼器结构更简单、容量更大、寿命更长、吸能率更高等特点[1]。 近年来，众学者对胶泥阻尼器的工作机理和力学性能的研究层出不穷。早期的力学模型大多基于实验和经验。Pekcan[3]对单出杆阻尼器进行不同频率拉伸压缩实验，提出了阻尼力与速度的幂函数呈正比。而后Terenzi[4]对单出杆粘弹型阻尼器和双出杆粘滞型阻尼器进行了实验，并在Pekcan公式的基础上进行了分析。Meram[5]为聚氨酯粘弹性阻尼器建立了Maxwell模型，Kelvin模型和Oscillation模型，并和实验结果进行比较。后来出现了将环形间隙的流动简化为平板间缝隙流动的思想。贾久红[6]和王虹琴[8]用微分的思想分析缝隙中的流体，推导出阻尼力正比于速度的幂函数，且阻尼系数与冲击实验环境、结构尺寸、粘度相关，但待定系数只能由实验确定。刘洪亮[9]、徐忠四[10]、仝志辉[11]等人对孔隙-缝隙式阻尼器进行分析，将流体运动分解为压差流、剪切流和孔隙流，推导出阻尼力的具体表达式，其中阻尼系数可由阻尼器结构尺寸直接求得。齐梦晓[12]在此基础上用AMESim液压系统仿真软件对胶泥阻尼器进行分析。随着流体动力学的发展，CFD技术开始应用于胶泥阻尼器的分析。郭强[13]和丁行武[14]建立双出杆孔隙式阻尼器的湍流CFD模型，并与正弦运动实验数据进行比较。孙靖雅[15]建立双出杆缝隙式阻尼器的层流CFD模型，并与落锤冲击实验数据比较。狄长春[16]建立了胶泥反后坐装置的三维流固耦合模型，并用射击实验进行验证。Wael Elsaady[17]建立磁流变阻尼器的湍流CFD模型，并进行了磁场和流场的数值计算。 以上研究是目前出现的针对各型胶泥阻尼器的常见方法。其中速度相关力学模型经过较大简化，参数确定较困难，准确性有限。而计算流体力学模型多数为正弦运动数值模拟，未进行速度冲击的数值模拟。且研究针对的阻尼器的结构和类型各不相同，无法横向对比。本文将以同型胶泥阻尼器为研究对象，基于计算流体力学和多体系统动力学，应用三种不同力学模型，经过静压实验的模型验证后，对阻尼器进行速度冲击仿真计算。分析和总结各方法的优缺点，为胶泥阻尼器的相关研究提供理论参考。 2 理论原理 2.1 阻尼器工作原理 目前出现的胶泥阻尼器有多种类型，按活塞杆结构可以分为单出杆式、等径双出杆式、非等径双出杆式。按活塞结构可以分为缝隙式、孔隙式、缝隙-孔隙式。本文选取较有代表性的非等径双出杆缝隙式胶泥阻尼器作为研究对象。如图1所示，该阻尼器主要由缸体、活塞、活塞杆、限位块和胶泥材料组成。相对等径双出杆式，其拥有自动弹性恢复的能力；相对单出杆式，同等力学性能下的活塞杆更粗，结构强度更大。胶泥在注入阻尼器时会被赋予一定的初始压力，这和限位结构能够保证阻尼器不受较小外力的扰动。当活塞和活塞杆向压缩方向运动时，缸体和活塞杆之间的空间变小，胶泥被压缩。弹性势能被储存，同时胶泥会流经环形间隙，产生阻尼力；当活塞和活塞杆向复原方向运动时，弹性势能被释放，胶泥向反方向流经环形间隙，继续耗散能量。 图1 胶泥阻尼器原理图 2.2 流体粘弹性模型 胶泥材料因其分子构型的独特性，在力学上既表现出液体的粘性，还有固体的弹性，且粘度比一般流体要大得多。工程上一般将胶泥材料的力学性能简化为kelvin模型，如式(1)所示，阻抗力为弹性力Fx和阻尼力Fv的叠加。 F=Fx+Fv (1) 液体的可压缩性，本文均采用简化的Tait可压缩液体状态方程，如式(2)所示。液体的体积弹性模量用来衡量液体的可压缩性，体积弹性模量值越小液体越容易被压缩。在工程应用上体积弹性模量是液体的固有属性，由液体本身确定，其定义为在一定温度下，液体承受的压强改变单位大小时，液体体积的相对变化值的倒数，如式(3)所示 (2) (3) 式中：m为液体密度指数；E为体积模量，且E=E0+nΔp；E0为参考体积模量；ρ为液体密度；ρ0为液体参考密度；p为液体压力；V为液体体积。 由于该阻尼器的最大压缩率不到1.5%，弹性模型的液体密度指数m可以视为1。由式(2)可以推导出弹性力Fk的表达式(4) (4) 式中：ΔS为活塞两侧的面积差，也等于两端的活塞杆截面积差；p0为胶泥的初始压力；V0为胶泥的初始体积；x为阻尼器的压缩行程。 胶泥属于典型的非牛顿流体，不满足牛顿内摩擦定律，其粘度值会随着剪切力和剪切率的变化而变化。本粘度模型使用广义非牛顿流体的幂律本构方