华工结构动力学编程

一.所编程序及编译:

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相关推荐

  1. 高中python编程例题-青少年Python实例编程入门

    在大数据.云计算的时代,Python程序设计语言具有先天的优势,目前高中信息技术新课程标准也选定Python做为计算思维训练的语言.提前掌握一些Python的基础概念与算法的知识,对老师.学生有重大的 ...

  2. AI:一个20年程序猿的学习资料大全—结构分析软件/办公软件/电气制造控制/高级语言编程/平面三维设计/视频编辑/FQ格式转换软件——只有你不想要的,没有你找不到的

    AI:一个20年程序猿的学习资料大全-结构分析软件/办公软件/电气制造控制/高级语言编程/平面三维设计/视频编辑/FQ格式转换软件--只有你不想要的,没有你找不到的 目录 (有偿提供,替朋友转载,扫描 ...

  3. Go 语言编程 — 高级数据类型 — 结构体

    目录 文章目录 目录 结构体 访问结构体成员 向函数传递结构体 结构体指针 结构体标签(Struct Tag) 结构体 Golang 中,结构体是由一系列具有相同类型或不同类型的数据构成的数据集合.与 ...

  4. C 语言编程 — 高级数据类型 — 结构体与位域

    目录 文章目录 目录 前文列表 结构体 定义结构体 初始化结构体变量 访问结构体成员 结构体的内存分布 将结构体作为实参传入函数 指向结构体变量的指针 位域 定义位域 使用位域结构体的成员 前文列表 ...

  5. 【MATLAB】流程控制 ( 循环结构 | for 循环 | while 循环 | 分支结构 | if end 分支结构 | if else end 分支结构 | switch case 分支结构 )

    文章目录 一.MATLAB 流程控制结构 二.for 循环 1.for 循环 2.嵌套 for 循环 三.while 循环 二.分支结构 1.if end 分支结构 2.if else end 分支结 ...

  6. labview实例_手把手以实例教你学LabVIEW编程,事件结构框架

    LabVIEW又称为G语言,简单易学.形象直观,采用图形化的编程方式,是专为测试.测量和控制应用而设计的系统工程软件. 前面几篇文章中依次讲解了LabVIEW中顺序顺序结构.条件结构及循环结构的编程方 ...

  7. 阿里巴巴对Java编程【应用结构】的规约

    转载自 阿里巴巴对Java编程[应用结构]的规约 应用分层 1. [推荐]图中默认上层依赖于下层,箭头关系表示可直接依赖,如:开放接口层可以依赖于Web 层,也可以直接依赖于 Service 层,依此 ...

  8. matlab s函数_matlab结构体 rmfield,arrayfun,structfun,struct2cell,cell2struct

    觉得有用的话,请点击右下角 推荐给更多小伙伴 neoken_xuAsurada2015Evacloud 本文为matlab自学笔记的一部分,之所以学习matlab是因为其真的是人工智能无论是神经网络还 ...

  9. [matlab]使用struct创建结构体

    [matlab]使用struct创建结构体 语法: struct(域名1,属性值1,域名2,属性值2,.....) 我们具体看看结构体怎么使用: >> person(1)=struct(' ...

最新文章

  1. c语言设计通讯录管理程序,C语言程序的设计学生通讯录管理系统方案.docx
  2. 用户体验很棒的Git学习网站
  3. 提前防止Non-PIE错误,检测app是否包含PIE标志
  4. 电脑今天关机,第二天重启 解决方案
  5. html面试要带电脑吗,前端面试巧妙回答浏览器兼容问题
  6. QTP教程03 - 手工参数化
  7. UnityShader 浮雕凹凸贴图BumpMap与法线贴图NormalMap的原理及其区别
  8. Oracle 中文转拼音
  9. echart 三维可视化地图_Echarts三维坐标系
  10. 通过 blast 结果查看 测序数据fastq是否被污染,以及污染reads所属物种、所占比例
  11. 计算机组成原理实验内存读数,计算机组成原理实验八内存系统实验(3页)-原创力文档...
  12. 丝滑!CVPR 2021 视频插帧论文+开源代码汇总
  13. chrome总是提示\喔唷,崩溃啦\的解决办法
  14. JS验证电话和传真号码格式
  15. 〖Python全栈白宝书-免费版⑲〗- 字符串的编码格式
  16. 计算机网络基本知识(一):分类
  17. grid布局浏览器兼容_grid布局适配IE
  18. 35岁以后找工作难,只能熬到65岁才能退休,太惨 了!
  19. 【翻译】压缩Windows XP中的ServicePackFiles/i386文件夹
  20. 手写基于Spring Cloud的TCC分布式事务框架

热门文章

  1. IT界的兄弟姐妹们,大家好。
  2. 织梦插件-站长必备免费织梦插件快速收录
  3. Milking Order(拓扑+二分+优先队列)
  4. 江苏旅游职业学院计算机专业,江苏旅游职业学院有哪些专业 附好的重点专业名单...
  5. uni-app聊天功能输入框删除emoji表情
  6. 深圳荣耀Java后端一面
  7. “金九银十”求职指南:这些城市和行业最热门!
  8. 如何学好C语言,一个成功人士的心得!
  9. ios kb转m_超圣iOS下载
  10. Solaris IPMP两种配置方法