A CAS Approach to Handle the Anisotropic Hooke’s Law for Cancellous Bone and Wood

Algorithm 2

  In[4]∶=MakeName[myexp__]  ∶=

   ToExpression  //@  (  StringJoin  @@  Map[  ToString[  #]  &,  List[  myexp]])

  MakeTensor[  mystring_,  2,  dim_  ]  ∶=

 Apply[(MakeName[mystring,  ##]  &),  SymIndex[2,  dim],  {2}]

MakeTensor[  mystring_,  4  ]  ∶=

 Apply[(MakeName[mystring,  ##]  &),  SymIndex[4],  {4}]

 In[5]∶=(  stress  =  MakeTensor["sigma",  2,  3])  //  MatrixForm

   In[6]∶=(  strain  =  MakeTensor["epsilon",  2,  3])  //  MatrixForm

 In[7]∶=(  C4  =  MakeTensor["C",  4]  )  //  MatrixForm