我怎么能修复”错误信谊/ subsref(第898行)R_tilde =装入的(“subsref”, L_tilde Idx);吗?

16个视图(30天)

显示旧的评论

麦特尔斯 2022年10月2日

0
链接

这个问题直接联系

https://ch.mathworks.com/matlabcentral/answers/1816150——如何-我-修复错误sym - subsref -线- 898 r_tilde装入的subsref - l_tilde idx

回答: 明星黾 2022年10月2日

错误也表示在43行;

错误Assignment2CodeQuadratic2(43行)

KgQ (indiciesQ indiciesQ) = KgQ (indiciesQ indiciesQ) +潜艇(KeQ, [x1Q x2Q x3Q], [xQ(我* 2 - 1)

xQ(我* 2)xQ(我* 2 + 1)]);%克成公斤的贡献

% ENME302-22S2, 2022年9月21日,线性有限元素示范

%来自元素使用线性形状函数方程。

                          清晰,近所有clc;
                         
                          % %二次形状元素:
                         
                          信谊xQ x1Q x2Q x3Q u1Q u2Q u3Q a b c
                         
                          W = 0.01;%宽度[m]
                         
                          tau0 = 1;%初始刚度(Pa)
                         
                          τ= (1 + 0.9 * sin (sqrt(2) *π* xQ) / W)) * tau0;%刚度方程
                         
                          f = 100;%的力密度(N / m ^ 3)
                         
                          uQ =一个* xQ ^ 2 + b * xQ + c;%插值函数(二次)从方程8.16
                         
                          xQ, eqn1Q =潜艇(uQ x1Q) = = u1Q;%评估左节点,方程8.18
                         
                          xQ, eqn2Q =潜艇(uQ x2Q) = = u2Q;%评估左节点,方程8.18 b
                         
                          xQ, eqn3Q =潜艇(uQ x3Q) = = u3Q;%评估左节点,方程8.18摄氏度
                         
                          consQ =解决([eqn1Q eqn2Q eqn3Q], [a b c]);%解联立方程的常数a, b和c,方程8.19
                         
                          uQ =潜艇(uQ [a b c], [consQ。consQ。b consQ.c]);%替代/插入插值函数,方程8.20
                         
                          [NQ ~] =多项式系数(uQ [u1Q u2Q u3Q]);%收集系数定义形状函数,方程8.20 b
                         
                          N1Q = NQ (1);习% N1 = 2 * ^ 2 - 3 * 11 + 1;
                         
                          N2Q = NQ (2);% N2 = 4 *ξ^ 2 + 4 * 11;
                         
                          N3Q = NQ (3);习% N3 = 2 * ^ 2 - 11;
                         
                          %元素方程,柯* Te = Fe(象征性的):
                         
                          KeQ =符号(0 (3));%元素刚度矩阵
                         
                          天文=符号(0 (3,1));%元素rhs迫使向量
                         
                          为i = 1:3% 3自由度的二次元素
                         
                          eleqQ = int(τ* diff (uQ xQ) * diff (NQ (i), xQ), xQ, x1Q, x3Q)——int (f * NQ (i), xQ, x1Q, x3Q);
                         
                          [coefQ ~] =多项式系数(eleqQ [u1Q u2Q u3Q]);%收集系数
                         
                          KeQ(我:)= coefQ (1:3);% u1, u2, u3的系数
                         
                          天文(i) = -coefQ (4);%剩余条款(不是系数u1, u2或u3)
                         
                          结束
                         
                          为kk = 2:4:22
                         
                          * Tg = Fg %全球方程,公斤(数字):
                         
                          L = W;%长度(毫米)
                         
                          一个=乐;%的元素数量
                         
                          景深=一个+ 1;%的自由度
                         
                          x = linspace (0 L景深);% x坐标的景深
                         
                          KgQ =符号(0(景深));%全球刚度矩阵
                         
                          FgQ =符号(0(景深,1));%全球rhs迫使向量
                         
                          为我= 1:一个%组装每一组元素的方程
                         
                          indiciesQ =我* 2 - 1:我* 2 + 1;
                         
                          KgQ (indiciesQ indiciesQ) = KgQ (indiciesQ indiciesQ) +潜艇(KeQ, [x1Q x2Q x3Q], [xQ(我* 2)xQ(我* 2)xQ(我* 2 + 1)));%克成公斤的贡献
                         
                          FgQ (indiciesQ) = FgQ (indiciesQ) +潜艇(天文,[x1Q x2Q x3Q], [xQ(我* 2)xQ(我* 2)xQ(我* 2 + 1)));% Fe Fg的贡献
                         
                          结束
                         
                          %的最终方程与边界条件,K *确认= F(数字):
                         
                          Ta = 0;%对lh狄利克雷边界条件
                         
                          FgQ (1) = FgQ (1) - KgQ(1,1) *助教;
                         
                          KgQ (1, - 1) = 1;
                         
                          FgQ (2) = FgQ (2) - KgQ(2, 1) *助教;
                         
                          KgQ (2, 1) = 0;
                         
                          KgQ (3,1) = 0;
                         
                          结核病= 0;% rhs狄利克雷边界条件
                         
                          FgQ(景深)= FgQ(景深)——KgQ(景深,景深)*结核病;
                         
                          KgQ(景深,景深)= 1;
                         
                          FgQ (dof-1) = FgQ (dof-1)——KgQ (dof-1景深)*结核病;
                         
                          KgQ (dof-1,景深)= 0;
                         
                          KgQ (dof-2,景深)= 0;
                         
                          KgQ =双(KgQ);
                         
                          FgQ =双(FgQ);
                         
                          %解决和情节的解决方案:
                         
                          TbcQ = KgQ \ FgQ;%未知变量,在这种情况下,我们计算dT / dx x1和dT / dx | | x5
                         
                          TgQ =[助教;TbcQ (2: dof-1);结核病);%全局位移矢量值,即取代狄利克雷公元前
                         
                          %在中间找到位移为每个元素的数量
                         
                          loc =找到(x = = 0.005);
                         
                          位移(kk) = Tg (loc) * 10 ^ 3;
                         
                          结束
                         
                             指数超过数组元素的数量。指标不得超过1。

错误索引(第1079行)
R_tilde =内装式(subsref, L_tilde Idx);

                          %位移数组
                         
                          centreDisplacement =位移(2:4:22)%去除零的数组
                         
                          %相对误差
                         
                          relError (kk) = abs (centreDisplacement(结束)-centreDisplacement (kk)) / centreDisplacement (kk)) * 100;
                         
                          情节(x * 10 ^ 3、Tg、“k -”);包含(“x(毫米));ylabel (“偏转(m)”);
                         
                          标题(“线性形状函数替代”)
                         
                          传奇(“线性”);

0评论
显示1年长的评论藏1年长的评论

登录置评。

在回答这个问题。

答案(2)

Torsten 2022年10月2日

0
链接

直接链接到这个答案

https://ch.mathworks.com/matlabcentral/answers/1816150——如何-我-修复错误sym - subsref -线- 898 r_tilde内装式- subsref l_tilde idx # answer_1065210

                              KgQ (indiciesQ indiciesQ) = KgQ (indiciesQ indiciesQ) +潜艇(KeQ, [x1Q x2Q x3Q], [xQ(我* 2)xQ(我* 2)xQ(我* 2 + 1)));%克成公斤的贡献
                             
                              FgQ (indiciesQ) = FgQ (indiciesQ) +潜艇(天文,[x1Q x2Q x3Q], [xQ(我* 2)xQ(我* 2)xQ(我* 2 + 1)));% Fe Fg的贡献