基于MATLAB梁模型振型动画程序设计

最新推荐文章于 2021-04-11 04:19:30 发布

原创

最新推荐文章于 2021-04-11 04:19:30 发布 · 6k 阅读

42 ·

CC 4.0 BY-SA版权

文章标签：

#MATLAB梁模型编程 #振型动画

基于MATLAB梁模型振型动画程序设计

梁模型振型动画

%该模型为一端固定梁模型振型动画
clear; close all; clc

% 系统参数
E = 1e7;
A = 1.5;
rho = 2.6e-4;%密度

% 节点坐标矩阵第一行为x坐标,对应列第二行为y坐标.既有坐标点(0 0)，(0 40)，(40 0)，(40 40)，(80 0)，(80 40)
p = [0 0 40 40 80 80;
0 40 0 40 0 40];
numberOfNodes = size(p, 2);% 取出系统节点个数即p.

% 对各节点进行建模如第一行表示：连接1 3节点形成一个单元，第二行1 4表示连接1 4节点形成一个单元
t = [1 3;
1 4;
2 4;
3 4;
3 5;
4 5;
4 6;
5 6];
numberOfElements = size(t, 1);% 取出系统总的单元个数即t.

% c
c = A * E;

% Initialization of K and F
K = zeros(2 * numberOfNodes);% 初始刚度阵
M = zeros(2 * numberOfNodes);
F = zeros(2 * numberOfNodes, 1);

% 以下for循环为计算整体刚度阵，质量阵为计算固有频率做准备，经过试验下列for循环可以使用自己写成的有限元整体刚度
%阵，质量阵来代替。但要注意
%在amp = a(:, mode)’ * exp(1i * omega(mode, mode) * timeStep); amp = reshape(amp, 2, numberOfNodes);中输入
%numberOfNodes要与自己实际模型一致
%
for e = 1 : numberOfElements
nodes = t(e, ?;
dofs = reshape([2 * nod