diff --git a/HyperElasticity/Integration_Lib.pro b/HyperElasticity/Integration_Lib.pro deleted file mode 100644 index 99599cf472e523368faa38885b860885bf1a97b9..0000000000000000000000000000000000000000 --- a/HyperElasticity/Integration_Lib.pro +++ /dev/null @@ -1,37 +0,0 @@ -/* - Integration method - GradGrad - CurlCurl -*/ - - -Integration { - { Name GradGrad ; - Case { {Type Gauss ; - Case { - { GeoElement Point ; NumberOfPoints 1 ; } - { GeoElement Line ; NumberOfPoints 4 ; } - { GeoElement Triangle ; NumberOfPoints 4 ; } - { GeoElement Quadrangle ; NumberOfPoints 4 ; } - { GeoElement Tetrahedron ; NumberOfPoints 4 ; } - { GeoElement Hexahedron ; NumberOfPoints 6 ; } - { GeoElement Prism ; NumberOfPoints 9 ; } } - } - } - } - { Name CurlCurl ; - Case { {Type Gauss ; - Case { - { GeoElement Point ; NumberOfPoints 1 ; } - { GeoElement Line ; NumberOfPoints 4 ; } - { GeoElement Triangle ; NumberOfPoints 4 ; } - { GeoElement Quadrangle ; NumberOfPoints 4 ; } - { GeoElement Tetrahedron ; NumberOfPoints 4 ; } - { GeoElement Hexahedron ; NumberOfPoints 6 ; } - { GeoElement Prism ; NumberOfPoints 9 ; } } - } - } - } -} - -/* --------------------------------------------------------------------------*/ diff --git a/HyperElasticity/Jacobian_Lib.pro b/HyperElasticity/Jacobian_Lib.pro deleted file mode 100644 index 71f7b0841d08149e47d76a4baa0918f0dd3387f9..0000000000000000000000000000000000000000 --- a/HyperElasticity/Jacobian_Lib.pro +++ /dev/null @@ -1,59 +0,0 @@ -/* - Jacobian methods - Vol -*/ - -/* I N P U T - --------- - - GlobalGroup : - ----------- - DomainInf Regions with Spherical Shell Transformation - - Parameters : - ---------- - Val_Rint, Val_Rext Inner and outer radius of the Spherical Shell - of DomainInf - -*/ - -/* --------------------------------------------------------------------------*/ - - - - - -Group { - DefineGroup[ DomainInf ] ; - DefineVariable[ Val_Rint, Val_Rext ] ; -} - - - -/* --------------------------------------------------------------------------*/ - - - - -Jacobian { - { Name Vol ; - Case { { Region DomainInf ; - Jacobian VolSphShell {Val_Rint, Val_Rext} ; } - { Region All ; Jacobian Vol ; } - } - } - { Name JSur ; - Case { - { Region All ; Jacobian Sur ; } - } - } - { Name JLin ; - Case { - { Region All ; Jacobian Lin ; } - } - } -} - - - -/* --------------------------------------------------------------------------*/ diff --git a/HyperElasticity/Total_Lagrangian.pro b/HyperElasticity/Total_Lagrangian.pro deleted file mode 100644 index 5f42170368c355c7a4c1391235b6cf3090d0dbe9..0000000000000000000000000000000000000000 --- a/HyperElasticity/Total_Lagrangian.pro +++ /dev/null @@ -1,104 +0,0 @@ -//===================================================================== - -FunctionSpace { - For i In {1:3} - { Name H_u~{i} ; Type Form0 ; - BasisFunction { - { Name sn~{i} ; NameOfCoef un~{i} ; Function BF_Node ; - Support Region[ { Domain_Mecha} ] ; Entity NodesOf[ All ] ; } - } - Constraint { - { NameOfCoef un~{i} ; - EntityType NodesOf ; NameOfConstraint Displacement~{i} ; } - } - } - EndFor -} - -Formulation { - { Name Total_Lagrangian ; Type FemEquation ; - Quantity { - For i In {1:3} - { Name u~{i} ; Type Local ; NameOfSpace H_u~{i} ; } - EndFor - } - Equation { - // The right-hand side - //=================== - For i In {1:3} - Galerkin { [ P~{i}[ (TensorDiag[1., 1., 1.] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] , {d u~{i}} ] ; - In Domain_Mecha; Jacobian Vol; Integration GradGrad; } - For j In {1:3} - Galerkin { [ C_Tot~{j}~{i}[ (TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] * Dof{d u~{i}}, {d u~{j}} ] ; - In Domain_Mecha ; Jacobian Vol ; Integration GradGrad ; } - Galerkin { [-C_Tot~{j}~{i}[ (TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] * {d u~{i}}, {d u~{j}} ] ; - In Domain_Mecha ; Jacobian Vol ; Integration GradGrad ; } - EndFor - EndFor - - Galerkin { [ - CompX[ -force_ext[] ] , {u~{1}} ] ; - In Domain_Mecha; Jacobian Vol; Integration GradGrad; } - Galerkin { [ - CompY[ -force_ext[] ] , {u~{2}} ] ; - In Domain_Mecha; Jacobian Vol; Integration GradGrad; } - Galerkin { [ - CompZ[ -force_ext[] ] , {u~{3}} ] ; - In Domain_Mecha; Jacobian Vol; Integration GradGrad; } - - } - } -} - -PostProcessing { - { Name Total_Lagrangian ; NameOfFormulation Total_Lagrangian ; - NameOfSystem Sys_Mec; - PostQuantity { - { Name u ; Value{ Term{ [ Vector[ {u~{1}}, {u~{2}}, {u~{3}} ]]; - In Domain_Mecha ; Jacobian Vol ; } } } - { Name um ; Value{ Term{ [ Norm[ Vector[ {u~{1}}, {u~{2}}, {u~{3}} ] ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name ux ; Value{ Term{ [ {u~{1}} ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name uy ; Value{ Term{ [ {u~{2}} ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name uz ; Value{ Term{ [ {u~{3}} ] ; - In Domain_Mecha ; Jacobian Vol; } } } - - { Name PK2_xx; Value{ Term { [ CompXX[PK2[ (TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name PK2_xy; Value{ Term { [ CompXY[PK2[ (TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name PK2_xz; Value{ Term { [ CompXZ[PK2[ (TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name PK2_yx; Value{ Term { [ CompYX[PK2[ (TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name PK2_yy; Value{ Term { [ CompYY[PK2[ (TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name PK2_yz; Value{ Term { [ CompYZ[PK2[ (TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name PK2_zx; Value{ Term { [ CompZX[PK2[ (TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name PK2_zy; Value{ Term { [ CompZY[PK2[ (TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name PK2_zz; Value{ Term { [ CompZZ[PK2[ (TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - - { Name F_tot_xx ; Value { Term { [ CompXX[(TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name F_tot_xy ; Value { Term { [ CompXY[(TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name F_tot_xz ; Value { Term { [ CompXZ[(TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name F_tot_yx ; Value { Term { [ CompYX[(TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name F_tot_yy ; Value { Term { [ CompYY[(TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name F_tot_yz ; Value { Term { [ CompYZ[(TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name F_tot_zx ; Value { Term { [ CompZX[(TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name F_tot_zy ; Value { Term { [ CompZY[(TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - { Name F_tot_zz ; Value { Term { [ CompZZ[(TensorDiag[1.0, 1.0, 1.0] + TensorV[{d u~{1}}, {d u~{2}}, {d u~{3}} ] ) ] ] ; - In Domain_Mecha ; Jacobian Vol; } } } - } - } -}