From Import From Import From Import From Import From Import From Import From Import From Import From Import From Import From Import From Import From Import From Import
From Import From Import From Import From Import From Import From Import From Import From Import From Import From Import From Import From Import From Import From Import
####BRICK GEOMETRY#####
# Brick Numbers & Size
x = 5 # number of bricks in x-dir.
y = 8 # number of bricks in y-dir.
brick_length = 0.024 # (metre)
brick_height = 0.012 # brick height (metre)
brick_width = 0.012 # (metre)
#Property##
E = 4e9 #material's young modulus (newton/sq.metre)
v = 0.25 #material's poissons ratio
density = 750 #(kilograms/cu.metre)
#Load##
velocity = 0.001 #"v1" displacement velocity (metre/second)
#Step##
step_self = 'Step-1'
timeperiod_self = 1 #Step time (seconds)
step_push = 'Step-2'
timeperiod = 10 #Step time (seconds)
nonlinear_geom = ON #"nlgeom" non-linear geometry(ON, OFF)
max_in_num = 2000 #"maxNumInc" maximum number of increments
ini_inc = 0.001 #"initialInc" initial increment size
min_inc = 1e-10 #"minInc" minimum increment size
sup_calc = ON #"nohaf" Supress calculation
amp = STEP #"amplitude" Instantaneous (STEP, RAMP)
##Solution Controls##
propagation = OFF #"allowPropogation" propagate from previous(ON, OFF)
discont = OFF #"discontinuous" (ON, OFF)
def_val = OFF #"resetDefaultValues" (ON, OFF)
time_inc=(8.0, 10.0, 9.0, 16.0, 10.0, 4.0, 12.0, 5.0, 6.0, 3.0, 50.0)
#"timeIncrementation"(Io, Ir, Ip, Ic, Il, Ig, Is, Ia, Ij, It, Ics)
line_search=(0.0, 1.0, 0.0001, 0.25, 0.1) #"lineSearch"(Nls, Slsmax, Slsmin, flss, nls)
#Interaction Property
##Friction
'Friction' #name
fric_form = PENALTY #"formulation"
direction = ISOTROPIC #"directionality"
stress_limit = None #"shear stress limit"
fric_coeff = 0.01 #"table" friction coefficient
frac = 0.005 #"fraction" maximum elastic slip
##Normal###
'Normal' #name
method = PENALTY #"constraintEnforcementMethod"
separate = ON #"allowSeparation"
##Cohesive##
'Cohesive' #name
Knn = 1e10 #"table" stiffness-normal (newton/metre)
Kss = 1e10 #"table" stiffness-shear (newton/metre)
Ktt = 1e10 #"table" stiffness-tangential (newton/metre)
coh_eligibility = INITIAL_NODES #"eligibility" slave-nodes initially in contact
##Damage##
'Damage' #name
crit = QUAD_TRACTION #"criterion"
init_normal = 1e6 #"iniTable" damage initiating stress(newton/sq.metre)
init_shear = 1e6 #"iniTable" damage initiating stress(newton/sq.metre)
init_tangential = 1e6 #"iniTable" damage initiating stress(newton/sq.metre)
evol_value = 0.002 #"evolTable" (metre)
stab_viscosity = 0.002 #viscosityCoeff" stabilization vicosity coeff
#Interaction###
inter_brick = ['Cohesive+Friction+Damage','Cohesive+Friction+50%Damage']
sliding_val_1 = SMALL #"sliding"
wall_to_support = 'Normal'
sliding_val_2 = FINITE #"sliding"
#Mesh###
meshsize = 0.012 #element size (metre)
meshtype = CPS4 #"elemCode" (quad plane stress full integration)
elementshape = QUAD #"elemShape" Quadrilateral(QUAD, QUAD_DOMINATED, TRI)