一种利用虚拟材料等效栓接结合部动态特性的方法技术

本发明专利技术公开一种利用虚拟材料等效栓接结合部动态特性的方法，将栓接结合部看作多种虚拟材料层，利用有限元法模拟栓接结合部在不同螺栓预紧力下，不同接触区域的受力情况，将栓接结合部的接触压力分布等效为关于接触半径的4阶多项式函数；利用材料力学知识建立栓接结合部等效虚拟材料模型，确定有效接触区域内接触半径与虚拟材料属性(弹性模量、剪切模量、泊松比、密度、虚拟层厚度)间的函数关系；通过有限元法将二层、三层或者多层虚拟材料与各子结构耦合建立整体结构的有限元模型，并进行动态响应分析，将获得的频响函数与模态试验数据比较，验证所建模型的正确性，为研究栓接结合部动态特性提供新的思路和方法。

A method of using the equivalent bolt of virtual material to connect the dynamic characteristics of the joint

The invention discloses a method for using the virtual equivalent dynamic characteristics of bolted joints, the bolted joints as multiple virtual material layer, simulation of bolted joints in different bolt pretightening force using the finite element method, stress condition of different contact area, the contact pressure of the bolt joint distribution of equivalent 4 order polynomial function of the contact radius; according to the theory of material mechanics of bolted joints virtual material model, determine the effective contact area and contact radius of virtual material properties (elastic modulus, shear modulus, Poisson's ratio, density, thickness of the virtual function relationship between); finite element model by finite element method the two layer, three layer or multi-layer virtual material and sub structure to establish the overall structure, and dynamic response analysis, the frequency response function and modal test data than The correctness of the model is verified, which provides a new idea and method for the study of the dynamic characteristics of the bolt joint joint.

【技术实现步骤摘要】
一种利用虚拟材料等效栓接结合部动态特性的方法
本专利技术涉及结合面理论研究领域中动力学特性领域，尤其涉及一种利用虚拟材料法等效栓接结合部动力学特性的方法。
技术介绍
数控机床中存在大量结合部，使得结构的连续性被破坏，从而表现出复杂的动力学特性，极大程度地影响了整机动态性能。研究表明，机床中结合部的接触刚度约占机床总刚度的60％～80％，结合部的接触阻尼占机床全部阻尼的90％以上，结合部引起的变形量约占机床总变形量的40％～60％。因此，开展结合部动态特性研究是机械系统由单件分析走向整机分析的关键，是整机性能预测的基础。通常采用弹簧-阻尼单元模拟结合部的动态特性，但是利用有限元仿真时难以将刚度和阻尼模型与各子结构之间进行耦合，整体结构的动态特性难以通过仿真技术实现。利用虚拟材料模拟栓接结合部的动态特性可以大大降低利用有限元仿真时带来的困难，但是目前的研究是将栓接结合部等效为一层虚拟材料，难以提高整体结构的仿真精度。
技术实现思路
本专利技术的目的在于提供一种利用虚拟材料等效栓接结合部动态特性的方法，以解决上述技术问题。本专利技术利用多层虚拟材料模拟栓接结合部的动态特性，降低栓接结合部与各子结构间的耦合难度，提高整体结构的仿真精度。为了实现上述目的，本专利技术采用如下技术方案：一种利用虚拟材料等效栓接结合部动态特性的方法，包括以下步骤：步骤一、确定栓接结合部接触压应力：假设垂直栓接结合部方向上压应力均匀分布，大小随螺栓作用力的影响半径变化而变化，且压应力分布呈截头圆锥体，利用有限元法模拟栓接结合部在不同螺栓预紧力下，不同接触区域的受力情况，并将接触压力等效为关于接触半径的4阶多项式函数；步骤二、确定栓接结合部虚拟材料的相关参数：建立栓接结合部等效虚拟材料模型，确定有效接触区域内接触半径与虚拟材料属性间的函数关系；虚拟材料属性包括弹性模量E、剪切模量G、泊松比υ、密度和虚拟层厚度；其中虚拟材料的弹性模量E、剪切模量G、泊松比υ均是关于螺栓作用力的影响半径的函数，而平均密度和虚拟材料厚度与接触半径无关，为确定值；步骤三、有限元法确定整体结构的频率响应：根据虚拟材料模型和接触压力的变化曲线，将栓接结合部等效为二层、三层或者多层虚拟材料，利用子结构耦合法建立整体结构的有限元模型，并进行动态响应分析，获得栓接结合部整体结构的频响函数。进一步的，垂直栓接结合部方向上压应力大小随螺栓作用力的影响半径变化而变化，且压应力分布呈截头圆锥体，将栓接结合部接触压力Pn(r)等效为关于(r/a)的4阶多项式函数：式中，r表示螺栓作用力的影响半径，a表示螺栓孔半径，a1，a2，a3和a4为4阶多项式函数的相关系数；根据边界条件：由于，螺栓接触载荷为：式中，积分上限c表示接触区域的半径：c＝b+htanα(4)其中，b表示连接件所受载荷半径，h表示连接件的厚度，α表示有效作用区域内的圆锥体半锥角；由式(1)、(2)和(3)求出系数式中，λ＝c/a。进一步的，步骤二中，利用有限元法和材料力学知识建立栓接结合部等效虚拟材料模型，获得有效接触区域内接触半径与虚拟材料属性间的函数关系；设虚拟材料厚度为l，名义接触面积为Aa，法向力Fn的作用下的变形量为Δl，则虚拟材料的应力为σ＝F/Aa，应变ε＝Δl/l；根据胡克定律σ＝Eε(6)得式中，kn表示虚拟材料的法向刚度其中，αn和βn表示结合部法向特性参数，该属性通过试验测定；将式(8)带入式(7)中，得虚拟材料的弹性模量为对于虚拟材料的泊松比根据Mindlin理论，得结合面法向刚度与切向刚度比为式中，υ表示虚拟材料的泊松比；kt表示虚拟材料的切向刚度其中，αt和βt表示结合部切向特性参数，该属性通过试验测定；根据式(8)、(10)和(11)得虚拟材料的泊松比为虚拟材料的剪切模量为由式(9)、(12)和(13)求出虚拟材料的剪切模量；虚拟材料的平均密度为式中，m1和m2分别表示上下试件微凸体的质量，ρ1和ρ2分别表示上下试件微凸体的密度，l1和l2分别表示上下试件微凸体的厚度；l＝1mm，l1≈l2＝0.5mm，式(14)转化为：相对于现有技术，本专利技术的优点在于：(1)利用有限元法仿真栓接结合部在不同螺栓预紧力下，不同接触区域的受力情况，建立栓接结合部的接触压应力分布模型。(2)根据螺栓预紧力的作用范围和螺栓接触压应力的规律，将接触区域划分为若干区域，分别确定各区域的虚拟材料属性，从而提高分析栓接结合部动态特性的精度。(3)利用材料力学建立各区域虚拟材料模型，该模型包含弹性模量、剪切模量、泊松比、密度和虚拟材料的厚度，其中虚拟材料的弹性模量E、剪切模量G、泊松比υ均是关于螺栓作用力的影响半径的函数，而平均密度和虚拟材料厚度与接触半径无关，为确定值。(4)通过模态试验和有限元仿真分析，对比频响函数曲线的一致性，验证虚拟材料模型的正确性，从而为研究栓接结合部动态特性提供新的思路和方法。附图说明图1是本专利技术确定栓接结合部动态特性的原理图。图2是通过确定螺栓预紧力规律获得多层虚拟材料层图。图3是试件图。图4和图5是仿真与试验频响函数比较图。具体实施方式结合图1、图2和图3，本专利技术是一种利用虚拟材料等效栓接结合部动态特性的方法，包含以下步骤：步骤一、确定栓接结合部接触压应力。假设垂直栓接结合部方向上压应力均匀分布，大小随r向位置变化而变化，且压应力分布呈截头圆锥体；利用有限元法模拟栓接结合部在不同螺栓预紧力下，不同接触区域的受力情况，将栓接结合部接触压力Pn(r)等效为关于(r/a)的4阶多项式函数式中，r表示螺栓作用力的影响半径，a表示螺栓孔半径，a0、a1、a2、a3和a4为4阶多项式函数的相关系数。根据边界条件：由于，螺栓接触载荷F为：式中，积分上限c表示接触区域的半径，根据图1可确定为c＝b+htanα(4)其中，b表示连接件所受载荷半径，h表示连接件的厚度，α表示有效作用区域内的圆锥体半锥角。所以，由式(1)、(2)和(3)可求出系数式中，λ＝c/a。步骤二、确定栓接结合部虚拟材料的参数值。设虚拟材料厚度为l，名义接触面积为Aa，法向力Fn的作用下的变形量为Δl，则虚拟材料的应力为σ＝F/Aa，应变ε＝Δl/l。根据胡克定律σ＝Eε(6)可得式中，kn表示虚拟材料的法向刚度其中，αn和βn表示结合部法向特性参数，该属性通过试验测定。将式(8)带入式(7)中，可得虚拟材料的弹性模量为对于虚拟材料的泊松比根据Mindlin理论，可得结合面法向刚度与切向刚度比为：式中，υ表示虚拟材料的泊松比；kt表示虚拟材料的切向刚度其中，αt和βt表示结合部切向特性参数，该属性通过试验测定。根据式(8)、(10)和(11)可得虚拟材料的泊松比为根据材料力学可知，虚拟材料的剪切模量为由式(9)、(12)和(13)可求出虚拟材料的剪切模量。虚拟材料的平均密度为式中，m1和m2分别表示上下试件微凸体的质量，ρ1和ρ2分别表示上下试件微凸体的密度，l1和l2分别表示上下试件微凸体的厚度。一般假设l＝1mm，l1≈l2＝0.5mm，则式(14)可转化为综上，根据所建立的结合部虚拟材料模型(9)、(12)、(13)和(15)可知，除了虚拟材料的密度以外，其他参数均与接触半径r相关。所以，理论上接触半径取得越密越接近真值。步骤三、有本文档来自技高网...

【技术保护点】
一种利用虚拟材料等效栓接结合部动态特性的方法，其特征在于，包括以下步骤：步骤一、确定栓接结合部接触压应力：假设垂直栓接结合部方向上压应力均匀分布，大小随螺栓作用力的影响半径变化而变化，且压应力分布呈截头圆锥体，利用有限元法模拟栓接结合部在不同螺栓预紧力下，不同接触区域的受力情况，并将接触压力等效为关于接触半径的4阶多项式函数；步骤二、确定栓接结合部虚拟材料的相关参数：建立栓接结合部等效虚拟材料模型，确定有效接触区域内接触半径与虚拟材料属性间的函数关系；虚拟材料属性包括弹性模量E、剪切模量G、泊松比υ、密度和虚拟层厚度；其中虚拟材料的弹性模量E、剪切模量G、泊松比υ均是关于螺栓作用力的影响半径的函数，而平均密度和虚拟材料厚度与接触半径无关，为确定值；步骤三、有限元法确定整体结构的频率响应：根据虚拟材料模型和接触压力的变化曲线，将栓接结合部等效为二层、三层或者多层虚拟材料，利用子结构耦合法建立整体结构的有限元模型，并进行动态响应分析，获得栓接结合部整体结构的频响函数。

【技术特征摘要】
1.一种利用虚拟材料等效栓接结合部动态特性的方法，其特征在于，包括以下步骤：步骤一、确定栓接结合部接触压应力：假设垂直栓接结合部方向上压应力均匀分布，大小随螺栓作用力的影响半径变化而变化，且压应力分布呈截头圆锥体，利用有限元法模拟栓接结合部在不同螺栓预紧力下，不同接触区域的受力情况，并将接触压力等效为关于接触半径的4阶多项式函数；步骤二、确定栓接结合部虚拟材料的相关参数：建立栓接结合部等效虚拟材料模型，确定有效接触区域内接触半径与虚拟材料属性间的函数关系；虚拟材料属性包括弹性模量E、剪切模量G、泊松比υ、密度和虚拟层厚度；其中虚拟材料的弹性模量E、剪切模量G、泊松比υ均是关于螺栓作用力的影响半径的函数，而平均密度和虚拟材料厚度与接触半径无关，为确定值；步骤三、有限元法确定整体结构的频率响应：根据虚拟材料模型和接触压力的变化曲线，将栓接结合部等效为二层、三层或者多层虚拟材料，利用子结构耦合法建立整体结构的有限元模型，并进行动态响应分析，获得栓接结合部整体结构的频响函数。2.根据权利要求1所述的一种利用虚拟材料等效栓接结合部动态特性的方法，其特征在于，垂直栓接结合部方向上压应力大小随螺栓作用力的影响半径变化而变化，且压应力分布呈截头圆锥体，将栓接结合部接触压力Pn(r)等效为关于(r/a)的4阶多项式函数：式中，r表示螺栓作用力的影响半径，a表示螺栓孔半径，a0、a1、a2、a3和a4为4阶多项式函数的相关系数；根据边界条件：由于，螺栓接触载荷为：式中，积分上限c表示接触区域的半径：c＝b+htanα(4)其中，b表示连接件所受载荷半径，h表示连接件的厚度，α表示有效作用区域内的圆锥体半锥角；由式(1)、(2)和(3)求出系数1式中，λ＝...

【专利技术属性】
技术研发人员：李玲，蔡安江，褚崴，
申请(专利权)人：西安建筑科技大学，
类型：发明
国别省市：陕西,61