JP6736108B2 - しわを抑制する空間薄膜構造治具形状を最適化する設計方法 - Google Patents

しわを抑制する空間薄膜構造治具形状を最適化する設計方法 Download PDF

Info

Publication number
JP6736108B2
JP6736108B2 JP2019525755A JP2019525755A JP6736108B2 JP 6736108 B2 JP6736108 B2 JP 6736108B2 JP 2019525755 A JP2019525755 A JP 2019525755A JP 2019525755 A JP2019525755 A JP 2019525755A JP 6736108 B2 JP6736108 B2 JP 6736108B2
Authority
JP
Japan
Prior art keywords
thin film
jig
shape
film structure
design
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
JP2019525755A
Other languages
English (en)
Other versions
JP2020504859A (ja
JP2020504859A5 (ja
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.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
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 Dalian University of Technology filed Critical Dalian University of Technology
Publication of JP2020504859A publication Critical patent/JP2020504859A/ja
Publication of JP2020504859A5 publication Critical patent/JP2020504859A5/ja
Application granted granted Critical
Publication of JP6736108B2 publication Critical patent/JP6736108B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, 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)

Description

本発明は、宇宙薄膜構造設計の分野に属し、空間薄膜構造の治具形状を最適化する設計方法に関する。
フレキシブル薄膜は軽量で、大きな変形に耐えることができ、折り畳み/展開が容易であるなどの利点があり、宇宙空間構造に広く使われている。薄膜と展開機構との間は剛性治具で接続され、宇宙で展開されると膜面内に引張応力が発生することによって、その特定の機能を実現する。しかし、薄膜は面内圧縮応力にほとんど耐えられないため、外荷重の作用で面外屈曲、すなわち、しわ現象が起こりやすい。従来の設計で与えられる構造形式に対して物理試験と有限要素分析を行うと、空間薄膜構造は展開後に大量のしわが発生しやすく、空間構造の表面精度と使用性能に大きく影響する。この問題に対して、薄膜の内部やエッジに穴を開けることで、しわを抑制することはできるが、必然的に薄膜の面積を減少させてしまう。高精度な表面の要求を保証するばかりではなく、十分に広い薄膜の作業面積を有する有効な構造形式を得るための非常に有効な方法は、設計を最適化する手段を用いて治具の構造形式を再設計し、薄膜の変位荷重境界条件を変更することで薄膜内部領域の応力分布状況を調整し、薄膜全域の最小主応力をプラスの値に増大させ、予期される完全な引張状態にすることである。
空間薄膜の高精度な表面の要求及び広い面積の作業膜面の要求に対して、本発明は、支持治具の形状を最適化する設計方法を提供しており、この方法により薄膜の最小主応力を高め、薄膜内のしわの発生を抑制することができると同時に、薄膜が十分に広い作業面積を有するように保証できる。本発明は空間アンテナ、ソーラーセイルなどの空間薄膜構造の治具設計に適しており、薄膜のしわを抑制し、構造の作業性能を保証するのに有利であり、しかも何ら製造、発射及び運転のコストを増加させることもない。
上記の目的を達成するため、本発明に係る技術案は以下のとおりである。
しわを抑制する薄膜構造治具形状を最適化する設計方法であって、主に治具部材の形状の最適化及び空間薄膜構造の数値検証という2つの部分を含み、具体的なステップは以下のとおりである。
第1のステップであって、治具部材に対して形状を最適化する。
従来の設計で与えられる薄膜治具のエッジラインは一般的に直線であり、治具で変位引張荷重を加えた後、薄膜の局部的な最小主応力がゼロかマイナスの値になるため、しわ現象が発生してしまい、表面の要求を満足できないが、膜面精度設計の要求と面積の要求を満足する空間薄膜構造を得るため、本発明は、薄膜最小主応力がプラスの値になるように、形状最適化の手段で治具形式を再設計し、しわの発生を抑制する。
1.1)構造の寸法要求と実際のロード状況に基づいて設計領域を決定し、有限要素セルメッシュを分割し、治具付き薄膜構造の有限要素モデルを作成する。前記の薄膜構造有限要素モデルでは、剛性治具部材は弾性率が十分大きい材料でシミュレーションを行い、薄膜と治具の連結線上でいくつかの設計点を選択し、Bスプライン関数で補間して治具と薄膜との境界を生成させ、前記治具材料の弾性率は薄膜の1000倍以上である。
1.2)剛性治具に変位荷重を加え、非線形有限要素分析方法を用いて前記薄膜構造有限要素モデルを分析する。
1.3)薄膜領域内の最小主応力の最大化を最終的な設計目標とし、治具のエッジラインを設計し、形状最適化モデルを構築することであって:
(a)目標については、薄膜領域の最小主応力を最大化すること、すなわち
ここで、eは有限要素ユニット番号、Sは最小主応力、Ωは薄膜領域である。
(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)で作成した形状最適化モデルに基づいて、グローバル最適化アルゴリズムを用いて解を求め、治具形状最適化の問題のグローバル最適解を得る。前記最適化アルゴリズムは代理モデルアルゴリズム、遺伝的アルゴリズム又は勾配に基づいた最適化アルゴリズムである。
第2のステップであって、空間薄膜構造に対して数値検証を行う。
第1のステップの1.5)で得られた最適化された治具形式に基づいて、薄膜のランダムな欠陥を導入することにより、非線形後座屈分析を行い、最適化された空間薄膜構造に対して効果検証を行う。
本発明の有益な効果は以下のとおりである。最適化される前には、治具は直線形式であり、薄膜領域は引張荷重の作用で最小主応力がゼロ又はマイナスの値になる場合があり、しわが発生してしまい、高精度な表面の要求を満足できない。本発明に係る方法で得られた「曲辺」形式の治具を採用した後は、薄膜は荷重作用でしわが発生しなくなり、最小主応力はプラスの値であり、薄膜面積を保証することもできる。この治具の配置は比較的簡単で、加工製造が容易で、組み立てと空間展開に有利で、有限要素分析と地上静力試験で検証され、この構造は性能要求を満足する。
本発明の実施例で提供される空間薄膜アンテナ構造の設計領域である。 空間薄膜アンテナ構造の単一の治具の最適化の設計図である。 本発明の方法を用いて設計された空間薄膜アンテナ構造の効果図である。 本発明の実施例で提供されるソーラーセイル構造の設計領域である。 ソーラーセイル構造の単一治具の最適化の設計図である。 本発明の方法を用いて設計されたソーラーセイル構造の効果図である。
以下、技術案と図面を組み合わせて本発明の具体的な実施例を説明する。図中、
は治具に加える変位荷重を示す。
第1のステップで、治具部材に対して形状を最適化する。
1.1)構造の寸法要求と実際のロード状況に基づいて設計領域を決定し、有限セルメッシュを分割し、治具付き薄膜構造の有限要素モデルを作成する。剛性治具部材は1000倍の薄膜の弾性率でシミュレーションする。薄膜と治具の連結線上でいくつかの設計点を選択し、Bスプライン関数で補間して治具と薄膜との境界を生成する。図1は空間薄膜アンテナ構造の設計領域であり、構造の上下対称性に基づいて、治具毎に11の設計点を選択し、図3はソーラーセイル構造の設計領域であり、その対称性に基づいて、治具毎に9つの設計点を選択する。2種の初期構造は引張荷重の作用で、ともに明らかなしわ行為が存在する。
1.2)剛性治具に変位荷重を加え、非線形有限要素分析方法を用いて薄膜構造有限要素モデルを分析し、各ユニットの最小主応力を得る。
1.3)薄膜領域内の最小主応力の最大化を最終的な設計目標とし、治具のエッジラインを設計し、形状最適化モデルを構築することであって、:
(a)目標については、薄膜領域の最小主応力を最大化すること、すなわち
であり、ここで、eは有限要素ユニット番号、S2は最小主応力、Ωは薄膜領域である。
(b)拘束については、薄膜面積の使用量を拘束下限として決定し、前記の面積の使用量は初期薄膜面積の95%である。
(c)設計変数については、治具エッジライン上の設計点の座標である。
1.4)ステップ1.3)で作成した形状最適化モデルに基づいて、極小極大最適化目標を凝集変換し、前記凝集変換は
であり、S*は達成しようとする最小主応力であり、0.5を取り、Pは凝集パラメータであり、20を取る。
1.5)ステップ1.3)と1.4)で作成した形状最適化の問いに基づいて、グローバル最適化アルゴリズム(例えば、Kriging代理モデル方法)を用いて解を求め、治具形状最適化の問いのグローバル最適解を得、最終的に「弓形」、「凸形」境界の治具形式を得る。空間薄膜アンテナとソーラーセイルの治具形状については、それぞれ図2(a)と図4(a)を参照されたい。
第2のステップで、空間薄膜構造に対して数値検証を行う。
第1のステップの1.5)で得られた最適化された治具形式に基づいて、薄膜のランダムな欠陥を導入することにより、非線形後座屈分析を行い、最適化された空間薄膜構造に対して効果検証を行う。その結果、最適化された「弓形」治具(図2(b)に示すように、薄膜面外の変位は0に近く、局所的なしわ現象はなく、図中の両側の色が濃いものは治具であり、薄い灰色の領域が薄膜を示す)と「凸形」治具(図4(b)に示すように、薄膜面外の変位は0に近く、局所的なしわ現象はなく、図中の四隅の黒色領域は治具を示し、薄い灰色の領域は薄膜を示す)により、2種の構造の薄膜の最小主応力をそれぞれ高めるのに有利で、薄膜を裁断せず、運転コストを増加させずに局部的なしわ現象を回避し、高精度な表面の要求を満足することがわかる。

