CN108595769B - 一种基于优化算法的减振器刚度模拟方法 - Google Patents

一种基于优化算法的减振器刚度模拟方法 Download PDF

Info

Publication number
CN108595769B
CN108595769B CN201810261714.6A CN201810261714A CN108595769B CN 108595769 B CN108595769 B CN 108595769B CN 201810261714 A CN201810261714 A CN 201810261714A CN 108595769 B CN108595769 B CN 108595769B
Authority
CN
China
Prior art keywords
shock absorber
finite element
rigidity
anisotropic
actual
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
Application number
CN201810261714.6A
Other languages
English (en)
Other versions
CN108595769A (zh
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.)
Beihang University
Original Assignee
Beihang University
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 Beihang University filed Critical Beihang University
Priority to CN201810261714.6A priority Critical patent/CN108595769B/zh
Publication of CN108595769A publication Critical patent/CN108595769A/zh
Application granted granted Critical
Publication of CN108595769B publication Critical patent/CN108595769B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于优化算法的减振器刚度模拟方法。首先，根据实际减振器的几何形状建立相应的物理模型；其次，对建立的模型进行有限元网格的划分；第三，将建立的减振器模型赋予正交各向异性材料属性；第四，对减振器一端施加固定约束，另外一端施加单位载荷，进行有限元的求解，从而得到减振器的各向刚度；最后，以正交各向异性材料的各向弹性模量为设计变量，以有限元模拟得到的减振器各向刚度与实际减振器各向刚度最接近为优化目标，进行优化，从而精确模拟实际减振器的刚度。将优化得到的材料属性赋予减振器有限元，进行整体工程结构的有限元计算，可以得到更精确的解。

Description

