CN113360972B - 一种带屋面板的凯威特网壳冲击仿真方法 - Google Patents
一种带屋面板的凯威特网壳冲击仿真方法 Download PDFInfo
- Publication number
- CN113360972B CN113360972B CN202110536419.9A CN202110536419A CN113360972B CN 113360972 B CN113360972 B CN 113360972B CN 202110536419 A CN202110536419 A CN 202110536419A CN 113360972 B CN113360972 B CN 113360972B
- Authority
- CN
- China
- Prior art keywords
- impact
- shell
- net shell
- roof
- net
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000004088 simulation Methods 0.000 title claims abstract description 30
- 239000011257 shell material Substances 0.000 claims abstract description 142
- 230000005484 gravity Effects 0.000 claims abstract description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000003116 impacting effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Civil Engineering (AREA)
- Computational Mathematics (AREA)
- Structural Engineering (AREA)
- Architecture (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Instructional Devices (AREA)
Abstract
一种带屋面板的凯威特球面网壳冲击仿真方法。是为了解决带屋面板凯威特球面网壳难以建立、网壳冲击仿真周期长效率低下的问题。本发明包括如下步骤：建立凯威特网壳各部件模型；将模型导入ABAQUS中，生成有限元模型各部件，建立冲击物模型；装配各部件；连接模型各部件，设置边界条件；对屋面板施加荷载并布置重力场，设置冲击物速度；设置分析步，创建冲击物与网壳接触关系，选择接触算法；赋予网壳材料，设置冲击物质量；选择各部件单元类型，部署网格划分策略，记录各位置网格布种数；判断冲击区网格是否过度畸变，运算周期是否正常。本发明用于凯威特网壳冲击仿真。
Description
技术领域
本发明属于结构冲击领域，涉及一种带屋面板的凯威特网壳冲击仿真方法。
背景技术
随着人类社会复杂程度日益加深，结构在使用阶段遭受炸弹爆炸和飞行物冲击等极端偶然荷载作用的风险也在加大。大跨空间结构因其造型优美以及受力性能良好，往往是一个地区的地标性建筑物，凯威特球面网壳是最具代表性的形式之一，因其可灵活地覆盖较大空间且受力均匀而被广泛应用，而这类建筑同时具有人口阶段性密集的特点，更易成为恐怖袭击的目标，一旦在极端偶然事件下发生破坏甚至倒塌，会造成无法估量的严重后果。在此背景下，研究大跨空间结构尤其是凯威特球面网壳的冲击动力响应具有极其重要的政治经济意义。
对大跨空间结构抗冲击研究，当前主要采用数值模拟为主、试验为辅的方法，通过有限元模拟分析结构冲击动力响应，此外，目前针对凯威特球面网壳冲击仿真多数未考虑屋面板，主要原因是带屋面板凯威特球面网壳部件众多，结构复杂，较难建立模型，而且带屋面板凯威特球面网壳形式复杂必然导致运算周期加长效率低下。
发明内容
本发明的目的是提供一种带屋面板的凯威特网壳冲击仿真方法，此方法可以快速建立带屋面板的凯威特球面网壳，精确快速地模拟凯威特球面网壳冲击动力响应，解决带屋面板凯威特球面网壳难以建立、网壳冲击仿真周期长效率低下的问题。
本发明提供一种带屋面板的凯威特网壳冲击仿真方法，包括以下步骤：
步骤1：建立凯威特球面网壳主杆、檩条、立柱、铆钉、屋面板模型；
步骤2：使用ABAQUS中“Import”功能，将步骤1建立的K8型网壳主杆、檩条、立柱、铆钉和屋面板导入ABAQUS中，生成有限元网壳模型各部件，同时建立圆柱体冲击物；
步骤3：将步骤2各部件按网壳主杆、立柱、檩条、铆钉、屋面板、冲击物顺序由下向上放置，完成有限元模型装配；
步骤4：连接网壳各部件，设置网壳支撑位三向不动铰支座，完成有限元模型边界条件设置；
步骤5：对屋面板施加荷载并布置重力场，设置冲击物冲击速度；
步骤6：设置分析步，创建冲击物与网壳之间接触关系，选择接触算法；
步骤7：赋予网壳材料参数，设置冲击物质量；
步骤8：选择各部件单元类型，部署网壳各部件网格划分策略，划分网壳各部件与冲击物网格数量，记录网壳各部件各位置网格布种数；
步骤9：进行冲击物不同速度不同质量时网壳冲击模拟，观察网壳云图冲击区是否有网格过度畸变现象发生以及模型运算是否高效，
若冲击区存在网格过度畸变，回到步骤7，加密冲击区网格数量，
若冲击区网格表现正常，但运算周期过长，回到步骤7，减小冲击区网格数量，
若冲击区网格表现正常且运算周期在可接受范围内，输出网壳冲击动力响应并进行分析。
所述的一种带屋面板凯威特球面网壳冲击仿真方法，所述步骤1中，建立凯威特球面网壳主杆、檩条、立柱、铆钉、屋面板模型，具体过程是，
1.1）确定凯威特球面网壳基本参数，根据网壳主杆与檩条、屋面板几何关系，确定檩条和屋面板的跨度、矢高，使用大跨空间结构软件MST建立网壳主杆、檩条、屋面板模型，
1.2）将网壳主杆—檩条、檩条—屋面板导入AUTO-CAD中，连接网壳主杆与檩条生成主杆，连接檩条与屋面板生成铆钉，
1.3）把屋面板导入Rhino软件中，生成屋面板蒙皮。
所述的一种带屋面板凯威特球面网壳冲击仿真方法，所述步骤1.1）中，网壳主杆与檩条、屋面板几何尺寸关系的确定，具体过程是，
檩条或屋面跨度L 1、矢高f 1，通过下式确定：
式中，L为网壳主杆跨度，f为网壳主杆矢高，a为立柱及铆钉高度。
所述的一种带屋面板凯威特球面网壳冲击仿真方法，所述步骤4中，网壳模型各部件连接方式，具体过程是，对网壳所有节点进行编号，将网壳主杆、立柱、檩条、铆钉、屋面板依次连接，除主杆与立柱采用MPC-Beam刚接外，其余部件均采用MPC-Pin铰接方式连接，设置网壳支撑位三向不动铰支座，完成有限元模型边界条件设置。
所述的一种带屋面板凯威特球面网壳冲击仿真方法，所述步骤6中，冲击物与网壳接触方式为GeneralContact，接触算法采用罚函数法。
所述的一种带屋面板凯威特球面网壳冲击仿真方法，所述步骤8中，网壳主杆、檩条、立柱、铆钉采用显示梁单元B31，屋面板采用显示壳单元S4R，冲击物采用离散刚体单元R3D4。
所述的一种带屋面板凯威特球面网壳冲击仿真方法，所述步骤8中，部署网壳各部件网格划分策略，具体过程是保证冲击物与网壳接触区域网壳网格尺寸小于冲击物尺寸，对冲击物与网壳接触区域网格加密，对非接触区域粗略划分。
有益效果：
1.本发明实现了带屋面板的凯威特球面网壳冲击仿真分析。
2.运用大跨空间结构软件MST建立网壳主杆、檩条、屋面板模型，利用AUTO-CAD建立网壳立柱，铆钉模型，使用Rhino软件生成屋面板蒙皮，解决了带屋面板模型难以建立的问题。
3.选用带屋面板网壳各部件单元，合理简化网壳模型，减少仿真运算周期，提高数值模拟运算效率。
4.合理部署网壳各部件网格划分策略，对网壳与冲击物接触区域网格进行加密，选择合理合适的接触类型和接触算法，以此方式提高冲击仿真精度与可信度，为网壳结构设计提供依据。
附图说明
图1为本发明采用的一种带屋面板的凯威特球面网壳冲击仿真方法流程图；
图2为网壳主杆与檩条、屋面板几何关系示意图；
图3为带屋面板的凯威特球面网壳各部件示意图；
图4为带屋面板的凯威特球面网壳各部件连接方式示意图；
图5为带屋面板的凯威特球面网壳整体示意图；
图6为带屋面板的凯威特球面网壳采用单元示意图；
图7为带屋面板的凯威特球面网壳网格划分策略。
具体实施方式
实施例1：
参照图1，为本发明采用的一种带屋面板的凯威特球面网壳冲击仿真方法流程图，按照一下步骤实施：
步骤1：建立凯威特球面网壳主杆、檩条、立柱、铆钉、屋面板模型；
步骤2：使用有限元分析软件ABAQUS中“Import”功能，将步骤1建立的K8型网壳主杆、檩条、立柱、铆钉和屋面板导入ABAQUS中，生成有限元网壳模型各部件，同时建立圆柱体冲击物；
步骤3：将步骤2各部件按网壳主杆、立柱、檩条、铆钉、屋面板、冲击物顺序由下向上放置，完成有限元模型装配；
步骤4：连接网壳各部件，设置网壳支撑位三向不动铰支座，完成有限元模型边界条件设置；
步骤5：对屋面板施加荷载并布置重力场，设置冲击物冲击速度；
步骤6：设置分析步，创建冲击物与网壳之间接触关系，选择接触算法；
步骤7：赋予网壳材料参数，设置冲击物质量；
步骤8：选择各部件单元类型，部署网壳各部件网格划分策略，划分网壳各部件与冲击物网格数量，记录网壳各部件各位置网格布种数；
步骤9：进行冲击物不同速度不同质量时网壳冲击模拟，观察网壳云图冲击区是否有网格过度畸变现象发生以及模型运算是否高效，
若冲击区存在网格过度畸变，回到步骤7，加密冲击区网格数量，
若冲击区网格表现正常，但运算周期过长，回到步骤7，减小冲击区网格数量，
若冲击区网格表现正常且运算周期在可接受范围内，输出网壳冲击动力响应并进行分析。
实施例2：
根据实施例1所述的一种带屋面板的凯威特网壳冲击仿真方法，所述步骤1中，建立凯威特球面网壳主杆、檩条、立柱、铆钉、屋面板模型，具体过程是，
1.1）确定凯威特球面网壳基本参数，根据网壳主杆与檩条、屋面板几何关系，确定檩条和屋面板的跨度、矢高，使用大跨空间结构软件MST建立网壳主杆、檩条、屋面板模型；
1.2）将网壳主杆—檩条、檩条—屋面板导入AUTO-CAD中，连接网壳主杆与檩条生成主杆，连接檩条与屋面板生成铆钉；
1.3）把屋面板导入Rhino软件中，生成屋面板蒙皮。
实施例3：
根据实施例1或2所述的一种带屋面板的凯威特网壳冲击仿真方法，所述步骤1.1）中，网壳主杆与檩条、屋面板几何尺寸关系的确定，具体过程是，
檩条或屋面跨度L 1、矢高f 1，通过下式确定：
式中，L为网壳主杆跨度，f为网壳主杆矢高，a为立柱及铆钉高度。
实施例4：
根据实施例1或2或3所述的一种带屋面板的凯威特网壳冲击仿真方法，所述步骤4中，网壳模型各部件连接方式，具体过程是，对网壳所有节点进行编号，将网壳主杆、立柱、檩条、铆钉、屋面板依次连接，除主杆与立柱采用MPC-Beam刚接外，其余部件均采用MPC-Pin铰接方式连接，设置网壳支撑位三向不动铰支座，完成有限元模型边界条件设置。
实施例5：
根据实施例1或2或3或4所述的一种带屋面板的凯威特网壳冲击仿真方法，所述步骤6中，冲击物与网壳接触方式为GeneralContact，接触算法采用罚函数法。
实施例6：
根据实施例1或2或3或4或5所述的一种带屋面板的凯威特网壳冲击仿真方法，所述步骤8中，网壳主杆、檩条、立柱、铆钉采用显示梁单元B31，屋面板采用显示壳单元S4R，冲击物采用离散刚体单元R3D4。
实施例7：
根据实施例1或2或3或4或5或6所述的一种带屋面板的凯威特网壳冲击仿真方法，所述步骤8中，部署网壳各部件网格划分策略，具体过程是保证冲击物与网壳接触区域网壳网格尺寸小于冲击物尺寸，对冲击物与网壳接触区域网格加密，对非接触区域粗略划分。
Claims (5)
1.一种带屋面板的凯威特球面网壳冲击仿真方法，其特征是：包括以下步骤：
步骤1：建立凯威特球面网壳主杆、檩条、立柱、铆钉、屋面板模型；
步骤2：使用有限元分析软件ABAQUS中Import功能，将步骤1建立的K8型网壳主杆、檩条、立柱、铆钉和屋面板导入ABAQUS中，生成有限元网壳模型各部件，同时建立圆柱体冲击物；
步骤3：将步骤2各部件按网壳主杆、立柱、檩条、铆钉、屋面板、冲击物顺序由下向上放置，完成有限元模型装配；
步骤4：连接网壳各部件，设置网壳支撑位三向不动铰支座，完成有限元模型边界条件设置；
步骤5：对屋面板施加荷载并布置重力场，设置冲击物冲击速度；
步骤6：设置分析步，创建冲击物与网壳之间接触关系，选择接触算法；
步骤7：赋予网壳材料参数，设置冲击物质量；
步骤8：选择各部件单元类型，部署网壳各部件网格划分策略，划分网壳各部件与冲击物网格数量，记录网壳各部件各位置网格布种数；
步骤9：进行冲击物不同速度不同质量时网壳冲击模拟，观察网壳云图冲击区是否有网格过度畸变现象发生以及模型运算是否高效，
若冲击区存在网格过度畸变，回到步骤7，加密冲击区网格数量，
若冲击区网格表现正常，但运算周期过长，回到步骤7，减小冲击区网格数量，
若冲击区网格表现正常且运算周期在可接受范围内，输出网壳冲击动力响应并进行分析；
所述步骤1中，建立凯威特球面网壳主杆、檩条、立柱、铆钉、屋面板模型，具体过程是，
1.1）确定凯威特球面网壳基本参数，根据网壳主杆与檩条、屋面板几何关系，确定檩条和屋面板的跨度、矢高，使用大跨空间结构软件MST建立网壳主杆、檩条、屋面板模型；
1.2）将网壳主杆—檩条、檩条—屋面板导入AUTO-CAD中，连接网壳主杆与檩条生成主杆，连接檩条与屋面板生成铆钉；
1.3）把屋面板导入Rhino软件中，生成屋面板蒙皮；
所述步骤1.1）中，网壳主杆与檩条、屋面板几何尺寸关系的确定，具体过程是，
檩条或屋面跨度L 1、矢高f 1，通过下式确定：
；
式中，L为网壳主杆跨度，f为网壳主杆矢高，a为立柱及铆钉高度。
2.根据权利要求1所述的一种带屋面板的凯威特球面网壳冲击仿真方法，其特征是：所述步骤4中，网壳模型各部件连接方式，具体过程是，对网壳所有节点进行编号，将网壳主杆、立柱、檩条、铆钉、屋面板依次连接，除主杆与立柱采用MPC-Beam刚接外，其余部件均采用MPC-Pin铰接方式连接，设置网壳支撑位三向不动铰支座，完成有限元模型边界条件设置。
3.根据权利要求1所述的一种带屋面板的凯威特球面网壳冲击仿真方法，其特征是：所述步骤6中，冲击物与网壳接触方式为GeneralContact，接触算法采用罚函数法。
4.根据权利要求1所述的一种带屋面板的凯威特球面网壳冲击仿真方法，其特征是：所述步骤8中，网壳主杆、檩条、立柱、铆钉采用显示梁单元B31，屋面板采用显示壳单元S4R，冲击物采用离散刚体单元R3D4。
5.根据权利要求1所述的一种带屋面板的凯威特球面网壳冲击仿真方法，其特征是：所述步骤8中，部署网壳各部件网格划分策略，具体过程是保证冲击物与网壳接触区域网壳网格尺寸小于冲击物尺寸，对冲击物与网壳接触区域网格加密，对非接触区域粗略划分。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110536419.9A CN113360972B (zh) | 2021-05-17 | 2021-05-17 | 一种带屋面板的凯威特网壳冲击仿真方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110536419.9A CN113360972B (zh) | 2021-05-17 | 2021-05-17 | 一种带屋面板的凯威特网壳冲击仿真方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113360972A CN113360972A (zh) | 2021-09-07 |
| CN113360972B true CN113360972B (zh) | 2023-12-26 |
Family
ID=77526430
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110536419.9A Active CN113360972B (zh) | 2021-05-17 | 2021-05-17 | 一种带屋面板的凯威特网壳冲击仿真方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113360972B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108875209A (zh) * | 2018-06-18 | 2018-11-23 | 西北工业大学 | 一种高速冲击载荷下聚碳酸酯薄板大变形韧性破坏问题仿真方法 |
| KR101959383B1 (ko) * | 2017-10-11 | 2019-03-18 | 단국대학교 산학협력단 | 파라메트릭 기법을 이용한 개폐식 대공간 구조물의 해석 모델링 시스템 |
| CN109877512A (zh) * | 2019-03-18 | 2019-06-14 | 兰州理工大学 | K6型单层球面网壳结构缩尺模型的搭建胎架及施工方法 |
| CN110705142A (zh) * | 2019-08-31 | 2020-01-17 | 南京理工大学 | 一种冲击载荷下Ti/CFRP层合板结构的数值模拟方法 |
| CN111611633A (zh) * | 2020-01-09 | 2020-09-01 | 中南建筑设计院股份有限公司 | 一种基于Revit与Dynamo参数化网壳结构模型快速建模方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10858818B2 (en) * | 2017-07-18 | 2020-12-08 | Kurt Evan Haberman | Interlocking building system using one-piece skin-and-frame panels, vacuum-insulation, vertical slide-locks, multi-story slides, and snap-locks |
-
2021
- 2021-05-17 CN CN202110536419.9A patent/CN113360972B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101959383B1 (ko) * | 2017-10-11 | 2019-03-18 | 단국대학교 산학협력단 | 파라메트릭 기법을 이용한 개폐식 대공간 구조물의 해석 모델링 시스템 |
| CN108875209A (zh) * | 2018-06-18 | 2018-11-23 | 西北工业大学 | 一种高速冲击载荷下聚碳酸酯薄板大变形韧性破坏问题仿真方法 |
| CN109877512A (zh) * | 2019-03-18 | 2019-06-14 | 兰州理工大学 | K6型单层球面网壳结构缩尺模型的搭建胎架及施工方法 |
| CN110705142A (zh) * | 2019-08-31 | 2020-01-17 | 南京理工大学 | 一种冲击载荷下Ti/CFRP层合板结构的数值模拟方法 |
| CN111611633A (zh) * | 2020-01-09 | 2020-09-01 | 中南建筑设计院股份有限公司 | 一种基于Revit与Dynamo参数化网壳结构模型快速建模方法 |
Non-Patent Citations (2)
| Title |
|---|
| 凯威特—联方型单层网壳结构冲击相似律及抗冲击性能研究;钟渝楷;工程科技Ⅱ辑(第2017年第02期期);全文 * |
| 飞行器冲击作用下凯威特型单层球面网壳动力响应研究;孙运慧;工程科技Ⅱ辑(第2018年第03期期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113360972A (zh) | 2021-09-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110555229B (zh) | 一种无网格固体力学仿真方法、电子设备及存储介质 | |
| CN102622485B (zh) | 确定索杆梁系空间结构放样态的形态分析逆迭代法 | |
| CN103093022B (zh) | 电功率分布仿真器 | |
| CN106934136A (zh) | 一种细长型飞行器的振动试验仿真系统及方法 | |
| CN110889166A (zh) | 基于热力耦合约束的航空用轴承支架轻量化设计方法 | |
| Mayes et al. | Design studies for the transmission simulator method of experimental dynamic substructuring | |
| CN113360972B (zh) | 一种带屋面板的凯威特网壳冲击仿真方法 | |
| CN108038277A (zh) | 一种航天器有限元模型的二次缩聚方法 | |
| CN115795725A (zh) | 一种基于拓扑优化的液压支架设计流程及优化方法 | |
| Glinsky et al. | DEVStone: a benchmarking technique for studying performance of DEVS modeling and simulation environments | |
| CN117436313A (zh) | 光伏组件动载测试方法、装置、计算机设备和存储介质 | |
| Caliò et al. | A macro-element approach for modeling the nonlinear behaviour of monumental buildings under static and seismic loadings | |
| Zhang et al. | Study on the collapse performances of fabricated single-layer spatial network structures based on discrete beam elements | |
| CN102930086A (zh) | 基于抽柱法的带断裂钢框架连续倒塌分析方法 | |
| Borgart | New challenges for the structural morphology group | |
| Wang et al. | A rapid design optimization framework | |
| CN110795871A (zh) | 一种水陆两栖飞机船底框强度计算和试验方法 | |
| Lochner-Aldinger et al. | Biomimetic approach to structural morphology | |
| Boudjemai et al. | Small satellite computer-aided design and manufacturing | |
| CN113343529B (zh) | 一种整体壁板结构损伤断裂的全局控制方法和装置 | |
| From et al. | A small satellite system design process | |
| Necevska-Cvetanovska et al. | Methodology for seismic design of R/C building structures | |
| Schaffenberger et al. | Simulating convection with cellular automata. | |
| Akış | Lateral load analysis of shear wall-frame structures | |
| CN118296696A (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 |