一种确定薄壁结构高刚度连接区载荷传递的方法技术

本发明专利技术提出一种确定薄壁结构高刚度连接区载荷传递的方法，通过刚度矩阵缩聚技术，将薄壁结构的刚度特性聚缩到连接区，由此精确得到连接区的聚缩刚度，避免了人为对于加权平均约束单元局部刚度权值的不合理设定处理方法。得到聚缩后的刚度特性后，即可由载荷刚度分配原则，按节点刚度的承载能力进行载荷分配，进而得到合理有效的高精度载荷传递结果，在航空航天类薄壁结构工程分析中具有明显的实用价值及意义。

【技术实现步骤摘要】
一种确定薄壁结构高刚度连接区载荷传递的方法
本专利技术涉及大型工程薄壁结构有限元分析
，具体为一种确定薄壁结构高刚度连接区载荷传递的方法，具体涉及一种离散化处理单元，通过消除薄壁结构离散化有限元分析中高刚度连接区域存在的病态矩阵特性，可将结构高刚度连接构件所施加的惯性或外部载荷精确分配至薄壁主结构。
技术介绍
采用有限元方法建立实际工程结构的离散化数值分析模型，并实施结构设计分析与强度计算，已成为现代结构工程的典型先进设计技术与方法。在航空器薄壁结构初步设计阶段，对主体薄壁结构设计分析时必须计及发动机、起落架等的刚体惯性载荷及外部载荷作用，且这些载荷通常通过高刚度构件连接到主体结构的连接区域。针对薄壁结构这一典型特征，需要在薄壁结构建立离散化有限元模型时，利用特殊算法技术处理这些刚体惯性载荷及外部载荷，将这些刚体的惯性载荷及外部作用的载荷等效地转移到薄壁结构连接区域上(局部结构的部分离散化节点)，从而避免结构有限元建模时局部区域刚度差异过大而导致的数值病态，避免由此导致的计算精度变差乃至得到错误结果的技术问题。查阅大量相关多点约束处理的技术文献，均是在文献1“MSCNastranLinearStaticAnalysisUser’sGuide”基础上的工程应用分析，如文献2“飞机外蒙皮上非受力结构有限元模型设计浅析”介绍了采用文献1的单元技术在外蒙皮非受力口盖等非承力结构气动载荷传递分析中的应用方法。文献1公开了一种加权平均约束单元技术(RBE3)，也是MSC.Nastran软件在刚性连接结构载荷传递分析中的应用方法。该技术通过设置权值的方式完成高刚度连接区载荷分配至薄壁主结构的计算，其权值大小理论上表征薄壁结构主要承力构件在连接区的刚度特性，由此避免高刚度连接件引入结构刚度矩阵所导致的病态问题。但该方法具有以下几个缺点：1.假定薄壁结构连接区域节点在不同方向刚度特性相同，事实上薄壁结构在不同方向上的刚度通常差异很大。2.不能传递弯矩及扭矩，弯矩及扭矩在薄壁结构力学性能分析中起着非常重要的作用，忽略弯矩及扭矩对薄壁结构的影响会导致不合理的分析结果，乃至错误。3.权值的确定具有主观性，且不同权值则会产生不同的计算结果。由于难以给定合理的权值，进而得到不合理的载荷分配结果。这是由于连接区薄壁结构承力构件刚度特性难以直观给出，与薄壁结构的形状、材料以及边界约束条件等密切相关。
技术实现思路
为解决现有技术存在的问题，本专利技术提出了一种确定薄壁结构高刚度连接区载荷传递的方法，改进并修改了加权平均约束单元如
技术介绍
所述的技术缺陷，结果表明该算法技术更合理、更可靠，可满足工程薄壁结构高刚度连接区的载荷传递计算精度需求。本专利技术通过刚度矩阵缩聚技术，将薄壁结构的刚度特性聚缩到连接区，由此精确得到连接区的聚缩刚度，避免了人为对于加权平均约束单元局部刚度权值的不合理设定处理方法。得到聚缩后的刚度特性后，即可由载荷刚度分配原则，按节点刚度的承载能力进行载荷分配，进而得到合理有效的高精度载荷传递结果，在航空航天类薄壁结构工程分析中具有明显的实用价值及意义。本专利技术的技术方案为：所述一种确定薄壁结构高刚度连接区载荷传递的方法，其特征在于：包括以下步骤：步骤1：确定薄壁主结构连接节点gi与刚体载荷节点r之间的线位移及角位移约束关系：薄壁主结构连接节点gi与刚体载荷节点r的线位移分别为及角位移分别为及根据刚体位移约束关系，gi点的线位移和角位移由r点线位移和角位移线性唯一表示为：其中i＝1,....,n，n为主结构连接区节点总数，R为刚体载荷节点r到薄壁结构连接节点gi的向量；并得到矩阵形式表示的薄壁主结构连接节点gi与刚体载荷节点r的位移约束关系：ugi＝Giur其中，步骤2：通过步骤1得到的薄壁主结构连接节点gi与刚体载荷节点r之间的线位移及角位移约束关系，组装得到薄壁主结构连接节点集g与刚体载荷节点r之间的位移转角约束矩阵：则薄壁主结构连接节点集g与刚体载荷节点r的位移及转角满足线性关系：ug＝Gur，其中为薄壁主结构连接节点集g的位移向量；步骤3：通过刚度矩阵的自由度缩聚方法得到聚缩刚度矩阵Kg：薄壁主结构中不包含连接结构的节点集为l集，连接点集为g集，薄壁主结构的整体刚度矩阵表示为：其中Kgg，Kgl，Klg，Kll为整体刚度矩阵K的子矩阵，薄壁主结构的整体刚度方程为：ul为薄壁主结构中不包含连接结构节点集l的位移向量；得到薄壁主结构连接节点集g的位移载荷关系为：(-KglKll-1Klg+Kgg)ug＝Fg从而得到聚缩刚度矩阵Kg为：Kg＝-KglKll-1Klg+Kgg步骤4：根据刚体载荷节点上的惯性载荷及外加结构载荷，得到完整的刚体载荷合力Fr：Fr＝Mrn+F0其中Mr为刚体载荷质量，F0为外加结构载荷，n为刚体载荷过载系数向量；步骤5：根据步骤3得到的聚缩刚度矩阵Kg，得到连接节点集的等效传递载荷向量Fg：Fg＝KgG(GTKgG)-1Fr。有益效果本专利技术的目的是针对薄壁结构刚性连接载荷分配问题中，采用传统大刚度连接件结构或采用加权平均约束单元技术出现的计算不合理缺陷，通过精确计算薄壁结构连接节点刚度特性，从而精确得到薄壁结构关于刚性连接的载荷分配关系。该技术计算过程简单，且得到精确解，计算结果为离散化模型的精确解。附图说明图1为薄壁结构中高刚度连接区载荷传递的有限元设计分析方法流程图；图2为机翼薄壁结构含发动机模型主视图；图3为机翼薄壁结构含发动机模型仰视图；图4为含简化的发动机受载示意图；图5为机翼薄壁结构在发动机载荷作用下应力云图。具体实施方式下面结合具体实施例描述本专利技术：本实施例为挂载发动机的飞机机翼薄壁结构，对于飞机类薄壁结构，由于发动机、起落架以及导弹的刚度通常远大于机翼薄壁结构的刚度，因此可作刚性假设，将这些刚体载荷传递至薄壁主结构。机翼的CAD数据采用CATIA软件建立，机翼主要由上下蒙皮、前后梁腹板和肋腹板构成。机翼主要参数如下：展长为1.196m，根弦长为0.814m，翼尖弦长为0.457m，前缘后掠角为30°，后缘后掠角为15.8°，无扭转角，蒙皮厚度为0.002m，翼肋厚度为0.0015m，梁腹板厚度为0.0015m。图2为机翼模型主视图，图3为机翼模型仰视图。机翼下方挂载由刚性假设得到的简化发动机模型，由一刚性梁代替，这是由于刚体的载荷分配只与刚体的连接方式有关，刚体结构形状并不影响刚体的载荷分配。采用Hypermesh软件对步骤1建立的机翼CAD模型进行网格划分。在机翼的有限元建模过程中，梁腹板及蒙皮等构件主要采用三节点三角形壳单元和四节点四边形壳单元进行有限元离散化，桁条、荆条及缘条均采用两节点的梁单元有限元离散化。为保证数值计算精度，在机翼薄壁结构的建模过程中，尽可能采用了四边形壳单元建模。将建立好以后的机翼有限元模型导入Patran有限元软件，并向Patran有限元软件输入机翼薄壁结构材料参数，机翼由铝合金及钛合金组成，具体材料参数如表1所示。表1机翼薄壁结构材料弹性常数对机翼前后梁施加位移约束工况。机翼的有限元模型和边界约束条件如图2所示。发动机载荷由本专利技术的精确刚性连接载荷传递单元连接到机翼连接节点处。在Patran中完成建立机翼薄壁结构的有限元模型后，生成BDF模型本文档来自技高网...

