|
|
|
####################################################################################
# HyperbolicGapUnloadReloadResponse
#
# Description:
#
# Determine response of hyperbolic gap material under unloading/reloading.
#
# Revision: $WCREV$
# Date: $WCDATE$
# URL: $WCURL$
# Date Created: 04/2008
# Author: Matthew Dryden (dryden@berkeley.edu)
####################################################################################
# define units
set in 1.;
set sec 1.;
set kip 1.;
set ksi [expr $kip/pow($in,2)];
set psi [expr $ksi/1000.];
set ft [expr $in*12.];
# define constants
set g [expr 32.2*$ft/pow($sec,2)];
set PI [expr 2*asin(1.0)];
set MatID 1;
set Kmax 20300;
set Kur 20300;
set Rf 0.7;
set Fult -326;
set gap -0.025;
model test 1
uniaxialMaterial HyperbolicGapMaterial $MatID $Kmax $Kur $Rf $Fult $gap
uniaxialTest 1
set count 0;
set maxCount 3229;
set factor 1.0;
set fileID [open HyperbolicGapUnloadReloadStrainPattern.txt r]
set outFileID [open HyperbolicGapUnloadReload.txt w]
foreach line [split [read $fileID] \n] {
incr count 1
if {$count > $maxCount} break;
set strain [lindex $line 0]
if {$strain != ""} {
set strain [expr $strain*$factor]
strainUniaxialTest $strain
set stress [stressUniaxialTest]
puts $outFileID "$strain $stress"
}
}
close $fileID
close $outFileID