结构动力学MATLAB编程例题,华工结构动力学编程 MATLAB
华工结构动力学编程
一.所编程序及编译:
clear%***********读入地震记录***********
fid=fopen('ei.txt');
[Accelerate,count]=fscanf(fid,'%g');%count读入的记录的量
time=0:0.02:(count-1)*0.02;
%***********线性加速度法计算各反应***********
%初始化各储存向量
Displace=zeros(1,count);%相对位移
Velocity=zeros(1,count);%相对速度
AbsAcce=zeros(1,count);%绝对加速度
Damp=0.05;%结构阻尼比取为0.05
Tc=0.0:0.05:10;%结构自振周期
Dt=0.02;%地震记录的步长
%记录计算得到的反应,MDis为最大相对位移,MVel为最大相对速度%MAcc某阻尼时最大绝对加速度,用于画图
MDis=zeros(1,length(Tc));
MVel=zeros(1,length(Tc));
MAcc=zeros(1,length(Tc));
t=1;%在下一个循环中控制不同的结构自振周期
forT=0.0:0.05:10
Frcy=2*pi/T;%结构自振频率
DamFrcy=Frcy*sqrt(1-Damp*Damp);%计算公式化简
e_t=exp(-Damp*Frcy*Dt);
s=sin(DamFrcy*Dt);
c=cos(DamFrcy*Dt);
A=zeros(2,2);
A(1,1)=e_t*(s*Damp/sqrt(1-Damp*Damp)+c);
A(1,2)=e_t*s/DamFrcy;
A(2,1)=-Frcy*e_t*s/sqrt(1-Damp*Damp);
A(2,2)=e_t*(-s*Damp/sqrt(1-Damp*Damp)+c);
d_f=(2*Damp^2-1)/(Frcy^2*Dt);%计算公式化简
d_3t=Damp/(Frcy^3*Dt);
B=zeros(2,2);
B(1,1)=e_t*((d_f+Damp/Frcy)*s/DamFrcy+(2*d_3t+1/Frcy^2)*c)-2*d_3t;B(1,2)=-e_t*(d_f*s/DamFrcy+2*d_3t*c)-1/Frcy^2+2*d_3t;
B(2,1)=-e_t*(((Damp/(Frcy*Dt)+1)*s/DamFrcy)+(1/(Frcy^2*Dt))*c)+1/(Frcy^2*Dt);B(2,2)=e_t*((Damp/(Frcy*Dt)*s/DamFrcy)+(1/(Frcy^2*Dt))*c)-1/(Frcy^2*Dt);fori=1:(count-1)%根据地震记录,计算不同的反应
结构动力学MATLAB编程例题,华工结构动力学编程 MATLAB相关推荐
- 高中python编程例题-青少年Python实例编程入门
在大数据.云计算的时代,Python程序设计语言具有先天的优势,目前高中信息技术新课程标准也选定Python做为计算思维训练的语言.提前掌握一些Python的基础概念与算法的知识,对老师.学生有重大的 ...
- AI:一个20年程序猿的学习资料大全—结构分析软件/办公软件/电气制造控制/高级语言编程/平面三维设计/视频编辑/FQ格式转换软件——只有你不想要的,没有你找不到的
AI:一个20年程序猿的学习资料大全-结构分析软件/办公软件/电气制造控制/高级语言编程/平面三维设计/视频编辑/FQ格式转换软件--只有你不想要的,没有你找不到的 目录 (有偿提供,替朋友转载,扫描 ...
- Go 语言编程 — 高级数据类型 — 结构体
目录 文章目录 目录 结构体 访问结构体成员 向函数传递结构体 结构体指针 结构体标签(Struct Tag) 结构体 Golang 中,结构体是由一系列具有相同类型或不同类型的数据构成的数据集合.与 ...
- C 语言编程 — 高级数据类型 — 结构体与位域
目录 文章目录 目录 前文列表 结构体 定义结构体 初始化结构体变量 访问结构体成员 结构体的内存分布 将结构体作为实参传入函数 指向结构体变量的指针 位域 定义位域 使用位域结构体的成员 前文列表 ...
- 【MATLAB】流程控制 ( 循环结构 | for 循环 | while 循环 | 分支结构 | if end 分支结构 | if else end 分支结构 | switch case 分支结构 )
文章目录 一.MATLAB 流程控制结构 二.for 循环 1.for 循环 2.嵌套 for 循环 三.while 循环 二.分支结构 1.if end 分支结构 2.if else end 分支结 ...
- labview实例_手把手以实例教你学LabVIEW编程,事件结构框架
LabVIEW又称为G语言,简单易学.形象直观,采用图形化的编程方式,是专为测试.测量和控制应用而设计的系统工程软件. 前面几篇文章中依次讲解了LabVIEW中顺序顺序结构.条件结构及循环结构的编程方 ...
- 阿里巴巴对Java编程【应用结构】的规约
转载自 阿里巴巴对Java编程[应用结构]的规约 应用分层 1. [推荐]图中默认上层依赖于下层,箭头关系表示可直接依赖,如:开放接口层可以依赖于Web 层,也可以直接依赖于 Service 层,依此 ...
- matlab s函数_matlab结构体 rmfield,arrayfun,structfun,struct2cell,cell2struct
觉得有用的话,请点击右下角 推荐给更多小伙伴 neoken_xuAsurada2015Evacloud 本文为matlab自学笔记的一部分,之所以学习matlab是因为其真的是人工智能无论是神经网络还 ...
- [matlab]使用struct创建结构体
[matlab]使用struct创建结构体 语法: struct(域名1,属性值1,域名2,属性值2,.....) 我们具体看看结构体怎么使用: >> person(1)=struct(' ...
最新文章
- c语言设计通讯录管理程序,C语言程序的设计学生通讯录管理系统方案.docx
- 用户体验很棒的Git学习网站
- 提前防止Non-PIE错误,检测app是否包含PIE标志
- 电脑今天关机,第二天重启 解决方案
- html面试要带电脑吗,前端面试巧妙回答浏览器兼容问题
- QTP教程03 - 手工参数化
- UnityShader 浮雕凹凸贴图BumpMap与法线贴图NormalMap的原理及其区别
- Oracle 中文转拼音
- echart 三维可视化地图_Echarts三维坐标系
- 通过 blast 结果查看 测序数据fastq是否被污染,以及污染reads所属物种、所占比例
- 计算机组成原理实验内存读数,计算机组成原理实验八内存系统实验(3页)-原创力文档...
- 丝滑!CVPR 2021 视频插帧论文+开源代码汇总
- chrome总是提示\喔唷,崩溃啦\的解决办法
- JS验证电话和传真号码格式
- 〖Python全栈白宝书-免费版⑲〗- 字符串的编码格式
- 计算机网络基本知识(一):分类
- grid布局浏览器兼容_grid布局适配IE
- 35岁以后找工作难,只能熬到65岁才能退休,太惨 了!
- 【翻译】压缩Windows XP中的ServicePackFiles/i386文件夹
- 手写基于Spring Cloud的TCC分布式事务框架