/
inclined_bc_2d.i
110 lines (97 loc) · 1.68 KB
/
inclined_bc_2d.i
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
[GlobalParams]
displacements = 'disp_x disp_y'
[]
[Mesh]
type = GeneratedMesh
dim = 2
nx = 4
ny = 8
xmin = 0.0
xmax = 1.0
ymin = 0.0
ymax = 2.0
elem_type = QUAD4
[]
[MeshModifiers]
[./rotate]
type = Transform
transform = ROTATE
vector_value = '0 0 -60'
[../]
[]
[Modules/TensorMechanics/Master/All]
strain = FINITE
add_variables = true
[]
[BCs]
[./Pressure]
[./top]
boundary = top
function = '-1000*t'
[../]
[../]
[./right_x]
type = PenaltyInclinedBC
variable = disp_x
boundary = right
penalty = 1.0e8
component = 0
[../]
[./right_y]
type = PenaltyInclinedBC
variable = disp_y
boundary = right
penalty = 1.0e8
component = 1
[../]
[./bottom_x]
type = PenaltyInclinedBC
variable = disp_x
boundary = bottom
penalty = 1.0e8
component = 0
[../]
[./bottom_y]
type = PenaltyInclinedBC
variable = disp_y
boundary = bottom
penalty = 1.0e8
component = 1
[../]
[]
[Materials]
[./elasticity_tensor]
type = ComputeIsotropicElasticityTensor
youngs_modulus = 1e6
poissons_ratio = 0.3
[../]
[./stress]
type = ComputeFiniteStrainElasticStress
[../]
[]
[Executioner]
type = Transient
solve_type = 'PJFNK'
petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'
petsc_options_value = 'lu superlu_dist'
# controls for linear iterations
l_max_its = 10
l_tol = 1e-4
# controls for nonlinear iterations
nl_max_its = 100
nl_rel_tol = 1e-8
nl_abs_tol = 1e-8
# time control
start_time = 0.0
dt = 1
end_time = 5
[]
[Preconditioning]
[./smp]
type = SMP
full = true
[../]
[]
[Outputs]
exodus = true
[]