CN110321571B - 一种蜂窝板壳结构的力学参数数值提取方法 - Google Patents

一种蜂窝板壳结构的力学参数数值提取方法 Download PDF

Info

Publication number
CN110321571B
CN110321571B CN201810269152.XA CN201810269152A CN110321571B CN 110321571 B CN110321571 B CN 110321571B CN 201810269152 A CN201810269152 A CN 201810269152A CN 110321571 B CN110321571 B CN 110321571B
Authority
CN
China
Prior art keywords
model
honeycomb
shell structure
plate shell
honeycomb plate
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201810269152.XA
Other languages
English (en)
Other versions
CN110321571A (zh
Inventor
骆海涛
王巍
富佳
矫利闯
于长帅
刘广明
陈宁
陈士朋
武廷课
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenyang Institute of Automation of CAS
Original Assignee
Shenyang Institute of Automation of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenyang Institute of Automation of CAS filed Critical Shenyang Institute of Automation of CAS
Priority to CN201810269152.XA priority Critical patent/CN110321571B/zh
Publication of CN110321571A publication Critical patent/CN110321571A/zh
Application granted granted Critical
Publication of CN110321571B publication Critical patent/CN110321571B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Computer Graphics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种蜂窝板壳结构的力学参数数值提取方法，在三维建模软件中建立蜂窝板壳结构的三维模型，并进行模型简化；将模型转换格式，导入前处理软件Hypermesh中，修复蜂窝板壳模型；在Hypermesh中选用壳单元对简化后的三维模型进行离散，添加连接方式，对材料和属性进行设置，完成蜂窝板壳结构有限元模型的构建；在Hypermesh中Nastran求解器下进行分析设置，添加求解目标的控制卡片，求解得到支反力；根据应力应变公式，代入支反力计算出相应的弹性参数。本发明使得在蜂窝板壳结构的力学参数提取过程中能更加快捷方便地获取结果，相比于理论推导等需要繁复的力学运算和简化，大大提高了计算精度和运算速度，节省了时间。

Description

一种蜂窝板壳结构的力学参数数值提取方法
技术领域
本发明涉及蜂窝板壳结构的力学仿真分析领域，具体地说是一种蜂窝板壳结构的力学参数数值提取方法。
背景技术
蜂窝夹层板是目前应用较广的一种多用途结构材料，具有复杂几何结构的正六边形铝蜂窝夹层板是最为常见的蜂窝夹层板。由于具有较高的比强度、比刚度以及较好的隔热隔振和耐冲击、抗疲劳的优点，蜂窝夹层结构广泛应用于航空航天、船舶、交通运输、建筑等领域。例如：飞机和卫星上大量采用铝蜂窝夹层板结构，基于其良好的力学性能和便于减重的特点，有的卫星的主要受力结构几乎全部由蜂窝夹层板结构构成。一般情况下，蜂窝结构质量能占整个卫星板结构的80-90％。
由于在Nastran等当前常用航天产品通用工程软件中没有蜂窝结构单元，因此在对包含蜂窝结构的产品进行数值分析时只能采用三维实体模型或基于等效理论的等效板模型。但是三维模型建模较为复杂，并且由于蜂窝板形状原因导致计算量会很大，同时，蜂窝夹层板壳结构的内部构造使得其整体剪切模量、弹性模量等弹性参数无法直接通过构成材料来直接确定或者进行简单转换得到。在如今科学研究中，计算量小的等效板模型的应用更为普遍。等效理论的重点主要是针对蜂窝板芯层结构进行等效，等效过程主要是保持蜂窝夹层板芯层等效前后剪切模量、弹性模量等弹性参数不变，在等效过程中需要通过一定的方法求解各个方向的弹性参数。常用的方法一般是理论推导求解，这种办法不仅容易出错，而且过程繁复，会耗费大量的时间和精力。于是，蜂窝夹层板壳结构的弹性参数的获取便成为一大问题。
蜂窝板壳结构的力学参数传统的提取方法是运用力学知识进行理论推导，这种方法虽然比较可靠，但是计算过程复杂繁琐，而且最终的数值结果也会忽略某些较小的量，例如因为蜂窝胞壁的厚度相对高度而言很小，因此二者比值的高次幂常常被忽略，因此比较费时间的同时精度也不是很高，这些弊端使得我们必须迫切找到一种更加方便快捷且精确可靠的弹性参数获取方法。
发明内容
针对现有技术的不足，本发明提供一种蜂窝板壳结构的力学参数数值提取方法，可以方便快捷、简单明了地来提取蜂窝板壳结构的各力学参数。
本发明为实现上述目的所采用的技术方案是：
一种蜂窝板壳结构的力学参数数值提取方法，包括以下过程：
步骤1：在三维建模软件中建立蜂窝板壳结构的三维模型，并进行模型简化；
步骤2：将模型转换格式，导入前处理软件Hypermesh中，修复蜂窝板壳模型；
步骤3：在Hypermesh中选用壳单元对简化后的三维模型进行离散，添加连接方式，对材料和属性进行设置，完成蜂窝板壳结构有限元模型的构建；
步骤4：在Hypermesh中Nastran求解器下进行分析设置，添加求解目标的控制卡片，求解得到支反力；
步骤5：根据应力应变公式，代入支反力计算出相应的弹性参数。
所述蜂窝板壳结构的三维模型包括上下面板和由若干正六边形胞元构成的蜂窝芯层。
所述修复蜂窝板壳模型包括生成局部因模型导入而缺失的面。
所述对简化后的三维模型进行离散为壳单元在承载计算过程中，忽略沿壳体厚度方向应力的变化。
所述连接方式为胶粘连接。
所述分析设置包括对边界条件的设置和载荷的施加。
所述载荷的施加是根据目标力学参数计算公式施加给定的应变位移载荷到相应位置。
本发明具有以下有益效果及优点：
1.本发明在蜂窝夹层板壳结构有限元模型的构建过程中，选用壳单元对简化后的三维模型进行离散，忽略沿壳体厚度方向应力的变化，认为应力沿该方向均匀变化，不仅使得在仿真过程中实现简化计算，而且也避免了体网格划分不同导致结果变化的问题；
2.本发明对蜂窝夹层板壳结构的有限元模型施加应变位移载荷，使其对应方向产生特定的应变，精确地模拟出蜂窝夹层板壳结构实际工况下的应力应变情况，最终得到有限元建模分析的结果要比理论计算值更精确；
3.本发明通过有限元建模分析，最终得到相应方向的支反力，再根据给定的应变位移，便可以根据基础的应力应变力学公式得到该力学参数，该方法思路清晰，分析结果精确。
附图说明
图1为本发明的方法流程图；
图2为蜂窝夹层板壳结构的三维图；
图3为蜂窝夹层板壳结构在Hypermesh中的壳单元简化模型示意图；
图4为蜂窝夹层板壳结构的有限元模型示意图；
图5为蜂窝夹层板壳结构坐标系示意图；
图6为求解Ez边界条件及载荷示意图。
具体实施方式
下面结合附图及实施例对本发明做进一步的详细说明。
如图1所示为本发明的方法流程图。
本发明的数值提取方法的过程主要分以下三步：
步骤一，在三维建模软件中建立蜂窝夹层板壳结构的三维模型，完整地画出正六边形蜂窝立体结构。然后，进行相应的模型简化，主要对蜂窝夹层板边界位置由于建模时镜像导致的蜂窝胞壁不规整进行修改，以便于有限元模型的构建；
步骤二，将模型转换格式导入前处理软件Hypermesh中，然后修复蜂窝板壳模型。蜂窝板模型从三维建模软件导入到Hypermesh中，蜂窝壁厚度相对高度而言很小，模型导入后实体变为由外壳包络，并且会导致部分面的缺失情况，因此需要将缺失面补齐，为后续有限元建模做准备；
步骤三，在Hypermesh中选用壳单元对简化后的三维模型进行离散，进一步完成蜂窝板壳结构有限元模型的构建。壳单元忽略沿壳体厚度方向应力的变化，认为应力沿该方向均匀变化，使得在仿真过程中实现简化计算；在前处理过程中，相对于体网格，使用壳单元能够大大减小计算量，而且也避免了体网格划分不同导致结果变化的问题；
步骤四，在Hypermesh中Nastran求解器下进行分析设置，最终求解得到支反力。设置边界条件以及载荷来仿真真实工况下的应力应变进行静力学分析，进而得到整个蜂窝夹层板壳结构在应变位移载荷作用下的下的支反力。对蜂窝夹层板施加边界条件和载荷，以计算求解Ez的分析设置为例；
本发明建立了精确的蜂窝夹层板壳有限元模型，将施加应变位移载荷后仿真分析得到的支反力代入应力应变计算公式，最终得到真实反映蜂窝结构应力应变特性的力学参数，过程简单，结果精确。
本发明针对蜂窝夹层板壳结构的的力学参数数值提取方法来进行有限元建模与仿真分析。蜂窝夹层板壳结构简单，由正六边形蜂窝排列而成。
首先，在三维建模软件中建立了蜂窝夹层板壳的三维模型，并对其进行相应的简化：主要对蜂窝夹层板边界位置由于建模时镜像导致的蜂窝胞壁不规整进行修改，以便于有限元模型的构建，最终三维模型如图2所示。
然后，将模型转换格式导入前处理软件Hypermesh中，然后修复蜂窝板壳模型：由于蜂窝壁厚度相对高度而言很小，模型导入后实体变为由外壳包络，并且会导致部分面的缺失情况，因此需要将缺失面补齐，为后续有限元建模做准备。模型采用壳单元简化如图3所示。
其次，在Hypermesh软件中对蜂窝夹层板壳结构进行前处理，蜂窝夹层板壳结构的有限元建模采用的是PSHELL壳单元，选用壳单元对简化后的三维模型进行离散，进一步完成蜂窝板壳结构有限元模型的构建。壳单元忽略沿壳体厚度方向应力的变化，认为应力沿该方向均匀变化，使得在仿真过程中实现简化计算；在前处理过程中，相对于体网格，使用壳单元能够大大减小计算量，而且也避免了体网格划分不同导致结果变化的问题。该方法是根据蜂窝夹层板结构复杂、计算量大的点特点，用计算量小、处理简单方便的壳网格来完成有限元建模，它既能够保证蜂窝夹层板壳力学参数提取过程中的计算精度，又能够减小网格整体的规模。最终得到蜂窝夹层板壳结构的有限元模型，如图4所示。整个模型共包含4719个单元和4694个节点。
最后，将整个有限元模型在Hypermesh中Nastran求解器下进行分析设置，具体需要对其进行边界条件和载荷的设置。根据受力和应力应变关系，确定蜂窝夹层板壳结构的边界条件设置和载荷设定。边界条件的设置为根据所求的力学参数(以Z方向弹性模量Ez为例)来确定固定位置，在Hypermesh中，蜂窝夹层板壳结构坐标系如图5所示，运用Constraints来对蜂窝夹层板Z向下部整体位置固定。为了提取蜂窝夹层板壳结构的力学参数，对蜂窝夹层板结构的Z向向下施加特定的应变位移载荷εz＝0.001m，由此载荷产生支反力，其具体分析设置如图6所示。
进行静力学求解即可得到最终的结果，仿真分析计算结束后，得到载荷下的支反力。
最终将仿真支反力代入位移载荷和支反力以及蜂窝夹层板的尺寸参数到公式
Figure BDA0001612119080000051
中，其中，F为计算得到的支反力，W为蜂窝夹芯宽度，L为蜂窝夹芯长度，εz为已知施加的位移应变。计算得到Z向弹性模量，同理可计算得到其他力学参数的值，各力学参数提取时的边界条件设置如下：
测定弹性模量Ex时，将X向上边界固定，在X向下边界施加应变位移载荷εx＝0.001m；
测定弹性模量Ey时，将Y向上边界固定，在Y向下边界施加应变位移载荷εy＝0.001m；
测定弹性模量Ez时，将Z向上边界固定，在Z向下边界施加应变位移载荷εz＝0.001m；
测定剪切模量Gxy时，将X向上边界固定，在X向下边界施加应变位移载荷εy＝0.001m；
测定剪切模量Gxz时，将X向上边界固定，在X向下边界施加应变位移载荷εz＝0.001m；
测定剪切模量Gyz时，将Y向上边界固定，在Y向下边界施加应变位移载荷εz＝0.001m；
当施加应变位移载荷时，根据仿真计算得到的应力云图，读取支反力。
最终可以得到如下结论：
1)根据力学公式计算得到Z向弹性模量为0.649GPa，而根据解析法理论计算得到的Z向弹性模量为0.647GPa，通过进一步计算对比发现，弹性模量与剪切模量等结果数值在同一个数量级上，并且符合地较好，说明此方法的合理和准确。尽管两种方法得到的参数有一定的差别，但是考虑到解析法在理论计算中问题的简化和极小量的省略以及本发明涉及的有限元法用壳单元简化等情况存在，数值上较小的出入能够接受。
2)由所得各力学参数计算结果对比可知，由于参数提取方法自身理论的限制，导致解析法的结果并不完整，部分参数无法计算得到，这也从一定程度上说明了本发明的力学参数提取方法的优势所在。

