Skip to content

Commit

Permalink
Adding an example of a periodic beam structure
Browse files Browse the repository at this point in the history
  • Loading branch information
mmalinen committed Apr 11, 2024
1 parent 7a16645 commit d8c398b
Show file tree
Hide file tree
Showing 2 changed files with 105 additions and 1 deletion.
2 changes: 1 addition & 1 deletion fem/tests/Beam_3D_Cantilever/beam1d/mesh.boundary
@@ -1,2 +1,2 @@
1 1 0 1 101 1
2 1 10 0 101 2
2 2 10 0 101 2
104 changes: 104 additions & 0 deletions fem/tests/Beam_3D_Cantilever/periodic.sif
@@ -0,0 +1,104 @@
! An example of a periodic beam structure
!

Check Keywords "Warn"

Header
Mesh DB "." "beam1d"
End

Simulation
Max Output Level = 5
Coordinate System = Cartesian 3D
Simulation Type = Steady
Output Intervals = 1
Steady State Max Iterations = 1
Post File = "periodic.vtu"
End

Body 1
Equation = 1
Material = 1
Body Force = 1
End

Material 1
Youngs Modulus = Real 2.0e-1
Shear Modulus = Real 1.0
! Youngs Modulus = Real 70.0e+9
! Shear Modulus = Real 26.0e+9

! Specify basis where the cross section parameters are expressed:
Principal Direction 2(3) = Real 0.0 0.0 -1.0

! Specify I_k:
Second Moment of Area 2 = Real 1.0
Second Moment of Area 3 = Real 1.0
! Specify A:
Cross Section Area = Real 1.0
! Specify J_T:
Torsional Constant = Real 1.0
! The density is used only in transient cases:
Density = 2700.0
End

Body Force 1
Body Force 1 = 0.0
Body Force 2 = Variable Coordinate
Real MATC "1.0e-2*cos(2.0*pi*tx(0))"
Body Force 3 = 0.0
End

Equation 1 :: Active Solvers(1) = 1

Solver 1
Equation = "Timoshenko Beam Equations"
Procedure = "BeamSolver3D" "TimoshenkoSolver"

Nonlinear System Max Iterations = 1
Nonlinear System Convergence Tolerance = 1e-7

Linear System Solver = "Iterative"
Linear System Preconditioning = Diagonal
Linear System Max Iterations = 1000
Linear System Convergence Tolerance = 1e-9
Linear System Iterative Method = IDRS
! Linear System Abort Not Converged = False
Steady State Convergence Tolerance = 1e-9
Apply Conforming BCs = True
End

! -------------------------------------------------------------------------
! The following can be used to save the deflection and the rotation
! at the end:
! -------------------------------------------------------------------------
Solver 2
Equation = "Save Scalars"
! Exec Solver = After Timestep
Procedure = "SaveData" "SaveScalars"
! Filename = cantilever.dat
Variable 1 = U 1
Variable 2 = U 2
Variable 3 = U 3
Variable 4 = Theta 1
Variable 5 = Theta 2
Variable 6 = Theta 3
Save Points(1) = 2
End

Boundary Condition 1
Target Boundaries(1) = 1
U 1 = Real 0.0
U 3 = Real 0.0
Theta 1 = Real 0.0
Theta 2 = Real 0.0
Theta 3 = Real 0.0
End

Boundary Condition 2
Target Boundaries(1) = 2
Conforming BC = 1
End

Solver 1 :: Reference Norm = Real 9.73572657E-05
Solver 1 :: Reference Norm Tolerance = Real 1.0e-5

0 comments on commit d8c398b

Please sign in to comment.