Claims (3)

  1. しわを抑制する空間薄膜構造治具形状を最適化する設計方法であって、以下のステップ、すなわち:治具部材に対して形状を最適化する第1のステップであって、
    1.1)構造の寸法要求と実際のロード状況に基づいて設計領域を決定し、有限要素セルメッシュを分割し、治具付き薄膜構造の有限要素モデルを作成し、前記の薄膜構造有限要素モデルで、剛性治具部材は薄膜の1000倍以上の弾性率を有する材料でシミュレーションを行い、薄膜と治具の連結線上でいくつかの設計点を選択し、Bスプライン関数で補間して治具と薄膜との境界を生成させ;
    1.2)剛性治具に変位荷重を加え、非線形有限要素分析方法を用いて前記薄膜構造有限要素モデルを分析し;
    1.3)薄膜領域内の最小主応力の最大化を最終的な設計目標として、治具のエッジラインを設計し、形状最適化モデルを構築し;
    (a)目標については、薄膜領域の最小主応力を最大化すること、すなわち

    ここで、eは有限要素ユニット番号、S2は最小主応力、Ωは薄膜領域であり;
    (b)拘束については、薄膜面積の使用量を拘束下限として決定し、前記の面積の使用量は初期薄膜面積の95%以上であり;
    (c)設計変数について、治具エッジライン上の設計点の座標であり;
    1.4)ステップ1.3)で作成した形状最適化モデルに基づいて、最適化目標を凝集変換し;
    1.5)ステップ1.3)と1.4)で作成した形状最適化の問題に基づいて、グローバル最適化アルゴリズムを用いて解を求め、治具形状最適化の問題のグローバル最適解を得るステップと、
    空間薄膜構造に対して数値検証を行う第2のステップであって、
    第1のステップの1.5)で得られた最適化された治具形式に基づいて、薄膜のランダムな欠陥を導入することにより、非線形後座屈分析を行い、最適化された空間薄膜構造に対して効果検証を行うステップと、を有し、
    ステップ1.4)に記載の凝集変換はp-norm方法を含み、ここで、p-norm凝集関数の式は

    であり、S*は達成しようとする最小主応力であって、0.1〜1.0を取り、Pは凝集パラメータであって、20〜50を取る
    ことを特徴とするしわを抑制する空間薄膜構造治具形状を最適化する設計方法。
  2. ステップ1.5)に記載の最適化アルゴリズムは代理モデルアルゴリズム、遺伝的アルゴリズム又は勾配に基づいた最適化アルゴリズムである
    ことを特徴とする請求項1に記載のしわを抑制する空間薄膜構造治具形状を最適化する設計方法。
  3. ステップ1.5)に記載の最適化アルゴリズムは代理モデルアルゴリズム、遺伝的アルゴリズム又は勾配に基づいた最適化アルゴリズムである
    ことを特徴とする請求項に記載のしわを抑制する空間薄膜構造治具形状を最適化する設計方法。
