机械阻抗法与频响分析
目录
- 1. 什么是频响分析
- 2. 机械阻抗和导纳的概念
- 3. 集中参数元件的阻抗和导纳
- 3.1 阻尼元件的阻抗和导纳
- 3.2 机械阻抗网络图的建立
- 3.3 机械阻抗的串并联计算
- 4. 单自由度紫铜的频响特性和导纳曲线
- 4.1 无阻尼系统
- 4.1.1 自由系统
- 4.2 约束系统
- 4.2 有阻尼系统
- 4.2.1 自由系统
- 4.2.2 约束系统
- 5.结构阻尼
1. 什么是频响分析
用一特定已知的激励力,以可控的方式来激励结构,同时测量输入和输出信号,据此分析并找出结构的动态参数,以评价其动态特性或建立数学模型,这就是频响分析。
频响分析的结果一般表示为频率——响应曲线:包括幅频曲线,相频曲线,实频曲线,虚频曲线,或是响应向量的矢端轨迹图。
对于单自由度系统,频响分析的一般步骤是:固有频率→阻尼→刚度→质量。
2. 机械阻抗和导纳的概念
位移阻抗,速度阻抗和加速度阻抗的计算公式
Z D = F X = k − ω 2 m + j ω c Z V = F V = c + j ω m + k j ω Z A = F A = m − k ω 2 + c j ω \begin{aligned} Z_{D} & =\frac{F}{X}=k-\omega^{2} m+\mathrm{j} \omega c \\ Z_{V} & =\frac{F}{V}=c+\mathrm{j} \omega m+\frac{k}{\mathrm{j} \omega} \\ Z_{A} & =\frac{F}{A}=m-\frac{k}{\omega^{2}}+\frac{c}{\mathrm{j} \omega} \end{aligned} ZDZVZA=XF=k−ω2m+jωc=VF=c+jωm+jωk=AF=m−ω2k+jωc
位移导纳,速度导纳和加速度导纳的计算公式
Y D = 1 / Z D = X / F Y V = 1 / Z V = V / F Y A = 1 / Z A = A / F \begin{array}{l} Y_{D}=1 / Z_{D}=X / F \\ Y_{V}=1 / Z_{V}=V / F \\ Y_{A}=1 / Z_{A}=A / F \end{array} YD=1/ZD=X/FYV=1/ZV=V/FYA=1/ZA=A/F
从阻抗以及导纳的表达式来看,阻抗的物理意义是产生单位响应需要施加的力,导纳的物理意义是单位力能够产生的响应。
3. 集中参数元件的阻抗和导纳
3.1 阻尼元件的阻抗和导纳
3.2 机械阻抗网络图的建立
根据节点法,确定各元件的连接点与串并联关系。之后把力施加到节点上。另外需要注意的是,代表实际质量的元件总有一个接地端。
3.3 机械阻抗的串并联计算
- 机械阻抗的并联计算
若有n个阻抗并联,则其等效阻抗 Z p Z_p Zp为
Z p = ∑ i = 1 n Z i Z_{\mathrm{p}}=\sum_{i=1}^{n} Z_{i} Zp=i=1∑nZi - 机械阻抗的串联情况
串联系统的阻抗的倒数是参与串联的各阻抗倒数之和
1 Z s = ∑ i = 1 n 1 Z i \frac{1}{Z_{\mathrm{s}}}=\sum_{i=1}^{n} \frac{1}{Z_{i}} Zs1=i=1∑nZi1
4. 单自由度紫铜的频响特性和导纳曲线
4.1 无阻尼系统
4.1.1 自由系统
单自由度系统的力学模型与机械阻抗网络为:
系统的导纳为:
Y = − 1 ω 2 m + 1 k = ω 2 m − k ω 2 m k = − 1 ω 2 m ( 1 − ω 2 ω A 2 ) = 1 k ( 1 − ω A 2 ω 2 ) Y=-\frac{1}{\omega^{2} m}+\frac{1}{k}=\frac{\omega^{2} m-k}{\omega^{2} m k}=\frac{-1}{\omega^{2} m}\left(1-\frac{\omega^{2}}{\omega_{A}^{2}}\right)=\frac{1}{k}\left(1-\frac{\omega_{A}^{2}}{\omega^{2}}\right) Y=−ω2m1+k1=ω2mkω2m−k=ω2m−1(1−ωA2ω2)=k1(1−ω2ωA2)
其中, ω A = k m \omega_A=\sqrt{\frac{k}{m}} ωA=mk 。
可知,系统的导纳幅频特性曲线为:
4.2 约束系统
带约束单自由系统的力学模型与机械阻抗网络为:
该系统的位移阻抗为:
Z = − ω 2 m + k = m ( k m − ω 2 ) = m ( Ω n 2 − ω 2 ) = k ( 1 − ω 2 Ω n 2 ) Z=-\omega^{2} m+k=m\left(\frac{k}{m}-\omega^{2}\right)=m\left(\Omega_{n}^{2}-\omega^{2}\right)=k\left(1-\frac{\omega^{2}}{\Omega_{n}^{2}}\right) Z=−ω2m+k=m(mk−ω2)=m(Ωn2−ω2)=k(1−Ωn2ω2)
其中, Ω n 2 = k m \Omega_{n}^{2}=\frac{k}{m} Ωn2=mk,该系统阻抗的幅频特性曲线与导纳的幅频特性曲线分别为:
4.2 有阻尼系统
4.2.1 自由系统
单自由度黏性自由系统的力学模型与机械阻抗示意图为:
该系统的总导纳为:
Y p − G = 1 Z p − G = 1 Z p − a + 1 Z m = 1 k + j ω c + 1 − ω 2 m \begin{aligned} Y_{p-G} & =\frac{1}{Z_{p-G}}=\frac{1}{Z_{p-a}}+\frac{1}{Z_{m}} \\ & =\frac{1}{k+\mathrm{j} \omega c}+\frac{1}{-\omega^{2} m} \end{aligned} Yp−G=Zp−G1=Zp−a1+Zm1=k+jωc1+−ω2m1
利用导纳函数计算出 p p p点和 a a a点的位移,定义位移之比为振动传递系数
∣ T ∣ = 1 + 4 ζ 2 ( ω ω A ) 2 [ 1 − ( ω ω A ) 2 ] 2 + 4 ζ 2 ( ω ω A ) 2 |T|=\sqrt{\frac{1+4 \zeta^{2}\left(\frac{\omega}{\omega_{A}}\right)^{2}}{\left[1-\left(\frac{\omega}{\omega_{A}}\right)^{2}\right]^{2}+4 \zeta^{2}\left(\frac{\omega}{\omega_{A}}\right)^{2}}} ∣T∣=[1−(ωAω)2]2+4ζ2(ωAω)21+4ζ2(ωAω)2
利用该函数便可以对系统进行隔振设计。从传递函数图像可以看出,只有当 ω / ω A > 2 时 \omega/\omega_A>\sqrt{2}时 ω/ωA>2 时,才能起到隔振的目的,所以隔振弹簧的刚度小些比较好。另外,大的阻尼对隔振对起到不好的效果,但是实际上还是要有一点阻尼的,可以让振动快速衰减下来。
在进行上述形式的隔振设计时,通常的步骤为:
- step1:选定 ∣ T ∣ |T| ∣T∣
- step2:根据 ∣ T ∣ |T| ∣T∣值确定 ω / ω A \omega/\omega_A ω/ωA,这一步可通过查图确定
- step3:根据外界振动频率 ω \omega ω确定 ω A \omega_A ωA;
- step4:根据设备质量 m m m定 k k k
- step5:校核 k k k元件的刚度
4.2.2 约束系统
对于一单自由度黏性约束系统,其物理模型及机械阻抗网络图为:
该系统的总阻抗为:
Z = Z m + Z k + Z c = − m ω 2 + k + j ω c = ( k − m ω 2 ) 2 + ( c ω ) 2 e j Φ ′ = k [ 1 − ( ω Ω n ) 2 ] 2 + ( 2 ζ ω Ω n ) 2 e j Φ ′ Φ ′ = arctan c ω k − m ω 2 = arctan 2 ζ ω Ω n 1 − ( ω Ω n ) 2 \begin{aligned} Z & =Z_{m}+Z_{k}+Z_{c}=-m \omega^{2}+k+\mathrm{j} \omega c \\ & =\sqrt{\left(k-m \omega^{2}\right)^{2}+(c \omega)^{2}} \mathrm{e}^{\mathrm{j} \Phi^{\prime}}=k \sqrt{\left[1-\left(\frac{\omega}{\Omega_{n}}\right)^{2}\right]^{2}+\left(2 \zeta \frac{\omega}{\Omega_{n}}\right)^{2}} \mathrm{e}^{\mathrm{j} \Phi^{\prime}} \\ \Phi^{\prime} & =\arctan \frac{c \omega}{k-m \omega^{2}}=\arctan \frac{2 \zeta \frac{\omega}{\Omega_{n}}}{1-\left(\frac{\omega}{\Omega_{n}}\right)^{2}} \end{aligned} ZΦ′=Zm+Zk+Zc=−mω2+k+jωc=(k−mω2)2+(cω)2 ejΦ′=k[1−(Ωnω)2]2+(2ζΩnω)2 ejΦ′=arctank−mω2cω=arctan1−(Ωnω)22ζΩnω
换成导纳表达式,为:
Y = 1 / k [ 1 − ( ω Ω n ) 2 ] 2 + ( 2 ζ ω Ω n ) 2 e j Φ = ∣ Y ∣ e j Φ Y=\frac{1 / k}{\sqrt{\left[1-\left(\frac{\omega}{\Omega_{n}}\right)^{2}\right]^{2}+\left(2 \zeta \frac{\omega}{\Omega_{n}}\right)^{2}}} \mathrm{e}^{\mathrm{j} \Phi}=|Y| \mathrm{e}^{\mathrm{j} \Phi} Y=[1−(Ωnω)2]2+(2ζΩnω)2 1/kejΦ=∣Y∣ejΦ
其中:
∣ Y ∣ = 1 / k [ 1 − ( ω Ω n ) 2 ] 2 + ( 2 ζ ω Ω n ) 2 |Y|=\frac{1 / k}{\sqrt{\left[1-\left(\frac{\omega}{\Omega_{n}}\right)^{2}\right]^{2}+\left(2 \zeta \frac{\omega}{\Omega_{n}}\right)^{2}}} ∣Y∣=[1−(Ωnω)2]2+(2ζΩnω)2 1/k
Φ = − Φ ′ = arctan − 2 ζ ω Ω n 1 − ( ω Ω n ) 2 \Phi=-\Phi^{\prime}=\arctan \frac{-2 \zeta \frac{\omega}{\Omega_{n}}}{1-\left(\frac{\omega}{\Omega_{n}}\right)^{2}} Φ=−Φ′=arctan1−(Ωnω)2−2ζΩnω
其阻抗和导纳的幅频与相频曲线分别为:
5.结构阻尼
在之前的推导中,我们一直假定系统的阻尼力为 F D = c x ˙ F_{D}=c \dot{x} FD=cx˙,但是在实际情况下,阻尼力和试验频率几何是无关的。结构阻尼的假定提供了一种与频率无关的数学模型,它假定阻尼力的大小与位移大小成正比,但是与速度同相。
F D = h ∣ x ∣ x ˙ ∣ x ˙ ∣ F_{D}=h|x| \frac{\dot{x}}{|\dot{x}|} FD=h∣x∣∣x˙∣x˙
新的阻尼模型下,系统的动力学方程为:
[ − ω 2 m + k ( 1 + j h k ) ] X = F \left[-\omega^{2} m+k\left(1+\mathrm{j} \frac{h}{k}\right)\right] X=F [−ω2m+k(1+jkh)]X=F
令 h = g k h=gk h=gk, g g g称为结构阻尼因子,并令 k ∗ = k ( 1 + j g ) k^{*}=k(1+\mathrm{j} g) k∗=k(1+jg),称 k ∗ k^{*} k∗为复刚度。但在实际操作中,我们通常将阻尼修正为与结构阻尼等效的数值,等效的方法是基于能量消耗相等的原则。修正后的阻尼为:
C e = k g ω C_{\mathrm{e}}=\frac{k g}{\omega} Ce=ωkg
此时,系统的动力学方程修正为:
m x ¨ + C e x ˙ + k x = f ( t ) m \ddot{x}+C_{\mathrm{e}} \dot{x}+k x=f(t) mx¨+Cex˙+kx=f(t)
机械阻抗法与频响分析相关推荐
- 汽车振动响应分析-频响函数法(附程序)
文章目录 汽车的四自由度振动模型 频响函数法 MATLAB程序 汽车的四自由度振动模型 m:车身的等效质量: m_1:前轴的等效质量(包含电机的质量在内): m_2:后轴的等效质量: m_3:人体和座 ...
- matlab编程刀尖频响,用半理论法预测主轴系统刀尖点频响函数
1问题的提出主轴系统刀尖点频响函数是衡量和评价数控机床切削稳定性的重要指标,是优化机床加工工艺.评价机床动态特性以及综合性能的依据.其中,机床的自激振动(又称再生颤振)对切削稳定的影响极大.目前,主要 ...
- LOTO示波器如何测试阻抗的频响曲线
LOTO示波器如何测试阻抗的频响曲线 模块的输入输出端口,在电路分析上,一般简单表征为电阻来进行计算和分析.但多数情况下,这些端口并不是纯电阻的特性,更精确一些,它可能是电阻电容以及电感的组合,表现为 ...
- 用步进正弦激励法测量计算机声卡输出输入幅值频响函数
一. 实验目的 (1) 通过对计算机声卡系统的实测加深对频响函数的理解: (2) 对声音振动及测量有直观了解: (3) 对计算机声卡系统的频响特性有直观的认识. 二. 实验内容 (1) 用声卡发出的步 ...
- 运放基本结构及频响(一)
运算放大器基本结构及频响 绪论 基本结构 运放设计 频响 增益 开环 闭环 回到开环 绪论 运放设计{结构设计增益,级数,五管差分基本够用了参数设计MOS的(W/L),补偿电容C运放设计\begin{ ...
- 伺服系统带宽、伺服电机频响带宽的测试
伺服系统带宽 原文:https://wenku.baidu.com/view/bde80b1b650e52ea551898ab.html 电机驱动器如何测试? https://zhidao.baidu ...
- matlab导纳圆,分析分量法及导纳圆分析法
2.3.2分析分量法 由于试验条件的关系,本文采用频域识别的方法对机床的模态参数进行识别, 下面简单介绍在频域中的一些模态参数识别方法. 分量分析法就是将频响函数分成实部和虚部分量进行分析,式(2-1 ...
- 频响测试低12dB问题
find . -name audio_device.xml | xargs ls -l audio_device.xml snd_kcontrol_new iface:control类型,通常是SND ...
- 位移传递率matlab编程,各种谱计算,频响函数,传递率
A.信号与谱的分类 由于时域信号有不同的分类, 变换后对应的频域也有不同的谱 信号可分为模拟(连续)信号和数字(离散)信号, 连续信号变换后称为谱密度, 离散信号变换 后称为谱. 连续信号又可分为绝对 ...
最新文章
- 清华自动驾驶前沿报告!解密六大关键技术,全球人才分布【附下载】| 智东西内参...
- SQL对象名无效怎么办法 华骨龙
- 阿里DataV可视化大屏介绍
- windows消息队列C语言,08年计算机二级C语言辅导:vc操作微软消息队列的方法
- mysql索引创建及使用注意事项
- 关于MATLAB fplot 出现matlab 警告: 函数处理数组输入时行为异常。要改善性能,请将您的函数正确向量化,以返回大小和形状与输入参数相同的输出。
- DataGrid显示字段过长用...代替
- 天行数据API智能机器人接口
- fastboot简介(android烧录指令)
- Django 创建项目app
- Cortex-M3 Bit-Banding
- SSR门户项目爬坑之路(三)
- NameNode HA的部署方法
- wpa_supplicant wpa_cli 无线网络配置
- 创业维艰---分享这三年的揪心历程
- 基于eBox旋转编码器
- virtual方法和abstract方法的使用(轉載)
- 【时间和空间复杂度】
- 农村土地确权之例会纪要—— 新蔡县土地确权第七次例会会议纪要
- osg 三维gis开发_三维GIS与游戏引擎的跨界融合,打造数字化孪生的平行世界