CN108133097B - 一种抑制褶皱的空间薄膜结构夹具形状优化设计方法 - Google Patents
一种抑制褶皱的空间薄膜结构夹具形状优化设计方法 Download PDFInfo
- Publication number
- CN108133097B CN108133097B CN201711388553.9A CN201711388553A CN108133097B CN 108133097 B CN108133097 B CN 108133097B CN 201711388553 A CN201711388553 A CN 201711388553A CN 108133097 B CN108133097 B CN 108133097B
- Authority
- CN
- China
- Prior art keywords
- film
- clamp
- film structure
- fixture
- shape
- 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
- 238000005457 optimization Methods 0.000 title claims abstract description 35
- 238000000034 method Methods 0.000 title claims abstract description 16
- 230000037303 wrinkles Effects 0.000 title claims abstract description 16
- 230000002401 inhibitory effect Effects 0.000 title claims abstract description 6
- 239000010408 film Substances 0.000 claims abstract description 82
- 239000010409 thin film Substances 0.000 claims abstract description 16
- 238000004458 analytical method Methods 0.000 claims abstract description 9
- 238000006073 displacement reaction Methods 0.000 claims description 8
- 238000012795 verification Methods 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 5
- 238000005054 agglomeration Methods 0.000 claims description 4
- 230000002776 aggregation Effects 0.000 claims description 4
- 238000009833 condensation Methods 0.000 claims description 4
- 230000005494 condensation Effects 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 4
- 230000007547 defect Effects 0.000 claims description 3
- 230000002068 genetic effect Effects 0.000 claims description 2
- 230000015271 coagulation Effects 0.000 claims 1
- 238000005345 coagulation Methods 0.000 claims 1
- 230000009471 action Effects 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 4
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- 238000009826 distribution Methods 0.000 abstract description 2
- 238000002360 preparation method Methods 0.000 abstract description 2
- 230000001105 regulatory effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- 239000003795 chemical substances by application Substances 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Aerials With Secondary Devices (AREA)
- Photovoltaic Devices (AREA)
Abstract
本发明属于航天薄膜结构设计领域，提供一种抑制褶皱的空间薄膜结构夹具形状优化设计方法，解决空间薄膜结构在传统夹具拉伸作用下容易产生褶皱现象的问题。在非线性有限元分析基础上，通过优化夹具的形状改变加载边界条件，最大化薄膜区域内单元的最小主应力，调控薄膜的主应力分布，采用全局优化算法寻找全局最优设计，进而获得“弓形”和“凸形”边界的新型夹具形式，以达到完全拟制褶皱的目的。本发明不仅拟制薄膜中褶皱的产生，而且不对薄膜进行裁剪，能够保证薄膜具有足够大的工作面积；适合于空间天线、太阳帆等空间薄膜结构的夹具设计，有利于拟制薄膜褶皱，保证结构工作性能，并且不增加任何制造、发射及运行成本。
Description
技术领域
本发明属于航天薄膜结构设计领域，涉及一种空间薄膜结构的夹具形状优化设计方法。
背景技术
柔性薄膜具备重量轻、能承受大变形、易于折叠/展开等优点，广泛用于航天空间结构。薄膜和展开机构之间通过刚性夹具连接，在太空中展开后使得膜面内产生张拉应力，用以实现其特定功能。然而，由于薄膜几乎不能承受面内压缩应力，在外荷载作用下容易出现面外屈曲，即褶皱现象。对传统设计给出的结构形式进行物理试验及有限元分析表明：空间薄膜结构在展开后容易产生大量褶皱，会严重影响空间结构的形面精度和使用性能。针对此问题，通过在薄膜内部或边缘进行挖孔，虽然能够抑制褶皱，但必然会减少薄膜面积。为得到既能保证高精度形面要求，又具有足够大薄膜工作面积的有效结构形式，一个非常有效的方法是采用优化设计的手段重新设计夹具结构形式，通过改变薄膜的位移加载边界条件来调整薄膜内部区域的应力分布情况，使整个薄膜区域的最小主应力增大为正值，达到预期的完全张拉状态。
发明内容
针对空间薄膜的高精度形面要求及大面积工作膜面要求，本发明提供一种支撑夹具的形状优化设计方法，该方法能够提高薄膜的最小主应力，拟制薄膜中褶皱的产生，同时能够保证薄膜具有足够大的工作面积。本发明适合于空间天线、太阳帆等空间薄膜结构的夹具设计，有利于拟制薄膜褶皱，保证结构工作性能，并且不增加任何制造、发射及运行成本。
为了达到上述目的，本发明的技术方案为：
一种抑制褶皱的薄膜结构夹具形状优化设计方法，主要包括夹具部件形状优化和空间薄膜结构数值验证两部分，具体步骤如下：
第一步，对夹具部件进行形状优化
传统设计给出的薄膜夹具边线一般为直线，在通过夹具施加位移拉伸载荷作用后，薄膜局部的最小主应力为零或负值，从而发生褶皱现象，不满足形面要求。为得到满足膜面精度设计要求和面积要求的空间薄膜结构，本发明通过形状优化的手段重新设计夹具形式，使得薄膜最小主应力为正值，拟制褶皱产生。
1.1)根据结构的尺寸要求和实际加载情况确定设计域，划分有限元单元网格，建立带夹具的薄膜结构有限元模型。所述的薄膜结构有限元模型中，刚性夹具部件用弹性模量足够大的材料进行模拟；在薄膜与夹具的连接线上选取若干设计点，用B样条函数插值生成夹具和薄膜边界；所述夹具材料的弹性模量不小于薄膜的1000倍。
1.2)对刚性夹具施加位移荷载，采用非线性有限元分析方法对所述薄膜结构有限元模型进行分析。
1.3)以薄膜区域内最小主应力最大化为最终设计目标，对夹具的边线进行设计，构造形状优化模型：
(a)目标：使薄膜区域的最小主应力最大化，即
其中，e为有限元单元编号，S2为最小主应力，Ωm为薄膜区域。
(b)约束：确定薄膜面积用量，作为约束下限。所述的面积用量不小于初始薄膜面积的95％。
(c)设计变量：夹具边线上设计点的坐标。
1.4)根据步骤1.3)建立的形状优化模型，对极小极大优化目标进行凝聚变换，得到等效的优化目标函数。所述凝聚变换包括p-norm方法，其中p-norm凝聚函数的表达式为
其中S*为期望达到的最小主应力，取0.1-1.0，p为凝聚参数，取为20-50。
1.5)根据步骤1.3)和1.4)建立的形状优化模型，采用全局优化算法进行求解，得到夹具形状优化问题的全局最优解。所述优化算法为代理模型算法、遗传算法或基于梯度的优化算法。
第二步，对空间薄膜结构进行数值验证
在第一步1.5)得到的优化后夹具形式的基础上，通过引入薄膜随机缺陷，进行非线性后屈曲分析，对优化后的空间薄膜结构进行效果验证。
本发明的有益效果为：在优化前，夹具为直线形式，薄膜区域在拉伸载荷作用下存在最小主应力为零或负值的情况，从而产生褶皱，不满足高精度形面要求。采用本发明方法得到的“曲边”形式的夹具后，薄膜在荷载作用下不再产生褶皱，最小主应力为正值，也能够保证薄膜面积。该夹具构型比较简单，易于加工制造，且有利于装配以及空间展开，经有限元分析和地面静力试验验证，该结构满足性能要求。
附图说明
图1为本发明实施例提供的一种空间薄膜天线结构设计域。图中：
表示施加在夹具上的位移荷载。
图2(a)为空间薄膜天线结构的单个夹具最优设计图。
图2(b)为采用本发明方法设计得到的空间薄膜天线结构效果图。
图3为本发明实施例提供的一种太阳帆结构设计域。
图4(a)为太阳帆结构的单个夹具最优设计图。
图4(b)为采用本发明方法设计得到的太阳帆结构效果图。
具体实施方式
以下结合技术方案和附图详细叙述本发明的具体实施例。
第一步，对夹具部件进行形状优化
1.1)根据结构的尺寸要求和实际加载情况确定设计域，并划分有限元单元网格，建立带夹具的薄膜结构有限元模型。刚性夹具部件用1000倍薄膜的弹性模量进行模拟。在薄膜与夹具的连接线上选取若干设计点，用B样条函数插值生成夹具和薄膜边界。图1为空间薄膜天线结构的设计域，根据结构上下对称性，每个夹具选取11个设计点，图3为太阳帆结构设计域，根据其对称性，每个夹具选取9个设计点。两种初始结构在拉伸荷载作用下均存在明显的褶皱行为。
1.2)对刚性夹具施加位移荷载，采用非线性有限元分析方法对薄膜结构有限元模型进行分析，获得每个单元的最小主应力。
1.3)以薄膜区域内最小主应力最大化为最终设计目标，对夹具的边线进行设计，构造形状优化模型：
(a)目标：使薄膜区域的最小主应力最大化，即
其中，e为有限元单元编号，S2为最小主应力，Ωm为薄膜区域；
(b)约束：确定薄膜面积用量，作为约束下限。所述的面积用量为初始薄膜面积的95％。
(c)设计变量：夹具边线上设计点的坐标。
1.4)根据步骤1.3)建立的形状优化模型，对极小极大优化目标进行凝聚变换，所述凝聚变换为
其中S*为期望达到的最小主应力，取0.5，p为凝聚参数，取为20。
1.5)根据步骤1.3)和1.4)建立的形状优化问题，采用全局优化算法(如Kriging代理模型方法)进行求解，得到夹具形状优化问题的全局最优解，最终得到“弓形”、“凸形”边界的夹具形式。空间薄膜天线和太阳帆的夹具形状分别见图2(a)和图4(a)所示。
第二步，对空间薄膜结构进行数值验证
在第一步1.5)得到的优化后夹具形式的基础上，通过引入薄膜随机缺陷，进行非线性后屈曲分析，对优化后的空间薄膜结构进行效果验证。结果表明优化得到的“弓形”夹具(如图2(b)所示，薄膜面外位移接近于0，没有局部褶皱现象，图中两侧颜色较深的为夹具，浅灰色区域表示为薄膜)和“凸形”夹具(如图4(b)所示，薄膜面外位移接近于0，没有局部褶皱现象，图中四个角黑色区域表示为夹具，浅灰色区域表示为薄膜)分别有利于提高两种结构的薄膜最小主应力，在不裁剪薄膜、不增加运行成本的情况下避免了局部褶皱现象，满足高精度形面要求。
Claims (3)
1.一种抑制褶皱的空间薄膜结构夹具形状优化设计方法，其特征在于以下步骤：
第一步，对夹具部件进行形状优化
1.1)根据结构的尺寸要求和实际加载情况确定设计域，划分有限元单元网格，建立带夹具的薄膜结构有限元模型；所述的薄膜结构有限元模型中，刚性夹具部件用弹性模量大的材料进行模拟；在薄膜与夹具的连接线上选取若干设计点，用B样条函数插值生成夹具和薄膜边界；
1.2)对刚性夹具施加位移荷载，采用非线性有限元分析方法对所述薄膜结构有限元模型进行分析；
1.3)以薄膜区域内最小主应力最大化为最终设计目标，对夹具的边线进行设计，构造形状优化模型：
(a)目标：使薄膜区域的最小主应力最大化，即
其中，e为有限元单元编号，S2为最小主应力，Ωm为薄膜区域；
(b)约束：确定薄膜面积用量，作为约束下限；所述的面积用量不小于初始薄膜面积的95％；
(c)设计变量：夹具边线上设计点的坐标；
1.4)根据步骤1.3)建立的形状优化模型，对极小极大优化目标进行凝聚变换；
所述凝聚变换包括p-norm方法，其中p-norm凝聚函数的表达式为
其中S*为期望达到的最小主应力，取0.1-1.0，p为凝聚参数，取为20-50；
1.5)根据步骤1.3)和1.4)建立的形状优化问题，采用全局优化算法进行求解，得到夹具形状优化问题的全局最优解；
第二步，对空间薄膜结构进行数值验证
在第一步1.5)得到的优化后夹具形式的基础上，通过引入薄膜随机缺陷，进行非线性后屈曲分析，对优化后的空间薄膜结构进行效果验证。
2.根据权利要求1所述的一种抑制褶皱的空间薄膜结构夹具形状优化设计方法，其特征在于，步骤1.1)中所述材料的弹性模量不小于薄膜的1000倍。
3.根据权利要求1或2所述的一种抑制褶皱的空间薄膜结构夹具形状优化设计方法，其特征在于，步骤1.5)所述优化算法为代理模型算法、遗传算法或基于梯度的优化算法。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711388553.9A CN108133097B (zh) | 2017-12-21 | 2017-12-21 | 一种抑制褶皱的空间薄膜结构夹具形状优化设计方法 |
| US16/346,513 US20200065447A1 (en) | 2017-12-21 | 2018-12-20 | Method for fixture shape optimization design of space membrane structure for inhibiting wrinkling |
| JP2019525755A JP6736108B2 (ja) | 2017-12-21 | 2018-12-20 | しわを抑制する空間薄膜構造治具形状を最適化する設計方法 |
| PCT/CN2018/122264 WO2019120239A1 (zh) | 2017-12-21 | 2018-12-20 | 一种抑制褶皱的空间薄膜结构夹具形状优化设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711388553.9A CN108133097B (zh) | 2017-12-21 | 2017-12-21 | 一种抑制褶皱的空间薄膜结构夹具形状优化设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108133097A CN108133097A (zh) | 2018-06-08 |
| CN108133097B true CN108133097B (zh) | 2021-01-19 |
Family
ID=62390976
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711388553.9A Active CN108133097B (zh) | 2017-12-21 | 2017-12-21 | 一种抑制褶皱的空间薄膜结构夹具形状优化设计方法 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US20200065447A1 (zh) |
| JP (1) | JP6736108B2 (zh) |
| CN (1) | CN108133097B (zh) |
| WO (1) | WO2019120239A1 (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108133097B (zh) * | 2017-12-21 | 2021-01-19 | 大连理工大学 | 一种抑制褶皱的空间薄膜结构夹具形状优化设计方法 |
| CN109033705B (zh) * | 2018-08-24 | 2022-11-04 | 南华大学 | 一种空间褶皱薄膜弹性应变能和主应力有限元分析方法 |
| CN112800642B (zh) * | 2020-12-29 | 2022-07-19 | 中国人民解放军国防科技大学 | 一种基于非线性有限元的薄膜褶皱强度评估方法及*** |
| CN115374542B (zh) * | 2022-10-24 | 2023-03-24 | 中汽研（天津）汽车工程研究院有限公司 | 基于三角形网格位置的网格调整方法、设备和存储介质 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2011137793A (ja) * | 2009-12-31 | 2011-07-14 | Korea Inst Of Science & Technology | フレキシブル基板上の金属層厚さ測定方法 |
| CN105868476A (zh) * | 2016-03-30 | 2016-08-17 | 西安电子科技大学 | 一种空间正多边形平面薄膜结构的最优弧边确定方法 |
| CN105975649A (zh) * | 2016-03-30 | 2016-09-28 | 西安电子科技大学 | 一种方形太阳帆表面褶皱形态的建模方法 |
| CN106777768A (zh) * | 2017-01-09 | 2017-05-31 | 大连理工大学 | 一种用于消除薄膜结构拉伸褶皱的优化设计方法 |
| CN106897491A (zh) * | 2017-01-19 | 2017-06-27 | 大连理工大学 | 一种抑制矩形膜受拉产生褶皱的结构设计方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6876956B1 (en) * | 1999-08-31 | 2005-04-05 | California Institute Of Technology | Method and system for thin-shell finite-element analysis |
| CN108133097B (zh) * | 2017-12-21 | 2021-01-19 | 大连理工大学 | 一种抑制褶皱的空间薄膜结构夹具形状优化设计方法 |
-
2017
- 2017-12-21 CN CN201711388553.9A patent/CN108133097B/zh active Active
-
2018
- 2018-12-20 WO PCT/CN2018/122264 patent/WO2019120239A1/zh active Application Filing
- 2018-12-20 US US16/346,513 patent/US20200065447A1/en not_active Abandoned
- 2018-12-20 JP JP2019525755A patent/JP6736108B2/ja not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2011137793A (ja) * | 2009-12-31 | 2011-07-14 | Korea Inst Of Science & Technology | フレキシブル基板上の金属層厚さ測定方法 |
| CN105868476A (zh) * | 2016-03-30 | 2016-08-17 | 西安电子科技大学 | 一种空间正多边形平面薄膜结构的最优弧边确定方法 |
| CN105975649A (zh) * | 2016-03-30 | 2016-09-28 | 西安电子科技大学 | 一种方形太阳帆表面褶皱形态的建模方法 |
| CN106777768A (zh) * | 2017-01-09 | 2017-05-31 | 大连理工大学 | 一种用于消除薄膜结构拉伸褶皱的优化设计方法 |
| CN106897491A (zh) * | 2017-01-19 | 2017-06-27 | 大连理工大学 | 一种抑制矩形膜受拉产生褶皱的结构设计方法 |
Non-Patent Citations (2)
| Title |
|---|
| An enhanced aggregation method for topology optimization with local stress constraints;yangjun luo 等;《Computer Methods in Applied Mechanics and Engineering》;20130228;第254卷(第2013期);31-41 * |
| 空间薄膜结构皱曲行为与特性研究;王长国;《中国博士学位论文全文数据库 基础科学辑》;20080516(第5期);A005-31 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108133097A (zh) | 2018-06-08 |
| US20200065447A1 (en) | 2020-02-27 |
| WO2019120239A1 (zh) | 2019-06-27 |
| JP6736108B2 (ja) | 2020-08-05 |
| JP2020504859A (ja) | 2020-02-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108133097B (zh) | 一种抑制褶皱的空间薄膜结构夹具形状优化设计方法 | |
| Pandey et al. | A new C0 higher-order layerwise finite element formulation for the analysis of laminated and sandwich plates | |
| CN109241562B (zh) | 基于多尺度有限元方法的微结构材料弹性性能测定方法 | |
| CN104123421A (zh) | 一种基于机电场耦合的静电成形薄膜反射面形态设计方法 | |
| WO2018126465A1 (zh) | 一种用于消除薄膜结构拉伸褶皱的优化设计方法 | |
| CN110162826B (zh) | 薄壁结构热气动弹性动响应分析方法 | |
| CN110245376A (zh) | 一种复合材料层裂二维计算方法 | |
| CN106777768B (zh) | 一种用于消除薄膜结构拉伸褶皱的优化设计方法 | |
| CN106021764A (zh) | 面向机电耦合的有源相控阵天线性能仿真置信度的计算方法 | |
| CN110955941A (zh) | 基于向量场的复合材料结构优化设计方法及设备 | |
| CN102902845A (zh) | 一种直升机旋翼桨叶剖面设计方法 | |
| CN103455670A (zh) | 基于多点约束的多组件结构***布局优化设计方法 | |
| CN105046021A (zh) | 非定常气动力最小状态有理近似的非线性优化算法 | |
| CN104408286A (zh) | 基于正交分解理论的平面构件变形分解与振型识别方法 | |
| CN106250605A (zh) | 基于精确变形控制的薄板压电智能结构协同拓扑优化设计方法 | |
| CN110704953A (zh) | 一种大展弦比机翼静气弹性能设计敏度的分析方法 | |
| Wang et al. | Voronoi polygonal hybrid finite elements with boundary integrals for plane isotropic elastic problems | |
| CN105740547A (zh) | 一种基于弹簧-质量块串子模型的设备和开口加筋板耦合振动基频预报方法 | |
| CN111079327A (zh) | 一种基于各向异性过滤技术的异形曲面加筋拓扑优化方法 | |
| Gerhold et al. | The parallel mesh deformation of the DLR TAU-code | |
| CN111551343B (zh) | 带栅格舵火箭子级全速域气动特性风洞试验设计方法 | |
| CN105426640A (zh) | 基于罚函数的多组件结构***布局优化设计方法 | |
| Shi et al. | Multibody dynamic analysis using a rotation-free shell element with corotational frame | |
| CN107992709A (zh) | 基于中间函数的热结构模型修正方法 | |
| CN107169174A (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 |