一种基于优化算法的减振器刚度模拟方法
技术领域
本发明涉及减振器刚度模拟方法，特别涉及基于优化算法的减振器刚度模拟方法。
背景技术
为提高设备的运行寿命和可靠性，减振器设计已成为力学设计中关注的重点。减振器的设计需要重点考虑其刚度特性。结构在外部力作用下抵抗变形的能力即为刚度，即引起单位位移所需要的外力。
在实际工程应用中，往往是选用现有的减振器型号，能够给出减震器的各向刚度。然而再进行结构的静力分析时，却无法直接建模得到与实际减振器各向刚度相同的有限元减振器模型。基于此，本发明提出基于优化算法的减振器刚度模拟方法，采用正交各向异性材料模拟减振器模型，并通过优化得到与实际减振器各向刚度相同的有限元减振器模型。
发明内容
本发明要解决的技术问题是：提出基于优化算法的减振器刚度模拟方法，采用正交各向异性材料模拟减振器模型，并通过优化得到与实际减振器各向刚度相同的有限元减振器模型，为有限元静力计算中减振器的模拟提供一种可选方法。
本发明采用的技术方案为：一种基于优化算法的减振器刚度模拟方法，其实现步骤如下：
步骤一：根据实际减振器的几何形状建立相应的物理模型；
步骤二：对建立的模型进行有限元网格的划分；
步骤三：将建立的减振器模型赋予正交各向异性材料属性；
步骤四：对减振器一端施加固定约束，另外一端施加单位载荷，进行有限元的求解，从而得到减振器的各向刚度；
步骤五：以正交各向异性材料的各向弹性模量为设计变量，以有限元模拟得到的减振器各向刚度与实际减振器各向刚度最接近为优化目标，进行优化，从而精确模拟实际减振器的刚度；
步骤六：将优化得到的正交各向异性材料的各向弹性模量赋予整体结构，进行整体结构的有限元分析，从而实现整体结构的有限元求解。
进一步的，所述步骤一中根据实际减振器的几何形状建立相应的物理模型。
进一步的，所述步骤二中对建立的模型进行有限元网格的划分。
进一步的，所述步骤三中将建立的减振器模型赋予正交各向异性材料属性。
进一步的，所述步骤四中对减振器一端施加固定约束，另外一端施加单位载荷，进行有限元的求解，从而得到减振器的各向刚度。
进一步的，步骤五中所述步骤五中优化列式为：
Figure BDA0001610396630000021
其中，E1,E2,E3,v12,v23,v13,G12,G23,G13为减振器正交各向异性材料的材料参数；K'x,K'y,K'z为有限元求解所得的各向刚度；Kx,Ky,Kz为实际减振器的各向刚度。
进一步的，所述步骤六中将优化得到的正交各向异性材料的各向弹性模量赋予整体结构，进行整体结构的有限元分析，从而实现整体结构的有限元求解。
本发明与现有技术相比的优点在于：本发明提出了基于优化算法的减振器刚度模拟方法，采用正交各向异性材料模拟减振器模型，并通过优化得到与实际减振器各向刚度相同的有限元减振器模型，为有限元静力计算中减振器的模拟提供一种可选方法。所提基于优化算法的减振器刚度模拟方法简单，容易实施，运算速度快，节省硬件资源和成本。
附图说明
图1是本发明针对基于优化算法的减振器刚度模拟方法的总体流程图；
图2是本发明中优化实施流程图；
图3是本发明中优化目标的迭代历程图；
图4是本发明中减振器x向刚度的迭代历程图；
图5是本发明中减振器y向刚度的迭代历程图；
图6是本发明中减振器z向刚度的迭代历程图。
具体实施方式
下面结合附图以及具体实施例进一步说明本发明。
本发明提出了一种基于优化算法的减振器刚度模拟方法，具体步骤如下：
步骤一：根据实际减振器的几何形状建立相应的物理模型。
步骤二：对建立的模型进行有限元网格的划分。
步骤三：将建立的减振器模型赋予正交各向异性材料属性。
步骤四：对减振器一端施加固定约束，另外一端施加单位载荷，进行有限元的求解，从而得到减振器的各向刚度。
步骤五：以正交各向异性材料的各向弹性模量为设计变量，以有限元模拟得到的减振器各向刚度与实际减振器各向刚度最接近为优化目标，进行优化，从而精确模拟实际减振器的刚度。优化列式为：
Figure BDA0001610396630000031
其中，E1,E2,E3,v12,v23,v13,G12,G23,G13为减振器正交各向异性材料的材料参数；K'x,K'y,K'z为有限元求解所得的各向刚度；Kx,Ky,Kz为实际减振器的各向刚度。
步骤六：将优化得到的正交各向异性材料的各向弹性模量赋予整体结构，进行整体结构的有限元分析，从而实现整体结构的有限元求解。
实施例1：
为了更充分地了解该发明的特点及其对工程实际的适用性，本发明针对某飞机上的液冷源***减振器的刚度进行模拟。为保证液冷源的稳定性，液冷源上装有4个减振器。然而在进行有限元静力仿真时，必须要准确模拟出减振器的刚度。已知减振器的各向刚度为Kx＝22N/mm,Ky＝22N/mm,Kz＝88N/mm。采用优化方法对减振器刚度进行模拟，优化列示为：
Figure BDA0001610396630000032
其中，E1,E2,E3,v12,v23,v13,G12,G23,G13为减振器正交各向异性材料的材料参数；K'x,K'y,K'z为有限元求解所得的各向刚度。
采用有限元网格划分软件ANSA对减振器划分网格，导入有限元前处理软件Patran中，将其赋予正交各向异性材料属性。施加载荷和边界条件，一端固支，另一端施加单位载荷。采用有限元求解软件Nastran进行有限元求解从而得到减振器的各向刚度。最后，采用优化软件iSIGHT搭建如图2所示平台，进行优化，图中采用Patran软件进行有限元前处理和弹性参数修改，采用Nastran软件进行减振器位移计算，使用Patran软件提取位移；采用Matlab软件计算减振器刚度。
在建立的优化平台上进行优化，目标函数迭代曲线如图3所示。各向材料参数优化结果为：
E1＝1.37MPa,E2＝5.00MPa,E3＝14.14MPa,
G12＝1050.28MPa,G23＝1083.70MPa,G13＝75.73MPa,
v12＝0.388,v23＝0.353,v13＝0.175
x向，y向和z向刚度优化迭代曲线分别如图4，图5，图6所示，基本接近真实值。
综上所述，本发明提出了一种基于优化算法的减振器刚度模拟方法。首先，根据实际减振器的几何形状建立相应的物理模型；其次，对建立的模型进行有限元网格的划分；第三，将建立的减振器模型赋予正交各向异性材料属性；第四，对减振器一端施加固定约束，另外一端施加单位载荷，进行有限元的求解，从而得到减振器的各向刚度；最后，以正交各向异性材料的各向弹性模量为设计变量，以有限元模拟得到的减振器各向刚度与实际减振器各向刚度最接近为优化目标，进行优化，从而精确模拟实际减振器的刚度。
以上仅是本发明的具体步骤，对本发明的保护范围不构成任何限制；凡采用等同变换或者等效替换而形成的技术方案，均落在本发明权利保护范围之内。
本发明未详细阐述部分属于本领域技术人员的公知技术。

Claims (1)

1.一种基于优化算法的减振器刚度模拟方法，其特征在于，实现步骤如下：
步骤一：根据实际减振器的几何形状建立相应的物理模型；
步骤二：对建立的模型进行有限元网格的划分；
步骤三：将建立的减振器模型赋予正交各向异性材料属性；
步骤四：对减振器一端施加固定约束，另外一端施加单位载荷，进行有限元的求解，从而得到减振器的各向刚度；
步骤五：以正交各向异性材料的各向弹性模量为设计变量，以有限元模拟得到的减振器各向刚度与实际减振器各向刚度最接近为优化目标，进行优化，从而精确模拟实际减振器的刚度；
步骤六：将优化得到的正交各向异性材料的各向弹性模量赋予整体结构，进行整体结构的有限元分析，从而实现整体结构的有限元求解；
所述步骤五中优化列式为：
Figure FDA0003168982940000011
其中，E1,E2,E3,v12,v23,v13,G12,G23,G13为减振器正交各向异性材料的材料参数；K'x,K'y,K'z为有限元求解所得的各向刚度；Kx,Ky,Kz为实际减振器的各向刚度。
CN201810261714.6A 2018-03-28 2018-03-28 一种基于优化算法的减振器刚度模拟方法 Active CN108595769B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810261714.6A CN108595769B (zh) 2018-03-28 2018-03-28 一种基于优化算法的减振器刚度模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810261714.6A CN108595769B (zh) 2018-03-28 2018-03-28 一种基于优化算法的减振器刚度模拟方法

