CN109299512B - 一种基于质量影响的快速灵敏度分析方法 - Google Patents
一种基于质量影响的快速灵敏度分析方法 Download PDFInfo
- Publication number
- CN109299512B CN109299512B CN201810981550.4A CN201810981550A CN109299512B CN 109299512 B CN109299512 B CN 109299512B CN 201810981550 A CN201810981550 A CN 201810981550A CN 109299512 B CN109299512 B CN 109299512B
- Authority
- CN
- China
- Prior art keywords
- quality
- sensitivity
- matrix
- admittance
- velocity
- 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
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)
- Complex Calculations (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提供了一种基于质量影响的快速灵敏度分析方法,构造速度导纳矩阵,并获得前m阶模态频率,从结构第一个节点开始添加质量摄动项,将速度导纳矩阵代入矩阵修正公式获得摄动后的速度导纳矩阵,辨识结构的频率,获得结构模态频率对质量的灵敏度,按照节点顺序改变质量摄动点位置,重复前述步骤获得对应得灵敏度,从而获得整个结构模态频率对质量的灵敏度,绘制灵敏度曲线。本发明当结构的质量发生变化时,利用矩阵变换公式无需进行有限元再次计算,只需要初始的速度导纳信息进行数值计算即可获得摄动后的速度导纳,简化计算效率,更加方便,实现了基于速度导纳对质量的灵敏度快速计算方法,具有实际工程意义。
Description
技术领域
本发明涉及一种灵敏度分析方法,具体涉及一种质量影响的灵敏度分析方法。
背景技术
灵敏度分析是一种评价因设计变量或参数的改变而引起结构响应特性变化率的方法。针对结构系统灵敏度的研究也越来重要,它是当前计算力学和结构工程领域的主要研究方向之一。
实际应用中,结构的灵敏度分析在结构优化、可靠性评估和参数识别等占据重要作用,方法按照策略可以分为离散法和变分法。传统的灵敏度分析方法一般需要在摄动下,重新进行有限元计算,特别当结构具有柔性特征时,有限元数量较多时,计算量大,效率低。有些工况需要测量速度,从而基于速度导纳信息进行如对质量参数等的灵敏度分析,如果解决传统灵敏度计算效率低,已成为亟待解决的实际工程问题。
发明内容
发明目的:本发明的目的在于针对现有技术的不足,提供一种基于质量影响的快速灵敏度分析方法。
技术方案:本发明提供了一种基于质量影响的快速灵敏度分析方法,包括以下步骤:
(1)基于Matlab软件进行有限元分析,获得速度导纳信息,构造速度导纳矩阵,并获得前m阶模态频率,从结构第一个节点开始添加质量摄动项;
(2)基于步骤(1)中速度导纳矩阵代入矩阵修正公式获得摄动后的速度导纳矩阵;
(3)辨识结构的频率,获得结构模态频率对质量的灵敏度;
(4)按照节点顺序改变质量摄动点位置,重复步骤(2)(3)获得对应的灵敏度,从而获得整个结构模态频率对质量的灵敏度,绘制灵敏度曲线。
进一步,步骤(1)包括以下步骤:
(11)结构的自由度为N,结构的速度导纳矩阵为:
其中,vij表示在结构节点j作用单位脉冲下节点i的速度导纳信息,i=1,2…N,j=1,2…N;
(12)采用有理正交多项式拟合算法对vij进行辨识,获得前m阶模态频率,分别为qr(r=1,2…m),将qr元素组成m维列向量q。
进一步,步骤(2)包括以下步骤:
(21)在结构i节点处添加质量摄动量Δmi(i=1,2…N),此时结构的速度导纳矩阵为V*,根据速度导纳矩阵的逆矩阵与动态刚度矩阵的关系,推导获得:
其中,ui∈RN×1,表示在列向量第i个元素为e为虚数单位,即e2=-1,ω为圆频率;
(22)由矩阵修正公式可知:
即可以建立摄动前后结构的速度导纳的联系;
(23)由公式(1)、(2)、(3)联列,化解可以获得质量摄动后速度导纳矩阵与初始导纳矩阵的关系:
进一步,步骤(3)包括以下步骤:
(31)基于摄动后的速度导纳矩阵V*,重复步骤(12)提取摄动后结构的频率(r=1, 2…m),将元素组成m维列向量q*;
(32)定义获取灵敏度函数gi,表示在节点i处发生质量摄动后的灵敏度,函数输入量为摄动前后的频率,具体公式如下:
其中,等式左边括号量即为输入函数gi的输入量,等式右边即为计算输出结果,具体形式为:
进一步确定灵敏度最大处的阶次信息,调用MATALB中find函数:
有益效果:本发明方法首先通过有限元计算获得结构的速度导纳信息,当结构的质量发生变化时,利用矩阵变换公式无需进行有限元再次计算,只需要初始的速度导纳信息进行数值计算即可获得摄动后的速度导纳,简化计算效率,更加方便,实现了基于速度导纳对质量的灵敏度快速计算方法,具有实际工程意义。
附图说明
图1为实施例中七自由度弹簧-阻尼-质量系统示意图;
图2为结构灵敏度曲线。
具体实施方式
下面对本发明技术方案进行详细说明,但是本发明的保护范围不局限于所述实施例。
本实施例采用一个简单的7自由度弹簧-阻尼-质量系统来验证,如图1所示,系统的参数分别为:mi=1kg(i=1,2…7),弹簧ki=10N/m(i=1,2…7),阻尼器ci=0.02N·m/s(i=1, 2…7),添加的质量摄动量为Δm=0.5kg。包括以下步骤:
步骤1,基于Matlab软件进行有限元分析,获得速度导纳信息,构造速度导纳矩阵,并获得前7阶模态频率,从结构第一个节点开始添加质量摄动项:
1.1)基于有限元计算获得结构的速度导纳矩阵为:
其中,vij表示在结构节点j作用单位脉冲下节点i的速度导纳信息,i=1,2…8,j=1,2…8;
1.2)采用有理正交多项式拟合算法对v11进行辨识,获得前7阶频率,分别为qr(r=1, 2…7),将qr元素组成7维列向量q,具体为:
q1=0.105,q2=0.311,q3=0.503,
q4=0.674,q5=0.814,q6=0.921,q7=0.988
q=[0.105 0.311 0.503 0.674 0.814 0.921 0.988]T
在结构节点1处添加质量摄动项Δm1=0.5。
步骤2,基于步骤(1)中速度导纳矩阵代入矩阵变换公式获得摄动后的速度导纳矩阵:
2.1)摄动后此时结构的速度导纳矩阵为V*,根据速度导纳矩阵的逆矩阵与动态刚度矩阵的关系,推导获得:
η=-1,Δmi>0
e为虚数单位,即e2=-1,ω为圆频率;
2.2)由矩阵修正公式可知:
即可以建立其摄动前后结构的速度导纳的联系;
2.3)由公式(1)、(2)、(3)联列化解可以获得质量摄动后速度导纳矩阵与初始导纳矩阵的关系:
步骤3,辨识结构的频率,获得结构模态频率对质量的灵敏度:
3.1)基于摄动后的速度导纳矩阵V*,重复步骤(12)提取摄动后结构的频率(r=1, 2…7),将元素组成7维列向量q*;
q*=[0.105 0.303 0.476 0.636 0.784 0.905 0.984]T
3.2)定义获取灵敏度函数g1,表示在节点i处发生质量摄动后的灵敏度,函数输入量为摄动前后的频率,具体公式如下:
其中,等式左边括号量即为输入函数g1的输入量,等式右边即为计算输出结果,具体形式为
进一步确定灵敏度最大处的阶次信息,调用MATALB中find函数:
步骤4,按照节点顺序改变质量摄动点位置,重复步骤(2)(3)获得对应的灵敏度,从而获得整个结构模态频率对质量的灵敏度,绘制灵敏度曲线,见图2。
图2结果表明,当结构的质量发生改变时,在节点6处灵敏度的绝对值最大,故在结构节点6处修改质量影响最大,而在节点7处影响最小。本发明的快速灵敏度分析方法,突破了传统灵敏度需要进行多次的计算的局限性,只需要进行一次有限元计算,利用初始结构的速度导纳进行数值计算获得质量摄动后的结构的响应,更加快捷。
Claims (1)
1.一种基于质量影响的快速灵敏度分析方法,其特征在于:包括以下步骤:
(1)基于Matlab软件进行有限元分析,获得速度导纳信息,构造速度导纳矩阵,并获得前m阶模态频率,从结构第一个节点开始添加质量摄动项;
(2)基于步骤(1)中速度导纳矩阵代入矩阵修正公式获得摄动后的速度导纳矩阵;
(3)辨识结构的频率,获得结构模态频率对质量的灵敏度;
(4)按照节点顺序改变质量摄动点位置,重复步骤(2)(3)获得对应的灵敏度,从而获得整个结构模态频率对质量的灵敏度,绘制灵敏度曲线;
具体的,步骤(1)包括以下步骤:
(11)结构的自由度为N,结构的速度导纳矩阵为:
其中,vij表示在结构节点j作用单位脉冲下节点i的速度导纳信息,i=1,2…N,j=1,2…N;
(12)采用有理正交多项式拟合算法对vij进行辨识,获得前m阶模态频率,分别为qr,r=1,2…m,将qr元素组成m维列向量q;
步骤(2)包括以下步骤:
(21)在结构i节点处添加质量摄动量Δmi,i=1,2…N,此时结构的速度导纳矩阵为V*,根据速度导纳矩阵的逆矩阵与动态刚度矩阵的关系,推导获得:
其中,ui∈RN×1,表示在列向量第i个元素为e为虚数单位,即e2=-1,ω为圆频率;
(22)由矩阵修正公式可知:
即可以建立摄动前后结构的速度导纳的联系;
(23)由公式(1)、(2)、(3)联列,化解可以获得质量摄动后速度导纳矩阵与初始导纳矩阵的关系:
步骤(3)包括以下步骤:
(31)基于摄动后的速度导纳矩阵V*,重复步骤(12)提取摄动后结构的频率r=1,2…m,将元素组成m维列向量q*;
(32)定义获取灵敏度函数gi,表示在节点i处发生质量摄动后的灵敏度,函数输入量为摄动前后的频率,具体公式如下:
其中,等式左边括号量即为输入函数gi的输入量,等式右边即为计算输出结果,具体形式为:
进一步确定灵敏度最大处的阶次信息,调用MATALB中find函数:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810981550.4A CN109299512B (zh) | 2018-08-27 | 2018-08-27 | 一种基于质量影响的快速灵敏度分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810981550.4A CN109299512B (zh) | 2018-08-27 | 2018-08-27 | 一种基于质量影响的快速灵敏度分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109299512A CN109299512A (zh) | 2019-02-01 |
CN109299512B true CN109299512B (zh) | 2019-06-18 |
Family
ID=65165440
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810981550.4A Active CN109299512B (zh) | 2018-08-27 | 2018-08-27 | 一种基于质量影响的快速灵敏度分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109299512B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107356387A (zh) * | 2017-07-21 | 2017-11-17 | 东南大学 | 一种模态试验中多传感器附加质量消除方法 |
CN108256256A (zh) * | 2018-01-31 | 2018-07-06 | 东南大学 | 一种保证阻尼系统特定模态频率不变的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10061878B2 (en) * | 2015-12-22 | 2018-08-28 | Dassault Systemes Simulia Corp. | Effectively solving structural dynamics problems with modal damping in physical coordinates |
-
2018
- 2018-08-27 CN CN201810981550.4A patent/CN109299512B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107356387A (zh) * | 2017-07-21 | 2017-11-17 | 东南大学 | 一种模态试验中多传感器附加质量消除方法 |
CN108256256A (zh) * | 2018-01-31 | 2018-07-06 | 东南大学 | 一种保证阻尼系统特定模态频率不变的方法 |
Non-Patent Citations (3)
Title |
---|
On the use of fourier descriptors for the assessment of frequency coupling matrices of power electronic devices;Roberto Langella 等;《2018 18th International Conference on Harmonics and Quality of Power (ICHQP)》;20180611;第1-6页 |
基于灵敏度分析的网络固有结构理论的应用研究;林志清 等;《四川电力技术》;20051230;第15-18页 |
复合地基动力特性及其应用研究;宋春雨;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20100915(第09期);第C038-35页 |
Also Published As
Publication number | Publication date |
---|---|
CN109299512A (zh) | 2019-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102081359B (zh) | 基于DSP Builder的变时滞超混沌数字电路设计方法及电路 | |
CN108984976B (zh) | 一种基于加速度响应结构灵敏度计算方法 | |
TW200630819A (en) | Method of using intelligent theory to design heat dissipation module and device thereof | |
CN108108559B (zh) | 一种基于子结构的结构响应获取方法及灵敏度获取方法 | |
CN109344433B (zh) | 基于响应信号的灵敏度数值计算方法 | |
CN105978725A (zh) | 一种基于传感器网络的非脆弱性分布式故障估计方法 | |
CN110889509A (zh) | 一种基于梯度动量加速的联合学习方法及装置 | |
CN109145446B (zh) | 一种基于模态应变能和卷积神经网络的结构损伤识别方法 | |
CN104075402A (zh) | 智能空调控制方法和系统 | |
CN113722966B (zh) | 一种集成电路板仿真多级分布式并行计算方法 | |
CN103729564A (zh) | 一种基于粒子图像测速技术的压力场计算方法和装置 | |
CN109147081A (zh) | 一种设备运行稳定性分析方法及系统 | |
CN109299512B (zh) | 一种基于质量影响的快速灵敏度分析方法 | |
CN105893448A (zh) | 一种生成应用热度榜单的方法及系统 | |
CN108493936A (zh) | 基于子空间辨识方法的电力系统低频振荡估计的改进方法 | |
CN113221475A (zh) | 一种用于高精度流场分析的网格自适应方法 | |
Zhu et al. | Analysis of a diffusive SIS epidemic model with spontaneous infection and a linear source in spatially heterogeneous environment. | |
CN110171426A (zh) | 一种基于智能设备的驾驶员危险驾驶行为识别方法及系统 | |
CN109299513B (zh) | 一种模态频率对质量的灵敏度分析方法 | |
JP4612989B2 (ja) | メカトロニクスシステムのシミュレーション方法 | |
CN108984480A (zh) | 一种多个加速度传感器附加质量影响的一步消除方法 | |
CN104124908B (zh) | 在线识别惯量比的系统及方法 | |
CN206133535U (zh) | 基于气体传感器的真随机数发生器及真随机数生成装置 | |
CN102682174A (zh) | 一种基于传感数据的柔性有限元自动分析方法 | |
CN105512359B (zh) | 一种处理cae仿真结果的方法及系统 |
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 |