Code: Select all
Header
Mesh DB "." "."
End
Constants
End
Simulation
Max Output Level = 3
Coordinate System = Cartesian 2D
Simulation Type = Steady State
Steady State Max Iterations = 5
Steady State Min Iterations = 1
Output Intervals = 1
Output File = "piezo.result"
Post File = "piezo.ep"
End
Body 1
Equation = 1
Material = 1
Body Force = 1
End
Body 2
Equation = 1
Material = 2
Body Force = 2
End
Body Force 1
Stress BodyForce 2 = -1e8
Piezo Material = Logical True
! Stress(3) = Variable Electric Field 2
! Real
! 0 -1e8 0 0
! 1 -1e8 0 0
! End
End
Body Force 2
Stress BodyForce 2 = -1e8
Piezo Material = Logical True
! Stress(3) = Variable Electric Field 2
! Real
! 0 1 0 0
! 1 1 0 0
! End
End
Equation 1
Active Solvers(2) = 1 2
End
Solver 2
Equation = Electrostatics
Variable = Potential
Variable DOFs = 1
Procedure = "StatElecSolve" "StatElecSolver"
Linear System Solver = Direct
Linear System Direct Method = UMFPack
Linear System Iterative Method = BiCGStab
Linear System Preconditioning = ILU1
Linear System Max Iterations = 500
Linear System Convergence Tolerance = 1.0e-8
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Newton After Iterations = 20
Nonlinear System Max Iterations = 1000
Nonlinear System Convergence Tolerance = 1.0e-5
Nonlinear System Relaxation Factor = 1.0
Steady State Convergence Tolerance = 1.0e-4
End
Solver 1
! Exec SOlver = Never
Equation = Stress Analysis
Variable = Displacement
Variable DOFs = 2
Linear System Solver = Direct
Linear System Direct Method = UMFPack
Linear System Iterative Method = BiCGStab
Linear System Preconditioning = ILU1
Linear System Max Iterations = 500
Linear System Convergence Tolerance = 1.0e-8
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Newton After Iterations = 20
Nonlinear System Max Iterations = 1
Nonlinear System Convergence Tolerance = 1.0e-5
Nonlinear System Relaxation Factor = 1.0
Steady State Convergence Tolerance = 1.0e-4
End
Material 1
Density = 2e3
Youngs Modulus = 70e9
Poisson Ratio = 0.3
Relative Permittivity = 1
Piezo Material Coefficients(2,4) = Real \
0 10 0 0 \
10 0 0 0
End
Material 2
Density = 2e3
Youngs Modulus = 70e9
Poisson Ratio = 0.3
Relative Permittivity = 1
Piezo Material Coefficients(2,4) = Real \
0 -10 0 0 \
-10 0 0 0
End
Boundary Condition 1
Target Boundaries(2) = 1 2 ! 5 6
Displacement 1 = 0
Displacement 2 = 0
End
Boundary Condition 2
Target Boundaries = 3
Potential = 0
End
Boundary Condition 3
Target Boundaries = 4
Potential = 1
End
$fprintf( stderr, "TEST CASE 1\n");
RUN
$fprintf( stderr, "END TEST CASE 1: Target NRM=0.595127795559, EPS=1.0e-3\n" );