使用ode45对特定节点的时间依赖强制来确定响应

1次观看(过去30天)

显示旧的注释

Shubham 2021年9月22日

0
链接

直接链接到这个问题

https://ch.mathworks.com/matlabcentral/answers/1458429-to-determine-response-using-ode45-for-a-time-dependent-forcing-at-specific-nodes

评论道: 1月 2021年9月22日

                          函数[qdot] = deriv(t,q,A,B)
                         
                          F(9,1) = 100*sin(100*t);共有7个节点(每个节点4个自由度)。我想用这两个力
                         
                          F(10,1) = 100*cos(100*t);%在x和y方向上的第3个节点
                         
                          qdot = A*q +B*F;%状态空间eqn - A - 56x56 | q - 56x1 | B - 56x28 | F - 28x1(除第9 /10行外均为零)
                         
                          如何写出这个随时间变化的力向量
                         
                          %写的不正确。在简单
                         
                          %系统，我们用直接写
                         
                          %力eqn。但是在这里
                         
                          组装系统我不知道怎么做
                         
                          按照我所写的去做
                         
                          %没有显示任何错误，但它是
                         
                          不正确，因为我越来越奇怪了
                         
                          %的回应。
                         
                          所有矩阵A和B都是预先定义的常量。力向量是a
                         
                          时间的%函数(28x1)，第9行和第10行只有两个值。如何把它合并到状态空间eqn中?
                         
                          选项= odesset;
                         
                          (时间,q) =数值(@deriv Time0, 0 (2 * ndof, 1),选择,A, B);
                         
                          时间=时间';
                         
                          response = C*q';
                         
                          返回