Therefore I made a simple Cylinder mesh and used the StressSolver and the MeshSolver. The StressSolver works well but the MeshSolver does not seem to change the Mesh at all.

Here I will post a part of my Sif file:

Solver 1 is the StresSolver which does not do anything at the moment, because I want to resolve the Issue with the MeshSolver first.

Solver 3 is the ResultOutputSolver so i can see the result in paraview.

I have therefore not inculded them.

- Code: Select all
`Simulation`

!-----------------OutPut-----------------

Max Output Level = 20

Coordinate System = Cartesian

Simulation Type = Transient

!-----------------Iterations-----------------

Steady State Max Iterations = 10

Steady State Min Iterations = 1

!-----------------Time Intervals-----------------

Timestep Intervals = 100

Timestep Size = real 1

!adaptive_in_time=Logical False !Logical True

!-----------------Numerical Parameters-----------------

tolerance=real 1.0e-07 !tolerance for checking if to points are diffrent must be significantly smaller then the samllest Material

!-----------------Jellyroll Parameters-----------------

End

Equation 1

Active Solvers (3) = 1 2 3

Stress Analysis = Logical true

End

! MATERIALS

Material 1

Name = "Copper (generic)"

Heat expansion Coefficient = 16.5e-6

Heat Conductivity = 401.0

Sound speed = 3810.0

Heat Capacity = 385.0

Mesh Poisson ratio = 0.34

Density = 8960.0

Poisson ratio = 0.34

Youngs modulus = 115.0e9

End

Solver 2

Equation = Mesh Update

Variable = -dofs 3 Mesh Update

Procedure = "MeshSolve" "MeshSolver"

Exec Solver = Always

Stabilize = True

Bubbles = False

Lumped Mass Matrix = False

Optimize Bandwidth = True

Steady State Convergence Tolerance = 1.0e-5

Nonlinear System Convergence Tolerance = 1.0e-7

Nonlinear System Max Iterations = 20

Nonlinear System Newton After Iterations = 3

Nonlinear System Newton After Tolerance = 1.0e-3

Nonlinear System Relaxation Factor = 1

Linear System Solver = Iterative

Linear System Iterative Method = BiCGStab

Linear System Max Iterations = 500

Linear System Convergence Tolerance = 1.0e-10

BiCGstabl polynomial degree = 2

Linear System Preconditioning = Diagonal

Linear System ILUT Tolerance = 1.0e-3

Linear System Abort Not Converged = False

Linear System Residual Output = 1

Linear System Precondition Recompute = 1

End

Boundary Condition 1

!====== Displacement on -y =====

Target Boundaries(1) = 1

Displacement 1=0

Displacement 2=0

Displacement 3=0

Mesh Update 1=0

Mesh Update 2=0

Mesh Update 3=0

End

Boundary Condition 2

!====== Displacement on Zylinder Mantel =====

Target Boundaries(1) = 3

Displacement 1=0

Displacement 2=0

Displacement 3=0

Mesh Update 1=0

Mesh Update 2=0

Mesh Update 3=0

End

Boundary Condition 3

!====== Displacement on y+ =====

Target Boundaries(1) = 2

Normal Force = real 1e6

Mesh Update 1= real 0.0

Mesh Update 2= real 1

Mesh Update 3= real 0.0

!Mesh Update 1= Equals Displacement 1

!Mesh Update 2= Equals Displacement 2

!Mesh Update 3= Equals Displacement 3

End

One problem which might be part of this might be that the DisplaceMesh function (which should be defined in MeshUtils) is not defined anywhere.

I would really appreciate some help. I Do not know what i am doing wrong.

Kaputo