单自由度系统中质量、阻尼和刚度变化对频率响应函数(FRF)影响图的绘制

作者:赵兵

日期:2020-02-17



目录

1.     背景

2.     VISIO绘制

3.     Matlab绘制

  (1)     M变化时

  (2)     K变化时

  (3)     C变化时

4.     参考文章


1. 背景

写文章时需要用到几张图,下面是从PDF上截图截出来的,用来表示单自由度系统在冲击激励下的频率响应曲线,当K(刚度),C(阻尼),M(质量)变化时,频率响应曲线的变化情况。

     

图 1 单自由度系统刚度,阻尼,质量影响曲线

用图1放在文章中,不太美观,一看就是影印的,要是一般的文章还好,如果用来发表的就拉下档次了。所以就尝试了下重绘,重绘有两个方案,

(1)用PPT、visio或者Adobe Illustrator等绘图工具自己绘制

(2)用matlab先得到曲线的函数,进一步把函数显示出来

2. Visio绘制

尝试了一下,质量变化影响曲线图相对还好画一些,但阻尼变化的曲线就不是那么容易画了,画出来总觉得差点意思。很难保证各曲线之间的间隔均匀;当然肯定是能做到,但自己非此方面的熟手,所以有两个解决方案,一个是尝试用Adobe Illustrator进行一定时间的专门学习,但时间成本太高;或者是淘宝上寻求供外包解决,也是能解决的;但最终我还是决定先尝试下第二种方案,先得到曲线的函数,然后用软件把函数绘出来。

  

图 2 Visio绘制曲线图

3. Matlab绘制

单自由度系统的物理模型如图 3 所示,

 

图 3单自由度系统模型

它的动力学方程为

其中

M:质量;

C:阻尼;

K:刚度

:加速度;

:速度;

x:位移;

f:外力;

t:时间。

这里的f(t)为脉冲激励:

% % 激励采用脉冲激励,脉冲激励为Hanning函数
function f = hanning_imp(t, Tc, A)
    f = zeros(size(t));
    f(t < Tc) = A / 2 * (1 - cos(2*pi * t(t < Tc) / Tc));
end

可以根据这个函数得到一个脉冲激励

画图此图:

dt = 0.00001;
t = 0:dt:200;
Tc = 0.001;
A = 10;
u = hanning_imp(t, Tc, A);
plot(t,u,‘LineWidth‘,1.5 );
axis([-5 200 -0.1 10.5 ]);
xlabel(‘t/s‘)
ylabel(‘Amp/N‘)
text(75,8,‘Impact Force‘);

图 4 冲击力

建立系统方程,求解频率响应函数(FRF)

function [freq_x , amp_y]=frf_bing(m,k,c)

% m 质量
% k 刚度
% c 阻尼

num = 1;
den = [m c k];
sys = tf(num, den);

%采样频率(Hz) 100Hz 实际并不需要这么高的采样频率,但是如果采样时间太小,hanning脉冲不完整
% 为了得到准确的响应dt一定要小,否则做出的相位可能不对
dt = 0.00001;
fs = 1/dt;  

t = 0:dt:200;
Tc = 0.001;
A = 10;
u = hanning_imp(t, Tc, A);
y = lsim(sys, u, t);
y = y‘;

N = length(u);
fy = fft(y);
fu = fft(u);
ft = fy ./ fu;
f = (0:N-1) * fs ./ N;
ft_r = real(ft);
ft_i = imag(ft);

part = (f < 30);
freq_x=f(part);
amp_y= abs(ft(part));

End

(1) M变化时

取m=100,120,140,160,180,200分别绘制FRF的响应曲线。

clc;
clear;
close all;

%%
k = 1000;      %初始化k
c=100;         %初始化c
M=100:20:200;  %初始化m, 取m=100,120,140,160,180,200分别绘图

f1= figure(1);
hold on
    for i= 1:length(M)
    [a(:,i) , b(:,i)]=frf_bing(M(i),k,c);
end

plot(a,log(b),‘b‘);
axis([0 1 -8.2 -5.3]);
title(‘Log(FRF)‘);
xlabel(‘Frequency‘)
ylabel(‘Amplitude‘)
set(gca,‘XTick‘,[],‘YTick‘,[]);
text(0.3,-7,‘M↑‘);
annotation(‘arrow‘,[0.7 0.3],[0.6 0.8]) ;
f1.Position=([ 0 0 400 300])

