CN105718684B - 一种主动控制结构动刚度的记忆合金驱动件优化布局方法 - Google Patents
一种主动控制结构动刚度的记忆合金驱动件优化布局方法 Download PDFInfo
- Publication number
- CN105718684B CN105718684B CN201610049017.5A CN201610049017A CN105718684B CN 105718684 B CN105718684 B CN 105718684B CN 201610049017 A CN201610049017 A CN 201610049017A CN 105718684 B CN105718684 B CN 105718684B
- Authority
- CN
- China
- Prior art keywords
- alloy wire
- memory alloy
- matrix
- based structures
- stiffness matrix
- 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
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]
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)
- Feedback Control In General (AREA)
Abstract
一种主动控制结构动刚度的记忆合金驱动件优化布局方法,利用刚度等效法计算出受到记忆合金丝拉力作用下的结构模态特性,然后利用MMA方法对记忆合金丝布局位置进行迭代优化,从而得到记忆合金丝布局的最优解,使用本方法进行设计时,记忆合金丝的布局不必受限于主观设定有限元网格,能够自由布置,从而比传统设计的结果更接近最优解。
Description
技术领域
本发明属于结构的模态特性、动刚度控制技术领域,具体涉及一种主动控制结构动刚度的记忆合金驱动件优化布局方法。
技术背景
结构的模态特性由结构本身的质量分布以及动刚度共同决定,而结构的动刚度由自身属性与所受应力的大小、位置共同决定,所以通过合理布置加载力的位置和大小可以实现对结构动刚度及模态特性的主动控制。工程中通常利用记忆合金丝作为驱动元件控制加载力的大小,并使用遗传算法或其改进算法确定力的加载位置,但是使用这种方法时,加载的力必须位于人为设定的有限元网格节点上,因此寻优的范围被大大缩减,很难得到最优结果。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供一种主动控制结构动刚度的记忆合金驱动件优化布局方法,能够确定记忆合金丝的最优布局,从而得到最好的刚度控制效果。
为达到上述目标,本发明采取的技术方案为:
一种主动控制结构动刚度的记忆合金驱动件优化布局方法,包括以下步骤:
1)、确定目标函数及设计变量:
设定优化目标函数为
其中C为设计变量,它是由所有记忆合金丝端点坐标构成的向量;为第b次迭代后待优化结构的a阶模态固有频率,为记忆合金丝初始布局时的待优化结构a阶模态固有频率;
2)、建立模型:
在有限元软件中使用shell63单元建立不包含记忆合金丝的待优化结构的有限元模型,该模型称为基结构模型;导出基结构模型的整体刚度矩阵KKG与整体质量矩阵M,假设整个基结构有n个节点,那么矩阵KKG和M为6*n阶矩阵;
3)、布置记忆合金丝的初始位置:
在基结构模型中均匀地布置等长的记忆合金丝,以此作为初始的记忆合金丝布局;初始布局达到均匀稀疏的涵盖整个基结构表面;将记忆合金丝的端点坐标记录到向量C中,并在下面的迭代优化中不断更新向量C;
4)、对记忆合金丝进行设置:
设定每根记忆合金丝只对基结构提供恒定的拉力,且拉力相同,同时不依靠记忆合金丝自身的刚度为基结构模型提供加强作用,其实现方式为:将记忆合金丝通过螺栓连接在基结构模型表面;
5)、求解各个记忆合金丝的应力刚化刚度矩阵并将其组装:
编辑程序得出所有记忆合金丝在基结构模型坐标系下的应力刚化刚度矩阵,这里得到用变量C表示出的单元应力刚化刚度矩阵KK;假定总共有m根记忆合金丝,定义每一根记忆合金丝都是一个独立单元,每个单元有2个节点12个自由度,那么能够得到m个单元应力刚化刚度矩阵KKt(t=1、2…...m)每个矩阵均为12阶方阵;将各个记忆合金丝的单元应力刚化刚度矩阵KKt(t=1、2…...m)依次沿对角线拼接在一起,形成一个12*m阶的总应力刚化刚度矩阵K,其中没有被定义到的元素均置为零;
6)、使用刚度等效法将总应力刚化刚度矩阵K与基结构模型整体刚度矩阵KKG进行融合,从而得到整个结构的总刚度矩阵KBB,其具体方法为:
6.1)求解转换矩阵T:
求解H:
其中:
上式取基结构n个节点中的一个节点i,定义该点坐标为xi、yi、zi;取m个记忆合金丝上所有2*m个端点中的一个节点j,定义该点坐标为xj、yj、zj;dsp值为基结构模型两个相邻节点间距离的2-10倍;
求解H0:
其中:
上式取基结构模型n个节点中的两个节点p和q,定义这两点坐标分别为xp、yp、zp和xq、yq、zq;
求解转换矩阵T:
6.2)、求出总刚度矩阵KBB:KBB=TT×K×T+KKG;
7)、计算待优化结构的模态特性:
利用(KBB-ω2M){φ}=0计算出第b次迭代后的a阶模态的固有频率以及初始记忆合金丝布局时的a阶模态的固有频率
8)、求解目标函数对设计变量C中每一个变量的偏导数;
9)、更新设计变量C:
将目标函数f(C)及其对各个变量的偏导数导入MMA算法,MMA算法将设计变量C更新一次;
10)、迭代更新设计变量C
重复进行步骤5)-步骤9),每次迭代都将所有设计变量更新一次,这组新的设计变量将被用于下一轮迭代之中;最终当目标函数低于初始值且收敛时,停止迭代,此时便得到了记忆合金丝的初步优化布局;
11)后处理:参照初步优化布局进行圆整处理,再根据加工工艺要求以及制造装配要求进一步修改布局,得到最终设计。
本发明的有益效果为:
由于本方法使用了刚度等效法,所以记忆合金丝的布局摆脱了人为建立的有限元网格的束缚,实现了在整个结构任意位置的自由布置,因此能够得到真正的最优布局方案,而非传统优化方法中的受约束的布局。
附图说明
图1为本发明实施例中基结构模型上记忆合金丝的初始布局图。
图2为传统的布局优化方法与应用刚度等效法的布局优化方法之间的对比;其中图a为传统优化方法得到的布局示意图,图b为应用刚度等效法后得到的布局示意图。
具体实施方式
下面结合在500mm*500mm铝板上布置记忆合金丝的实施例对本发明作进一步说明。
一种主动控制结构动刚度的记忆合金驱动件优化布局方法,包括以下步骤:
1)、确定目标函数及设计变量:
本实施例以尽可能提高待优化结构的1阶模态固有频率为目标,所以设定优化目标函数为:
其中C为设计变量,它是由所有记忆合金丝端点坐标构成的向量,本实施例中布置有4根记忆合金丝,所以C为16维向量;为记忆合金丝初始布局时的待优化结构1阶模态固有频率,为第b次迭代后结构的1阶模态固有频率;
2)、建立模型:
在有限元软件中使用shell63单元均匀建立11*11个节点的铝板有限元模型,该模型称为基结构模型,这一基结构模型不包括记忆合金丝;导出基结构模型的整体刚度矩阵KKG与整体质量矩阵M,此时矩阵KKG和M为726阶矩阵;
3)、布置记忆合金丝的初始位置:
如图1所示,在基结构模型中均匀地布置等长的记忆合金丝;将记忆合金丝的端点坐标记录到向量C中,并在下面的迭代优化中不断更新向量C;
4)、对记忆合金丝进行设置:
设定每根记忆合金丝只对基结构提供恒定的拉力,且拉力相同,同时不依靠记忆合金丝自身的刚度为基结构模型提供加强作用,其实现方式为:将记忆合金丝通过螺栓连接在基结构模型表面;
5)、求解各个记忆合金丝的应力刚化刚度矩阵并将其组装:
编辑程序得出所有记忆合金丝在基结构模型坐标系下的单元应力刚化刚度矩阵,这里得到用变量C表示出的单元应力刚化刚度矩阵KK;由于本实施例有4根记忆合金丝,定义每一根记忆合金丝都是一个单元,每个单元有2个节点12个自由度,那么可以得到4个单元应力刚化刚度矩阵KKt(t=1、2、3、4)每个矩阵均为12阶方阵;将各个记忆合金丝的单元刚度矩阵KKt(t=1、2、3、4)依次沿对角线拼接在一起,形成一个48阶的总应力刚化刚度矩阵K,其中没有被定义到的元素均置为零;
6)、使用刚度等效法将总应力刚化刚度矩阵K与基结构模型整体刚度矩阵KKG进行融合,从而得到整个结构的总刚度矩阵KBB;使用这种方法时,尽管记忆合金丝已经脱离了基结构模型的节点,如图2中b图所示(传统方法如图2中a图所示),但其产生的应力刚化刚度矩阵仍能等效到基结构模型节点上,从而满足有限元方法的要求;其具体方法为:
6.1)求解转换矩阵T:
求解H:
其中:
上式取基结构121个节点中的一个节点i,定义该点坐标为xi、yi、zi;取4个记忆合金丝上所有8个端点中的一个节点j,定义该点坐标为xj、yj、zj;dsp值为基结构模型两个相邻节点间距离的3倍;
求解H0:
其中:
上式取基结构模型121个节点中的两个节点p和q,定义这两点坐标分别为xp、yp、zp和xq、yq、zq;
求解转换矩阵T:
6.2)、求出总刚度矩阵KBB:KBB=TT×K×T+KKG;
7)、计算待优化结构的模态特性:
利用(KBB-ω2M){φ}=0计算出第b次迭代后的1阶模态的固有频率以及初始记忆合金丝布局时的1阶模态的固有频率
8)、求解目标函数对设计变量C中每一个变量的偏导数;
9)、更新设计变量C:
将目标函数f(C)及其对各个变量的偏导数导入MMA算法,MMA算法将设计变量C更新一次;
10)、迭代更新设计变量C
重复进行步骤5)-步骤9),每次迭代都将所有设计变量更新一次,这组新的设计变量将被用于下一轮迭代之中;迭代到200次后停止,此时便得到了记忆合金丝的初步优化布局;
11)、后处理:参照初步优化布局进行圆整处理,再根据加工工艺要求以及制造装配要求进一步修改布局,得到最终设计。
Claims (1)
1.一种主动控制结构动刚度的记忆合金驱动件优化布局方法,其特征在于,包括以下步骤:
1)、确定目标函数及设计变量:
设定优化目标函数为
其中C为设计变量,它是由所有记忆合金丝端点坐标构成的向量;为第b次迭代后待优化结构的a阶模态固有频率,为记忆合金丝初始布局时的待优化结构a阶模态固有频率;
2)、建立模型:
在有限元软件中使用shell63单元建立不包含记忆合金丝的待优化结构的有限元模型,该模型称为基结构模型;导出基结构模型的整体刚度矩阵KKG与整体质量矩阵M,假设整个基结构有n个节点,那么矩阵KKG和M为6*n阶矩阵;
3)、布置记忆合金丝的初始位置:
在基结构模型中均匀地布置等长的记忆合金丝,以此作为初始的记忆合金丝布局;初始布局达到均匀稀疏的涵盖整个基结构表面;将记忆合金丝的端点坐标记录到向量C中,并在下面的迭代优化中不断更新向量C;
4)、对记忆合金丝进行设置:
设定每根记忆合金丝只对基结构提供恒定的拉力,且拉力相同,同时不依靠记忆合金丝自身的刚度为基结构模型提供加强作用,其实现方式为:将记忆合金丝通过螺栓连接在基结构模型表面;
5)、求解各个记忆合金丝的应力刚化刚度矩阵并将其组装:
编辑程序得出所有记忆合金丝在基结构模型坐标系下的应力刚化刚度矩阵,这里得到用变量C表示出的单元应力刚化刚度矩阵KK;假定总共有m根记忆合金丝,定义每一根记忆合金丝都是一个独立单元,每个单元有2个节点12个自由度,那么能够得到m个单元应力刚化刚度矩阵KKt,其中t=1、2…...、m,每个矩阵均为12阶方阵;将各个记忆合金丝的单元应力刚化刚度矩阵KKt依次沿对角线拼接在一起,形成一个12*m阶的总应力刚化刚度矩阵K,其中没有被定义到的元素均置为零;
6)、使用刚度等效法将总应力刚化刚度矩阵K与基结构模型整体刚度矩阵KKG进行融合,从而得到整个结构的总刚度矩阵KBB,其具体方法为:
6.1)求解转换矩阵T:
转换矩阵T由两个矩阵点乘得到,分别记为矩阵H和矩阵H0,即
求解H:
其中:
其中i=1、2、......、n;j=1、2、......、2×m;
上式取基结构n个节点中的一个节点i,定义该点坐标为xi、yi、zi;取m个记忆合金丝上所有2*m个端点中的一个节点j,定义该点坐标为xj、yj、zj;dsp值为基结构模型两个相邻节点间距离的2-10倍;
求解H0:
其中:
其中p=1、2、......、n;q=1、2、......、n;
上式取基结构模型n个节点中的两个节点p和q,定义这两点坐标分别为xp、yp、zp和xq、yq、zq;
求解转换矩阵T:
6.2)、求出总刚度矩阵KBB:KBB=TT×K×T+KKG;
7)、计算待优化结构的模态特性:
利用(KBB-ω2M){φ}=0计算出第b次迭代后的a阶模态的固有频率以及初始记忆合金丝布局时的a阶模态的固有频率
8)、求解目标函数对设计变量C中每一个变量的偏导数;
9)、更新设计变量C:
将目标函数f(C)及其对各个变量的偏导数导入MMA算法,MMA算法将设计变量C更新一次;
10)、迭代更新设计变量C
重复进行步骤5)-步骤9),每次迭代都将所有设计变量更新一次,这组新的设计变量将被用于下一轮迭代之中;最终当目标函数低于初始值且收敛时,停止迭代,此时便得到了记忆合金丝的初步优化布局;
11)后处理:参照初步优化布局进行圆整处理,再根据加工工艺要求以及制造装配要求进一步修改布局,得到最终设计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610049017.5A CN105718684B (zh) | 2016-01-25 | 2016-01-25 | 一种主动控制结构动刚度的记忆合金驱动件优化布局方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610049017.5A CN105718684B (zh) | 2016-01-25 | 2016-01-25 | 一种主动控制结构动刚度的记忆合金驱动件优化布局方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105718684A CN105718684A (zh) | 2016-06-29 |
CN105718684B true CN105718684B (zh) | 2018-08-07 |
Family
ID=56154766
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610049017.5A Active CN105718684B (zh) | 2016-01-25 | 2016-01-25 | 一种主动控制结构动刚度的记忆合金驱动件优化布局方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105718684B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106650155B (zh) * | 2016-12-30 | 2020-07-10 | 中国机械工业集团有限公司 | 一种大规模振动工程控制装置优化布置方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104850698A (zh) * | 2015-05-15 | 2015-08-19 | 西安交通大学 | 一种考虑测量与调整工艺的精密机床公差设计方法 |
CN104899376A (zh) * | 2015-06-09 | 2015-09-09 | 西安交通大学 | 一种面向低碳制造的焊接箱型梁智能优化设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0507618D0 (en) * | 2005-04-15 | 2005-05-25 | Lms Internat Nv | Method and system for dynamic analysis of complex systems |
-
2016
- 2016-01-25 CN CN201610049017.5A patent/CN105718684B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104850698A (zh) * | 2015-05-15 | 2015-08-19 | 西安交通大学 | 一种考虑测量与调整工艺的精密机床公差设计方法 |
CN104899376A (zh) * | 2015-06-09 | 2015-09-09 | 西安交通大学 | 一种面向低碳制造的焊接箱型梁智能优化设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105718684A (zh) | 2016-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105354388A (zh) | 一种加强筋的生长式拓扑优化设计方法 | |
CN108846472A (zh) | 一种自适应遗传粒子群混合算法的优化方法 | |
CN107391824B (zh) | 增材制造中自支撑结构的拓扑优化设计方法 | |
CN106934147A (zh) | 一种空间桁架索网天线等几何找形方法 | |
CN108629140B (zh) | 一种基于测地线距离的带孔复合材料结构设计优化方法 | |
CN106026187B (zh) | 一种含分布式电源的配电网重构的方法及系统 | |
CN107844632A (zh) | 基于和声搜索算法的稀布线阵栅瓣抑制方法 | |
CN110232212B (zh) | 基于蝙蝠算法的天线设计优化方法及系统 | |
CN105718684B (zh) | 一种主动控制结构动刚度的记忆合金驱动件优化布局方法 | |
CN101242103A (zh) | 电力系统稳定器参数的智能优化方法 | |
CN110909496B (zh) | 约束阻尼结构两级协同优化设计方法 | |
CN109460577A (zh) | 一种保证制造质量的变刚度复合材料结构优化设计方法 | |
CN106712019B (zh) | 一种有源滤波器优化配置方法 | |
CN106651011B (zh) | 一种基于粒子群算法的渠系优化配水方法 | |
CN116244988A (zh) | 板料旋压多网格法仿真的高质量四边形网格保形构造方法 | |
CN102779211B (zh) | 大型电除尘器烟箱结构优化设计方法 | |
CN113158510B (zh) | 一种索网结构的优化找形方法 | |
CN117035106B (zh) | 量子潮流计算修正方程的分阶计算方法、系统和存储介质 | |
CN107301290B (zh) | 一种电磁带隙电源板的电源分配网络生长式拓扑优化方法 | |
CN110174842A (zh) | 空间太阳能发电卫星在轨组装的分布式振动控制器设计方法 | |
CN108021758A (zh) | 基于有限元的索网天线等张力找形方法 | |
CN113890042B (zh) | 基于有源配电网的潮流计算方法 | |
CN113595528B (zh) | 一种基于幂函数的自适应变步长lms滤波器及其实现方法 | |
CN110188498B (zh) | 一种基于拓扑优化变密度法的最优非设计空间划分方法 | |
CN103258098B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |