使用Matlab 1 d纽马克

32视图(30天)

显示旧的评论

努尔尤 2016年1月21日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/264582-1d-newmark-using-matlab

评论道: Chaudhary P·帕特尔2022年11月29日

你好,

我试图解决动态结构书中例子5.1 Anil K。Chopra 168页。但是,我有一个问题与Matlab代码。以下是我的代码:-

                         clc
                        
                         m = 0.2533;
                        
                         k = 10;
                        
                         c = 0.1592;
                        
                         gam = 1/2;β= 1/6;
                        
                         %初始计算
                        
                         阿宝= 0;
                        
                         情况= 0;
                        
                         udot0 = 0;
                        
                         u2dot1 = 1 / m * (po-k * u0-c * udot0);
                        
                         %时间步
                        
                         ti = 0。;
                        
                         tf = 4。;
                        
                         dt = 0.1;
                        
                         t = ti: dt: tf;
                        
                         元=修复((tf-ti) / dt);
                        
                         kcap = k + gam /(β* dt) * c + m / dt(β* ^ 2);
                        
                         a = m /(β* dt) + gam * c /β;
                        
                         b = 0.5 * m /β+ dt * (0.5 * gam /β- 1)* c;
                        
                         dp =差异(p);
                        
                         为我= 1:nt
                        
                         deltaP = dp(我)+ * udot(我)+ b * u2dot(我);
                        
                         du_i = deltaP / kcap;
                        
                         dudot_i = gam /(β* dt) * du_i——gam /β* udot (i) + dt * (1 - 0.5 * gam /β)* u2dot(我);
                        
                         du2dot_i = 1 /(β* dt ^ 2) * du_i - 1 /(β* dt) * udot (i) - 0.5 /β* u2dot(我);
                        
                         u (i + 1) = du_i + u(我);
                        
                         udot (i + 1) = dudot_i + udot(我);
                        
                         u2dot (i + 1) = du2dot_i + u2dot(我);
                        
                         结束

当我运行编码这个警告出现

                         未定义的函数或变量“p”。
                        
                         错误(第25行)
                        
                         dp =差异(p);

为什么它发生,如何解决这个问题呢?

谢谢你！