CN103970968B - 一种解决离散杆件结构优化中对称性约束的分组方法 - Google Patents
一种解决离散杆件结构优化中对称性约束的分组方法 Download PDFInfo
- Publication number
- CN103970968B CN103970968B CN201410230946.7A CN201410230946A CN103970968B CN 103970968 B CN103970968 B CN 103970968B CN 201410230946 A CN201410230946 A CN 201410230946A CN 103970968 B CN103970968 B CN 103970968B
- Authority
- CN
- China
- Prior art keywords
- node
- rod member
- symmetry
- memory element
- discrete
- 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
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 title claims abstract description 11
- 238000003860 storage Methods 0.000 claims abstract description 17
- 238000005516 engineering process Methods 0.000 claims description 10
- 238000013461 design Methods 0.000 description 11
- 238000004364 calculation method Methods 0.000 description 3
- 238000012217 deletion Methods 0.000 description 3
- 230000037430 deletion Effects 0.000 description 3
- 238000005303 weighing Methods 0.000 description 3
- 238000003672 processing method Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种解决离散杆件结构优化中对称性约束的分组方法,包括以下步骤:1)初始化对称面/轴的系数;2)计算节点到对称面/轴的距离d;3)将d<0、d>0和d=0的节点分别储存在第一、第二和第三存储单元中;4)计算在第一存储单元中的节点关于对称面/轴对称点的坐标;5)寻找第二存储单元中与对称点的距离小于给定极小值的节点;6)将节点对储存在第四存储单元中;7)判断第一存储单元中的所有节点是否均已经历步骤4)~6);8)将杆件编号储存到第五存储单元中,计算第五存储单元中包含的杆件数目,并初始化索引参数;9)在第五存储单元中搜索排序第一的杆件两端的节点,并从第四存储单元中搜索节点的对称节点;10)杆件连接的两个节点是否为对称节点;11)判断索引参数是否为杆件数目;12)判断第五存储单元是否为空集。
Description
技术领域
本发明涉及一种根据对称性约束要求进行分组的方法,特别是关于一种解决离散杆件结构优化中对称性约束的分组方法。
背景技术
在实际工程以及生活中,桁架、梁架等离散杆件结构的应用非常广泛,为提高离散杆件结构的综合性能,相关学者将优化方法引入到离散杆件设计中。目前,离散杆件结构优化是结构优化领域非常重要的研究方向之一。根据优化的内容,离散杆件优化可分为拓扑优化、尺寸优化和形状优化,其中拓扑优化的拓扑设计变量取值对应杆件的存在或删除,尺寸优化的尺寸设计变量对应杆件横截面的尺寸参数,而形状优化的形状设计变量决定节点的位置。
为了降低制造、安装成本,提高优化效率,离散杆件结构进行优化时,通常添加对称性约束,为了满足对称性约束,对称杆件的拓扑设计变量、尺寸设计变量取值需要相同,对称节点的形状设计变量取值也需要相同。常见的处理方法是:离散杆件结构优化前人为地将对称的杆件或者节点分为一组,同一组中杆件的拓扑设计变量或尺寸设计变量取相同的值,同一组中节点的形状设计变量取相同的值。对于实际工程以及生活中大量存在的复杂的离散杆件结构而言,上述处理方法存在工作量大、分组效率低和容易出错等问题。
发明内容
针对上述问题,本发明的目的是提供一种能够减少工作量、分组准确全面地解决离散杆件结构优化中对称性约束的分组方法。
为实现上述目的,本发明采取以下技术方案:一种解决离散杆件结构优化中对称性约束的分组方法,包括以下步骤:1)确定离散杆件结构的坐标系并在该坐标系内建立对称面或对称轴的方程f(x,y,z);2)计算离散杆件结构中具有对称性约束的节点i到对称面或对称轴的距离d;3)将d<0的节点储存在第一存储单元中,将d>0的节点储存在第二存储单元中,将d=0的节点储存在第三存储单元中;4)计算在第一存储单元中的节点i关于对称面或对称轴对称的点i'的坐标xi'、yi'、zi';5)在第二存储单元中,与点i'的距离小于给定极小值ε的节点k是离散杆件结构中的节点i关于对称面或对称轴的对称节点,而第三存储单元中的节点关于对称面或对称轴的对称节点是节点本身;6)两两对称的节点构成节点对,将节点对储存在第四存储单元中,更新第四存储单元,使其包含对称节点对(i,k);7)判断第一存储单元中的所有节点是否均已经历所述步骤4)~6),如果否,则返回所述步骤4)重新循环,否则,则继续下一步;8)将离散杆件结构中所有的杆件编号储存到第五存储单元中,计算第五存储单元中包含的杆件数目Ne,并初始化索引参数j,即令j=2;9)第五存储单元中排序第一的杆件为Q(1),杆件Q(1)两端的节点为i1和i2,从第四存储单元中搜索i1和i2的对称节点i′1和i′2;10)杆件Q(j)连接的两个节点是否为i′1和i′2,如果是,则杆件Q(j)为杆件Q(1)的对称杆,杆件Q(1)和Q(j)为一组,杆件Q(1)和Q(j)为一组,将其储存到第六存储单元中,并从第五存储单元中删除杆件Q(1)和Q(j),转至所述步骤12);否则,则转至所述步骤11);11)判断j是否为Ne,如果是,则杆件Q(1)不存在对称杆件,单独为一组,放入第六存储单元中,并将杆件Q(1)从第五存储单元中删除,转至所述步骤12);否则j=j+1,转至所述步骤10);12)判断第五存储单元是否为空集,如果第五存储单元为空集则所有离散杆件结构均已完成分组,整个过程结束;否则,令j=2,更新第五存储单元包含的杆件数目Ne,转至所述步骤9),重新循环。
所述步骤1)中,对称面或对称轴的方程f(x,y,z)为:
f(x,y,z)=Ax+By+Cz+D=0
式中,如果离散杆件结构是二维的,则C=0。
所述步骤2)中,距离d的计算公式为:
式中,xi,yi,zi为节点i的坐标值。
所述步骤4)中,节点i关于对称面或对称轴对称的点i'的坐标xi'、yi'、zi'的计算公式为:
本发明由于采取以上技术方案,其具有以下优点:1、由于本发明是根据用户提供的对称面或对称轴的信息,采用解析法寻找离散杆件结构中对称的杆件和节点,不需要人为地将具有对称性约束的杆件或者节点进行分组,从而使得分组简单且准确。2、采用本发明方法,可以非常简单方便地寻找到离散杆件结构中每一组具有对称性约束的杆件和节点,因此有效地保证了分组的全面性。3、本发明由于不是对具体的设计变量进行分组,而是将杆件和节点作为分组对象,因此可以适用拓扑优化、尺寸优化和形状优化。本发明可以广泛用于解决离散杆件结构优化过程中具有对称性约束杆件或者节点的分组。
附图说明
图1是本发明的流程示意图
图2是平面六节点杆件的结构示意图
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明提供的解决离散杆件结构优化中对称性约束的分组方法,包括以下步骤:
1)确定离散杆件结构的坐标系并在该坐标系内建立对称面或对称轴的方程:
f(x,y,z)=Ax+By+Cz+D=0 (1)
式中,如果离散杆件结构是二维的,则C=0。
2)计算离散杆件结构中具有对称性约束的节点i到对称面或对称轴的距离d:
式中,xi,yi,zi为节点i的坐标值。
3)将d<0的节点储存在节点集合NodeSet1中,将d>0的节点储存在节点集合NodeSet2中,将d=0的节点储存在节点集合NodeSet3。
4)计算在NodeSet1中的节点i关于对称面或对称轴对称的点i'的坐标xi'、yi'、zi':
5)在NodeSet2中,与点i'的距离小于给定极小值ε的节点k是离散杆件结构中的节点i关于对称面或对称轴的对称节点,而NodeSet3中的节点关于对称面或对称轴的对称节点是节点本身。
6)两两对称的节点构成一节点对,将节点对储存在节点对集合SymNode中,更新SymNode,使其包含节点对(i,k)。
7)判断NodeSet1中的所有节点i是否均已经历步骤4)~6),如果否,则返回步骤4)重新循环,否则,则继续下一步。
8)将离散杆件结构中所有的杆件编号储存到集合Q中,计算集合Q中包含的杆件数目Ne,并初始化索引参数j,即令j=2。
9)集合Q中排序第一的杆件为Q(1),杆件Q(1)两端的节点为i1和i2,从SymNode中搜索i1和i2的对称节点i′1和i′2。
10)判断杆件Q(j)连接的两个节点是否为i′1和i′2,如果是,则杆件Q(j)为杆件Q(1)的对称杆,杆件Q(1)和Q(j)为一组,将其储存到集合SymElm中,并从集合Q中删除杆件Q(1)和Q(j),转至步骤12);否则,则转至步骤11)。
11)判断j是否为Ne,如果是,则杆件Q(1)不存在对称杆件,单独为一组,放入集合SymElm中,并将杆件Q(1)从集合Q中删除,转至步骤12);否则j=j+1,转至步骤10)。
12)判断集合Q是否为空集,如果集合Q为空集则所有离散杆件结构均已完成分组,整个过程结束;否则,令j=2,更新集合Q包含的杆件数目Ne,转至步骤9),重新循环。
下面列举具体实施例,以对本发明有进一步的了解。
实施例:如图2所示,以平面六节点杆件结构为例,平面六节点杆件包含6个节点1、2、3、4、5、6,11根杆件①、②、③、④、⑤、⑥、⑦、⑧、⑨、⑩、若以节点4作为原点,杆件⑥所在的直线作为X轴,杆件①所在的直线作为Y轴,建立坐标系,则节点4的坐标为(0,0),其它各节点的坐标分别为(0,1)、(1,0)、(1,1)、(2,0)、(2,1)。要求该结构优化过程中所有杆件和节点均关于直线x=1对称,使用本发明方法对平面六节点杆件结构进行分组,包括以下步骤:
1)建立对称轴方程:
f(x)=x-1=0 (4)
2)计算平面六节点杆件结构中具有对称性要求的节点i(i=1,2,3,4,5,6)到对称轴f(x)的距离d:
d=xi-1 (5)
式中,xi为节点i在X轴上的坐标,计算结果见表1:
表1
节点编号 | 到对称轴的距离d |
1 | -1 |
2 | 0 |
3 | 1 |
4 | -1 |
5 | 0 |
6 | 1 |
3)将d<0的节点储存在NodeSet1中,NodeSet1={1,4};将d>0的节点储存在NodeSet2中,NodeSet2={3,6};将d=0的节点储存在NodeSet3中,NodeSet3={2,5}。
4)计算NodeSet1中的节点i=1,4关于对称轴对称的点i'的坐标xi'、yi':
式中,xi、yi为节点i的坐标值,计算结果如下:
表2
节点编号 | 关于对称轴对称的点坐标 |
1 | (2,1) |
4 | (2,0) |
5)在NodeSet2中,与节点i'距离小于给定极小值ε=0.001的节点k是平面六节点杆件结构中的节点i关于对称轴的对称节点,因此节点1与3对称,节点4与6对称,认为NodeSet3中的节点关于对称轴的对称节点为该节点本身,因此节点2和节点5的对称节点为本身。
6)两两对称的节点构成节点对,将节点对储存在节点对集合SymNode中,更新SymNode,SymNode={1,3;4,6;2,2;5,5}。
7)判断NodeSet1中的所有节点是否均已经历步骤4)~6),如果否,则返回步骤4)重新循环,否则,则继续下一步。
8)将平面六节点杆件结构中的所有杆件编号添加至集合Q,Q={1,2,3,4,5,6,7,8,9,10,11},计算集合Q中包含的杆件数目Ne,Ne=11。并初始化索引参数j,即令j=2。
下面以第一次循环为例说明步骤9)~12):
9)集合Q中排序第一的杆件为Q(1),杆件Q(1)两端的节点为i1和i2,从SymNode中搜索i1和i2的对称节点i′1和i′2,第一次循环时,Q(1)=1,(i1,i2)=(1,4),(i′1,i′2)=(3,6)。
10)判断杆件Q(j)连接的两个节点是否为i′1和i′2,如果是,则杆件Q(j)为杆件Q(1)的对称杆,杆件Q(1)和Q(j)为一组,将其储存到集合SymElm中,并从集合Q中删除杆件Q(1)和Q(j),转至步骤12);否则,则转至步骤11)。
第一次循环的第一次判断中j=2,Q(j)=Q(2)=2,连接的节点为{2,5}≠{3,6},转至步骤11)。第一次循环的第二次判断中j=3,Q(j)=Q(3)=3,连接的节点为{3,6}={3,6},因此杆件③为杆件①的对称杆件,杆件③和杆件①为一组,将其放入集合SymElm中,SymElm={1,3},并从集合Q中删除杆件③和杆件①,更新集合Q={2,4,5,6,7,8,9,10,11},转至步骤12)。
11)判断j是否为Ne,如果是,则杆件Q(1)不存在对称杆件,单独为一组,放入集合SymElm中,并将杆件Q(1)从集合Q中删除,转至步骤12);否则j=j+1,转至步骤10)。
第一次循环第一次判断时,j=2≠Ne,令j=j+1=3,转至步骤10),进行第二次判断。
12)判断集合Q是否为空集,如果Q为空集则所有杆件已完成分组,结束;否则令j=2,更新集合Q包含的杆件数目Ne,转至步骤9)。
第一次循环时,Q非空,令j=2,计算得Ne=9,第一次循环结束,转至步骤9)开始第二次循环。
采用本发明方法,平面六节点杆件结构分组情况如表3所示:
表3
组号 | 包含的杆件 |
1 | 1,3 |
2 | 2 |
3 | 4,5 |
4 | 6,7 |
5 | 8,9 |
6 | 10,11 |
上述各实施例仅用于说明本发明,其中各参数的表示、各集合的命名等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (4)
1.一种解决离散杆件结构优化中对称性约束的分组方法,包括以下步骤:
1)确定离散杆件结构的坐标系并在该坐标系内建立对称面或对称轴的方程f(x,y,z);
2)计算离散杆件结构中具有对称性约束的节点i到对称面或对称轴的距离d;
3)将d<0的节点储存在第一存储单元中,将d>0的节点储存在第二存储单元中,将d=0的节点储存在第三存储单元中;
4)计算在第一存储单元中的节点i关于对称面或对称轴对称的点i′的坐标xi′、yi′、zi′;
5)在第二存储单元中,与点i′的距离小于给定极小值ε的节点k是离散杆件结构中的节点i关于对称面或对称轴的对称节点,而第三存储单元中的节点关于对称面或对称轴的对称节点是节点本身;
6)两两对称的节点构成节点对,将节点对储存在第四存储单元中,更新第四存储单元,使其包含对称节点对(i,k);
7)判断第一存储单元中的所有节点是否均已经历所述步骤4)~6),如果否,则返回所述步骤4)重新循环,否则,则继续下一步;
8)将离散杆件结构中所有的杆件编号储存到第五存储单元中,计算第五存储单元中包含的杆件数目Ne,并初始化索引参数j,即令j=2;
9)第五存储单元中排序第一的杆件为Q(1),杆件Q(1)两端的节点为i1和i2,从第四存储单元中搜索i1和i2的对称节点i′1和i′2;
10)杆件Q(j)连接的两个节点是否为i′1和i′2,如果是,则杆件Q(j)为杆件Q(1)的对称杆,杆件Q(1)和Q(j)为一组,杆件Q(1)和Q(j)为一组,将其储存到第六存储单元中,并从第五存储单元中删除杆件Q(1)和Q(j),转至所述步骤12);否则,则转至所述步骤11);
11)判断j是否为Ne,如果是,则杆件Q(1)不存在对称杆件,单独为一组,放入第六存储单元中,并将杆件Q(1)从第五存储单元中删除,转至所述步骤12);否则j=j+1,转至所述步骤10);
12)判断第五存储单元是否为空集,如果第五存储单元为空集则所有离散杆件结构均已完成分组,整个过程结束;否则,令j=2,更新第五存储单元包含的杆件数目Ne,转至所述步骤9),重新循环。
2.如权利要求1所述的一种解决离散杆件结构优化中对称性约束的分组方法,其特征在于:所述步骤1)中,对称面或对称轴的方程f(x,y,z)为:
f(x,y,z)=Ax+By+Cz+D=0
式中,如果离散杆件结构是二维的,则C=0。
3.如权利要求2所述的一种解决离散杆件结构优化中对称性约束的分组方法,其特征在于:所述步骤2)中,距离d的计算公式为:
式中,xi,yi,zi为节点i的坐标值。
4.如权利要求3所述的一种解决离散杆件结构优化中对称性约束的分组方法,其特征在于:所述步骤4)中,节点i关于对称面或对称轴对称的点i′的坐标xi′、yi′、zi′的计算公式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410230946.7A CN103970968B (zh) | 2014-05-28 | 2014-05-28 | 一种解决离散杆件结构优化中对称性约束的分组方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410230946.7A CN103970968B (zh) | 2014-05-28 | 2014-05-28 | 一种解决离散杆件结构优化中对称性约束的分组方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103970968A CN103970968A (zh) | 2014-08-06 |
CN103970968B true CN103970968B (zh) | 2017-01-25 |
Family
ID=51240458
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410230946.7A Expired - Fee Related CN103970968B (zh) | 2014-05-28 | 2014-05-28 | 一种解决离散杆件结构优化中对称性约束的分组方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103970968B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108520137B (zh) * | 2018-04-02 | 2020-08-04 | 华中科技大学 | 一种适用于空间桁架结构的平面节点耦合处理方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101409439A (zh) * | 2008-11-19 | 2009-04-15 | 中国电力工程顾问集团公司 | 一种确定重冰区输电线路杆塔导地线布置的方法 |
CN101582130A (zh) * | 2009-05-27 | 2009-11-18 | 清华大学 | 一种改进遗传算法结构优化效率的方法 |
CN101582131A (zh) * | 2009-05-27 | 2009-11-18 | 清华大学 | 一种离散结构遗传优化的二维编码及解码方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10009195A1 (de) * | 2000-02-26 | 2001-08-30 | Topack Verpacktech Gmbh | Verfahren und Einrichtung zum Verkleben von Verpackungsteilen |
DE10044577A1 (de) * | 2000-09-08 | 2002-03-21 | Hauni Maschinenbau Ag | Verfahren und Vorrichtung zum Trennen einer doppeltbreiten Ausgangs-Hüllmaterialbahn in zwei Teilbahnen |
WO2008115339A1 (en) * | 2007-03-15 | 2008-09-25 | White Nicholas R | Open-ended electromagnetic corrector assembly and method for deflecting, focusing, and controlling the uniformity of a traveling ion beam |
-
2014
- 2014-05-28 CN CN201410230946.7A patent/CN103970968B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101409439A (zh) * | 2008-11-19 | 2009-04-15 | 中国电力工程顾问集团公司 | 一种确定重冰区输电线路杆塔导地线布置的方法 |
CN101582130A (zh) * | 2009-05-27 | 2009-11-18 | 清华大学 | 一种改进遗传算法结构优化效率的方法 |
CN101582131A (zh) * | 2009-05-27 | 2009-11-18 | 清华大学 | 一种离散结构遗传优化的二维编码及解码方法 |
Non-Patent Citations (3)
Title |
---|
基于离散变量网架结构优化设计研究;蔡建良;《中国优秀博硕士学位论文全文数据库 (硕士) (月刊)》;20061115;第2006年卷(第11期);全文 * |
大客车侧翻安全性多目标优化;范子杰等;《***仿真学报》;20120531;第24卷(第5期);全文 * |
采用分段样条插值的半离散方法分析薄壁杆件;李华煜等;《工程力学》;19920930;第9卷(第3期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103970968A (zh) | 2014-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102163218B (zh) | 基于图索引的图数据库关键词邻近搜索方法 | |
CN107844676A (zh) | 一种基于多性能约束的结构拓扑优化设计方法 | |
CN103645725B (zh) | 一种机器人示教轨迹规划方法和*** | |
CN103810244A (zh) | 一种基于数据分布的分布式数据存储***的扩容方法 | |
CN105844364A (zh) | 基于启发函数的服务机器人最优路径规划方法 | |
CN107862617A (zh) | 一种基于用户综合相似度的微博社区划分方法 | |
CN105550368A (zh) | 一种高维数据的近似最近邻检索方法及检索*** | |
CN104933624A (zh) | 复杂网络的社团发现方法及社团重要节点发现方法 | |
CN107330589A (zh) | 卫星网络协调风险的定量化评估方法及*** | |
CN104268629A (zh) | 一种基于先验信息和网络固有信息的复杂网络社区检测方法 | |
CN108388603A (zh) | 基于Spark框架的分布式概要数据结构的构建方法及查询方法 | |
CN107066663A (zh) | 一种基于满应力约束准则的桁架结构非概率可靠性拓扑优化方法 | |
CN103970968B (zh) | 一种解决离散杆件结构优化中对称性约束的分组方法 | |
CN106919998B (zh) | 一种基于aoe网络的物流模型及方法 | |
CN103559209A (zh) | 一种Voronoi Diagram与虚拟网格结合的高效空间最近邻查询方法 | |
CN103346556A (zh) | 一种配电网环路快速定位方法 | |
CN105045967A (zh) | 基于集团度的复杂网络重要节点排序方法及模型演化方法 | |
CN103337856B (zh) | 一种配电网网络重构的开关分组方法 | |
CN102819581B (zh) | 一种地理信息***拓扑集中的多边形链生成方法 | |
CN107038212A (zh) | 一种基于蒙特卡洛方法逆向求解PageRank的算法 | |
CN103761270A (zh) | 一种串数据词典的有序构造及检索方法 | |
CN106991195A (zh) | 一种分布式的子图枚举方法 | |
CN103346802B (zh) | Qc-ldpc码的构造方法 | |
CN101556604A (zh) | 一种面向复杂数据仓库环境的优化策略自动生成方法 | |
Flaherty et al. | Distributed octree data structures and local refinement method for the parallel solution of three-dimensional conservation laws |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170125 Termination date: 20190528 |