JP2019525755A 2017-12-21 2018-12-20 しわを抑制する空間薄膜構造治具形状を最適化する設計方法 Expired - Fee Related JP6736108B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
CN201711388553.9 2017-12-21
CN201711388553.9A CN108133097B (zh) 2017-12-21 2017-12-21 一种抑制褶皱的空间薄膜结构夹具形状优化设计方法
PCT/CN2018/122264 WO2019120239A1 (zh) 2017-12-21 2018-12-20 一种抑制褶皱的空间薄膜结构夹具形状优化设计方法

Publications (3)

Publication Number Publication Date
JP2020504859A JP2020504859A (ja) 2020-02-13
JP2020504859A5 JP2020504859A5 (ja) 2020-05-21
JP6736108B2 true JP6736108B2 (ja) 2020-08-05

Family

ID=62390976

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019525755A Expired - Fee Related JP6736108B2 (ja) 2017-12-21 2018-12-20 しわを抑制する空間薄膜構造治具形状を最適化する設計方法

Country Status (4)

Country Link
US (1) US20200065447A1 (ja)
JP (1) JP6736108B2 (ja)
CN (1) CN108133097B (ja)
WO (1) WO2019120239A1 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
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 中汽研(天津)汽车工程研究院有限公司 基于三角形网格位置的网格调整方法、设备和存储介质