Claims (5)

1.一种蜂窝板壳结构的力学参数数值提取方法，其特征在于：包括以下过程：
步骤1：在三维建模软件中建立蜂窝板壳结构的三维模型，并进行模型简化；
步骤2：将模型转换格式，导入前处理软件Hypermesh中，修复蜂窝板壳模型；
步骤3：在Hypermesh中选用壳单元对简化后的三维模型进行离散，添加连接方式，对材料和属性进行设置，完成蜂窝板壳结构有限元模型的构建；
步骤4：在Hypermesh中Nastran求解器下进行分析设置，添加求解目标的控制卡片，求解得到支反力；所述分析设置包括对边界条件的设置和载荷的施加；
所述载荷的施加是根据目标力学参数计算公式施加给定的应变位移载荷到相应位置；
步骤5：根据应力应变公式，代入支反力计算出相应的弹性参数；
将仿真支反力代入位移载荷和支反力以及蜂窝夹层板的尺寸参数到公式
Figure FDA0003162067350000011
中，其中，F为计算得到的支反力，W为蜂窝夹芯宽度，L为蜂窝夹芯长度，εz为已知施加的位移应变，得到相应的弹性参数。
2.根据权利要求1所述的蜂窝板壳结构的力学参数数值提取方法，其特征在于：所述蜂窝板壳结构的三维模型包括上下面板和由若干正六边形胞元构成的蜂窝芯层。
3.根据权利要求1所述的蜂窝板壳结构的力学参数数值提取方法，其特征在于：所述修复蜂窝板壳模型包括生成局部因模型导入而缺失的面。
4.根据权利要求1所述的蜂窝板壳结构的力学参数数值提取方法，其特征在于：所述对简化后的三维模型进行离散为壳单元在承载计算过程中，忽略沿壳体厚度方向应力的变化。
5.根据权利要求1所述的蜂窝板壳结构的力学参数数值提取方法，其特征在于：所述连接方式为胶粘连接。
CN201810269152.XA 2018-03-29 2018-03-29 一种蜂窝板壳结构的力学参数数值提取方法 Expired - Fee Related CN110321571B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810269152.XA CN110321571B (zh) 2018-03-29 2018-03-29 一种蜂窝板壳结构的力学参数数值提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810269152.XA CN110321571B (zh) 2018-03-29 2018-03-29 一种蜂窝板壳结构的力学参数数值提取方法

