应用之动力学与振动.pptVIP

上一页 目录 下一页 返回 * 优选应用之动力学与振动 当前第1页\共有32页\编于星期二\7点 目录 6.1 轨迹 6.2 单自由度系统 6.3 多自由度系统 习题 当前第2页\共有32页\编于星期二\7点 6.1 轨迹 举例说明:重力场中有两个物体,其中质量为m2的物体固定,而质量为m1的物体绕m2做平面圆周运动.做圆周运动的m1物体的轨道半径用变量r表示,角度用变量a表示. m2 m1 a r 当前第3页\共有32页\编于星期二\7点 6.1 轨迹 例6.1：卫星绕地球转动时，m2等于地球的质量，m1等于卫星的质量，r为卫星球心与地球球心间的距离。其运动轨迹由下列方程组决定： 式中： ，其中t是时间变量，p为物体在地球表面做圆周运动的周期。在地球表面，r=6.373x106 m。 当前第4页\共有32页\编于星期二\7点 6.1 轨迹 用龙格—库塔法可以实现求解： 引入新状态变量： 当前第5页\共有32页\编于星期二\7点 建立函数文件Orbit.m function xd=Orbit(t,x) xd=[x(2) x(1)*x(4)^2-4.0*pi^2/x(1)^2 x(4) -2.0*x(2)*x(4)/x(1)]; 6.1 轨迹 组 X1初始 X2初始 X3初始 X4初始 轨迹类型 1 2 0 0 1.5 椭圆 2 1 0 0 2pi 圆 3 2 0 0 4 双曲线 三组初始条件(t=0)： 当前第6页\共有32页\编于星期二\7点 由初始条件建立执行文件execute_61.m initcond=[2 0 0 1.5;1 0 0 2*pi;2 0 0 4]; tspan=linspace(0,5,1000); options=odeset(RelTol,1e-6,AbsTol,[1e-6 1e-6 1e-6 1e-6]); lintype=[‘k- ‘b-. ‘r--]; for i=1:3 [t,x]=ode45(‘Orbit,tspan,[initcond(i,:)],options); polar(x(:,3),x(:,1),lintype(2*(i-1)+1:2*i)); hold on end text(0.5,-1.2,椭圆轨迹); text(-1.2,1,圆轨迹); text(1.75,2,双曲线轨迹); 6.1 轨迹 常微分方程的数值求解函数 当前第7页\共有32页\编于星期二\7点 程序运行结果 6.1 轨迹 当前第8页\共有32页\编于星期二\7点 6.2 单自由度系统 6.2.1 概述 一.力学模型 m c K，a X（t） F（t）=X（0）kf（t） 弹簧—质量—阻尼系统 其中：振体质量为m，弹簧的线性系数为k，非线性系数为a，阻尼系数为c，外力F（t）。 当前第9页\共有32页\编于星期二\7点 6.2 单自由度系统 二. 运动微分方程 用x表示系统的位移，则运动微分方程为： 式中： 固有频率: 非线性系数: 阻尼因子: 当前第10页\共有32页\编于星期二\7点 6.2 单自由度系统 引入新变量转化状态空间方程形式： 当前第11页\共有32页\编于星期二\7点 6.2 单自由度系统 6.2.2 线性系统的自由振动 一.运动微分方程 当 时，得到线性振动系统的自由振动方程。 当前第12页\共有32页\编于星期二\7点 6.2 单自由度系统 二.MATLAB求解 编写方程对应的函数文件FreeOscillation.m 0 三种阻尼系数( )（1）阻尼系数为0.1时是欠阻尼情况（2）阻尼系数为1时是临界阻尼情况（3）阻尼系数为5时是过阻尼情况 function xdot=FreeOscillation(t,x,zeta,Alpha) xdot=[x(2);-2.0*zeta*x(2)-x(1)-Alpha*x(1)^3]; end 当前第13页\共有32页\编于星期二\7点 6.2 单自由度系统 由初始条件（位移和速度均为1时, ）建立执行文件(execute_62.m) zeta=[0.1 1.0 5.0]; Alpha=[0.0,0.0,0.0]; tspan=linspace(0,40,400); %生成0-40的四百个线性点 lintype=char(-k,--k,-.k); for i=1:3 [t,x]=ode45(@FreeOscillation,tspan,[1 1],[],zeta(i),Alpha(