【技术实现步骤摘要】
基于非协调实体单元的管路系统参数化有限元建模方法
[0001]本专利技术涉及管路系统动力学建模领域,具体而言,尤其涉及一种基于非协调实体单元的管路系统参数化有限元建模方法。
技术介绍
[0002]航空发动机管路系统是发动机燃油、滑油、空气等工作介质的运输通道,是发动机能量供给系统的重要组成部分,因而管路系统工作可靠性直接影响到发动机的安全性及寿命。在实际工作中,管路系统通过卡箍、支架等支撑结构与机匣相连,这样由航空发动机转子产生的激振力会通过机匣传递给管路,从而引起管路系统振动。为了提高管路系统的稳定性和可靠性,通常要对管路系统的动力学进行分析和设计,而有限元法由于具有较好的适用性一直是管路系统动力学建模的最主要的方法。
[0003]有限元法需要构建管路系统的有限元模型,而现有的文献大多基于Hypermesh或ANSYS平台构建管路系统的有限元模型,这种方式在分析单管路方面效果较好。但是,针对航空发动机管路系统中的复杂管路进行系统敷设或拓扑优化时,需要不断更新管路系统的有限元模型,此时基于Hypermesh或ANSYS平台重新创建管路系统的有限元模型,需要实现不同软件之间的数据传输和读取,建模效率极低。
技术实现思路
[0004]根据上述提出的航空发动机管路系统中的复杂管路进行系统敷设或拓扑优化时管路系统建模效率低下的技术问题,而提供一种基于非协调实体单元的管路系统参数化有限元建模方法。本专利技术能够适应任意管形的管路系统,并能够进行快速的管路系统参数化有限元建模。
[0005]本专利
【技术保护点】
【技术特征摘要】 【专利技术属性】
1.一种基于非协调实体单元的管路系统参数化有限元建模方法,其特征在于,包括:获取主节点的形函数以及附加节点的形函数,通过静力凝聚的方式消除附加节点的自由度;基于所述附加节点的自由度对协调实体单元进行修正,从而获取非协调实体单元的刚度矩阵和质量矩阵;对管路系统进行网格划分,获取管路直线段的单元节点坐标;根据管路折弯中线处的折弯半径和空间变换关系,得到了管路弯曲圆弧段的节点对应的折弯半径;获取管路弯曲圆弧段的单元节点坐标;根据管路直线段的单元节点坐标和管路弯曲圆弧段的单元节点坐标获取单元的刚度矩阵和质量矩阵并组集得到总体的刚度矩阵和质量矩阵,从而完成管路系统建模。2.根据权利要求1所述的一种基于非协调实体单元的管路系统参数化有限元建模方法,其特征在于,所述主节点的形函数为:其中,是母单元节点处的坐标值,所述附加节点的形函数为:N9(ξ)=1
‑
ξ2,N
10
(η)=1
‑
η2,N
11
(ζ)=1
‑
ζ23.根据权利要求1所述的一种基于非协调实体单元的管路系统参数化有限元建模方法,其特征在于,单元的刚度矩阵和质量矩阵根据以下计算获取:法,其特征在于,单元的刚度矩阵和质量矩阵根据以下计算获取:法,其特征在于,单元的刚度矩阵和质量矩阵根据以下计算获取:法,其特征在于,单元的刚度矩阵和质量矩阵根据以下计算获取:法,其特征在于,单元的刚度矩阵和质量矩阵根据以下计算获取:法,其特征在于,单元的刚度矩阵和质量矩阵根据以下计算获取:其中,K
e
为单元的刚度矩阵,B为附加应变矩阵,B为主自由度对应的应变矩阵,D为弹性矩阵;M
e
为单元的质量矩阵,ρ为密度,N为形函数矩阵,V为体积。4.根据权利要求1所述的一种基于非协调实体单元的管路系统参数化有限元建模方法,其特征在于,对管路系统进行网格划分,获取管路直线段的单元节点坐标,包括:对应单元的8个主节点分别编号为i,l,k,j,m,p,o,n,其中i,l,p,m组成一个扇区,j,k,o,n组成另一个扇区,针对任意单元,都可以看作其中的初始扇区节点坐标为已知的,选择i,l,p,m扇区的节点坐标为已知的,则j,k,o,n的节点坐标可以表示为x
j
=x
i
+S
L
d
x y
j
=y
i
+S
L
d
y z
j
=z
i
+S
L
d
z
x
n
=x
m
+S
L
d
x y
n
=y
m
+S
L
d
y z
n
=z
m
+S
L
d
z
x
k
=x
l
+S
L
d
x y
k
=y
技术研发人员:季文豪,孙伟,王鑫,张秉杰,李晖,汪博,马辉,林君哲,韩清凯,
申请(专利权)人:中国航发沈阳发动机研究所,
类型:发明
国别省市:
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。