I am trying to find the truss member stiffness matrix and I am inputting
k_member _ 1=1/L1*[c1 s1 0 0; 0 0 c1 s1]'*[1 -1; -1 1]*… [c1 s1 0 0; 0 0 c1 s1];
I have defined c1, s1, and L1 previously.
It says "Unexpected matlab expression"
MATLAB: What should be adjusted to make the script run
stiffness matrix
Related Solutions
c1 = 0, c2 = 0, c4 = 0, c5 = -1000, c6 = 0, d3 = -500, s1 = 1/1000, s2 = -1, s4 = -1/s6, s5 = 0, s6 = s6c1 = 0, c2 = 0, c4 = 0, c5 = -1000, c6 = 0, d3 = -500, s1 = 1/1000, s2 = -1, s4 = -1/s6, s5 = 0, s6 = s6c1 = 0, c2 = 0, c4 = 0, c5 = 1000, c6 = 0, d3 = 500, s1 = 1/1000, s2 = 1, s4 = 1/s6, s5 = 0, s6 = s6c1 = 0, c2 = 0, c4 = 0, c5 = 1000, c6 = 0, d3 = 500, s1 = 1/1000, s2 = 1, s4 = 1/s6, s5 = 0, s6 = s6c1 = 250/77, c2 = -1, c4 = -5929000/62500005929, c5 = 0, c6 = 1, d3 = 0, s1 = 1/1000, s2 = 0, s4 = 19250000000/62500005929, s5 = 1, s6 = 0c1 = 250/77, c2 = -1, c4 = -5929000/62500005929, c5 = 0, c6 = 1, d3 = 0, s1 = 1/1000, s2 = 0, s4 = 19250000000/62500005929, s5 = 1, s6 = 0c1 = 250/77, c2 = -1, c4 = 5929000/62500005929, c5 = 0, c6 = -1, d3 = 0, s1 = 1/1000, s2 = 0, s4 = -19250000000/62500005929, s5 = -1, s6 = 0c1 = 250/77, c2 = -1, c4 = 5929000/62500005929, c5 = 0, c6 = -1, d3 = 0, s1 = 1/1000, s2 = 0, s4 = -19250000000/62500005929, s5 = -1, s6 = 0c1 = 250/77, c2 = 0, c4 = 0, c5 = 0, c6 = 0, d3 = -5929/125000*s6, s1 = 0, s2 = 1/s6, s4 = 1, s5 = 77/250, s6 = s6c1 = 250/77, c2 = 0, c4 = 0, c5 = 0, c6 = 0, d3 = -5929/125000*s6, s1 = 0, s2 = 1/s6, s4 = 1, s5 = 77/250, s6 = s6c1 = 250/77, c2 = 0, c4 = 0, c5 = 0, c6 = 0, d3 = 5929/125000*s6, s1 = 0, s2 = -1/s6, s4 = -1, s5 = -77/250, s6 = s6c1 = 250/77, c2 = 0, c4 = 0, c5 = 0, c6 = 0, d3 = 5929/125000*s6, s1 = 0, s2 = -1/s6, s4 = -1, s5 = -77/250, s6 = s6c1 = 250/77, c2 = 1, c4 = -5929000/62500005929, c5 = 0, c6 = 1, d3 = 0, s1 = 1/1000, s2 = 0, s4 = -19250000000/62500005929, s5 = -1, s6 = 0c1 = 250/77, c2 = 1, c4 = -5929000/62500005929, c5 = 0, c6 = 1, d3 = 0, s1 = 1/1000, s2 = 0, s4 = -19250000000/62500005929, s5 = -1, s6 = 0c1 = 250/77, c2 = 1, c4 = 5929000/62500005929, c5 = 0, c6 = -1, d3 = 0, s1 = 1/1000, s2 = 0, s4 = 19250000000/62500005929, s5 = 1, s6 = 0c1 = 250/77, c2 = 1, c4 = 5929000/62500005929, c5 = 0, c6 = -1, d3 = 0, s1 = 1/1000, s2 = 0, s4 = 19250000000/62500005929, s5 = 1, s6 = 0c1 = 125/77-1/77*(-5929000*s1^2+5929*s1+15625)^(1/2), c2 = 0, c4 = 0, c5 = -5929* s1/(5923071*s1^2+250*(-5929000*s1^2+5929*s1+15625)^(1/2)-5929*s1-31250), c6 = 0, d3 = 1/500*(125+(-5929000*s1^2+5929*s1+15625)^(1/2))/s1, s1 = s1, s2 = 1, s4 = -1, s5 = (9625-77*(-5929000*s1^2+5929*s1+15625)^(1/2))/(5923071*s1^2+250*(-5929000*s1^2+5929*s1+15625)^(1/2)-5929*s1-31250), s6 = -1c1 = 125/77-1/77*(-5929000*s1^2+5929*s1+15625)^(1/2), c2 = 0, c4 = 0, c5 = -5929*s1/(5923071*s1^2+250*(-5929000*s1^2+5929*s1+15625)^(1/2)-5929*s1-31250), c6 = 0, d3 = 1/500*(125+(-5929000*s1^2+5929*s1+15625)^(1/2))/s1, s1 = s1, s2 = 1, s4 = 1, s5 = (-9625+77*(-5929000*s1^2+5929*s1+15625)^(1/2))/(5923071*s1^2+250*(-5929000*s1^2+5929*s1+15625)^(1/2)-5929*s1-31250), s6 = 1c1 = 125/77-1/77*(-5929000*s1^2+5929*s1+15625)^(1/2), c2 = 0, c4 = 0, c5 = 5929*s1/(5923071*s1^2+250*(-5929000*s1^2+5929*s1+15625)^(1/2)-5929*s1-31250), c6 = 0, d3 = 1/500*(-125-(-5929000*s1^2+5929*s1+15625)^(1/2))/s1, s1 = s1, s2 = -1, s4 = -1, s5 = (9625-77*(-5929000*s1^2+5929*s1+15625)^(1/2))/(5923071*s1^2+250*(-5929000*s1^2+5929*s1+15625)^(1/2)-5929*s1-31250), s6 = 1c1 = 125/77-1/77*(-5929000*s1^2+5929*s1+15625)^(1/2), c2 = 0, c4 = 0, c5 = 5929*s1/(5923071*s1^2+250*(-5929000*s1^2+5929*s1+15625)^(1/2)-5929*s1-31250), c6 = 0, d3 = 1/500*(-125-(-5929000*s1^2+5929*s1+15625)^(1/2))/s1, s1 = s1, s2 = -1, s4 = 1, s5 = (-9625+77*(-5929000*s1^2+5929*s1+15625)^(1/2))/(5923071*s1^2+250*(-5929000*s1^2+5929*s1+15625)^(1/2)-5929*s1-31250), s6 = -1c1 = 125/77+1/77*(-5929000*s1^2+5929*s1+15625)^(1/2), c2 = 0, c4 = 0, c5 = -5929*s1/(-5923071*s1^2+250*(-5929000*s1^2+5929*s1+15625)^(1/2)+5929*s1+31250), c6 = 0, d3 = 1/500*(-125+(-5929000*s1^2+5929*s1+15625)^(1/2))/s1, s1 = s1, s2 = -1, s4 = -1, s5 = (-9625-77*(-5929000*s1^2+5929*s1+15625)^(1/2))/(-5923071*s1^2+250*(-5929000*s1^2+5929*s1+15625)^(1/2)+5929*s1+31250), s6 = 1c1 = 125/77+1/77*(-5929000*s1^2+5929*s1+15625)^(1/2), c2 = 0, c4 = 0, c5 = -5929*s1/(-5923071*s1^2+250*(-5929000*s1^2+5929*s1+15625)^(1/2)+5929*s1+31250), c6 = 0, d3 = 1/500*(-125+(-5929000*s1^2+5929*s1+15625)^(1/2))/s1, s1 = s1, s2 = -1, s4 = 1, s5 = (9625+77*(-5929000*s1^2+5929*s1+15625)^(1/2))/(-5923071*s1^2+250*(-5929000*s1^2+5929*s1+15625)^(1/2)+5929*s1+31250), s6 = -1c1 = 125/77+1/77*(-5929000*s1^2+5929*s1+15625)^(1/2), c2 = 0, c4 = 0, c5 = 5929*s1/(-5923071*s1^2+250*(-5929000*s1^2+5929*s1+15625)^(1/2)+5929*s1+31250), c6 = 0, d3 = 1/500*(125-(-5929000*s1^2+5929*s1+15625)^(1/2))/s1, s1 = s1, s2 = 1, s4 = -1, s5 = (-9625-77*(-5929000*s1^2+5929*s1+15625)^(1/2))/(-5923071*s1^2+250*(-5929000*s1^2+5929*s1+15625)^(1/2)+5929*s1+31250), s6 = -1c1 = 125/77+1/77*(-5929000*s1^2+5929*s1+15625)^(1/2), c2 = 0, c4 = 0, c5 = 5929*s1/(-5923071*s1^2+250*(-5929000*s1^2+5929*s1+15625)^(1/2)+5929*s1+31250), c6 = 0, d3 = 1/500*(125-(-5929000*s1^2+5929*s1+15625)^(1/2))/s1, s1 = s1, s2 = 1, s4 = 1, s5 = (9625+77*(-5929000*s1^2+5929*s1+15625)^(1/2))/(-5923071*s1^2+250*(-5929000*s1^2+5929*s1+15625)^(1/2)+5929*s1+31250), s6 = 1Hi,
In you code, first you define L as followed:
L=[2*a 2*a a];Then, you redifine it to:
L = L(i); % where i is 1
After that L is a scalar and not a vector. SO Next time MATLAB reached the line
L = L(i); % this time i = 2
Since L isn't a vector of 3 number but just 2*a it doesn't works.
Don't overwrite L and it will be fine
Related Question
- How I can get Px,Py and Pz value
- Is this assignment displaying an error
- I want to find the determinant for 5.5 matrix
- Tried solving the following complex equation for W varying values of k from 0.1 to 100 and plot the positive imaginary part of the solution (W) against k. Got wrong roots using roots command : S.*((W.^2)-(C1+C2).*W.*k+C1.*C2.*(k.^2))+i.*(W-C0.
- How to find correlation between the row values of a matrix
Best Answer