matlab多自由度弹簧响应,用matlab编程实现法计算多自由度体系的动力响应..doc
用matlab编程实现法计算多自由度体系的动力响应.
用matlab编程实现法计算多自由度体系的动力响应
姓名:王 辉
学号 班级:硕1205班
专业:桥梁与隧道工程
2012年11月27日
用matlab编程实现法
计算多自由度体系的动力响应
一、法的基本原理
Newmark-?法是一种逐步积分的方法,避免了任何叠加的应用,能很好的适应非线性的反应分析。
Newmark-?法假定:
(1-1)
(1-2)
式中,?和?是按积分的精度和稳定性要求进行调整的参数。当?=0.5,?=0.25时,为常平均加速度法,即假定从t到t+?t时刻的速度不变,取为常数。研究表明,当?≥0.5,??≥0.25(0.5+?)2时,Newmark-?法是一种无条件稳定的格式。
由式(2-141)和式(2-142)可得到用及,,表示的,表达式,即有
(1-3)
(1-4)
考虑t+?t时刻的振动微分方程为:
(1-5)
将式(2-143)、式(2-144) 代入(2-145),得到关于ut+?t的方程
(1-6)
式中
求解式(2-146)可得,然后由式(2-143)和式(2-144)可解出和。
由此,Newmark-?法的计算步骤如下:
1.初始计算:
(1)形成刚度矩阵[K]、质量矩阵[M]和阻尼矩阵[C];
(2)给定初始值, 和;
(3)选择积分步长?t、参数?、?,并计算积分常数,,,,
,,,;
(4)形成有效刚度矩阵;
2.对每个时间步的计算:
(1)计算t+?t时刻的有效荷载:
(2)求解t+?t时刻的位移:
(3)计算t+?t时刻的速度和加速度:
Newmark-?方?t的大小不影响解的稳定性,?t的选择主要根据解的精度确定。
本文用法计算的基本问题
四层框架结构在顶部受一个简谐荷载的作用,力的作用时间=5s,计算响应的时间为100s,分2000步完成。阻尼矩阵由Rayleigh阻尼构造。
具体数据如下图:
图一:结构基本计算简图
计算法的源程序
m=[1,2,3,4];
m=diag(m);
k= [800 -800 0 0;
-800 2400 -1600 0;
0 -1600 4800 -3200;
0 0 -3200 8000];
c=0.05*m+0.02*k;
f0=100;
t1=5;
nt=2000;
dt=0.01;
alfa=0.25;
beta=0.5;
a0=1/alfa/dt/dt;
a1=beta/alfa/dt;
a2=1/alfa/dt;
a3=1/2/alfa-1;
a4=beta/alfa-1;
a5=dt/2*(beta/alfa-2);
a6=dt*(1-beta);
a7=dt*beta;
d=zeros(4,nt);
v=zeros(4,nt);
a=zeros(4,nt);
for i=2:nt
t=(i-1)*dt;
if (t
ke=k+a0*m+a1*c;
fe=f+m*(a0*d(:,i-1)+a2*v(:,i-1)+a3*a(:,i-1))+c*(a1*d(:,i-1)+a4*v(:,i-1)+a5*a(:,i-1));
d(:,i)=inv(ke)*fe;
a(:,i)=a0*(d(:,i)-d(:,i-1))-a2*v(:,i-1)-a3*a(:,i-1);
v(:,i)=v(:,i-1)+a6*a(:,i-1)+a7*a(:,i);
end
计算结果截图
最后程序分别计算出四个质点的位移、速度、加速度响应。
现将部分截图如下:
1、位移响应:
图二:1质点的位移响应
图三:4质点的位移响应
2、速度响应
图四:1质点的速度响应
图五:4质点的速度响应
3、加速度响应
图六:1质点的加速度响应
图七:4质点的加速度响应
硕1205班 王辉 1
matlab多自由度弹簧响应,用matlab编程实现法计算多自由度体系的动力响应..doc相关推荐
- 网络拓扑 代码 matlab,matlabdaima 复杂网络拓扑特征统计的编程代码,可以计算 中各种 值 Other systems 其他 246万源代码下载- www.pudn.com...
文件名称: matlabdaima下载 收藏√ [ 5 4 3 2 1 ] 开发工具: matlab 文件大小: 85 KB 上传时间: 2015-11-05 下载次数: 0 提 供 者: ...
- 用matlab分析时间响应教程,基于Matlab的多自由度耦合滑移模型的动力响应可靠度分析...
Hi~我是mania 参考资料: (1)基于耦合滑移模型的岩体地震可靠性分析方法 (3)结构可靠度分析:方法与程序 简化模型及运动平衡方程 如下图1所示的多层隔震建筑,当受到地震作用时每层结构体都会产 ...
- 主振型 matlab 振动,基于MATLAB语言的多自由度振动系统的固有频率及主振型计算分析...
基于 MATLAB 语言的多自由度振动系统的固有频率及主振型计算分析 文 涛 ,胡青春 (华南理工大学 机械工程学院 ,广东 广州 510640) 摘要 :多自由度振动系统固有频率及主振型计算分析是研 ...
- 位移传递率matlab编程,机械振动设计计算与VB编程实例
1 编写机械振动计算程序的一般步骤1 1.1 明确编写程序的目的和了解清楚相关的计算公式.物理量纲1 1.2 在简单的界面上用BASIC编写计算程序并完成调试2 1.3 设计一个较为完备的用户界面4 ...
- matlab二阶阻尼震荡衰减,二阶欠阻尼电路的零响应输入-MATLAB课程设计.doc
二阶欠阻尼电路的零响应输入-MATLAB课程设计.doc 摘要 MATLAB具有强大的运算和作图功能,为电路分析引入计算机数值方法提供了技术保证.MATLAB中的变量与常量都是矩阵,其元素可以使复数和 ...
- 中心差分法 matlab,(最新整理)结构动力学使用中心差分法计算单自由度体系动力反应的MATLAB程序...
<(最新整理)结构动力学使用中心差分法计算单自由度体系动力反应的MATLAB程序>由会员分享,可在线阅读,更多相关<(最新整理)结构动力学使用中心差分法计算单自由度体系动力反应的MA ...
- 杜哈梅 matlab,应用MatLab软件探讨结构动力响应时域和频域数值模拟教学
中图分类号:G642.0 文献标志码:A 文章编号:1005-2909(2017)03-0119-05 随着计算机技术的发展,数值模拟已经成为结构科学研究的基本方法,也是重要的教学内容.结构动力分析是 ...
- matlab麦克斯韦速率,MATLAB编程首先建立计算麦克斯韦分布律的子程序mxwm-Read.PPT...
MATLAB编程首先建立计算麦克斯韦分布律的子程序mxwm-Read * * 第5章 综合实例 用 MATLAB 求解问题时,一般要经历建模和编程两个过程,只有在建模正确的前提下,方能得出正确的结果. ...
- matlab离散系统 响应,离散系统的频率响应和输出响应的matlab实现
离散系统的频率响应和输出响应的matlab实现 上海电力学院 信号与系统实验报告 题目:离散系统的频率响应和输出响应 班级: 2011023 专业:电气工程及其自动化 学号:20111257 2013 ...
最新文章
- 在写游戏时钟类时,应确保时钟的计算是以某个固定的CPU为标准的
- 人脸对齐--Robust face landmark estimation under occlusion
- mac 查看端口_如何重置mac上的系统管理控制器smc教程
- UVa 11107 (后缀数组 二分) Life Forms
- Redis Scan命令
- 两张超级大表join优化
- Java 里的HashMap(HashTable) 简介.
- java中的JDBC用户管理系统_Java+JDBC学生信息管理系统
- centos php redhat,RHEL / CentOS 安装 OPcache 提升 PHP 效能
- java date.from_java datefromat
- oracle触发器无效且,oracle创建触发器成功但是插值失败
- php session 为空,ThinkPHP_session问题_dump显示session为空
- tensorflow学习之常用函数总结:tensorflow.reduce_mean()函数
- PL/pgSQL的RETURN QUERY例子
- mpvue 中使用 iconfont
- IDEA中如何配置Tomcat和项目?
- 苹果电脑上几款不错的图片编辑工具
- 闫令琪图形学入门笔记(光栅化篇)
- 苹果电池显示维修_iFixit拆解苹果iPhone 12/Pro:显示屏和电池可互换
- 3分钟带你了解Web前端开发工程师需要掌握什么?
热门文章
- 命令行清除redis缓存
- 美国网红python微博_利用Python对微博网红进行对比
- 使用DWM API画自定义Aero效果窗口
- vmware linux centos 安装 vscode 打不开怎么办 No protocol specified Unable to open X display.
- [3月总结,4月计划]
- 小程序源码:炫酷手持滚动弹幕生成小工具微信小程序-多玩法安装简单
- poj2228(dp)
- PDF怎么编辑,如何去除PDF水印
- c# php 时间戳,C#实现时间戳的简单方法
- Spring——如何理解Spring,DI和IoC