The OpenSees Community

Code: Select all

#~~~~~ Material Parameters ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
			   #  1   2   3      4     5      6    7    8     9      10   11    12   13   14   15   16  17    18    19
set MatCons		"$rho $e0 $G0 $poisson $P_at $k_c $M_csl $c $lambda_c $xi $e_r $m_low $h0 $c_h $n_b $A0 $n_d $z_max $c_z" 
			
set intTens		"$alpha_T $z_T"

set MP			"MaterialConstant 19 $MatCons InternalTensor 2 $intTens"



#~~~~~ Elastic State ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
set ES	     	"Dafalias-Manzari 3 4 5 6 2 $initS"


#~~~~~ Yield Function ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
set YF	      	"Dafalias-Manzari 0 12 2 1"


#~~~~~ Plastic Flow ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
set PF	      	"Dafalias-Manzari 0 2 0 11 0 9 0 10 0 5 0 12 0 7 0 8 0 16 0 17 2 1 2 2"


#~~~~~ Tensor Evolution Law ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
set tel 		"Dafalias-Manzari 2 11 9 10 5 12 7 8 15 13 14 3 1 2"
set ftel		"Dafalias-Manzari-fabric 12 18 19 1 2"	

set TE	      	"$tel $ftel"


nDMaterial NewTemplate3Dep 1 	-MaterialParameter $MP \
								-ElasticState $ES \
								-YieldFunction $YF \
								-PlasticFlow $PF \
								-TensorEvolution $TE \
								-Algorithm 0

##############################################################

for {set i 1} {$i <= $Num_Elem} {incr i} {
	set N1 [expr 4*$i+1] ; set N2 [expr 4*$i+2]
	set N3 [expr 4*$i+3] ; set N4 [expr 4*$i+4]
	set N5 [expr 4*$i-3] ; set N6 [expr 4*$i-2]
	set N7 [expr 4*$i-1] ; set N8 [expr 4*$i-0]
	element Brick8N_u_p_U  $i $N1 $N2 $N3 $N4 $N5 $N6 $N7 $N8 1 0.0 0.0 $g $poro $alpha $rho_s $rho_f \
				     	   $perm_x $perm_y $perm_z $B_s $B_f $PRESS
}