Publications (2)

Publication Number Publication Date
CN110321571A CN110321571A (zh) 2019-10-11
CN110321571B true CN110321571B (zh) 2021-09-28

Family

ID=68110639

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810269152.XA Expired - Fee Related CN110321571B (zh) 2018-03-29 2018-03-29 一种蜂窝板壳结构的力学参数数值提取方法

Country Status (1)

Country Link
CN (1) CN110321571B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111027125B (zh) * 2019-12-26 2024-01-26 东南大学 一种蜂窝等效单元及基于该蜂窝等效单元的参数计算方法
JP7348127B2 (ja) 2020-03-30 2023-09-20 日本碍子株式会社 ハニカム構造体の破壊応力の解析方法、ハニカム構造体の目封止部の目封止材のヤング率の解析方法及びハニカム構造体の端面の応力の解析方法
CN112231783A (zh) * 2020-10-22 2021-01-15 上海卫星工程研究所 基于三维模型的蜂窝板精度设计控制和检测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102663153A (zh) * 2012-03-08 2012-09-12 北京航空航天大学 一种异型蜂窝结构的有限元建模方法
CN103984802A (zh) * 2014-04-17 2014-08-13 中国航空工业集团公司沈阳飞机设计研究所 一种蜂窝夹层结构的有限元建模简化方法
CN107515982A (zh) * 2017-08-22 2017-12-26 电子科技大学 一种三维力学有限元模态分析中的接触分析方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10318675B2 (en) * 2015-07-13 2019-06-11 The Boeing Company Post-processing system for finite element analysis
US10502719B2 (en) * 2015-08-21 2019-12-10 The Boeing Company Analysis of a structure modeled with inconsistencies mapped thereon

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102663153A (zh) * 2012-03-08 2012-09-12 北京航空航天大学 一种异型蜂窝结构的有限元建模方法
CN103984802A (zh) * 2014-04-17 2014-08-13 中国航空工业集团公司沈阳飞机设计研究所 一种蜂窝夹层结构的有限元建模简化方法
CN107515982A (zh) * 2017-08-22 2017-12-26 电子科技大学 一种三维力学有限元模态分析中的接触分析方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
菱形夹芯结构蜂窝纸板及其基本力学性能研究;周晟;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20120215;第23-25、27-28页，图3-1、3-3、表3-3 *