【技术保护点】
一种确定薄壁结构高刚度连接区载荷传递的方法，其特征在于：包括以下步骤：步骤1：确定薄壁主结构连接节点gi与刚体载荷节点r之间的线位移及角位移约束关系：薄壁主结构连接节点gi与刚体载荷节点r的线位移分别为及角位移分别为及根据刚体位移约束关系，gi点的线位移和角位移由r点线位移和角位移线性唯一表示为：u‾gi=u‾r+θ‾r×R]]>θ‾gi=θ‾r]]>其中i＝1,….,n，n为主结构连接区节点总数，R为刚体载荷节点r到薄壁结构连接节点gi的向量；并得到矩阵形式表示的薄壁主结构连接节点gi与刚体载荷节点r的位移约束关系：ugi＝Giur其中，ugi=u‾giθ‾gi,ur=u‾rθ‾r,]]>步骤2：通过步骤1得到的薄壁主结构连接节点gi与刚体载荷节点r之间的线位移及角位移约束关系，组装得到薄壁主结构连接节点集g与刚体载荷节点r之间的位移转角约束矩阵：G=G1G2...Gn;]]>则薄壁主结构连接节点集g与刚体载荷节点r的位移及转角满足线性关系：ug＝Gur，其中ug=ug1ug2...ugn]]>为薄壁主结构连接节点集g的位移向量；步骤3：通过刚度矩阵的自由度缩聚方法得到聚缩刚度矩阵Kg：薄壁主结构中不包含连接结构的节点集为l集，连接点集为g集，薄壁主结构的整体刚度矩阵表示为：K=KggKglKlgKll]]>其中Kgg，Kgl，Klg，Kll为整体刚度矩阵K的子矩阵，薄壁主结构的整体刚度方程为：KggKglKlgKllugul=Fg0]]>ul为薄壁主结构中不包含连接结构节点集l的位移向量；得到薄壁主结构连接节点集g的位移载荷关系为：(‑KglKll‑1Klg+Kgg)ug＝Fg从而得到聚缩刚度矩阵Kg为：Kg＝‑KglKll‑1Klg+Kgg步骤4：根据刚体载荷节点上的惯性载荷及外加结构载荷，得到完整的刚体载荷合力Fr：Fr＝Mrn+F0其中Mr为刚体载荷质量，F0为外加结构载荷，n为刚体载荷过载系数向量；步骤5：根据步骤3得到的聚缩刚度矩阵Kg，得到连接节点集的等效传递载荷向量Fg：Fg＝KgG(GTKgG)‑1Fr。...

【技术特征摘要】
1.一种确定薄壁结构高刚度连接区载荷传递的方法，其特征在于：包括以下步骤：步骤1：确定薄壁主结构连接节点gi与刚体载荷节点r之间的线位移及角位移约束关系：薄壁主结构连接节点gi与刚体载荷节点r的线位移分别为及角位移分别为及根据刚体位移约束关系，gi点的线位移和角位移由r点线位移和角位移线性唯一表示为：其中i＝1,....,n，n为主结构连接区节点总数，R为刚体载荷节点r到薄壁主结构连接节点gi的向量；并得到矩阵形式表示的薄壁主结构连接节点gi与刚体载荷节点r的位移约束关系：ugi＝Giur其中，步骤2：通过步骤1得到的薄壁主结构连接节点gi与刚体载荷节点r之间的线位移及角位移约束关系，组装得到薄壁主结构连接节点集g与刚体载荷节点r之间的位移转角约束矩阵：则薄壁主结构连接节...

【专利技术属性】
技术研发人员：何鹏秋，孙秦，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西;61