Family Cites Families (7)

* Cited by examiner, † Cited by third party
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
KR101181274B1 (ko) * 2009-12-31 2012-09-11 한국과학기술연구원 플렉서블 기판상의 금속층 두께 측정방법
CN105868476B (zh) * 2016-03-30 2018-10-30 西安电子科技大学 一种空间正多边形平面薄膜结构的最优弧边确定方法
CN105975649B (zh) * 2016-03-30 2019-03-19 西安电子科技大学 一种方形太阳帆表面褶皱形态的建模方法
CN106777768B (zh) * 2017-01-09 2022-09-30 大连理工大学 一种用于消除薄膜结构拉伸褶皱的优化设计方法
CN106897491B (zh) * 2017-01-19 2020-03-17 大连理工大学 一种抑制矩形膜受拉产生褶皱的结构设计方法
CN108133097B (zh) * 2017-12-21 2021-01-19 大连理工大学 一种抑制褶皱的空间薄膜结构夹具形状优化设计方法

Also Published As

Publication number Publication date
JP2020504859A (ja) 2020-02-13
CN108133097A (zh) 2018-06-08
CN108133097B (zh) 2021-01-19
WO2019120239A1 (zh) 2019-06-27
US20200065447A1 (en) 2020-02-27

Similar Documents

Publication Publication Date Title
JP6736108B2 (ja) しわを抑制する空間薄膜構造治具形状を最適化する設計方法
US11429765B2 (en) Meshless method for solid mechanics simulation, electronic device, and storage medium
US11003807B2 (en) Techniques for generating materials to satisfy design criteria
CN105183996B (zh) 面元修正与网格预先自适应计算方法
RU2733001C2 (ru) Способ изготовления прокладки и сборочный узел летательного аппарата, содержащий такую прокладку
Maurin et al. Isogeometric analysis for nonlinear planar pantographic lattice: discrete and continuum models
CN108763658B (zh) 基于等几何方法的组合薄壁结构固有频率设计方法
JP2020504859A5 (ja)
CN109858133B (zh) 一种基于应力映射的点阵结构设计与优化方法
CN104881514A (zh) 一种面向3d打印的壳状部件轻量化建模系统
CN111079327A (zh) 一种基于各向异性过滤技术的异形曲面加筋拓扑优化方法
CN111079326A (zh) 二维各向异性网格单元度量张量场光滑化方法
Zhang et al. A morphing wing with cellular structure of non-uniform density
CN110704950A (zh) 一种消除自由飞配平载荷下飞机变形中刚性位移的方法
CN109063401B (zh) 采用等效静态位移实现结构接触力优化的方法
Wei et al. A structured mesh boundary motion approach for simulating wind effects on bluff bodies with changing boundaries
Zhang et al. Hencky bar-net model for buckling analysis of plates under non-uniform stress distribution
O’Hara et al. A two-scale generalized finite element method for fatigue crack propagation simulations utilizing a fixed, coarse hexahedral mesh
CN116757026A (zh) 一种基于等几何分析的板架结构极限强度分析实现方法
CN115618682A (zh) 片状几何微结构的参数化生成和力学性能优化方法与系统
CN115859715A (zh) 一种周期性多孔结构材料等效性能的计算方法
Anagnostou et al. Optimized tooling design algorithm for sheet metal forming over reconfigurable compliant tooling
Pankonien et al. Multi-material printed trailing edge control surface for an aeroservoelastic wind tunnel model
CN111177855B (zh) 一种全局气动弹性优化中的气动结构求解方法及系统
Leicht et al. Design, Manufacturing, and Experiments of Additive Manufactured Stiffened Thin-walled Plates

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190514

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20190514

A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20200123

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200128

A524 Written submission of copy of amendment under article 19 pct

Free format text: JAPANESE INTERMEDIATE CODE: A524

Effective date: 20200409

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200410

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20200707

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200708

R150 Certificate of patent or registration of utility model

Ref document number: 6736108

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees