$dy = 0.3 $mu = 0.1 $fileid = "a" Header CHECK KEYWORDS Warn Mesh DB "." "plates" Include Path "" Results Directory "" End Simulation Max Output Level = 5 Coordinate System = Cartesian 3D Coordinate Mapping(3) = 1 2 3 ! Simulation Type = Scanning Simulation Type = Transient Steady State Max Iterations = 1 Mesh Levels = 1 ! This leads to maximum dynamic friction coefficient of 0.2 ! Convergence is reached until 0.3 Timestep Intervals = 10 Timestep Size = 0.1 ! The ElasticSolver does not really like the Dirichlet conditions at the start ! of the nonlinear iteration. Initialize Dirichlet Conditions = False Simulation Timing = True End Constants Gravity(4) = 0 -1 0 9.82 End Body 1 Name = "upper" Target Bodies(1) = 1 Equation = 1 Material = 1 End Body 2 Name = "lower" Target Bodies(1) = 2 Equation = 1 Material = 2 End Solver 1 Equation = "NonlinElast" Procedure = "ElasticSolve" "ElasticSolver" ! Procedure = "StressSolve" "StressSolver" Variable = -dofs 3 Displacement Nonlinear System Convergence Tolerance = 1.0e-6 Nonlinear System Max Iterations = 30 Nonlinear System Relaxation Factor = 1.0 Linear System Solver = "Iterative" Linear System Preconditioning = Ilu0 Linear System Residual Output = 10 Linear System Max Iterations = 5000 Linear System Iterative Method = GCR !Bicgstab Linear System Convergence Tolerance = 1.0e-5 Calculate Stresses = Logical True Apply Contact BCs = Logical True Linear System Residual Mode = Logical True Nonlinear System Convergence Without Constraints = Logical True ! Eliminate Linear Constraints = Logical True ! Eliminate Slave = Logical True !! unused tricks for testing etc. ! Elasticity Solver Linear = Logical True ! Apply Limiter Conservative Remove After Iterations = Integer 5 ! Apply Limiter Conservative Add After Iterations = Integer 5 ! Optimize Bandwidth = False ! Linear System Scaling = False End Solver 2 Equation = "SaveLine" Exec Solver = never Procedure = "SaveData" "SaveLine" Filename = l$fileid$.dat End Solver 3 Equation = "SaveScalars" Exec Solver = After timestep Procedure = "SaveData" "SaveScalars" Filename = f$fileid$.dat ! Compute the net force in normal and tangential direction Variable 1 = displacement contact normalload Operator 1 = boundary int Variable 2 = displacement contact slipload Operator 2 = boundary int End Solver 4 ! Exec Solver = never Equation = "result output" Procedure = "ResultOutputSolve" "ResultOutputSolver" Output File Name = case_$fileid$ Vtu Format = Logical True Displace Mesh = Logical False Single Precision = Logical False End Equation 1 Name = "Deformation" Active Solvers(1) = 1 End Material 1 Name = "upper" Youngs modulus = 2.0e3 Density = 0.0 Poisson ratio = 0.3 End Material 2 Name = "lower" Youngs modulus = 2.0e3 Density = 0.0 Poisson ratio = 0.3 End Boundary Condition 1 Name = "Contact" Target Boundaries(1) = 8 Mortar BC = 5 Mortar BC Nonlinear = Logical True Mortar BC Initial Contact Depth = Real .1 Slide Contact = Logical True Dynamic Friction Coefficient = Real 0.2 Static Friction Coefficient = Real 0.3 ! If we have a steady state slide contact then this keyword ! can be used to define the relative contact velocity. ! Only the distance of the tangential component is used. ! Contact Velocity(3) = Real 1.0 0.0 0.0 ! Initial Contact Friction = Logical False ! Contact Active Set Minimum = Integer 10 ! Create Dual Projector = Logical True ! Note: if this is activated then the elimination fails for some reason! ! Projector Max Distance = Real $0.7*dy ! Flatten the interfaces to enable 1D projection Plane Projector = Logical True Galerkin Projector = Logical True ! Use N-T conditions which allow friction for generic cases be set Normal-Tangential Displacement = Logical True Mass Consistent Normals = Logical True Save Scalars = True ! Use Biorthogonal Basis = Logical True ! Save projector = Logical True ! Save line = Logical True ! Save Contact = Logical True End Boundary Condition 2 Name = "Top" Target Boundaries(1) = 7 Displacement 1 = Variable Time Real 0.0, 0.0 1.0, .35 end Displacement 2 = 0.0 Displacement 3 = -0.05 ! Force 2 = -1000. End Boundary Condition 3 Name = "Sides" Target Boundaries(2) = 3 5 End Boundary Condition 4 Name = "BaseBottom" Target Boundaries(5) = 1 3 5 9 11 Displacement 1 = 0.0 Displacement 2 = 0.0 Displacement 3 = 0.0 End Boundary Condition 5 Name = "BaseTop" Target Boundaries(1) = 4 Normal-Tangential Displacement = Logical True Mass Consistent Normals = Logical True ! If these are not obtained they are fetched from the master ! Dynamic Friction Coefficient = Variable time ! Real MATC "mu*(tx-1)" ! Contact Velocity(3) = Real -1.0 0.0 0.0 End