Publications (2)

Publication Number Publication Date
CN108595769A CN108595769A (zh) 2018-09-28
CN108595769B true CN108595769B (zh) 2021-10-01

Family

ID=63624741

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810261714.6A Active CN108595769B (zh) 2018-03-28 2018-03-28 一种基于优化算法的减振器刚度模拟方法

Country Status (1)

Country Link
CN (1) CN108595769B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109296693A (zh) * 2018-11-23 2019-02-01 中国人民解放军国防科技大学 一种液电式馈能减振器及其参数确定方法
CN110381433A (zh) * 2019-06-13 2019-10-25 苏州上声电子股份有限公司 一种行人警示扬声器的减震分析方法
CN111460601B (zh) * 2019-12-06 2020-11-17 西南石油大学 基于岩石物理建模的正交各向异性地层地应力预测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6077302A (en) * 1998-02-12 2000-06-20 Egs, Inc. System and method for analyzing and designing vibration isolators
EP2778992A1 (en) * 2013-03-15 2014-09-17 Sikorsky Aircraft Corporation Additive typology optimized manufacturing for multi-functional components
CN104850696A (zh) * 2015-05-15 2015-08-19 燕山大学 基于等效弹性模量的大型机械结构静刚度优化方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI519987B (zh) * 2014-11-14 2016-02-01 財團法人工業技術研究院 結構拓樸最佳化設計方法
CN105701287B (zh) * 2016-01-11 2018-03-20 东南大学 一种平台式惯导***的三向等刚度橡胶减振器设计方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6077302A (en) * 1998-02-12 2000-06-20 Egs, Inc. System and method for analyzing and designing vibration isolators
EP2778992A1 (en) * 2013-03-15 2014-09-17 Sikorsky Aircraft Corporation Additive typology optimized manufacturing for multi-functional components
CN104850696A (zh) * 2015-05-15 2015-08-19 燕山大学 基于等效弹性模量的大型机械结构静刚度优化方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《桥梁结构刚度参数静力方法优化识别及有限元模型修正》;何志军;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20120515(第5期);第15-18页 *
《橡胶减振器刚度有限元数值预测方法研究》;黄德东;《机械科学与技术》;20130531;第32卷(第5期);第661-663页 *

Also Published As

Publication number Publication date
CN108595769A (zh) 2018-09-28

Similar Documents

Publication Publication Date Title
CN109543243B (zh) 基于adams的舰载机拦阻着舰动力学仿真分析方法
CN108595769B (zh) 一种基于优化算法的减振器刚度模拟方法
JP6608637B2 (ja) モード動的解析におけるラグランジュ乗数を回復するシステムおよび方法
CN106202794A (zh) 基于多场有限元仿真的继电器类单机输出特性获取方法
CN104850697B (zh) 基于ansys和adams的大型天线动力学建模方法
CN105260581A (zh) 舰船机电控制设备虚拟振动和冲击试验方法
CN103838913B (zh) 曲线箱梁弯桥的有限单元法
CN113094946B (zh) 一种用于模拟材料开裂的相场模型局部化自适应算法
CN109446731A (zh) 一种基于abaqus的岩土工程数值模拟方法
CN109740182A (zh) 一种基于再生核粒子的无网格物理变形仿真方法
JP6736108B2 (ja) しわを抑制する空間薄膜構造治具形状を最適化する設計方法
CN107480357B (zh) 一种多杆系结构有限元模型边界模拟方法
CN111125963A (zh) 基于拉格朗日积分点有限元的数值仿真***及方法
CN109271655B (zh) 一种基于非对称有限元算法的材料尺度效应分析方法
CN109063401B (zh) 采用等效静态位移实现结构接触力优化的方法
CN114330043A (zh) 一种多自由度多功能高负载液压机械臂刚柔耦合分析方法
CN111079327A (zh) 一种基于各向异性过滤技术的异形曲面加筋拓扑优化方法
CN109388833B (zh) 一种基于疲劳寿命的弹性元件结构优化设计方法
CN104834795B (zh) 包带连接结构接触摩擦非线性特性模拟方法及***
CN110989497B (zh) 基于迭代优化的多轴多激励振动控制组合选择方法及***
CN108595893A (zh) 一种基于三层预处理子的三维力学模态仿真模拟方法
Ren et al. Structural reanalysis based on FRFs using sherman–morrison–woodbury formula
CN115906333A (zh) 一种桁架结构的几何非线性等效板动力学建模与响应分析方法
CN109726477A (zh) 一种间隙转动副有限元建模及分析方法
CN111177855B (zh) 一种全局气动弹性优化中的气动结构求解方法及***

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant