case.sif

!---------------------------------------------------------
! Case for solving two displacement equations as coupled.
! This one with monolithic coupled system + direct solver.
! This case additionally with eigen modes.
!
! P.R. 8.9.2021
!------------------------------------------------------------

Header
  Mesh DB "." "beam"
  Include Path ""
  Results Directory ""
End 

Simulation
  Max Output Level = 5
  Coordinate System = Cartesian
  Simulation Type = Steady State
  Output Intervals = 1
  Steady State Max Iterations = 1

  Post File = case.vtu

! Create a union of "dispa" and "dispb"
!  vtu: Vector Field 1 = String "dispa"
!  vtu: Vector Field 1 Complement = String "dispb"
  vtu: Eigen Analysis = Logical True
End

Body 1
  Name = "Body1"
  Equation = 1
  Material = 1 
  Body Force = 1
End

Body 2
  Name = "Body2"
  Equation = 2
  Material = 1
  Body Force = 1
End

Equation 1
  Name = "ElasticityA"
  Active Solvers = 1 
End

Equation 2
  Name = "ElasticityB"
  Active Solvers = 2 
End


Material 1
  Name = "Ideal"
  Youngs Modulus = 1.0e9
  Density = 1.0e3
  Poisson Ratio = 0.3
  Shell Thickness = 0.001
End

Body Force 1
  Name = "Load"
!  Stress BodyForce 2 = -1.0e7
End 


Solver 1
  Equation = LinElastA
  Variable = -dofs 3 DispA
  Procedure = "StressSolve" "StressSolver"

  Nonlinear System Max Iterations = 1
  Nonlinear System Consistent Norm = Logical True

  Linear System Solver = direct
  Linear System Direct Method = MUMPS
  Linear System Scaling = False
  Optimize Bandwidth = False

  Block Monolithic = Logical True
  Block System Mass Lumping = Logical False


  Automated Structure-Structure Coupling = True
  Displace Mesh = False

  Eigen Analysis = True
  Eigen System Values = 10
  Eigen System Convergence Tolerance = Real 1.0e-6
  Eigen System Select = smallest magnitude

!  Eigen System Complex = Logical True
End


Solver 2
  Equation = LinElastB
  Variable = -dofs 3 DispB
  Procedure = "ElasticSolve" "ElasticSolver"
!  Procedure = "StressSolve" "StressSolver"

  Exec Solver = never
  Nonlinear System Max Iterations = 1
!  Nonlinear System Consistent Norm = Logical True

  Eigen Analysis = True
  Eigen System Values = 10
  Nonlinear System Max Iterations = 1
  Optimize Bandwidth = False
  Displace Mesh = Logical False
  Linear System Scaling = False
End

Solver 3
  Equation = SaveScalars
  Procedure = "SaveData" "SaveScalars"
  
!  Filename = f.dat
  Parallel Reduce = Logical True

  Save Eigenvalues = True
  Show Norm Index = 7
End


Boundary Condition 1
  Name = "LeftWall"
  Target Boundaries(1) = 1
  DispA 1 = Real 0.0
  DispA 2 = Real 0.0
  DispA 3 = Real 0.0
End

Boundary Condition 2
  Name = "RightWall"
  Target Boundaries(1) = 2
  DispB 1 = Real 0.0
  DispB 2 = Real 0.0
  DispB 3 = Real 0.0
End


Solver 1:: Reference Norm = 3.37092254E-02
Solver 3:: Reference Norm = 2.74498900E+06