!################################# !#### My MATC Variables ######### !############################### Header CHECK KEYWORDS Warn !Mesh DB "." "cube" !Mesh DB "." "3stabgrob" !Mesh DB "." "bulk_grob_meter_alle_fl" Mesh DB "." "MeshParkett3Staebchen" !Mesh DB "." "MeshParkett3Staebchen22000" !Mesh DB "." "MeshParkett3Staebchen80000" Include Path "MyLibrary" Results Directory "" End Simulation Max Output Level = 3 Coordinate System = Cartesian Coordinate Mapping(3) = 1 2 3 Simulation Type = Transient Steady State Max Iterations = 15 Steady State Min Iterations = 1 Timestepping Method = BDF BDF Order = 1 Output Intervals = 1 Timestep intervals = 120 Timestep Sizes = 10800 !Adaptive Timestepping = True !Adaptive Time Error = 1.0e-7 !Adaptive Min Timestep = 2500 Solver Input File = case3.sif !Post File = case3.ep Nonlinear Update Exported Variables = LOGICAL True End Body 1 Target Bodies(1) = 1 Name = "Boden" Equation = 1 Material = 1 Initial condition = 1 Body Force = 1 End Body 2 Target Bodies(1) = 2 Name = "Leim" Equation = 1 Material = 2 Initial condition = 1 Body Force = 2 End Body 3 Target Bodies(1) = 3 Name = "Deckschicht" Equation = 1 Material = 3 Initial condition = 1 Body Force = 3 End Solver 1 Equation = "Eriksson2006 Equation" Variable = STRING Fields[MoistureContent:1 Temperature:1] Variable DOFS = 2 Procedure = "MyLibrary/Eriksson2006Solve" "Eriksson2006Solve" Exec Solver = Always Steady State Convergence Tolerance = 1.0e-9 Nonlinear System Convergence Tolerance = 1.0e-11 Nonlinear System Convergence Measure = Residual Nonlinear System Max Iterations = 50 Nonlinear System Relaxation Factor = 1 Linear System Solver = Iterativ Linear System Iterative Method = BiCGStabl Linear System Max Iterations = 500 Linear System Convergence Tolerance = 1.0e-8 Linear System Preconditioning = ILU0 Linear System ILUT Tolerance = 1.0e-3 Linear System Abort Not Converged = TRUE Linear System Residual Output = 1 Linear System Precondition Recompute = 1 End !Solver 2 ! Procedure = File "SaveData" "SaveMaterials" ! Equation = STRING "SaveScalar" ! Exec Solver = before saving ! Parameter 1 = String "RH" ! Parameter 2 = String "Specific Heat capacity" ! Parameter 3 = String "dwdH" ! Parameter 4 = String "Eb" ! !Parameter 5 = String "Outer_RH" !END ! !Solver 3 ! Procedure = File "SaveData" "SaveBoundaryValues" ! Equation = String "SaveBoundaryValues" ! Exec Solver = before saving ! Body Force Parameters 1 = String "MoistureSupply" ! Body Force Parameters 2 = String "HeatSource" ! Parameter 1 = String "T_SURF" ! Parameter 2 = String "MC_SURF" ! Parameter 3 = String "HeatFlux" ! Parameter 4 = String "MoistureFlux" !END Solver 2 Procedure = File "ResultOutputSolve" "ResultOutputSolver" Exec Solver = String "before saving" exec interval = 1 Equation = String "ResultOutput" Output File Name = String "post/cubecase3." !or any other output file name of your choice Output Format = String "vtk" Vtk Format = Logical True End Equation 1 Name = "MyCoupledEquation" !Active Solvers(5) = 1 2 3 4 5 Active Solvers(1) = 1 End Material 1 Name = STRING "Pine" Density = 450.0 !try density in kg/m^3 !Specific Heat capacity = 1.8e03 !J/Kg/K Specific Heat capacity = VARIABLE MoistureContent REAL Procedure "MyLibrary/Materialfunktionen" "cDeliiski" !REAL Procedure "MyLibrary/Materialfunktionen" "c_Deliiski" !Conductivity(3,3) = REAL 0.15 0 0 \ ! 0 0.15 0 \ ! 0 0 0.15 Conductivity(3,3) = VARIABLE Temperature REAL Procedure "MyLibrary/Materialfunktionen" "D_T_Pine_IHTP" !Diffusivity(3,3) = REAL 5e-9 0 0\ ! 0 5e-9 0\ ! 0 0 5e-9 Diffusivity(3,3) = VARIABLE MoistureContent REAL Procedure "MyLibrary/Materialfunktionen" "D_w_Eriksson" ! Activation Energy of Bound Water Eb = VARIABLE MoistureContent !REAL 38 !290 !MCJ/mol REAL Procedure "MyLibrary/Materialfunktionen" "Eb" ! Universal gas Constant R = REAL 8.3144621 !J/(Mol K) !relative humidity in Equilibrium with wood moisture content omega !RH = 0.3 RH = VARIABLE MoistureContent REAL Procedure "MyLibrary/Materialfunktionen" "RH_Zurwitz" !delta MoistureContent/ delta relative humidy is part of the sorret effect !dwdH = 0.1 dwdH = VARIABLE MoistureContent REAL Procedure "MyLibrary/Materialfunktionen" "dwdH_Zurwitz" ! ########################################################## ! ############ mechanische Eigenschaften ! ########################################################## Youngs modulus(6,6) = 2500 -33 -600 0 0 0 \ -33 61.6 -27 0 0 0 \ -600 -27 1430 0 0 0 \ 0 0 0 1590 0 0 \ 0 0 0 0 27000 0 \ 0 0 0 0 0 1290 Rotate Elasticity Tensor = Logical True Material Coordinates Unit Vector 1(3) = REAL 0.0\ 1.0\ 0.0 Material Coordinates Unit Vector 2(3) = REAL 0.0\ 0.0\ 1.0 Material Coordinates Unit Vector 3(3) = REAL 1.0\ 0.0\ 0.0 Reference Temperature = REAL 293.15 Heat Expansion Coefficient(6) = REAL 50e-06\ 60e-06\ 20e-06\ 0\ 0\ 0 !per degree Reference Moisture Content = REAL 0.15 moisture Expansion Coefficient(6) =REAL 1.9e-1\ 3.6e-01\ 0.3e-1\ 0\ 0\ 0 End Material 2 Name = STRING "Leim" Density = REAL 360!1e-5 !Heat Conductivity = 3350 Specific Heat capacity = 1.75e03 !J/Kg/K !Specific Heat capacity = VARIABLE MoistureContent !REAL Procedure "MyLibrary/Materialfunktionen" "cDeliiski" Conductivity(3,3) = REAL 0.15 0 0 \ 0 0.15 0 \ 0 0 0.15 !Conductivity(3,3) = VARIABLE Temperature !REAL Procedure "MyLibrary/Materialfunktionen" "D_T_Pine_IHTP" Diffusivity(3,3) = REAL 5e-5 0 0\ 0 5e-5 0\ 0 0 5e-5 !Diffusivity(3,3) = VARIABLE MoistureContent !REAL Procedure "MyLibrary/Materialfunktionen" "D_w_Eriksson" ! Activation Energy of Bound Water Eb = VARIABLE MoistureContent !REAL 38 !290 !MCJ/mol REAL Procedure "MyLibrary/Materialfunktionen" "Eb" ! Universal gas Constant R = REAL 8.3144621 !J/(Mol K) !relative humidity in Equilibrium with wood moisture content omega !RH = 0.3 RH = VARIABLE MoistureContent REAL Procedure "MyLibrary/Materialfunktionen" "RH_Zurwitz" !delta MoistureContent/ delta relative humidy is part of the sorret effect !dwdH = 0.1 dwdH = VARIABLE MoistureContent REAL Procedure "MyLibrary/Materialfunktionen" "dwdH_Zurwitz" !ausgedachtes isotropes Material: E = 3000 v=0.48 G=1700 Youngs modulus(6,6) = 26351 24324 24324 0 0 0\ 24324 26351 24324 0 0 0\ 24324 24324 26351 0 0 0\ 0 0 0 26351 0 0 \ 0 0 0 0 26351 0 \ 0 0 0 0 0 26351 Rotate Elasticity Tensor = Logical false ReferenceTemperature = REAL 293.15 ReferenceMoistureContent = REAL 0.15 ! HeatExpansionCoefficient(6) = 34e-06\ ! 20e-06\ ! 50e-06\ ! 0\ ! 0\ ! 0 !per degree HeatExpansionCoefficient(6) = 67e-09\ 67e-09\ 67e-09\ 0\ 0\ 0 !per degree HygroExpansionCoefficient(6) = 5e-03\ 5e-3 \ 5e-3\ 0\ 0\ 0 End Material 3 Name = STRING "Rotbuche" Moisture Capacity = REAL 0.2 Diffusivity = REAL 5.0e-3 !7.4716e-10 ! valide for Temp. = 297K, see Mamlquist1991 Density = REAL 360!7.2e-7 !Specific Heat capacity = 1.8e03 !J/Kg/K Specific Heat capacity = VARIABLE MoistureContent REAL Procedure "MyLibrary/Materialfunktionen" "cDeliiski" !REAL Procedure "MyLibrary/Materialfunktionen" "c_Deliiski" !Conductivity(3,3) = REAL 0.15 0 0 \ ! 0 0.15 0 \ ! 0 0 0.15 Conductivity(3,3) = VARIABLE Temperature REAL Procedure "MyLibrary/Materialfunktionen" "D_T_Beech_IHTP" !Diffusivity(3,3) = REAL 5e-9 0 0\ ! 0 5e-9 0\ ! 0 0 5e-9 Diffusivity(3,3) = VARIABLE MoistureContent REAL Procedure "MyLibrary/Materialfunktionen" "D_w_Eriksson" ! Activation Energy of Bound Water Eb = VARIABLE MoistureContent !REAL 38 !290 !MCJ/mol REAL Procedure "MyLibrary/Materialfunktionen" "Eb" ! Universal gas Constant R = REAL 8.3144621 !J/(Mol K) !relative humidity in Equilibrium with wood moisture content omega !RH = 0.3 RH = VARIABLE MoistureContent REAL Procedure "MyLibrary/Materialfunktionen" "RH_Zurwitz" !delta MoistureContent/ delta relative humidy is part of the sorret effect !dwdH = 0.1 dwdH = VARIABLE MoistureContent REAL Procedure "MyLibrary/Materialfunktionen" "dwdH_Zurwitz" !Orthotrope Material Fichte from Grimsel1999 ! E_t = 400 E_l = 16200 E_r = 699 Youngs modulus(6,6) = 862 -37 -310 0 0 0 \ -37 71.4 -32 0 0 0 \ -310 -32 428 0 0 0 \ 0 0 0 610 0 0 \ 0 0 0 0 2150 0 \ 0 0 0 0 0 929 Rotate Elasticity Tensor = Logical False ReferenceTemperature = REAL 20.0 ! HeatExpansionCoefficient(6) = 34e-06\ ! 20e-06\ ! 50e-06\ ! 0\ ! 0\ ! 0 !per degree HeatExpansionCoefficient(6) = 67e-06\ 67e-06\ 67e-06\ 0\ 0\ 0 !per degree ReferenceMoistureContent = REAL 0.01 HygroExpansionCoefficient(6) = 3.5e-3 2e-3\ 3.5e-3 \ 0\ 0\ 0 End Initial Condition 1 Name = "InitialCondition 1" Temperature = 283.15 !Temperature in Kelvin MoistureContent = 0.15 End Body Force 1 Name = "Gravity" !BodyForce 3 = -9.81 !MoistureSupply = 0.0000004 !HeatSource = 1 End Body Force 2 Name = "Gravity" !BodyForce 3 = -9.81 !MoistureSupply = 0.0000004 !HeatSource = 2 End Body Force 3 Name = "Gravity" !BodyForce 3 = -9.81 !MoistureSupply = 0.0000004 !HeatSource = 3 End Boundary Condition 10 Target Nodes(1) = 1 !(0,0,0) displacement 1 = 0 displacement 2 = 0 displacement 3 = 0 END Boundary Condition 11 Target Nodes(1) = 2 !(0.1,0,0) displacement 2 = 0 displacement 3 = 0 END Boundary Condition 12 Target Nodes(1) = 3 !(0.1,0.1,0) displacement 3 = 0 END Boundary Condition 1 Target Boundaries(2) = 2 3 Temperature Boundary Type = "cauchy" S_T =REAL 20.0 T_SURF = Variable time REAL 0 283.15 129600 283.15 604800 343.15 1209600 343.15 1296000 283.15 END END Boundary Condition 2 Target Boundaries(1) = 1 Temperature Boundary Type = "cauchy" S_T = REAL 20.0 T_SURF = 283.15 Moisture Boundary Type = "cauchy" S_MC = REAL 2.8e-8 MC_SURF = Variable time REAL Procedure "MyLibrary/Materialfunktionen" "EMC_Zurwitz_Parkett" End !Boundary Condition 1 ! Target Boundaries(3) = 1 2 3 ! ! !! Define outer RH ! !Outer_RH = REAL 0.4 ! !!Outer_RH = Variable time ! !! REAL ! !! 0 0.9 ! !! 86400 0.9 ! !! 864001 0.8 ! !! 345600 0.8 ! !! 432000 0.6 ! !! 1209600 0.4 ! !! 1209601 0.55 ! !! 1296000 0.55 ! !! END ! Moisture Boundary Type = "cauchy" ! S_MC = REAL 2.8e-8 !m/s ! MC_SURF = Variable time !% kg/kg ! REAL Procedure "MyLibrary/Materialfunktionen" "EMC_Zurwitz" ! ! Heating ! ! 0 - 1 0.15 ! ! 1 - 129600s diese dann halten ! ! Drying ! ! 129600 - 1209600 Rampe auf 0.05 ! ! Conditioning ! ! 1209600 - 1209601 kurze rampe auf 0.08 ! ! 1209601 - 1296000 halten auf 0.08 ! Temperature Boundary Type = "cauchy" ! S_T = 20.0 !W/(m^2 K) ! T_SURF = VARIABLE time !313.15 !Kelvin ! REAL ! 0 333.15 ! 129600 333.15 ! 604800 343.15 ! 1209600 343.15 ! 1296000 343.15 ! END ! ! ! Heating ! ! 0 - 1 338.15 sofort auf 60 Grad ! ! 1 - 129600s 338.15 diese dann halten ! ! Drying ! ! 129600 - 604800 Rampe auf 343.15 ! ! 604800 - 1209600 halten auf 343.15 ! ! Conditioning ! ! 1209600 - 1296000 halten auf 343.15 !End