viewFourthOrderTensor

Below is a demonstration of the features of the viewFourthOrderTensor2D function

Contents

clear; close all; clc;

Syntax

viewFourthOrderTensor(C,numDigits,fontSizeIm,fontSize);

Description

This function creates the 9x9 and the 6x6 (aka Voigt) array mappings for 3D 4th order input tensor C. The function can take up to 4 inputs:

Examples

Viewing fourth-order stiffness tensors

Creating the stiffness tensor for Hooke's law of linear elasticity

%Constructing 4th order base tensor set
I=eye(3,3); %The 2nd order identity tensor
II1=dyadicProduct(I,I,1); %4th order base tensor 1
II3=dyadicProduct(I,I,3); %4th order base tensor 3

%Parameters for Hooke's law
mu=1; %The shear modulus
lambda=5; %The lambda lame parameter
C=lambda.*II1+2.*mu.*II3; %Construct 4th order stiffness tensor

Visualizing the tensor using viewFourthOrderTensor

viewFourthOrderTensor(C); %Visualize tensor C

Viewing fourth-order stiffness tensors with symbolic variables

syms mu lambda; %Create symbolic Lame parameters
C=lambda.*II1+2.*mu.*II3; %Construct 4th order stiffness tensor

Visualizing the tensor using viewFourthOrderTensor

viewFourthOrderTensor(C); %Visualize tensor C

GIBBON www.gibboncode.org

Kevin Mattheus Moerman, [email protected]