Also Published As

Publication number Publication date
CN110321571A (zh) 2019-10-11

Similar Documents

Publication Publication Date Title
US11429765B2 (en) Meshless method for solid mechanics simulation, electronic device, and storage medium
CN110321571B (zh) 一种蜂窝板壳结构的力学参数数值提取方法
CN103488850B (zh) 一种热振复合残余应力定位均化的方法
CN105184390A (zh) 一种壁板结构静强度、刚度、稳定性的综合优化方法
CN104077428A (zh) 一种服务产业联盟的远程有限元分析方法
CN101604353B (zh) 翼组件的设计方法
CN110631792B (zh) 基于卷积神经网络的抗震混合试验模型更新方法
CN114970260A (zh) 一种用于模拟复合材料破坏的格构相场方法
CN115083541A (zh) 基于icvefg法计算正交各向异性材料应力强度因子的方法
CN112580241A (zh) 一种基于结构降阶模型的非线性气动弹性动响应分析方法
CN107688691B (zh) 一种木结构用齿板连接节点受力性能的数值模拟方法
CN112115630A (zh) 一种预测炸药热损伤的热力耦合仿真方法
Chaput et al. Vehicle sketch pad structural analysis module enhancements for wing design
CN115859715A (zh) 一种周期性多孔结构材料等效性能的计算方法
CN114239369B (zh) 一种蜂窝夹芯结构成型过程的受力分析方法
CN103399974B (zh) 量化比较随机振动仿真数据与实验数据的方法
CN115146402A (zh) 高温合金的氧化过程模拟方法及装置
CN108334694B (zh) 结构-土-结构体系三维参数化数值建模分析方法及系统
Zhang et al. The on-orbit thermal-structural analysis of the spacecraft component using MSC/NASTRAN
CN113378321B (zh) 一种试验件的建模方法、装置、设备、存储介质及试验件
Scotti et al. Active control of three surface wind tunnel aeroelastic demonstrator: Modelling and correlation
CN113297709B (zh) 一种试验件的建模方法、装置、设备、存储介质及试验件
Gemma et al. Multi-disciplinary and multi-objective optimization of an unconventional aircraft concept
CN113361099B (zh) 一种裂隙岩体模拟方法及系统
CN117725859A (zh) 飞机疲劳损伤关键部位的疲劳裂纹扩展预测方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210928