[发明专利]一种柔性桁架结构基于刚度影响的重分析方法

说明书

本发明提供了一种柔性桁架结构基于刚度影响的重分析方法，首先基于有限元分析获得柔性平面桁架位移频响函数，构造位移频响矩阵，当某一单元的弹性模量发生改变时，确定全局总刚度矩阵变化量，基于矩阵修正公式，根据初始位移频响矩阵快速获得修正后的结构响应，完成频响动态重分析求解。因此，无需进行多次有限元计算，利用初始的频响动态响应信号及明确结构的局部刚度变化，即可完成刚度摄动后结构的动态分析，简化计算效率，更加方便，具有实际工程意义。

技术领域

本发明涉及一种重分析方法，具体涉及一种基于刚度影响的重分析方法。

背景技术

重分析方法作为一种能够根据初始计算结果快速估计修改后结构的快速计算方法，近几十年来得到了广泛研究，并取得了一系列具有理论价值和工程意义的成果，在机械、土木等结构设计领域得到了广泛的应用。

实际工程中，有时候需要对柔性桁架结构进行局部设计调整，如改变特定桁架结构的弹性模量，必然会导致局部结构刚度矩阵发生改变，从而影响整个结构的刚度矩阵分布，结构的动态特性也随之发生改变。如何利用初始响应的信号，避免重新进行有限元计算，从而快速有效地获得结构刚度修正后动态响应，已成为亟待解决的实际工程问题。

发明内容

发明目的：本发明的目的在于针对现有技术的不足，提供一种柔性桁架结构基于刚度影响的重分析方法。

技术方案：本发明提供了一种柔性桁架结构基于刚度影响的重分析方法，包括以下步骤：

(1)基于Matlab软件进行有限元分析，获得柔性平面桁架位移频响函数，构造位移频响矩阵，当结构某一单元的弹性模量发生改变时，计算总刚度矩阵变化量，并确定其与摄动后的位移频响矩阵的关系；

(2)基于矩阵修正公式，根据步骤(1)初始位移频响函数计算获得修正后的结构响应，完成刚度重分析求解。

进一步，步骤(1)包括以下步骤：

(11)结构的自由度为N，获得的柔性平面桁架位移频响矩阵为：

其中，h_pq表示在结构节点q作用单位脉冲下、结构节点p的位移响应函数，p＝1，2…N，q＝1，2…N；

(12)以柔性平面桁架某一结构单元t为例，对应节点编号分别为i、j，单元与x轴正方向逆时针倾斜角为θ，该单元的刚度矩阵为：

其中，E为弹性模量，A为单元横截面积，l为单元长度；

(13)由倾斜角为θ可知，转换矩阵为：

(14)当弹性模量E增加ΔE时，对应的单元刚度矩阵变化量为ΔK^t：

(15)根据公式(3)(4)获得全局坐标下的单元刚度矩阵变化量：

其中，S，G分别为

(16)由刚度矩阵组装的性质可知，全局坐标下的单元刚度矩阵变化量中的元素分别放置在总刚度矩阵变化量ΔK中第(2i-1),(2i),(2j-1),(2j)行与列交叉位置处，则总刚度矩阵变化量ΔK可表示为：

其中表示该矩阵中的第(2i-1),(2i),(2j-1),(2j)行列交叉处分别对应着矩阵S，G中的元素；

(17)根据位移频响矩阵与动态刚度矩阵的关系，推导可以获得：