<

Question: Ploting matrix elements and assigning/ solving matrices


 

restart

with(LinearAlgebra)

InvT := Matrix([[c^2, s^2, -2*s*c], [s^2, c^2, 2*s*c], [s*c, -s*c, c^2-s^2]])

Matrix(%id = 18446746411704779590)

(1)

T := Matrix([[c^2, s^2, 2*s*c], [s^2, c^2, -2*s*c], [-s*c, s*c, c^2-s^2]])

Matrix(%id = 18446746411704773678)

(2)

c := cos(p)

cos(p)

(3)

s = sin(p)

s = sin(p)

(4)

Q := Matrix([[Q11, Q12, 0], [Q12, Q22, 0], [0, 0, Q66]])

Matrix(%id = 18446746411704759470)

(5)

Q11 := E1/(-v12*v21+1); Q12 := E2/(-v12*v21+1); Q66 := G12

E1/(-v12*v21+1)

 

E2/(-v12*v21+1)

 

G12

(6)

E1 := 0.233e12; E2 := 0.231e11; v21 := 0.2e-1; v12 := .2; G12 := 0.717e10

0.233e12

 

0.231e11

 

0.2e-1

 

.2

 

0.717e10

(7)

R := Matrix([[1, 0, 0], [0, 1, 0], [0, 0, 2]])

Matrix(%id = 18446746411704747062)

(8)

Qbar := Matrix([[Qb11, Qb12, Qb16], [Qb12, Qb22, Qb26], [Qb16, Qb26, Qb66]])

Matrix(%id = 18446746411704742726)

(9)

InvR := MatrixInverse(R)

Matrix(%id = 18446746411704745726)

(10)

eq1 := Qbar = InvT.Q.R.T.InvR

Matrix(%id = 18446746411704742726) = Matrix(%id = 18446746411704732718)

(11)

plot(Qb11, p = 0 .. 9)

Warning, expecting only range variable p in expression Qb11 to be plotted but found name Qb11

 

 

NULL


 

Download plotting_elements_of_matrix.mw

Please Wait...