CN103455670A - 基于多点约束的多组件结构系统布局优化设计方法 - Google Patents
基于多点约束的多组件结构系统布局优化设计方法 Download PDFInfo
- Publication number
- CN103455670A CN103455670A CN2013103768210A CN201310376821A CN103455670A CN 103455670 A CN103455670 A CN 103455670A CN 2013103768210 A CN2013103768210 A CN 2013103768210A CN 201310376821 A CN201310376821 A CN 201310376821A CN 103455670 A CN103455670 A CN 103455670A
- Authority
- CN
- China
- Prior art keywords
- assembly
- finite element
- displacement
- design
- connected node
- 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.)
- Granted
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于多点约束的多组件结构系统布局优化设计方法,用于解决现有多组件结构系统布局优化设计方法效率低的技术问题。技术方案是建立组件和结构的有限元模型,使用多点约束方法建立组件和结构之间的连接关系,创建优化问题的数学模型,进行结构系统有限元分析并求解灵敏度信息,最终优化得到合理的结果。由于引入多点约束技术建立组件和支撑结构的连接关系,避免了组件附近的结构网格重新划分,不仅减少了有限元模型前处理所需要的时间,而且也减少了整体有限元模型的单元总数目,缩短了有限元分析计算所需的时间。对于二维支撑多组件结构系统布局优化设计问题,整个优化过程所需时间由背景技术的20分钟46秒缩短为5分13秒。
Description
技术领域
本发明涉及一种多组件结构系统布局优化设计方法。特别涉及一种基于多点约束的多组件结构系统布局优化设计方法。
背景技术
航空航天、船舶、汽车等领域存在着大量的多组件结构系统设计问题。由于复杂的服役状况和苛刻的性能要求,多组件结构系统的力学性能设计问题在航空航天飞行器结构设计领域尤为突出。由于组件和支撑结构都具有其自身的力学性能,因此组件的摆放和结构的构型这两方面布局因素从根本上决定了结构系统的综合力学性能。为了提高结构系统的总体刚度,需要对整个系统的支撑结构材料分布问题和组件装填布局问题同时进行优化设计。
参照图1。文献1“Zhu J.H.,Beckers P.Zhang W.H.,On the multi-component layoutdesign with inertial force.Journal of Computational and Applied Mathematics.2010,234(7):2222-2230”公开了一种多组件结构系统布局优化设计方法,这种方法结合了结构拓扑优化技术和装填布局优化技术,采用了网格重新划分技术,通过重新划分可移动组件周围的结构网格1,实现了组件网格2和结构网格1之间节点对节点的连接,将组件嵌入到结构网格中,并进行优化设计得到合理的结果。对于二维支撑多组件结构系统布局优化设计问题,整个优化过程所需时间为20分钟46秒。
文献1公开的方法虽然能够解决多组件结构系统的布局优化设计问题,但是该方法采用的网格重新划分技术不可避免的增加了有限元模型单元的数目,使得有限元分析不仅在前处理划分网格阶段时花费了更多时间,还增加了有限元分析计算的时间,降低了优化求解的效率,不利于处理大规模复杂问题。
发明内容
为了克服现有多组件结构系统布局优化设计方法效率低的不足,本发明提供一种基于多点约束的多组件结构系统布局优化设计方法。该方法通过建立组件和结构的有限元模型,使用多点约束技术建立组件和结构之间的连接关系,创建优化问题的数学模型,进行结构系统有限元分析并求解灵敏度信息,最终优化得到合理的结果。由于引入多点约束技术建立组件和支撑结构的连接关系,避免了组件附近的结构网格重新划分,不仅减少了有限元模型前处理所需要的时间,而且也减少了整体有限元模型的单元总数目,缩短了有限元分析计算所需的时间,可以提高多组件结构系统布局优化设计的效率。
本发明解决其技术问题所采用的技术方案是:一种基于多点约束的多组件结构系统布局优化设计方法,其特点是包括以下步骤:
步骤一、通过组件和支撑结构的CAD模型建立有限元模型;为组件和支撑结构的材料属性赋值;定义载荷和边界条件。
步骤二、使用多点约束方法,建立组件网格2的组件连接节点3和结构网格1的结构连接节点4的位移线性关系式:
当组件连接节点3的数目大于等于1时,则建立相应的组件网格2上的组件连接节点3与结构网格1上的节点之间的位移线性关系式,并且将这些位移线性关系式改写为矩阵形式,如公式2所示:
HU=0 (2)
式中,H为上述位移线性关系式中的系数项,U表示由结构网格1上所有的节点的位移列向量和组件网格2上所有的节点位移列向量组合形成的系统总节点位移列向量。
步骤三、建立优化问题的数学模型为:
findX=(x1,x2,,xn),ξ=(ξjx,ξjy,ξjθ),j=1,2,Nc
s.t.KU+HTλ=F
HU=0 (3)
V(X)-V0≤0
式中,X为设计域上的拓扑优化设计变量;n为拓扑优化设计变量个数;ξ为表征组件安装位置和安装角度的几何设计变量,Nc表示组件个数,ξjx、ξjy和ξjθ分别表示第j个组件在x、y的安装位置坐标和安装角度;Φ为目标函数,表示为结构的总体应变能;K为有限元模型总体刚度矩阵;F为节点等效载荷向量;U为节点整体位移向量;V(X)为支撑结构材料用量,V0为材料用量上限;Ωj、Ωj1和Ωj2分别是第j个、第j1和第j2个组件占据的区域,Ωd表示拓扑优化设计域。
步骤四、将有限元模型进行一次有限元分析;分别对几何设计变量和伪密度设计变量进行灵敏度分析,求得目标函数和约束条件的灵敏度,选取梯度优化算法进行优化设计,得到优化结果。
本发明的有益效果是:该方法通过建立组件和结构的有限元模型,使用多点约束技术建立组件和结构之间的连接关系,创建优化问题的数学模型,进行结构系统有限元分析并求解灵敏度信息,最终优化得到合理的结果。由于引入多点约束技术建立组件和支撑结构的连接关系,避免了组件附近的结构网格重新划分,不仅减少了有限元模型前处理所需要的时间,而且也减少了整体有限元模型的单元总数目,缩短了有限元分析计算所需的时间,提高了多组件结构系统布局优化设计的效率。对于二维支撑多组件结构系统布局优化设计问题,整个优化过程所需时间由背景技术的20分钟46秒缩短为5分13秒。
下面结合附图和实施例对本发明作详细说明。
附图说明
图1是背景技术方法采用网格重新划分技术的示意图。
图2是本发明方法多点约束的示意图。
图3是本发明方法实施例中的工况图。
图4是本发明方法实施例中的优化结果图。
图中,1-结构网格;2-组件网格;3-组件连接节点;4-结构连接节点;5-拓扑设计域;6-组件;7-连接位置;8-结构构型。
具体实施方式
参照图2-4。本发明基于多点约束的多组件结构系统布局优化设计方法具体包括以下步骤。
二维支撑结构拓扑设计域5的尺寸为长1.8m,高0.6m,厚度为0.01m,共划分了2700个壳单元,因此拓扑设计变量的个数为2700。两个相同的长方形组件6长均为0.24m,高0.12m,厚度为0.01m。支撑结构右上端和上端靠右三分之一处固定。两个组件6的x向初始位置坐标分别是0.6m和1.2m,y向初始位置坐标均为0.3m,初始角度都为0度,并且通过连接位置7处的节点与拓扑设计域5连接。同时优化设计支撑结构和组件位置布局,使得整体结构刚度最大,材料用量体积分数最大为50%。方法步骤如下:
(a)通过组件6和支撑结构的CAD模型建立有限元模型。支撑结构材料的杨氏模量设为70Gpa,泊松比为0.3,;两个相同长方形组件6的材料的杨氏模量为200Gpa,泊松比为0.3。拓扑设计域5的右上端和上端距离设计区域右边界三分之一处固定。拓扑设计域5的左下端和下端距离设计域左边界三分之一处分别施加了向左和向下的2000N的力。
(b)使用多点约束方法,建立组件网格2的组件连接节点3和结构网格1的结构连接节点4的位移线性关系式:
当组件连接节点3的数目大于等于1时,则建立相应的组件网格2上的组件连接节点3与结构网格1上的节点之间的位移线性关系式,并且将这些位移线性关系式改写为矩阵形式:
HU=0 (2)
式中,H表示位移线性关系式中的系数项,U表示由结构上的所有节点的位移列向量和组件6上的所有节点的位移列向量组合形成的系统节点位移总列向量。
(c)建立优化模型为:
findX=(x1,x2,,x2700),ξ=(ξjx,ξjy,ξjθ),j=1,2
s.t.KU+HTλ=F
HU=0 (3)
V(X)-0.5≤0
式中,X为设计域上的拓扑优化设计变量;ξ为表征组件安装位置和安装角度的几何设计变量,ξjx、ξjy和ξjθ分别表示第j个组件在x、y的安装位置坐标和安装角度;Φ为目标函数,表示为结构的总体应变能;K为有限元模型总体刚度矩阵;F为节点等效载荷向量;U为节点整体位移向量;V(X)为支撑结构材料用量,材料用量上限为0.5。Ω1和Ω2分别是第一个和第二个组件占据的区域,Ωd表示拓扑设计域10。
(d)用有限元软件Ansys将模型进行一次有限元分析;再通过结构优化平台Boss-Quattro进行优化灵敏度分析,求得目标函数和约束条件的灵敏度,选取梯度优化算法GCMMA(Globally Converged Method of Moving Asymptotes)进行优化设计,得到了最终的优化结果。
从图4可以看出,最终结果中结构构型8清晰合理,组件6也在适当的位置摆放。通过本发明方法进行多组件结构系统布局优化设计,在有限元分析阶段不需要重新划分网格,也减少了有限单元的总数目,因此减少了优化所需的时间。应用背景技术文献1的方法,整个优化过程所需时间为20分钟46秒,而应用本发明方法,整个优化过程只需要了5分13秒,提高了多组件结构系统布局优化设计的效率。
Claims (1)
1.一种基于多点约束的多组件结构系统布局优化设计方法,其特征在于包括以下步骤:
步骤一、通过组件和支撑结构的CAD模型建立有限元模型;为组件和支撑结构的材料属性赋值;定义载荷和边界条件;
步骤二、使用多点约束方法,建立组件网格(2)的组件连接节点(3)和结构网格1的结构连接节点(4)的位移线性关系式:
当组件连接节点(3)的数目大于等于1时,则建立相应的组件网格(2)上的组件连接节点(3)与结构网格(1)上的节点之间的位移线性关系式,并且将这些位移线性关系式改写为矩阵形式:
HU=0 (2)
式中,H为上述位移线性关系式中的系数项,U表示由结构网格(1)上所有的节点的位移列向量和组件网格(2)上所有的节点位移列向量组合形成的系统总节点位移列向量;
步骤三、建立优化问题的数学模型为:
findX=(x1,x2,,xn),ξ=(ξjx,ξjy,ξjθ),j=1,2,Nc
s.t.KU+HTλ=F
HU=0 (3)
V(X)-V0≤0
式中,X为设计域上的拓扑优化设计变量;n为拓扑优化设计变量个数;ξ为表征组件安装位置和安装角度的几何设计变量,Nc表示组件个数,ξjx、ξjy和ξjθ分别表示第j个组件在x、y的安装位置坐标和安装角度;Φ为目标函数,表示为结构的总体应变能;K为有限元模型总体刚度矩阵;F为节点等效载荷向量;U为节点整体位移向量;V(X)为支撑结构材料用量,V0为材料用量上限;Ωj、Ωj1和Ωj2分别是第j个、第j1和第j2个组件占据的区域,Ωd表示拓扑优化设计域;
步骤四、将有限元模型进行一次有限元分析;分别对几何设计变量和伪密度设计变量进行灵敏度分析,求得目标函数和约束条件的灵敏度,选取梯度优化算法进行优化设计,得到优化结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310376821.0A CN103455670B (zh) | 2013-08-27 | 2013-08-27 | 基于多点约束的多组件结构系统布局优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310376821.0A CN103455670B (zh) | 2013-08-27 | 2013-08-27 | 基于多点约束的多组件结构系统布局优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103455670A true CN103455670A (zh) | 2013-12-18 |
CN103455670B CN103455670B (zh) | 2016-02-24 |
Family
ID=49738026
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310376821.0A Active CN103455670B (zh) | 2013-08-27 | 2013-08-27 | 基于多点约束的多组件结构系统布局优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103455670B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104484527A (zh) * | 2014-12-17 | 2015-04-01 | 清华大学 | 一种离散结构拓扑优化过程中均布载荷自动动态修改方法 |
CN104732043A (zh) * | 2015-04-14 | 2015-06-24 | 河北省电力勘测设计研究院 | 变电架构的设计方法 |
CN104992012A (zh) * | 2015-06-26 | 2015-10-21 | 华晨汽车集团控股有限公司 | 汽车后门刚度分析方法 |
CN105426640A (zh) * | 2015-12-28 | 2016-03-23 | 西北工业大学 | 基于罚函数的多组件结构系统布局优化设计方法 |
CN106295028A (zh) * | 2016-08-16 | 2017-01-04 | 西北工业大学 | 一种局部结构动力学建模方法及装置 |
CN109145495A (zh) * | 2018-09-11 | 2019-01-04 | 吉林大学 | 固定分型线双模铸造件多组件的拓扑优化方法 |
CN109543207A (zh) * | 2018-09-11 | 2019-03-29 | 吉林大学 | 考虑变分型线实现双模铸造件多组件设计的方法 |
CN109766564A (zh) * | 2018-10-31 | 2019-05-17 | 中国飞机强度研究所 | 考虑组件保形约束的多组件结构系统布局优化设计方法 |
CN109977526A (zh) * | 2019-03-21 | 2019-07-05 | 北京航空航天大学 | 一种基于三维cst技术的调整机翼有限元模型的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008096478A1 (ja) * | 2007-02-07 | 2008-08-14 | Nsk Ltd. | Pdmシステムとの連携機能を有するfem解析処理システム |
CN101697176A (zh) * | 2009-10-29 | 2010-04-21 | 西北工业大学 | 多组件结构系统布局优化设计方法 |
US20100204965A1 (en) * | 2009-01-30 | 2010-08-12 | Aisin Aw Co., Ltd. | Apparatus, method, and program for acoustic characteristics optimization model analysis |
-
2013
- 2013-08-27 CN CN201310376821.0A patent/CN103455670B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008096478A1 (ja) * | 2007-02-07 | 2008-08-14 | Nsk Ltd. | Pdmシステムとの連携機能を有するfem解析処理システム |
US20100204965A1 (en) * | 2009-01-30 | 2010-08-12 | Aisin Aw Co., Ltd. | Apparatus, method, and program for acoustic characteristics optimization model analysis |
CN101697176A (zh) * | 2009-10-29 | 2010-04-21 | 西北工业大学 | 多组件结构系统布局优化设计方法 |
Non-Patent Citations (3)
Title |
---|
卢其辉: "复杂结构的动态响应有限元分析技术和应用", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
彭慧莲 等: "双面约束多点摩擦多体系统的建模和数值方法", 《力学学报》 * |
曲晓锐 等: "多点约束(MPC)法与换热器整体有限元分析", 《设计计算》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104484527A (zh) * | 2014-12-17 | 2015-04-01 | 清华大学 | 一种离散结构拓扑优化过程中均布载荷自动动态修改方法 |
CN104484527B (zh) * | 2014-12-17 | 2017-12-19 | 清华大学 | 一种离散结构拓扑优化过程中均布载荷自动动态修改方法 |
CN104732043A (zh) * | 2015-04-14 | 2015-06-24 | 河北省电力勘测设计研究院 | 变电架构的设计方法 |
CN104732043B (zh) * | 2015-04-14 | 2017-09-12 | 河北省电力勘测设计研究院 | 变电架构的设计方法 |
CN104992012B (zh) * | 2015-06-26 | 2017-11-03 | 华晨汽车集团控股有限公司 | 汽车后门刚度分析方法 |
CN104992012A (zh) * | 2015-06-26 | 2015-10-21 | 华晨汽车集团控股有限公司 | 汽车后门刚度分析方法 |
CN105426640A (zh) * | 2015-12-28 | 2016-03-23 | 西北工业大学 | 基于罚函数的多组件结构系统布局优化设计方法 |
CN105426640B (zh) * | 2015-12-28 | 2018-11-02 | 西北工业大学 | 基于罚函数的多组件结构系统布局优化设计方法 |
CN106295028A (zh) * | 2016-08-16 | 2017-01-04 | 西北工业大学 | 一种局部结构动力学建模方法及装置 |
CN106295028B (zh) * | 2016-08-16 | 2019-12-03 | 西北工业大学 | 一种局部结构动力学建模方法及装置 |
CN109145495A (zh) * | 2018-09-11 | 2019-01-04 | 吉林大学 | 固定分型线双模铸造件多组件的拓扑优化方法 |
CN109543207A (zh) * | 2018-09-11 | 2019-03-29 | 吉林大学 | 考虑变分型线实现双模铸造件多组件设计的方法 |
CN109766564A (zh) * | 2018-10-31 | 2019-05-17 | 中国飞机强度研究所 | 考虑组件保形约束的多组件结构系统布局优化设计方法 |
CN109977526A (zh) * | 2019-03-21 | 2019-07-05 | 北京航空航天大学 | 一种基于三维cst技术的调整机翼有限元模型的方法 |
CN109977526B (zh) * | 2019-03-21 | 2021-01-05 | 北京航空航天大学 | 一种基于三维cst技术的调整机翼有限元模型的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103455670B (zh) | 2016-02-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103455670B (zh) | 基于多点约束的多组件结构系统布局优化设计方法 | |
Wang et al. | Concurrent design of hierarchical structures with three-dimensional parameterized lattice microstructures for additive manufacturing | |
Wang et al. | Structure-material integrated design by level sets | |
Coenen et al. | Computational homogenization for heterogeneous thin sheets | |
Patel et al. | Crashworthiness design using topology optimization | |
Wang et al. | Shape equilibrium constraint: a strategy for stress-constrained structural topology optimization | |
Wang et al. | A coupled smoothed finite element method (S-FEM) for structural-acoustic analysis of shells | |
CN104765922A (zh) | 基于保形约束的悬臂梁结构拓扑优化设计方法 | |
Zhou et al. | Band gap characteristics of periodically stiffened-thin-plate based on center-finite-difference-method | |
Kang et al. | A method using successive iteration of analysis and design for large-scale topology optimization considering eigenfrequencies | |
Kulkarni et al. | A review of research and recent trends in analysis of composite plates | |
CN106940739A (zh) | 一种振动对机翼共形相控阵天线电性能影响的快速预测方法 | |
Zhao et al. | Concurrent topology optimization with uniform microstructure for minimizing dynamic response in the time domain | |
Changizi et al. | Robust topology optimization of frame structures under geometric or material properties uncertainties | |
CN103615054A (zh) | 一种基于区格剪切变形的屈曲约束支撑布置方法 | |
Xing et al. | An eigenelement method of periodical composite structures | |
Li et al. | A peridynamic model for the nonlinear static analysis of truss and tensegrity structures | |
Yan et al. | Structure/material concurrent optimization of lattice materials based on extended multiscale finite element method | |
Oktay et al. | Three-dimensional structural topology optimization of aerial vehicles under aerodynamic loads | |
Ma et al. | Generative design of stiffened plates based on homogenization method | |
Xiao et al. | Estimation of effective properties of composite sandwich panels with negative Poisson's ratio by using variational asymptotic multiscale method | |
Li et al. | Dynamic stiffness design of plate/shell structures using explicit topology optimization | |
Wu et al. | Effects of higher-order global–local shear deformations on bending, vibration and buckling of multilayered plates | |
Liu et al. | A hybrid direct FE2 method for modeling of multiscale materials and structures with strain localization | |
Awrejcewicz et al. | Design of composite structures with extremal elastic properties in the presence of technological constraints |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210304 Address after: Room 101, building 6, area B, civil military integration innovation park, No. 699, aerospace East Road, Xi'an national civil aerospace industry base, Shaanxi 710000 Patentee after: Shaanxi Aoxiang Tianyu New Material Technology Co.,Ltd. Address before: 710072 No. 127 Youyi West Road, Shaanxi, Xi'an Patentee before: Northwestern Polytechnical University |
|
TR01 | Transfer of patent right |