资源描述：
《振动控制 主动控制算法简介》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、一、主动控制简介1．概念：结构主动控制需要实时测量结构反应或环境干扰，采用现代控制理论的主动控制算法在精确的结构模型基础上运算和决策最优控制力，最后作动器在很大的外部能量输入下实现最优控制力。2．特点：主动控制需要实时测量结构反应或环境干扰，是一种需要额外能量的控制技术，它与被动控制的根本区别是有无额外能量的消耗。3．优缺点：主动控制具有提高建筑物的抵抗不确定性地面运动，减少输入的干扰力，以及在地震时候自动地调整结构动力特征等能力，特别是在处理结构的风振反应具有良好的控制效果，与被动控制相比，主动控制具有更好的控制效果。但是，主动控制实际应用价格昂贵，在实际应用过
2、程中也会存与其它控制理论相同的问题，控制技术复杂、造价昂贵、维护要求高。4．组成：传感器、控制器、作动器5．工作方式：开环、闭环、开闭环。二、简单回顾主动控制的应用与MATLAB应用1.主动变刚度ＡＶＳ控制装置工作原理：首先将结构的反应反馈至控制器，控制器按照事先设定好的控制算法并结合结构的响应，判断装置的刚度状态，然后将控制信号发送至电液伺服阀以操纵其开关状态，实现不同的变刚度状态。锁定状态（ON）：电液伺服阀阀门关闭，双出杆活塞与液压缸之间没有相对位移，斜撑的相对变形与结构层变形相同，此时结构附加一个刚度；打开状态（ＯＦＦ）：电液伺服阀阀门打开，双出杆活塞与液
3、压缸之间有相对位移，液压缸的压力差使得液体发生流动，此过程中产生粘滞阻尼，此时结构附加一个阻尼。示意图如下：2.主动变阻尼ＡＶＤ控制装置工作原理：变孔径阻尼器以传统的液压流体阻尼器为基础，利用控制阀的开孔率调整粘性油对活塞的运动阻力，并将这种阻力通过活塞传递给结构，从而实现为结构提供阻尼的目的。关闭状态（ＯＮ）：开孔率一定，液体的流动速度受限，流动速度越小，产生的粘滞阻尼力越大，开孔率最小时，提供最大阻尼力，此时成为ＯＮ状态；打开状态（ＯＦＦ）：控制阀完全打开，由于液体的粘滞性可提供最小阻尼力。示意图如下：133．振动实例已知多自由度有阻尼线性结构的参数：，，阻尼
4、矩阵采用瑞利阻尼，根据前两阶自振频率及阻尼比确定，阻尼比取0.05，该多自由度结构（参数同上）所受地震波数据见dzb.xls文件，文件第一列为时间，单位s，文件第2列为加速度，单位m/s2。方法采用中心差分法。3.1变刚度对比了刚度分别为K、10*K以及0.1*K时M1的响应时程曲线以及最大位移。ＭＡＴＬＡＢ程序如下：clearclcM=diag([276227602300]);%质量矩阵K=100000*[4.406-1.9210;-1.9213.443-1.522;0-1.5221.522];kk={K,10.*K,0.1.*K}%细胞矩阵-变刚度W=[4.1
5、041;10.4906;14.9514];%各阶频率zuni=0.05area=2*W(1)*W(2)*zuni/(W(1)+W(2));byta=2*zuni/(W(1)+W(2));C=area*M+byta*K;%阻尼矩阵num=xlsread('dzb.xls',1,'B1:B1501');P=M*ones(3,1)*num';%读入外荷载*********中心差分法**********h=0.02;%步长para=[1/h^2,1/(2*h),2/h^2,h^2/2];%参数向量Kx=para(1)*M+C*para(2);%x(i+1)前系数x(:,1
6、)=zeros(3,1);%初位移v(:,1)=zeros(3,1);%初速度a(:,1)=-0.00082*num(1)*ones(3,1);%初加速度forj=1:3fori=1:1:1501%差分迭代第一步ifi<2;x0=x(:,1)-h*v(:,1)+h^2/2*a(:,1);Px(:,i)=P(:,i)-(kk{j}-para(3)*M)*x(:,i)-(para(1)*M-para(2)*C)*x0;13x(:,i+1)=inv(Kx)*Px(:,i);a(:,i+1)=para(1)*(x0-2*x(:,i)+x(:,i+1));%加速度响应v(:
7、,1)=para(2)*(x(:,i+1)-x0);%速度响应else%差分迭代Px(:,i)=P(:,i)-(kk{j}-para(3)*M)*x(:,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(
8、abs(X
