Add files via upload

ECheynet · Mar 27, 2023 · 2e52df8 · 2e52df8
1 parent f63d090
commit 2e52df8
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/dynaRespTD.m b/dynaRespTD.m
@@ -183,7 +183,7 @@
         A = (1./aDT2.*M+beta/aDT*C+K);
         B1 = F+M.*(1./aDT2*x0+1./aDT*dx0+(1/(2*alphaCoeff)-1)*ddx0);
         B2 = C.*(beta/aDT*x0+(beta/alphaCoeff-1).*dx0);
-        B3 = (beta/alphaCoeff-2)*dt/2*ddx0;
+        B3 = C.*((beta/alphaCoeff-2)*dt/2*ddx0);
 
         x1 = A\(B1+B2+B3);
         ddx1= 1/aDT2.*(x1-x0)-1/aDT.*dx0-(1/(2*alphaCoeff)-1).*ddx0;

diff --git a/dynaRespTD_noCoupling.m b/dynaRespTD_noCoupling.m
@@ -159,7 +159,7 @@
         A = (1./aDT2.*M+beta/aDT*C+K);
         B1 = F+M.*(1./aDT2*x0+1./aDT*dx0+(1/(2*alphaCoeff)-1)*ddx0);
         B2 = C.*(beta/aDT*x0+(beta/alphaCoeff-1).*dx0);
-        B3 = (beta/alphaCoeff-2)*dt/2*ddx0;
+        B3 = C.*(beta/alphaCoeff-2)*dt/2*ddx0;
 
         x1 = A\(B1+B2+B3);
         ddx1= 1/aDT2.*(x1-x0)-1/aDT.*dx0-(1/(2*alphaCoeff)-1).*ddx0;