航空航天飞行器进行高超音速飞行的情况下,其表面受到高速气流的气动加热和对流换热联合作用,飞行器前缘最高表面温度甚至超高1800℃,会对飞行器结构产生不良影响,因此必须对飞行器前缘结构进行数值仿真,通过数值仿真确定出强迫对流换热系数。
由于飞行器前缘承受气动加热,气动加热作为施加热流与结构表面温度存在耦合关系,因此,建立有限元模型时,必须能够快速提取飞行器前缘体单元表面编号和飞行器前缘体内部圆孔体单元表面编号。
技术实现要素:
本发明的目的是:提供一种飞行器前缘结构有限元模型建立方法。
本发明的技术方案如下:
附图说明
图1是飞行器前缘结构示意图;
图2是单孔壳单元连接图;
图3是完整壳单元有限元网格;
图4是沿飞行器前缘母线拉伸生成的体单元有限元网格;
图5是沿通过镜像生成的完整体单元有限元网格;
图6是飞行器前缘外表面施加的热流载荷;
图7是在飞行器前缘的九个孔的内表面施加强迫对流边界条件;
图8是飞行器前缘某时刻的温度云图显示。
具体实施方式
发明原理:
一个完整的有限元模型包含网格划分、载荷、约束的施加以及边界条件施加,对于结构比较复杂的模型,根据载荷、约束、边界条件的不同,在建立模型时需要通盘考虑。
对于给定的飞行器前缘,其母线、施加对流的内部通道尺寸和数量已经完全确定,有限元数值分析时要求在飞行器前缘外表面施加热流,在飞行器前缘内部圆孔表面施加对流边界条件,由于结构形式的复杂性,在有限元模型建立时必须采取一些技巧才能快速、准确地对家里的模型施加载荷与边界条件。
下面结合附图和实施例对本发明做进一步说明:
步骤一、对飞行器前缘在对称中面实施切分,得到结构中心面(图1),对结构中心面的第一个圆孔区域进行平面网格划分,划分结果见图2,在有限元建模软件平台自动形成的单元连接方式见表1;
表1quad4网格连接方式
步骤二、表1满足右手螺旋定律模型建立规则,对当前壳单元沿着面法线方向拉伸为单层体单元(拉伸控制参数取为1),对体单元孔内表面施加对流边界条件,得知单元表面有两种编号,必须把两种编号改变为统一编号,则在全结构体单元孔内表面施加对流边界条件才能方便地选择;
步骤三、手动修改表1的壳单元cquad4编号为表2,由于单层体单元的建立仅仅是为了显示圆孔的单元内表面的编号,圆孔的单元内表面的编号确定后,删除单层体单元;
步骤四、对图2在当前坐标系进行平移复制,形成对称面上的全部壳单元(图3);
表2修改后的quad4网格连接方式
步骤五、控制参数取为30,将对称面上的全部壳单元作为基单元沿两条曲线滑动生成体单元,如图4所示。
步骤六、镜像(mirror)生成另一部分体单元网格,如图5所示;
步骤七、删除壳单元网格,合并重合点;
步骤八、在飞行器前缘外表面施加热流载荷,如图6所示;
步骤九、在飞行器前缘的九个孔的内表面施加强迫对流边界条件,如图7所示。