1.t=0 0=O)=0.,0)=6(0)=0,5(0)=00=0PO)=0/计算步骤 已知t时刻的状态向量及 f(0)+fD(0)+f(0)=P(O)j(0)=0 △w(t-1)、△(t1-1) y(O)0)k(O)=60,k(0)=60+947410074 求t时刻的状态向量及增量 j(O}={0}△P(O)=2△P()=25 (1)求t时刻的状态向量 (1)=y(t1)+△y(t ()(04y(O)=△P(O)/k(0)=000248 y(t1)=y(t1-1)+△y(t1-1) 2.t=0.1s △j(0)=0.0744 y()=P()-f0()-f(c y(0.1)=y(0)+△yO)=000248y(0.1)=y(0)+△(0)=0.074 (2)求AP(1)k(1) y(0.1)<0.05弹性阶段k(0.1)=60 k(t)=k()+2m+c() f(0.1)=k(0.1)×y(0.1)=0.1488 )=△P(1)+my(1)+3(t) f(0.1)=cxy(0.1)=00744P(0.1)=2.5 f(0.)=P()-f(0.1)-f(0.1)=22768,y(0.1)=14891 c()[3j(t)+xj(t) P(t(kN) (3)求A4)4(4) f 3kN 积分步长:M≤T/10 △P()=△P(1)+9474(t)+4.637y() 0.1 y(m △()=30△y()-3j(t)-0.05j(t)P(t)(kN) 0.1 0.8 t(s) 2.5 4 3.5 2.5 1.5 1 0.5 fs y(m) 3kN 0.05 ( ) ( ) 94.74 ( ) 4.637 ( )] ~ P t = P t + y  t +  y  t y (t) = 30y(t) −3y (t) − 0.05 y (t) 计算步骤 已知ti-1时刻的状态向量及 求ti时刻的状态向量及增量 ( ) ( )  i−1  i−1 y t 、 y  t （1）求ti时刻的状态向量 [ ( ) ( ) ( )] 1 ( ) ( ) ( ) ( ) ( ) ( ) ( ) 1 1 1 1 i i D i s i i i i i i i P t f t f t m y t y t y t y t y t y t y t = − − = +  = +  − − − −     （2）求 ) ~ ( ) ~ i i P t 、k（t （3）求 ( ) ) i i y t 、y （t 积分步长： t T /10 ( ) 3 ( ) 6 ( ) ( ) ~ 2 c t t m t k t k t  +  = + ( )] 2 ( )[3 ( ) ( ) 3 ( )] 6 ( ) ( ) [ ~ y t t c t y t y t y t t P t P t m      + + +   =  + 1. t=0 y(0) = 0; y(0) = 0; f (0) = k y(0) = 0; f (0) = cy(0) = 0;P(0) = 0 s D   f (0) + f (0) + f (0) = P(0); y(0) = 0 I D s            =           0 0 0 (0) (0) (0) y y y   (0) 60 947.4 1007.4 ~ k(0) = 60;k = + = (0) 2.5 ~ P(0) = 2.5;P = (0) 0.00248 ~ (0)/ ~ y(0) = P k = y (0) = 0.0744 2. t=0.1s y(0.1) = y(0) + y(0) = 0.00248 y(0.1)  0.05 f (0.1) = k(0.1) y(0.1) = 0.1488 s P(0.1) = 2.5 弹性阶段 k(0.1) = 60 f (0.1) = c y(0.1) = 0.0744 D  f (0.1) = P(0.1) − f (0.1) − f (0.1) = 2.2768; y(0.1) =1.4891 I D s  y (0.1) = y (0) + y (0) = 0.0744