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

摘要

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

主权项

1.一种柔性桁架结构基于刚度影响的重分析方法，其特征在于：包括以下步骤：(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)根据位移频响矩阵与动态刚度矩阵的关系，推导可以获得：其中，H^*为摄动以后的位移频响矩阵；步骤(2)基于矩阵修正公式，获得重分析后的位移频响矩阵：其中，I∈R^N×N为单位矩阵。