非线性振动 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日,山东大学数学与交叉科学博士(后)论坛会议在青岛校区华岗苑东楼举行,本次论坛由数学与交叉科学研究中心(以下简称交叉中心)与山东大学数学学院联合主办,会议由交叉中心副主任栾贻会主持. ...
最新文章
- 使用python完成的一个烟花小程序-人人都可以写的可视化Python小程序第二篇:旋转的烟花...
- flash builder 关联svn
- day13 Java学习(常见对象正则表达式)
- java intern_java String的intern方法
- php ci 处理图片 裁剪,Codeigniter实现智能裁剪图片的方法
- oracle查询排序asc/desc 多列 order by
- [转贴]使用jQuery自动缩图片 - [jQuery]
- [转载] Python 继承
- php 当前文件的文件夹大小,PHP获取文件夹大小函数用法实例
- 图解设计模式阅读(概述)
- 【杂项】VScode中C C++代码美化工具Uncrustify配置详解
- c#winform——Gobang五子棋简易版双人对战制作(基本结构+代码)
- Git将分支代码合并到主干/将主干代码合并到分支
- 金蝶云苍穹集成服务云,让企业数据自由流动
- 最新CCF会议|2022-2023顶会会议时间+投稿时间+官网链接(视觉+多媒体+数据挖掘+数据库+通用人工智能)
- 每日刷题记录 (十七)
- Oracle的hint
- Linux安装7z压缩软件记录
- 戴尔sc系列存储阵列柜服务器,Dell存储Compellent SC4020
- 漂流瓶php源码,微信小程序之仿微信漂流瓶