CN111475980A - 一种综合作动器质量影响的薄壁件动力学参数获取方法 - Google Patents
一种综合作动器质量影响的薄壁件动力学参数获取方法 Download PDFInfo
- Publication number
- CN111475980A CN111475980A CN202010271907.7A CN202010271907A CN111475980A CN 111475980 A CN111475980 A CN 111475980A CN 202010271907 A CN202010271907 A CN 202010271907A CN 111475980 A CN111475980 A CN 111475980A
- Authority
- CN
- China
- Prior art keywords
- mass
- thin
- wall part
- vibration control
- actuator
- 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
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
Abstract
本发明一种综合作动器质量影响的薄壁件动力学参数获取方法,属于结构振动控制领域;该方法首先对薄壁件进行有限元分析获得其总质量矩阵MT、总刚度矩阵KT、总阻尼矩阵CT。然后通过有限元模态分析获得各个极大变形的位置,按本发明提出的方法对极大变形点的质量、刚度和阻尼参数进行确定,累加作动器的随动件的质量和薄壁件极大变形点的质量获得新的质量参数,为结构振动控制提供控制参数。解决了薄壁件振动控制时作动器作用位置、振动控制动力学质量参数模糊的问题。使得振动控制对象更为合理,对于具体应用具有很强工程应用意义。
Description
技术领域
本发明属于结构振动控制领域,具体涉及一种综合作动器质量影响的薄壁件动力学参数获取方法,适用于作为壳体的薄壁件结构振动控制时,作动器的作用位置不确定、薄壁件在作用点的动力学参数不确定以及作动器随动件的质量对结构振动影响较大的情形。
背景技术
结构振动的主动控制是通过改变结构质量、阻尼与刚度等动力特性来控制结构响应。作动器作为主动控制的输出器件被广泛应用于振动控制领域取得了良好的振动控制效果。现有的振动控制技术仅将作动器与被控制件独立建模进行控制律的设计。然而在结构振动控制中存在两方面问题:
1.一般作动器的滑动件(例如:液压阻尼器、磁流变阻尼器的活塞杆及与活塞杆固接的部件)与被减振件为固定连接,现有的技术进行振动控制时将作动器视为力发生器件进行设计,忽略作动器对被减振结构的质量影响。在大型结构振动控制中,作动器随动件质量相对被控结构质量较小影响不大,但在薄壁件振动控制领域则较为严重。
2.在薄壁件振动控制领域,由于薄壁件往往为具有弱刚度特性的连续体,难以采用集中质量的建模方法进行建模。
因此,针对现有薄壁件结构振动控制时其动力学参数M、C、K不确定,作动器对薄壁件结构振动控制的影响被忽略的问题,提出了一种综合作动器质量影响的薄壁件振动控制动力学参数获取方法。
发明内容
要解决的技术问题:
为了避免现有技术的不足之处,本发明提出一种综合作动器质量影响的薄壁件动力学参数获取方法,既一种静力学变形对极大变形点周边区域的质量进行加权积分的方法,对作为壳体的薄壁件极大变形点的振动控制质量参数进行估计的方法。本发明通过将作动器随动件的质量累加到被减振薄壁件的振动控制质量参数来综合考虑作动器对薄壁件振动控制的影响,所述薄壁件又称被减振件或工件。
本发明的技术方案是:一种综合作动器质量影响的薄壁件动力学参数获取方法,其特征在于具体步骤如下:
步骤一:采用有限元方法对被减振件分析,获得其总质量矩阵MT、总刚度矩阵KT、总阻尼矩阵CT,维度为n×n;
步骤二:通过有限元模态分析,获得极大变形点的位置坐标,设定其中一个极大变形点坐标为(x0,y0);
步骤三:对被减振件进行静力学分析,所施加的静力载荷为被减振件实际受到激励的幅值,所施加位置为步骤二确定的极大变形点,从而获得该点的最大变形量δ0;
步骤五:对区域S所对应的质量采用加权积分方法进行质量估计,公式为:
其中ρ(x,y)面密度函数,MW为被减振件区域S的估计质量;
步骤六:确定对被减振件振动控制所用作动器,获得其随动件质量MA;
步骤七:将步骤五所得估计质量MW与步骤六所得随动件质量MA相加,为被减振件振动控制动力学质量参数MC,即MC=MA+MW;
步骤八:通过选择向量组获得总刚度矩阵在极大变形点的刚度值,作为被减振件振动控制动力学刚度参数KC,公式为:
步骤九:通过选择向量组获得总阻尼矩阵在极大变形点的阻尼值,作为被减振件振动控制动力学阻尼参数CC,公式为:
步骤十:基于上述得出的参数MC、KC、CC能够进行被减振件结构振动控制算法的设计。
本发明的进一步技术方案是:所述步骤六的作动器为具有活塞杆的阻尼器,所述随动件为随被减振件运动的部件,即阻尼器的活塞杆及与活塞杆固定的部件。
有益效果
本发明的有益效果在于:
2.通过对作动器的结构分析,将作动器中与薄壁件固接的随动件的质量累加到薄壁件结构振动控制质量参数上,使得薄壁件振动控制的动力学质量参数更为完善,从而提高了薄壁件在加工过程中的振动控制精度;通过所述加权积分式获得了薄壁件结构的动力学估计质量,但是当作动器的随动件质量相对较大时动力学质量参数任然不精确,本发明解决了薄壁件振动控制时结构振动控制的动力学质量参数任然不精确的问题。
附图说明
图1是本发明的总体流程图;
图2是本发明实施例所采用的薄壁件变形前示意图;
图3是本发明实施例对薄壁件模态分析获得法向振幅图;
图4是本发明实施例在一极大点施加一静力后变形分布图;
图5是本发明实施例相对变形分布图;
图6是本发明的实施例所采用作动器的剖视图。
具体实施方式
下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
以下通过参考附图详细描述本发明的较优秀实施例。在以下描述中,名词“上”,“下”对应于图中的上下方向。通过对一阶模态的分析进一步描述本发明具体技术路线。本实施例薄壁件为均匀材质,即面密度ρ(x,y)为常量。本实施例基于有限元方法,采用离散积分的得出薄壁件极大变形点的估计质量值。
首先,有限元方法对被减振件(也称为薄壁件、工件)分析获得总质量矩阵MT、总刚度矩阵KT、总阻尼矩阵CT,维度为n×n;采用模态分析获得各阶振型的多个振幅极大点,其中一个极大点的坐标设为(x0,y0),对振型的极大点(坐标为(x0,y0))进行静力学分析,施加的外力为薄壁件受到外界激励的幅值,确定在该点的变形δ0,按获得需要考虑变形质量的区域S,用于调节变形区域S的精度。δ(x,y)是标为(x,y)处的变形。之后通过公式
然后,对作动器进行结构分析确定其随动件的质量参数MA,则薄壁件结构振控制在该振型极大点的质量参数MC为:MC=MA+MW。
最后,由于结构总刚度矩阵KT和阻尼矩阵CT反映各节点的真实刚度阻尼值,对于一般连续结构在一确定节点处的刚度、阻尼与其周边节点的刚度、阻尼基本一致。因此只需通过选择向量对(n×n)维的总刚度MT、总阻尼CT矩阵中表征节点位置为(x0,y0)的刚度、阻尼值进行选取即可。选择向量组包含行向量和列向量,其中行向量: 其中列向量:其一般表达式为:
按图1所示流程图实施,本实施例对图2所示薄壁件进行具体方案的描述。
首先,对薄壁件进行模态分析,得到二阶下的振型分布如图3,得到两个振幅极大点为幅值方向相反,选其中一点进行静力学分析,其载荷为薄壁件使用场合可能受到激励的最大值。得到变形分布图,将其均匀分割为等面积的正方形,设每一小格具有质量为Δm,按图1流程图所示质量计算积分规则得离散求和公式:求和上下界按图1所示设定阈值ζ,通过不等式条件确定。
对作动器进行结构分析获得其随动件的质量参数,由图1流程图所述技术,可知图6作动器的活塞杆为随动件,获得作动器的质量为MA。
薄壁件振动控制动力学质量参数质量MC为作动器质量MA与薄壁件极大变形点估计质量MW之和。即MC=MW+MA。
通过有限元软件获得薄壁件的总刚度、总阻尼矩阵为(n×n)维,当二维矩阵元素坐标值与薄壁件节点的位置一致时可通过选择向量组获得节点处的刚度阻尼值分别为KC、CC。其一般表达式为:
本实施例中n为所划分等面积正方形在行、列方向个数的最大值。(x0,y0)为极大变形点的坐标。
通过上述技术方案的实施获得薄壁件振动控制动力学参数值分别为:MC、KC、CC。按该参数可进行后续薄壁件振动主动控制工作。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (2)
1.一种综合作动器质量影响的薄壁件动力学参数获取方法,其特征在于具体步骤如下:
步骤一:采用有限元方法对被减振件分析,获得其总质量矩阵MT、总刚度矩阵KT、总阻尼矩阵CT,维度为n×n;
步骤二:通过有限元模态分析,获得极大变形点的位置坐标,设定其中一个极大变形点坐标为(x0,y0);
步骤三:对被减振件进行静力学分析,所施加的静力载荷为被减振件实际受到激励的幅值,所施加位置为步骤二确定的极大变形点,从而获得该点的最大变形量δ0;
步骤五:对区域S所对应的质量采用加权积分方法进行质量估计,公式为:
其中ρ(x,y)面密度函数,MW为被减振件区域S的估计质量;
步骤六:确定对被减振件振动控制所用作动器,获得其随动件质量MA;
步骤七:将步骤五所得估计质量MW与步骤六所得随动件质量MA相加,为被减振件振动控制动力学质量参数MC,即MC=MA+MW;
步骤八:通过选择向量组获得总刚度矩阵在极大变形点的刚度值,作为被减振件振动控制动力学刚度参数KC,公式为:
步骤九:通过选择向量组获得总阻尼矩阵在极大变形点的阻尼值,作为被减振件振动控制动力学阻尼参数CC,公式为:
步骤十:基于上述得出的参数MC、KC、CC能够进行被减振件结构振动控制算法的设计。
2.根据权利要求1所述综合作动器质量影响的薄壁件动力学参数获取方法,其特征在于:所述步骤六的作动器为具有活塞杆的阻尼器,所述随动件为随被减振件运动的部件,即阻尼器的活塞杆及与活塞杆固定的部件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010271907.7A CN111475980B (zh) | 2020-04-09 | 2020-04-09 | 一种综合作动器质量影响的薄壁件动力学参数获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010271907.7A CN111475980B (zh) | 2020-04-09 | 2020-04-09 | 一种综合作动器质量影响的薄壁件动力学参数获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111475980A true CN111475980A (zh) | 2020-07-31 |
CN111475980B CN111475980B (zh) | 2022-07-26 |
Family
ID=71751205
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010271907.7A Active CN111475980B (zh) | 2020-04-09 | 2020-04-09 | 一种综合作动器质量影响的薄壁件动力学参数获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111475980B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112861353A (zh) * | 2021-02-08 | 2021-05-28 | 清华大学 | 薄壁结构阻尼涂层涂敷位置的确定方法 |
CN114278695A (zh) * | 2021-12-24 | 2022-04-05 | 西北工业大学 | 一种基于磁流变阻尼器的薄壁件加工振动半主动控制方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009053946A2 (en) * | 2007-10-26 | 2009-04-30 | Ecole Polytechnique Federale De Lausanne (Epfl) | Structures with adaptive stiffness and damping integrating shear thickening fluids |
US9073623B1 (en) * | 2013-03-15 | 2015-07-07 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | System and method for dynamic aeroelastic control |
CN107030522A (zh) * | 2017-05-24 | 2017-08-11 | 西北工业大学 | 面向复杂薄壁零件加工振动抑制的磁流变阻尼支撑器件 |
CN107423489A (zh) * | 2017-06-27 | 2017-12-01 | 西北工业大学 | 薄壁件铣削过程稳定性快速预测方法 |
CN107563005A (zh) * | 2017-08-04 | 2018-01-09 | 大连理工大学 | 一种拉压不同刚度结构振动瞬时最优控制方法 |
CN108897213A (zh) * | 2018-05-28 | 2018-11-27 | 宝鸡文理学院 | 一种基于面内调控的褶皱薄膜主动振动控制方法 |
-
2020
- 2020-04-09 CN CN202010271907.7A patent/CN111475980B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009053946A2 (en) * | 2007-10-26 | 2009-04-30 | Ecole Polytechnique Federale De Lausanne (Epfl) | Structures with adaptive stiffness and damping integrating shear thickening fluids |
US9073623B1 (en) * | 2013-03-15 | 2015-07-07 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | System and method for dynamic aeroelastic control |
CN107030522A (zh) * | 2017-05-24 | 2017-08-11 | 西北工业大学 | 面向复杂薄壁零件加工振动抑制的磁流变阻尼支撑器件 |
CN107423489A (zh) * | 2017-06-27 | 2017-12-01 | 西北工业大学 | 薄壁件铣削过程稳定性快速预测方法 |
CN107563005A (zh) * | 2017-08-04 | 2018-01-09 | 大连理工大学 | 一种拉压不同刚度结构振动瞬时最优控制方法 |
CN108897213A (zh) * | 2018-05-28 | 2018-11-27 | 宝鸡文理学院 | 一种基于面内调控的褶皱薄膜主动振动控制方法 |
Non-Patent Citations (4)
Title |
---|
EDUARDO DIEZ 等: "Feasibility study of in-process compensation of deformations in flexible milling", 《INTERNATIONAL JOURNAL OF MACHINE TOOLS & MANUFACTURE》 * |
ZHAO JU 等: "The Machining Error Compensation Method for Aeronautical Thin-walled Parts Machining Process", 《INTERNATIONAL CONFERENCE ON INTELLIGENT SYSTEMS RESEARCH AND MECHATRONICS ENGINEERING (ISRME 2015)》 * |
王晨晨 等: "滚动振动试验台垂向激振作动器的选型", 《机车电传动》 * |
陈冰 等: "面向航空发动机薄壁零件加工的自适应夹具设计现状与进展", 《航空制造技术》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112861353A (zh) * | 2021-02-08 | 2021-05-28 | 清华大学 | 薄壁结构阻尼涂层涂敷位置的确定方法 |
CN112861353B (zh) * | 2021-02-08 | 2022-05-03 | 清华大学 | 薄壁结构阻尼涂层涂敷位置的确定方法 |
CN114278695A (zh) * | 2021-12-24 | 2022-04-05 | 西北工业大学 | 一种基于磁流变阻尼器的薄壁件加工振动半主动控制方法 |
CN114278695B (zh) * | 2021-12-24 | 2023-06-27 | 西北工业大学 | 一种基于磁流变阻尼器的薄壁件加工振动半主动控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111475980B (zh) | 2022-07-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111475980B (zh) | 一种综合作动器质量影响的薄壁件动力学参数获取方法 | |
Lee et al. | Statistical optimization and assessment of a thermal error model for CNC machine tools | |
CN111475976B (zh) | 考虑混合不确定性的颗粒增强材料构件稳健拓扑优化方法 | |
CN113434921A (zh) | 一种考虑介纳观尺度效应的结构等几何拓扑优化方法 | |
Kitayama et al. | Optimization of variable blank holder force trajectory by sequential approximate optimization with RBF network | |
WO2020219544A1 (en) | Topology optimization of structure with multiple targets | |
CN105446264B (zh) | 基于特征的机床精度优化设计方法 | |
Lin et al. | Application of the fuzzy-based Taguchi method for the structural design of drawing dies | |
CN108763810A (zh) | 一种桥梁静载试验荷载加载布置与调整方法 | |
CN113221478A (zh) | 一种主动控制旋翼的数值分析方法及系统 | |
CN108897935A (zh) | 面向增材制造的点阵结构动力学响应优化设计方法 | |
KR100749850B1 (ko) | 회전 소성 가공의 수치 시뮬레이션 방법 및 기록 매체 | |
CN116560301A (zh) | 一种基于梯度优化的机床进给系统数理模型参数辨识方法 | |
Kayacan et al. | Velocity control of a spherical rolling robot using a grey-PID type fuzzy controller with an adaptive step size | |
CN109255141A (zh) | 一种汽车车身正向概念设计截面形状优化方法 | |
Burkart et al. | Optimization of a part holder design considering dynamic loads during return stroke of tool and ram | |
Burkart et al. | Compensation of elastic die and press deformations during sheet metal forming by optimizing blank holder design | |
CN109541939B (zh) | 高速电梯机械振动的多尺度近似显式模型预测控制方法 | |
CN113642129B (zh) | 一种工件校正载荷快速施加及变形有限元获取方法 | |
Gao et al. | A 6-sigma robust optimization method for stamping forming of automobile covering parts based on residual error and radial basis interpolation | |
JP2019128621A (ja) | 定常変形解析方法、定常変形解析システムおよびコンピュータ読み取り可能な記録媒体 | |
CN112131665A (zh) | 一种基于网格分类的火箭流固载荷混合转换方法 | |
CN108563856B (zh) | 一种基于自由节点b样条建模的自适应采样方法 | |
Winter et al. | NURBS-based shape and parameter optimization of structural components with an adaptive amount of control points | |
CN107463760B (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 |