文档介绍：一、主动控制简介
.概念：结构主动控制需要实时测量结构反应或环境干扰，采用现代控制理论的主动控制算法在精确的结构模型基础上运算和决策最优控制力，最后作动器在很大的外部能量输
入下实现最优控制力。
.特点：主动控制需要实时测量结构反应或(:,i)-(para(1)*M-para(2)*C)*x(:,i-1);x(:,i+1)=inv(Kx)*Px(:,i);
a(:,i+1)=para(1)*(x(:,i-1)-2*x(:,i)+x(:,i+1));%加速度响应
v(:,i)=para(2)*(x(:,i+1)-x(:,i-1));%速度响应
endend*************中心差分法*************
X=x(:,1:1501);Y=max(abs(X),[],2);Z(j)=max(Y);saveX%保存位移相应
subplot(3,1,j)%画图
plot(X(1,:))xlabel('时间t/')
ylabel('位移X1/m');
end运行结果如下：
4006008G01000120014001600
时间t/002s
-■.』o2OOJJEfx漆归
-3o
4006008001000120014001600
时间打002s
onu
2
4006008001000120014001600
时间VO02s
最大位移分别为：
.5C的情况下M1的响应时程曲线
依旧使用上述系统，对比了无阻尼，阻尼为C和。和最大位移。
MATLAB程序:
clearclc
M=diag([276227602300]);%质量矩阵
K=100000*[-;--;0-];%刚度矩阵
W=[;;];%各阶频率
zuni==2*W(1)*W(2)*zuni/(W(1)+W(2));byta=2*zuni/(W(1)+W(2));
C=area*M+byta*K;
%变阻尼
%读入外荷载
%步长
%参数向量
%x(i+1)前系数
%初位移
%初速度
%初加速度
%差分迭代第一步
cc={0*C,C,*C};
num=xlsread('',1,'B1:B1501');P=M*ones(3,1)*num';
**************中心差分法************
h=;
para=[1/hA2,1/(2*h),2/hA2,hA2/2];
Kx=para(1)*M+C*para(2);
x(:,1)=zeros(3,1);
v(:,1)=zeros(3,1);
a(:,1)=-*num(1)*ones(3,1);
forj=1:3
fori=1:1:1501
ifi<2;
x0=x(:,1)-h*v(:,1)+hA2/2*a(:,1);
Px(:,i)=P(:,i)-(K-para(3)*M)*x(:,i)-(para(1)*M-para(2)*cc{j})*x0;
x(:,i+1)=inv(Kx)*Px(:,i);
a(:,i+1)=para(1)*(x0-2*x(:,i)+x(:,i+1));%加速度响应
v(:,1)=para(2)*(x(:,i+1)-x0);%速度响应
else%差分迭代
Px(:,i)=P(:,i)-(K-para(3)*M)*x(:,i)-(para(1)*M-para(2)*cc{j})*x(:,i-1);
x(:,i+1)=inv(Kx)*Px(:,i);
a(:,i+1)=para(1)*(x(:,i-1)-2*x(:,i)+x(:,i+1));%加速度响应
v(:,i)=para(2)*(x(:,i+1)-x(:,i-1));%速度响应
end
end
**************
中心差分法
******************
X=x(:,1:1501);
Y=max(abs(X),[],2);
Z(j)=max(Y);
saveXsubplot(3,1,j)
%保存位移相应
%画图
plot(X(1,:))
xlabel('时间t/')
ylabel('位移X1/m');
end
运行结果是:
O
1-nWno
O.-O
4006008001000120014001600
时间tJQ02s
4006008001000120014001600
时间"0