!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! !! !! Test ForceToStress Solver !! !! !! !! !! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! check keywords warn echo on ! Define the parameter in MPa - a - m $yearinsec = 365.25*24*60*60 $rhoi = 900.0/(1.0e6*yearinsec^2) $rhow = 1000.0/(1.0e6*yearinsec^2) $gravity = -9.81*yearinsec^2 ! Prefactor from Cuffey and Paterson (2010) in MPa^{-3} a^{-1} $A1 = 2.89165e-13*yearinsec*1.0e18 $A2 = 2.42736e-02*yearinsec*1.0e18 $Q1 = 60.0e3 $Q2 = 115.0e3 Header Mesh DB "." "cube" End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Simulation Coordinate System = Cartesian 3D Simulation Type = Steady State Steady State Min Iterations = 1 Steady State Max Iterations = 1 Output File = "simple.result" Post File = "simple.vtu" max output level = 4 End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Body 1 !As body 1 does not contains Equation, ! there is a Warning in screen output Initial Condition = 1 End Body 2 Equation = 1 Material = 1 Initial Condition = 1 End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Initial Condition 1 StressAna = Variable Coordinate 1, Coordinate 2 Real MATC "100.0*(cos(6.0*pi*tx(0)))*(cos(6.0*pi*tx(1)))" End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Material 1 Density = Real $rhoi Viscosity Model = String "glen" Viscosity = 1.0 ! Dummy but avoid warning output Glen Exponent = Real 3.0 Limit Temperature = Real -10.0 Rate Factor 1 = Real $A1 Rate Factor 2 = Real $A2 Activation Energy 1 = Real $Q1 Activation Energy 2 = Real $Q2 Glen Enhancement Factor = Real 1.0 Critical Shear Rate = Real 1.0e-10 Constant Temperature = Real -1.0 Cauchy = Logical True End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !Solver 1 ! Equation = "ForceToStress" ! Procedure = "ElmerIceSolvers" "ForceToStress" ! Variable = String "Stress" ! Variable DOFs = 1 ! ! Variable = Fw[Fx:1 Fy:1 Fz:1] ! Variable DOFs = 3 ! ! Force Variable Name = String "Fz" ! ! Exported Variable 1 = StressAna ! Exported Variable 1 DOFS = Integer 1 ! ! Linear System Solver = Direct ! Linear System Direct Method = umfpack !End Solver 1 Equation = "Navier-Stokes" Stabilization Method = String Stabilized Flow Model = Stokes Linear System Solver = Iterative Linear System Iterative Method = BICGStab Linear System Preconditioning = ILU0 Linear System Convergence Tolerance = Real 1.0e-3 Linear System Max Iterations = 50 Linear System Residual Output = 500 Linear System Abort Not Converged = True Nonlinear System Max Iterations = 50 Nonlinear System Convergence Tolerance = 1.0e-3 Nonlinear System Newton After Iterations = 5 Nonlinear System Newton After Tolerance = 1.0e-02 Nonlinear System Relaxation Factor = 1.00 Steady State Convergence Tolerance = Real 1.0e-3 End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Equation 1 Active Solvers(1) = 1 End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Boundary Condition 1 Target Boundaries = 5 Body Id = 2 External Pressure = 20 End