非线性振动 matlab,数学与非线性科学 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断 - Powered by Discuz!...
17.png (2.32 KB, 下载次数: 4)
2018-2-28 16:07 上传 Matlab代码:
clear all
clc
tic
%定义各参数
syms t
w0=3;
epsR=0.001;
m=[1 0;0 1];
epsilon=0.24;
r=0.4;
delta=0.56;
k=[1+epsilon*r -epsilon*r;-epsilon*r 1+epsilon*delta];
Cs=[cos(t) cos(3*t) sin(t) sin(3*t)];
Cs1=diff(Cs,t,1);
S=[cos(t) cos(3*t) sin(t) sin(3*t) 0 0 0 0;0 0 0 0 cos(t) cos(3*t) sin(t) sin(3*t)];
A1=[1 1 1 1]';
A2=[1 1 1 1]';
A0=[A1;A2];
T1=[eye(4,4) zeros(4,4)];
T2=[zeros(4,4) eye(4,4)];
S2=diff(S,t,2);
fm=inline(S'*m*S2);
M=quadv(fm,0,2*pi);
fk=inline(S'*k*S);
K=quadv(fk,0,2*pi);
S1=diff(S,t,1);
c=[1 0;0 1];
fc=inline(S'*c*S1);
C=quadv(fc,0,2*pi);
c3=diag(S*A0).^2;
fc3=inline(S'*c3*S1);
C3=quadv(fc3,0,2*pi);
k2=diag(S*A0).*diag(S1*A0);
fk2=inline(S'*k2*S);
K2=quadv(fk2,0,2*pi);
%代入推导出的公式
Kmc=w0^2*M+epsilon*w0*(C3-C)+K+2*epsilon*w0*K2;
R=-(w0^2*M+epsilon*w0*(C3-C)+K)*A0;
Rmc=-(2*w0*M+epsilon*(C3-C))*A0;
%AA首元素已知a1=0.0,求ww
a1=0.0;
%变换矩阵,使ww变量代替a1
Kmc11=-Rmc(:,1);
Kmcr=[Kmc11 Kmc(:,2:size(Kmc,2))];
%求未知变量
AA=inv(Kmcr)*R;
%drtA1(1)
ww=AA(1);
%drtW(1)
%赋予新变量新值
A01=A0+[a1; AA(2:length(A0),1)];
%A(1)+drtA(1)
% Aw0=AA+A00;
%A1(0)+drtA1(1)=A1(1)
w01=w0+ww;
%W+drtW(1)
n=1;
tol=1;
while tol>epsR
A0=A01;
w0=w01;
c3=diag(S*A0).^2;
fc3=inline(S'*c3*S1);
C3=quadv(fc3,0,2*pi);
k2=diag(S*A0).*diag(S1*A0);
fk2=inline(S'*k2*S);
K2=quadv(fk2,0,2*pi);
%带入推导出的公式
Kmc=w0^2*M+epsilon*w0*(C3-C)+K+2*epsilon*w0*K2;
R=-(w0^2*M+epsilon*w0*(C3-C)+K)*A0;
Rmc=-(2*w0*M+epsilon*(C3-C))*A0;
%%%%%
tol=norm(R);
if(n>1000)
disp('迭代步数太多,可能不收敛')
return;
end
Kmc11=-Rmc(:,1);
Kmcr=[Kmc11 Kmc(:,2:size(Kmc,2))];
AA=inv(Kmcr)*R;
ww=AA(1);
%A00=[w0;A0(2:6,1)];
A01=A0+[a1;AA(2:length(A0),1)];
w01=w0+ww;
n=n+1;
end
X0=S*A0;
dX0=S1*A0;
%绘范德波图
tt=0:.1:10;
xo1=subs(X0(1),tt);
xo2=subs(X0(2),tt);
dxo1=subs(dX0(1),tt);
dxo2=subs(dX0(2),tt);
figure(1)
plot(xo1,dxo1,'b','linewidth',2)
hold on
plot(xo2,dxo2,'b','linewidth',2)
axis([-3 3 -3 3])
title('范德波极限环')
xlabel('x0')
ylabel('dx0')
toc
运行结果:范德波极限环
18.png (38.84 KB, 下载次数: 0)
2018-2-28 16:08 上传 ——以上代码由声振之家会员zhangwenjing分享,代码未经验证。
非线性振动 matlab,数学与非线性科学 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断 - Powered by Discuz!...相关推荐
- 计算机科学与技术syc,数学与非线性科学 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断 - Powered by Discuz!...
续上... ======================================================================== 多物理多尺度数值模拟第四届国际研讨会将在北 ...
- 理想点法matlab p趋向无穷大,Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断 - Powered by Discuz!...
7.png (183.66 KB, 下载次数: 0) 2018-1-25 16:06 上传 图3 实验所用的随机信号 采样点数N分别取128.256.512和1024,周期图法matlab代码如下: ...
- matlab vdp1000,Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断 - Powered by Discuz!...
搜狗截图20140423191008.png (21.33 KB, 下载次数: 1) 2014-4-23 19:10 上传 程序如下: syms t y; u=exp(-5*t)*cos(2*t-1) ...
- iso 2631 matlab,声学基础及其分析软件 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断 - Powered by Discuz!...
人体振动测试仪MAESTRO 特点 多种应用场所 4通道同时测量 数字滤波器 坚固 容易使用 完整 噪声选项 现场工作中的振动越来越受到关注,原因是逐渐增多的工伤是由高量级的振动引起的. 很多的产品, ...
- 悬置位移matlab,Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断 - Powered by Discuz!...
clear clc %动力总成质量 m=323.98;%kg %惯性积 Jx=25.706;Jy=13.607;Jz=23.014;Jxy=-3.029;Jyz=3.359;Jzx=-0.876;%k ...
- MATLAB里面的mghglobal函数,Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断 - Powered by Discuz!...
function exam8_2 % 本程序为第八章的第二个算例,采用平面梁单元计算两铰抛物线拱的在初始条件下 % 自由振动,并对时程曲线结果进行FFT变换,求得的频率可与exam8_1.m的结果进 ...
- matlab 除去振动,Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断 - Powered by Discuz!...
下面是徐荣桥老师书中exam8_2例子中的一段,不知对你有用否 step3. 计算时程响应(Newmark法) % step3.1 初始计算 gama = 0.5 ; beta = 0.25 ; C ...
- 没有shell63号单元_Ansys讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断 - Powered by Discuz!...
水下圆柱壳体的建模与声学分析 使用有限元软件ANSYS进行计算和分析时水下环肋圆柱壳体有限元模型的建立及结构声学分析主要分为以下一些步骤: 1.建立壳体的实体模型(包括有圆柱壳体的建立,给圆柱壳体加环 ...
- 生物数学计算机交叉科学论坛,数学与交叉科学博士(后)论坛在山东大学(青岛)举行...
[本站讯]11月8日,山东大学数学与交叉科学博士(后)论坛会议在青岛校区华岗苑东楼举行,本次论坛由数学与交叉科学研究中心(以下简称交叉中心)与山东大学数学学院联合主办,会议由交叉中心副主任栾贻会主持. ...
最新文章
- 关于神经网络模型的一些入门知识
- unity3d游戏开发第2版 pdf_学开发2个月,1个人做出个TapTap排行榜第4的游戏?
- Laravel_post访问【419】问题解决方法
- Java多线程:线程安全和非线程安全的集合对象
- PHP+jquery 树状菜单
- python3.6 websocket异步高并发_在Python3.6上的websocket客户端中侦听传入消息时出现问题...
- python数据分析函数大全_python中数据分析常用函数整理
- Kafka : 报错 KafkaController NoSuchElementException : : key not found : [xxx]
- LR 报错误: C interpreter run time error: Error -- Unresolved symbol : vuser_end解决方法
- 算法学习笔记:连通图详解
- 市场调研-全球与中国矢量网络分析仪(VNA)校准套件市场现状及未来发展趋势
- 计算机硬盘有坏道,硬盘有坏道怎么办 硬盘坏道分类
- 新计算机的word无法输入文字,【WORD为什么无法输入文字?】word如何转换excel
- identifier标签语句
- 华为v2服务器系统安装系统,华为服务器RH 2288H v2安装系统
- 怎样使用ApowerMirror实现将手机屏幕投屏到电脑
- 机器学习之三:降维技术
- WiFi穿墙完全指南
- 数据链路层---差错检测和纠正
- 电脑操作最忌讳18个小动作