Wolfram Function Repository
Instant-use add-on functions for the Wolfram Language
Function Repository Resource:
Transform components of tensors with arbitrary rank with regard to their transformation behavior under any given mapping
ResourceFunction["TensorCoordinateTransform"][tensor,mat] transform tensor components according to their transformation behavior matrix mat. All tensor slots are considered to be contravariant and not normalized. |
| "TransformationBehavior" | "AllContravariant" | transformation behavior of individual components |
| "Normalize" | False | whether the components are represented with respect to a normalized basis |
Apply a rotation transformation to a tensor:
| In[1]:= | ![]() |
| Out[1]= |
Rotate vector components about the three axes counterclockwise:
| In[2]:= | ![]() |
| Out[2]= |
Change the axis:
| In[3]:= | ![]() |
| Out[3]= |
Rotate the rank-2 tensor about the three axes counterclockwise:
| In[4]:= | ![]() |
| Out[4]= | ![]() |
Rotate a rank-4 tensor with three symmetries about the three axes counterclockwise:
| In[5]:= | ![]() |
| Out[6]= | ![]() |
The contravariant vector components with respect to covariant bases are a=vjbj=vi'bi'. The old system is Cartesian (aj=vj):
| In[7]:= |
Consider a transformation between new ei' and old ej covariant base vectors e1'=2e1+3e2+1e3, e2'=1e1+2e2+2e3 and e3'=1e1+2e3. Here are the contravariant components vi' of the new system:
| In[8]:= | ![]() |
| Out[8]= |
They should be the same vector:
| In[9]:= |
| In[10]:= |
| Out[10]= |
| In[11]:= |
| Out[11]= |
The covariant vector components with respect to contravariant bases are a=vjbj=vi'bj'
| In[12]:= | ![]() |
| Out[12]= |
Find related contravariant base vectors bi' by transforming the covariant base in contravariant sense (index juggling bi'→bi'):
| In[13]:= | ![]() |
| Out[13]= |
| In[14]:= |
| Out[14]= |
| In[15]:= |
| Out[15]= |
Note that covariant and contravariant base vectors are inverse to each other bi'bj'=δi'j':
| In[16]:= |
| Out[16]= |
Consider the aibi≡cijklskltjbi=ci'j'k'l'sk'l'tj'bi'≡ai'bi'-invariant tensors in Cartesian frame:
| In[17]:= | ![]() |
The transformation from Cartesian to local torus coordinates (i→i'):
| In[18]:= | ![]() |
The torus surface at α=β=0…2π, with R and r constant:
| In[19]:= |
| Out[19]= | ![]() |
Transform tensors' different ranks regarding the local torus base with respect to transformation behavior. This is faster than transformation of an entire assembled rank-7 tensor):
| In[20]:= | ![]() |
Local covariant base vector bi':
| In[21]:= |
Assemble the tensor term by contracting the right indices ci'j'k'l'sk'l'tj'bi≡ai'bi':
| In[22]:= |
| Out[22]= |
| In[23]:= |
| Out[23]= |
The result should be the same in the Cartesian system:
| In[24]:= |
| Out[24]= |
If the base vectors are not normalized ("Normalize"→False), i.e. if they have local-dependent lengths (and are not orthogonal in general reference systems either), the corresponding tensor components are transformed according to their co- and contravariant transformation behaviors. This leads to different representations of the same tensor object.
Here is an example of all four transforming possibilities of a rank-2 tensor to a local, non-normalized cylindrical system:
| In[25]:= | ![]() |
All coordinates transform the contravariant tij (default):
| In[26]:= | ![]() |
| Out[26]= | ![]() |
All coordinates transform the covariant tij:
| In[27]:= | ![]() |
| Out[27]= | ![]() |
A tensor object with coordinates of mixed-transformation behavior tij:
| In[28]:= | ![]() |
| Out[28]= | ![]() |
A tensor object with coordinates of opposite mixed-transformation behavior tji:
| In[29]:= | ![]() |
| Out[29]= | ![]() |
Nomenclature referencing all indices together or each index separately is interchangeable:
| In[30]:= | ![]() |
| Out[30]= |
This is analogous for the covariant transformation:
| In[31]:= | ![]() |
| Out[31]= |
Transformed contravariant vector components with respect to normalized base vectors:
| In[32]:= | ![]() |
| Out[32]= |
The same transformed contravariant vector components with respect to non-normalized base vectors (default):
| In[33]:= |
| Out[33]= |
Transformed covariant normalized base vectors:
| In[34]:= | ![]() |
| Out[34]= |
The contraction of the assigned tensors should be produce invariance:
| In[35]:= |
| Out[35]= |
Tensor components with respect to an orthonormal cylindrical system:
| In[36]:= | ![]() |
| Out[35]= | ![]() |
TransformedField also assumes a normalized base:
| In[37]:= | ![]() |
| Out[37]= | ![]() |
| In[38]:= | ![]() |
| Out[38]= |
With orthonormal reference systems (e.g. normalized cylindrical systems), no distinction between co- and contravariant transformation behavior is required:
| In[39]:= | ![]() |
| Out[39]= |
| In[40]:= |
Hook's general law describes a linear relationship between the components of the rank-2 stress tensor σ and the two-stage strain tensor ε using the rank-4 tensor C: σij=Cijklεkl with the symmetries σij=σji, εkl=εlk and Eijkl=Ejikl=Eijlk=Eklij.
Components of the general, fully anisotropic stiffness tensor with initial symmetries:
| In[41]:= | ![]() |
| Out[37]= | ![]() |
The number of independent components:
| In[42]:= |
| Out[42]= |
One plane of symmetry, rotated by 180 degrees, results in the same stiffness:
| In[43]:= |
| Out[44]= |
This results in a stiffness tensor with fewer independent components:
| In[45]:= | ![]() |
| Out[20]= | ![]() |
Determination and comparison of different material symmetries:
| In[46]:= | ![]() |
| In[47]:= | ![]() |
| In[48]:= | ![]() |
Summary:
| In[49]:= | ![]() |
| Out[49]= | ![]() |
This work is licensed under a Creative Commons Attribution 4.0 International License