% 保存为emf 矢量格式
set(gcf,‘unit‘,‘centimeters‘,‘position‘,[10 5 6.5 4.8]);
print(f1,‘-dmeta‘,‘M.emf‘)

可以生成一幅想要的曲线,并且可以保存为矢量格式,无论放大多少倍,图片还是很清晰。

  

图 5 质量变化的影响

(2)   K变化时

取K=1000,1200,1400,1600,1800,2000分别绘制FRF的响应曲线。

clc;
clear;
close all;

%%
k = 1000:200:2000;
c=100;
M=100;
f1= figure(1);
hold on
for i= 1:length(k)
[a(:,i) , b(:,i)]=frf_bing(M,k(i),c);
end

plot(a,log(b),‘b‘);
axis([0 1 -8.2 -5.3]);
title(‘Log(FRF)‘);
xlabel(‘Frequency‘)
ylabel(‘Amplitude‘)
set(gca,‘XTick‘,[],‘YTick‘,[]);
text(0.6,-7,‘K↑‘);
annotation(‘arrow‘,[0.4 0.8],[0.63 0.6]) ;
f1.Position=([ 0 0 400 300])

% 保存为emf 矢量格式
set(gcf,‘unit‘,‘centimeters‘,‘position‘,[10 5 6.5 4.8]);
print(f1,‘-dmeta‘,‘K.emf‘)


图 6刚度变化的影响

(3)   C变化时

取K=1000,1200,1400,1600,1800,2000分别绘制FRF的响应曲线。

clc;
clear;
close all;
%%
fontsizevalue=18;
c=100:20:200;
f1= figure(1);
hold on
for i= 1:length(c)
[a(:,i) , b(:,i)]=frf_bing(c(i));
end

plot(a,log(b),‘b‘);
axis([0 1 -8 -5.5]);
t1=title(‘Log(FRF)‘);
xl=xlabel(‘Frequency‘)
y1=ylabel(‘Amplitude‘)
t1.FontSize =fontsizevalue;
y1.FontSize =fontsizevalue;
xl.FontSize =fontsizevalue;
set(gca,‘XTick‘,[],‘YTick‘,[]);
text(0.5,-7,‘C↑‘);
annotation(‘arrow‘,[0.5 0.5],[0.9 0.5]) ;
f1.Position=([ 0 0 400 300])

% 保存为emf 矢量格式
set(gcf,‘unit‘,‘centimeters‘,‘position‘,[10 5 6.5 4.8]);
print(f1,‘-dmeta‘,‘K.emf‘)
 

图 7阻尼变化的影响

4.     参考文章

1. CSDN博主「whoispo」文 https://blog.csdn.net/WhoisPo/article/details/46865401

原文地址:https://www.cnblogs.com/Nicoooolas/p/12321596.html

时间: 2024-09-30 14:22:50

单自由度系统中质量、阻尼和刚度变化对频率响应函数(FRF)影响图的绘制的相关文章

多功能PCIE交换机之九:单NT系统中需要注意的问题

多功能PCIE交换机之九:单NT系统中需要注意的问题 1.单NT与双NT系统差异 NTB通常用在双控或者多控的系统上,用来实现跨节点的数据传输.在由多个节点构成的系统上,通常每个节点上都有一个NTB芯片,以实现全系统地址空间划分和路由.顾名思义,单NT是指用一个PCIE非透明桥来实现节点间数据传输的方案,通常用在两个节点之间.和多NT的系统相比,这种方法实现相对比较简单,成本较低,开发周期较短,能够缩短整个系统的关键功能和性能的验证周期. 2.单NT系统的实现 单NT的实现需要结合整个系统的设计

【转载】Ansys中的阻尼

原文地址:http://www.cnblogs.com/ylhome/archive/2009/08/26/1554195.html ANSYS动力学分析中提供了各种的阻尼形式,这些阻尼在分析中是如何计算,并对分析有什么影响呢?本文将就此做一些说明何介绍. 一. 首先要清楚,在完全方法和模态叠加法中定义的阻尼是不同.因为前者使用节点坐标,而后者使用总体坐标. 1. 在完全的模态分析.谐相应分析和瞬态分析中,振动方程为: 阻尼矩阵为下面的各阻尼形式之和: α为常值质量阻尼(α阻尼)(ALPHAD命

理论力学单自由度小阻尼受迫振动阻尼系数近似公式的推到

对于单自由度小阻尼受迫振动阻尼系数近似公式的推到有许多方法: 1.可以从基本表达式出发,一点一点的推到,也就是从振幅的表达式出发,求解出来最大振幅频率,然后求出最大振幅的表达式,让振幅等于最大振幅的一半,求解两个驱动频 率,由表达式一步步的变换就可以得到解析解,然后近似变换得到近似公式,但是过程有点繁琐,而且近似变换放在最后,没有发挥应有的功能 2.也可以由振动方程的形式,直接从能量的角度出发,应用积分方法与公式变换加上近似变换,舍去高阶小项得到 3.从振幅表达式出发,直接推导,中间过程近似舍去

工单在客服系统中的应用

从国内来看,目前企业的线上客服主要有 三种模式,在线客服.机器人客服.工单.前者是基于人工坐席会话的基础上,提供一对一的即时聊天会话,即时接起客户服务请求,后者则可适用于离线处理. 在线客服是国内目前应用非常广泛的服务模式,但也因为服务成本高,服务水平参差不齐而为企业诟病.机器人客服在很大程度上能够缓解客服人员的压力,同时标准化的客户会话语术,能有效减少因沟通过程中带来的理解误差以及因客服人员语气或态度差带来的客服投诉. 工单最常见于IT运维支持,如用于处理内部员工的设备报障问题.采购需求等.用

系统中生成编号/单号问题的实现方案讨论

应用场景 场景:对于大多数电商系统或财务系统来说,系统中的单号一般都不是从1开始的自增数字,而是一串有一定意义的字符串序列. 而往往这样的单号是要全局唯一的,不可重复. 那么,每次新增订单记录时,这个单号就要按照指定的规则来生成. 常见的订单号规则是字母前缀+日期+时间+定长的数字,如DD201610201559060001(format:XXYYYYMMDDHHmmssNNNN),PZ16102015000012(format:XXYYMMDDHHNNNNNN) 常见方案... 实现这样的场景

Android 系统中,那些能大幅提高工作效率的 API 汇总(持续更新中...)

前言 "条条大路通罗马."工作中,实现某个需求的方式往往不是唯一的,这些不同实现方式不仅表现在代码质量上,还影响着我们的工作效率.就像,在 Android 系统中,总有那么一些鲜为人知的 API 能够减少我们很多零碎的工作量.于是,就想凭着一些经验,整理一些常用的,找个地方归纳总结,也供日后翻阅. getResources().getIdentifier(String name, String defType, String defPackage) 根据资源名称获取资源 id.正常情况

幂等和高并发在电商系统中的使用

在Java web项目开发中,经常会听到在做订单系统中生成订单的时候,要做幂等性控制和并发控制,特对此部分内容作出总结,在高并发场景下,代码层面需要实现并发控制:但是幂等性,其实更多的是系统的接口对外的一种承诺,承诺一次请求和多次请求会返回同样的数据.关于幂等性将分别从高等代数中的幂等性.HTTP中的幂等性和订单生成系统中的幂等性阐述:并发性控制则提供了分布式锁等方式来对并发场景进行代码实现. 一.幂等性 idempotence  ['a?d?mpo?t?ns] 1.高等代数中关于幂等idemp

搭建属于你的在线实时采集系统 ——HTML5 在嵌入式系统中的应用

※已刊登在<无线电>04月刊上   搭建属于你的在线实时采集系统--HTML5 在嵌入式系统中的应用   作者:刘琛,徐洋   摘要: 本应用摆脱了以往嵌入式系统的数据采集方式,借助于最新的HTML5的Canvas API及WebSocket API两大特性,实现了数据的在线实时采集功能.提升了嵌入式采集系统的性能及体验.为嵌入式开发工作者提供参考.   关键字:HTML5:HTTP Server:Canvas:WebSocket:W5500:实时:采集系统:      当今信息社会,信息就是

管道设计CAD系统中重量重心计算

管道设计CAD系统中重量重心计算 [email protected] Abstract. 管道设计CAD系统中都有涉及到重量重心计算的功能,这个功能得到的重心数据主要用于托盘式造船时方便根据重心设置吊装配件.重量信息主要用于采购订货.本文主要介绍相关软件中重量重心功能,及重量重心计算实现原理.最后结合OpenCASCADE计算管道模型重量重心来验证. Key Words. CoG, CentreOfMass, Piping CAD, Piping Design 1.Introduction 船舶