CN102081696B - 层合板离散铺层角度设计方法 - Google Patents
层合板离散铺层角度设计方法 Download PDFInfo
- Publication number
- CN102081696B CN102081696B CN 201110045001 CN201110045001A CN102081696B CN 102081696 B CN102081696 B CN 102081696B CN 201110045001 CN201110045001 CN 201110045001 CN 201110045001 A CN201110045001 A CN 201110045001A CN 102081696 B CN102081696 B CN 102081696B
- Authority
- CN
- China
- Prior art keywords
- design
- centerdot
- laying angle
- value
- design variable
- 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
Images
Landscapes
- Tires In General (AREA)
Abstract
本发明公开了一种层合板离散铺层角度设计方法,用于解决现有的离散铺层角度设计方法中优化问题规模大的技术问题。技术方案是采用二值编码材料插值模型,建立离散铺层角度优化模型,以结构整体刚度最大即柔顺度C最小为设计目标,采用数学规划法进行优化设计。本方法适用于任意数目离散铺层角度的设计问题,所需设计变量数目和优化问题规模均不大于现有方法,且当铺层角度数目较多时能够大大缩减变量数目和优化问题规模。
Description
技术领域
本发明涉及一种层合板铺层设计方法,特别涉及一种层合板离散铺层角度设计方法。
背景技术
层合板由于其轻质、高性能、可设计性强,在航空航天领域得到越来越广泛的应用。层合板通常由多种不同纤维角的铺层复合而成,而每一层可能由具有不同纤维角的补片构成。工程应用中,通常纤维角度取一系列离散值。因此,层合板离散铺层角度设计能够实现结构的轻质高效的性能要求,尽可能的挖掘设计潜力,具有重要的理论和应用价值。
现有的层合板离散铺层角度设计方法主要有两大类。一类是采用遗传算法等作为优化算法。此类方法理论上具有全局寻优能力,并且能处理各类复杂性能要求;然而由于求解效率低下,此类方法无法用于多变量大规模问题的求解。另一类方法是将离散铺层角度设计问题转化为多相材料拓扑优化问题,采用灵敏度驱动的定量优化算法求解。此类方法求解效率高,因此适用于多变量大规模设计问题,但难以保证获得全局最优解,某些复杂设计要求也难以处理。
文献1“Stegmann J,Lund E.Discrete material optimization of general composite shellstructures.International Journal for Numerical Methods in Engineering.2005.62:2009-2027.”公开了一种层合板离散铺层角度优化设计方法;在每一独立设计区域内,该方法每一铺层角度对应一个设计变量;对具有n个独立设计区域问题,每一区域内设计变量数目mv与离散的铺层角度数目m相等;当铺层角度数目m较大时,该方法的设计变量数目也较大,从而在处理大规模问题时需要消耗大量存储空间和计算时间。
文献2“Bruyneel M.SFP-a new parameterization based on shape functions foroptimal material selection:Application to conventional composite plies.Structural andMultidisciplinary Optimization.2011.43(1):17-27.”公开了一种针对0°、90°和±45°四种铺层角度的层合板的优化设计方法;在每一独立设计区域内,该方法需要2个设计变量;对相同问题,该方法所需设计变量数是文献1设计方法的一半,但该方法只适用于4种铺层角度的层合板设计问题。
发明内容
为了克服现有的离散铺层角度设计方法中优化问题规模大的不足,本发明提供一种层合板离散铺层角度设计方法,该方法采用二值编码材料插值模型和相应的数学规划法,建立离散铺层角度优化模型,以结构整体刚度最大即柔顺度C最小为设计目标。该方法适用于任意数目离散铺层角度的设计问题,所需设计变量数目和优化问题规模均不大于现有方法,且当铺层角度数目较多时能够大大缩减变量数目和优化问题规模,可以满足实际工程需求。
本发明解决其技术问题所采用的技术方案:一种层合板离散铺层角度设计方法,其特点是包括以下步骤:
(a)对m个铺层角度设计问题,每个单元所需设计变量数目建立层合板设计空间有限元模型和设计变量初始值xik=0,其中,i表示单元编号,k表示该单元相应的第k个设计变量;给定第j种铺层角度下单元i的刚度矩阵为Ki (j);采用以下公式设置系数sjk的数值
其中
istep=2k
(b)根据当前设计变量xik,采用以下公式计算每一有限元单元的刚度矩阵Ki
式中,wij是每一铺层角度的加权系数,wij计算式是
wij=1表示第i个设计区域采用第j个铺层角度,wij=0表示不采用该铺层角度,p是给定的惩罚系数;
(e)重复步骤(b)至步骤(d),直至最近两次迭代计算得到结构整体柔顺度相对误差小于1%或达到预设的最大迭代次数。
本发明的有益效果是:由于采用二值编码材料插值模型和相应的数学规划法,建立离散铺层角度优化模型,以结构整体刚度最大即柔顺度C最小为设计目标。该方法适用于任意数目离散铺层角度的设计问题,所需设计变量数目和优化问题规模均不大于现有方法,且当铺层角度数目较多时能够大大缩减变量数目和优化问题规模,满足了实际工程需求。
下面结合附图和具体实施方式对本发明作详细说明。
附图说明
图1是实施例1的模型示意图。
图2是实施例1的设计结果图。
图3是背景技术文献1中设计方法的设计结果图。
图4是背景技术文献2中设计方法的设计结果图。
图5是实施例2的模型示意图。
图6是实施例2的设计结果图。
具体实施方式
以下实施例参照图1~6。
实施例1:(a)将长宽均为1m的正方形平面结构划分为16×16的正方形网格;每4×4个单元构成一个设计区域,其中所有单元具有相同的铺层角度;结构左侧固支,载荷为作用在结构右下角的集中载荷F=1000N;正交各相异性材料属性:Ex=146.86GPa,Ey=10.62GPa,Gxy=5.45GPa,vxy=0.33;铺层角度数目m=4,铺层角度分别为90°/45°/0°/-45°;给定设计变量初始值xik=0。
对m个铺层角度设计问题,每个单元所需设计变量数目建立层合板设计空间有限元模型和设计变量初始值xik=0,其中i表示单元编号,k表示该单元相应的第k个设计变量;给定第j种铺层角度下单元i的刚度矩阵为Ki (j);采用以下公式设置系数sjk的数值
其中
istep=2k
计算得到系数sjk的数值分别为s11=s12=s22=s31=-1,s21=s32=s41=s42=1。
(b)给定的惩罚系数p=3;根据当前设计变量值xik,采用以下公式计算每一有限元单元的刚度矩阵Ki
(c)从有限元分析结果中提取每一单元的弹性应变能Senei和单元节点位移向量ui,采用以下公式计算结构整体柔顺度C及其对每一设计变量的灵敏度
(d)根据当前设计变量值xik和灵敏度值以结构整体柔顺度C为目标函数,采用数学规划法对优化问题进行求解得到新的设计变量值xik;
(e)重复步骤(b)至步骤(d),直至最近两次迭代计算得到结构整体柔顺度相对误差小于1%或达到预设的最大迭代次数15。
每一设计区域需要2个设计变量,依据以上迭代步骤得到的铺层角度设计结果,其整体柔顺度为1.162×10-4W。
文献1中设计方法每一设计区域需要4个设计变量,设计结果的整体柔顺度为1.220×10-4W;文献2中设计方法每一设计区域需要2个设计变量,设计结果的整体柔顺度为1.182×10-4W。
实施例1每一设计区域所需设计变量数目与文献2相等,仅为文献1的一半,即优化问题规模仅为文献1的50%;由于结构整体刚度与柔顺度成倒数关系,即柔顺度越小结构整体刚度越大,因此,实施例1的结果整体刚度比文献1大约5%,比文献2大约2%。
实施例2:(a)将长宽分别为6m和1m的矩形形平面结构划分为240×40的正方形网格;每10×10个单元构成一个设计区域,其中所有单元具有相同的铺层角度;结构简支,载荷为作用在结构对称轴上的集中载荷F=1000N;正交各相异性材料属性正交各相异性材料属性:Ex=146.86GPa,Ey=10.62GPa,Gxy=5.45GPa,vxy=0.33;铺层角度数目m=36,铺层角度分别为
90°/85°/80°/75°/70°/65°/60°/55°/50°/45°/40°/35°/30°/25°/20°/15°/10°/5°/0°/-5°/-10°/-15°/-20°/-25°/-30°/-35°/-40°/-45°/-50°/-55°/-60°/-65°/-70°/-75°/-80°/-85°;给定设计变量初始值xik=0;采用以下公式设置系数sjk的数值
其中
istep=2k
计算得到系数sjk的数值分别为
(b)给定的惩罚系数p=3;根据当前设计变量值xik,采用以下公式计算每一有限元单元的刚度矩阵Ki
(e)重复步骤(b)至步骤(d),直至最近两次迭代计算得到结构整体柔顺度相对误差小于1%或达到预设的最大迭代次数30。
每一设计区域需要6个设计变量,依据以上迭代步骤得到的铺层角度设计结果,其整体柔顺度为9.679×10-3W。
文献1中设计方法每一设计区域需要36个设计变量;文献2无法处理本实施例。
实施例1每一设计区域所需设计变量数目仅为文献1的1/6,即优化问题规模仅为文献1的16.7%,文献2无法处理本实施例。
Claims (1)
1.一种层合板离散铺层角度设计方法,其特征在于包括下述步骤:
(a)对m个铺层角度设计问题,每个单元所需设计变量数目建立层合板设计空间有限元模型和设计变量初始值xik=0,其中,i表示单元编号,k表示该单元相应的第k个设计变量;给定第j种铺层角度下单元i的刚度矩阵为采用以下公式设置系数sjk的数值
其中
istep=2k
(b)根据当前设计变量值xik,采用以下公式计算每一有限元单元的刚度矩阵Ki
式中,wij是每一铺层角度的加权系数,wij计算式是
wij=1表示第i个设计区域采用第j个铺层角度,wij=0表示不采用该铺层角度,p是给定的惩罚系数;
(c)从有限元分析结果中提取每一单元的弹性应变能Senei和单元节点位移向量ui,计算结构整体柔顺度C及其对每一设计变量的灵敏度计算式分别是
(e)重复步骤(b)至步骤(d),直至最近两次迭代计算得到结构整体柔顺度相对误差小于1%或达到预设的最大迭代次数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110045001 CN102081696B (zh) | 2011-02-24 | 2011-02-24 | 层合板离散铺层角度设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110045001 CN102081696B (zh) | 2011-02-24 | 2011-02-24 | 层合板离散铺层角度设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102081696A CN102081696A (zh) | 2011-06-01 |
CN102081696B true CN102081696B (zh) | 2013-04-17 |
Family
ID=44087658
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110045001 Active CN102081696B (zh) | 2011-02-24 | 2011-02-24 | 层合板离散铺层角度设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102081696B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103400017B (zh) * | 2013-08-19 | 2016-02-17 | 中航沈飞民用飞机有限责任公司 | 复合材料铺层调整的工程优化方法 |
CN110705100B (zh) * | 2019-09-30 | 2020-10-23 | 西北工业大学 | 对称复合材料层合板弯曲刚度铺层序优化设计方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1673036A (zh) * | 2004-03-25 | 2005-09-28 | 北京航空航天大学 | 遗传算法结构优化网格系统 |
CN101794332A (zh) * | 2010-02-04 | 2010-08-04 | 南京航空航天大学 | 各向异性复材制件热压罐成形复材工装模板设计方法 |
-
2011
- 2011-02-24 CN CN 201110045001 patent/CN102081696B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1673036A (zh) * | 2004-03-25 | 2005-09-28 | 北京航空航天大学 | 遗传算法结构优化网格系统 |
CN101794332A (zh) * | 2010-02-04 | 2010-08-04 | 南京航空航天大学 | 各向异性复材制件热压罐成形复材工装模板设计方法 |
Non-Patent Citations (2)
Title |
---|
修英姝等.复合材料层合板稳定性的铺层优化设计.《工程力学》.2005,第22卷(第6期),212-216. * |
张二亮等.基于材料2结构协同设计的层合板多级优化方法.《机械科学与技术》.2006,第25卷(第1期),53-56. * |
Also Published As
Publication number | Publication date |
---|---|
CN102081696A (zh) | 2011-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110110413B (zh) | 一种基于材料场缩减级数展开的结构拓扑优化方法 | |
Setoodeh et al. | Combined topology and fiber path design of composite layers using cellular automata | |
Eidini | Zigzag-base folded sheet cellular mechanical metamaterials | |
CN103324850A (zh) | 基于多文件流的有限元两级分区两次缩聚并行方法 | |
Ngo-Cong et al. | Free vibration analysis of laminated composite plates based on FSDT using one-dimensional IRBFN method | |
CN107563010B (zh) | 基于形状特征的多尺度结构材料一体化设计方法 | |
Jameson | Airfoils admitting non-unique solutions of the Euler equations | |
CN108829914B (zh) | 一种frp结构件的结构与工艺一体化设计方法 | |
Belarbi et al. | An efficient eight‐node quadrilateral element for free vibration analysis of multilayer sandwich plates | |
CN103294857B (zh) | 以频响位移为优化目标的约束阻尼板的拓扑优化方法 | |
CN102514709B (zh) | 一种采用格栅结构的飞行器机翼盒段及设计方法 | |
CN106372347A (zh) | 改进双向渐进法的等效静载荷法动态响应拓扑优化方法 | |
CN102081696B (zh) | 层合板离散铺层角度设计方法 | |
Nguyen et al. | Experimental study on thrust and power of flapping-wing system based on rack-pinion mechanism | |
CN103455670A (zh) | 基于多点约束的多组件结构系统布局优化设计方法 | |
CN109446585B (zh) | 一种周期性多材料结构的拓扑优化设计方法 | |
Morris et al. | High‐fidelity aerodynamic shape optimization of modern transport wing using efficient hierarchical parameterization | |
Yu et al. | Numerical simulation of oscillating-wing based energy harvest mechanism using the high-order spectral difference method | |
CN104182598A (zh) | 基于水平集法的约束阻尼结构优化设计方法 | |
Lv et al. | Topology optimization of adaptive fluid-actuated cellular structures with arbitrary polygonal motor cells | |
CN205149038U (zh) | 组合式纸蜂窝夹心复合板 | |
CN107726035B (zh) | 一种手风琴式蜂窝结构 | |
Shabir et al. | Explore the chirality engineering of carbon nanotubes for solar cell | |
CN205822456U (zh) | 一种钢结构节点 | |
CN112922992A (zh) | 平面小半周界面型负泊松比-蜂窝复合吸能结构 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |