Skip to content

Commit

Permalink
Add content
Browse files Browse the repository at this point in the history
  • Loading branch information
@loganharbour
loganharbour committed Aug 30, 2023
1 parent cda0e7f commit f22477f
Show file tree
Hide file tree
Showing 33 changed files with 5,527 additions and 0 deletions.
3 changes: 3 additions & 0 deletions 64_Fiber_Failure/fiber_64.e
View file
Git LFS file not shown
402 changes: 402 additions & 0 deletions 64_Fiber_Failure/fiber_64_biaxis_out_0.csv
View file

Large diffs are not rendered by default.

3 changes: 3 additions & 0 deletions 64_Fiber_Failure/fiber_64_biaxis_out_0.e
View file
Git LFS file not shown
38 changes: 38 additions & 0 deletions 64_Fiber_Failure/large_p_fiber_64_biaxis_out_0.csv
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,38 @@
time,Reaction_Force_Along_X,Reaction_Force_Along_Y
0,0,0
0.01,0.00043903097578746,0.00043903565520959
0.02,0.00087806622776507,0.00087807560352611
0.03,0.0013171052189148,0.0013171192924236
0.04,0.0017561479369537,0.0017561667094654
0.05,0.0021951943811467,0.0021952178539052
0.06,0.0026342445508681,0.0026342727251077
0.07,0.0030732984491005,0.0030733313260914
0.08,0.0035123560737742,0.0035123936547594
0.09,0.0039514174242096,0.0039514597104217
0.1,0.0043904825166967,0.0043905295095701
0.11,0.004829551324882,0.0048296030255216
0.12,0.0052686238644877,0.0052686802741921
0.13,0.0057077001399254,0.0057077612600469
0.14,0.006146780133199,0.0061468459648651
0.15,0.0065858638654555,0.0065859344100544
0.16,0.0070249513215621,0.0070250265802932
0.17,0.0074640425130206,0.0074641224872237
0.18,0.007903137422541,0.0079032221133406
0.19,0.0083422360721631,0.0083423254809557
0.2,0.0087813384446266,0.0087814325725932
0.21,0.0092204445497844,0.0092205433982265
0.22,0.0096595543874654,0.0096596579576809
0.23,0.010098667951737,0.010098776244948
0.24,0.010537785242686,0.010537898260117
0.25,0.010976906268512,0.010977024011484
0.26,0.011416031017946,0.011416153487641
0.27,0.011855159502164,0.011855286699901
0.28,0.012294291717624,0.012294423644676
0.29,0.012733427661912,0.012733564319521
0.3,0.013172567336134,0.013172708725555
0.31,0.013451557315447,0.013502500273785
0.32,0.013042156570533,0.013072935860226
0.33,0.009525519206865,0.0092854427322531
0.335,0.007274843031071,0.0068844361778043
0.345,0.006554717427005,0.0060426422327763
0.355,0.0060176866015123,0.0055165584682952
3 changes: 3 additions & 0 deletions 64_Fiber_Failure/large_p_fiber_64_biaxis_out_0.e
View file
Git LFS file not shown
348 changes: 348 additions & 0 deletions 64_Fiber_Failure/model_fiber_64_biaxis.i
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,348 @@
[Mesh]
[fiber_model_mesh]
type = FileMeshGenerator
file = fiber_64.e
[]
[breakmesh]
type = BreakMeshByBlockGenerator
input = fiber_model_mesh
block_pars = 'Matrix Fiber'
interface_name = matrix_fiber
[]

[]

[GlobalParams]
displacements = 'disp_x disp_y disp_z'
[]

[Modules/TensorMechanics/Master]
[all]
strain = FINITE
add_variables = true
generate_output = 'stress_xx stress_yy stress_zz stress_xy stress_xz stress_yz hydrostatic_stress vonmises_stress'
extra_vector_tags = react
material_output_family = MONOMIAL
material_output_order = FIRST
[]
[]
[InterfaceKernels]
[ifk_x]
type = Rczm
component = 0
variable = 'disp_x'
neighbor_var = 'disp_x'
displacements = 'disp_x disp_y disp_z'
#boundary = 'interface'
boundary = 'matrix_fiber'
stabilized_para =2
[]
[ifk_y]
type = Rczm
component = 1
variable = 'disp_y'
neighbor_var = 'disp_y'
displacements = 'disp_x disp_y disp_z'
#boundary = 'interface'
boundary = 'matrix_fiber'
stabilized_para = 2
[]
[ifk_z]
type = Rczm
component = 2
variable = 'disp_z'
neighbor_var = 'disp_z'
displacements = 'disp_x disp_y disp_z'
#boundary = 'interface'
boundary = 'matrix_fiber'
stabilized_para = 2
[]
[]

[Materials]
[Elasticity_Matrix]
type = ComputeIsotropicElasticityTensor
youngs_modulus = 0.003
poissons_ratio = 0.3
block = 'Matrix'
[]
[Elasticity_Fiber]
type = ComputeIsotropicElasticityTensor
youngs_modulus = 0.072
poissons_ratio = 0.21
block = 'Fiber'
[]
[stress]
type = ComputeFiniteStrainElasticStress
block = 'Matrix Fiber'
[]
[./Rigid_czm_model]
type = OrtizPandolfiRczm
boundary = 'matrix_fiber'
maximum_effective_traction = 8e-6
maximum_effective_gap = 0.01
shear_weight = 0.707
displacements = 'disp_x disp_y disp_z'
outputs = all
[../]
[]

[Functions]
[topfunc]
type = ParsedFunction
value = 't'
#value = '0.1*sin(3.1415926/2.0*t)'
[]
[]

[BCs]
# Prescribed displacements at three boundary surfaces
# Enforce plain strain condition by constrainting
# zero-displacement on z_plus and z-minus

[roller_at_X_minus]
type = DirichletBC
preset = true
variable = disp_x
value = 0
boundary = 'X_Minus'
[]

[roller_at_Y_minus]
type = DirichletBC
preset = true
variable = disp_y
value = 0
boundary = 'Y_Minus'
[]

[roller_at_Z_minus]
type = DirichletBC
preset = true
variable = disp_z
value = 0
boundary = 'Z_Minus'
[]

[roller_at_Z_plus]
type = DirichletBC
preset = true
variable = disp_z
value = 0
boundary = 'Z_Plus'
[]

# restrict the rigid body motion after fully debonding
#[spring_fiber_x]
# type = PenaltyDirichletBC
# variable = disp_x
# value = 0
# penalty = 10
# boundary = 'fiber_Z_minus'
#[]
#[spring_fiber_y]
# type = PenaltyDirichletBC
# variable = disp_y
# value = 0
# penalty = 10
# boundary = 'fiber_Z_minus'
#[]
# [spring_fiber_z]
# type = PenaltyDirichletBC
# variable = disp_z
# value = 0
# penalty = 10
# boundary = 'fiber_Z_minus'
# []
#
#Apply displacement at X_plus
#
[U_at_X_plus]
type = FunctionDirichletBC
variable = disp_x
function = topfunc
preset = true
boundary = 'X_Plus'
[]
[U_at_Y_plus]
type = FunctionDirichletBC
variable = disp_y
function = topfunc
preset = true
boundary = 'Y_Plus'
[]
[]
[AuxVariables]
[stress_xx]
order = FIRST
family = MONOMIAL
[]
[stress_yy]
order = FIRST
family = MONOMIAL
[]
[stress_zz]
order = FIRST
family = MONOMIAL
[]
[stress_xy]
order = FIRST
family = MONOMIAL
[]
[stress_xz]
order = FIRST
family = MONOMIAL
[]
[stress_yz]
order = FIRST
family = MONOMIAL
[]
# [hydrostatic_stress]
# order = FIRST
# family = MONOMIAL
# []
# [vonmises_stress]
# order = FIRST
# family = MONOMIAL
# []
[RF_X]
order = FIRST
[]

[RF_Y]
order = FIRST
[]
[]

[AuxKernels]
[stress_xx]
type = RankTwoAux
rank_two_tensor = stress
variable = stress_xx
index_i = 0
index_j = 0
[]
[stress_yy]
type = RankTwoAux
rank_two_tensor = stress
variable = stress_yy
index_i = 1
index_j = 1
[]
[stress_zz]
type = RankTwoAux
rank_two_tensor = stress
variable = stress_zz
index_i = 2
index_j = 2
[]
[stress_xy]
type = RankTwoAux
rank_two_tensor = stress
variable = stress_xy
index_i = 0
index_j = 1
[]
[stress_xz]
type = RankTwoAux
rank_two_tensor = stress
variable = stress_xz
index_i = 0
index_j = 2
[]
[stress_yz]
type = RankTwoAux
rank_two_tensor = stress
variable = stress_yz
index_i = 1
index_j = 2
[]
#[hydrostatic_stress]
# type = RankTwoAux
# rank_two_tensor = stress
# variable = hydrostatic_stress
#[]
#[vonmises_stress]
# type = RankTwoAux
# rank_two_tensor = stress
# variable = vonmises_stress
#[]
[RF_X]
type = TagVectorAux
variable = RF_X
v = disp_x
vector_tag = react
[]
[RF_Y]
type = TagVectorAux
variable = RF_Y
v = disp_y
vector_tag = react
[]
[]

[Postprocessors]
[Reaction_Force_Along_X]
type = NodalSum
boundary = X_Plus
variable = RF_X
[]
[Reaction_Force_Along_Y]
type = NodalSum
boundary = Y_Plus
variable = RF_Y
[]
[]

[Problem]
extra_tag_vectors = react
[]

[Preconditioning]
[SMP]
type = SMP
full = true
[]
[]

[Executioner]
type = Transient
#[Quadrature]
# side_order = 'constant'
#[]
dt = 0.005
dtmax = 0.1
dtmin = 0.000001
end_time = 2

petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'
# petsc_options_value = 'lu'
petsc_options_value = 'lu superlu_dist'
solve_type = PJFNK

line_search = none

# end_time = 1.0
# dt = 0.1
l_max_its = 30
# l_tol = 1e-6
nl_max_its = 15

# nl_abs_tol = 1e-8
# nl_rel_tol = 1e-8

nl_abs_tol = 1e-8
[]

[Outputs]
file_base = fiber_64_biaxis_out_0
exodus = true
[csv]
type = CSV
[]
# checkpoint = true
[]

[Debug]
show_material_props = true
[]
Loading

0 comments on commit f22477f

Please sign in to comment.