strain rate与strain increment关系验证

11月前浏览4132

摘要：本文对FLAC3D中的strain rate与strain increment的关系进行验证。

1.问题来源

在Itasca专有技术交流群（289141433）中，Itasca官方的老师对FLAC3D中的strain increment与strain rate的关系进行了说明，现对该说法进行验证。

2.代码

model newmodel largestrain off
zone create brick size 1 1 1zone cmodel assign elasticzone property young 1e8 poisson 0.35 density 1800
zone face skinzone face apply velocity-normal 0 range group 'south' or 'north'zone face apply velocity-normal 0 range group 'east' or 'west'zone face apply velocity-normal 0 range group 'bottom'
history interval 1zone history name 'xxInc' strain-increment quantity xx position  1 1 1 method constantzone history name 'xxRate'  strain-rate quantity xx position  1 1 1 method constantzone history name 'yyInc' strain-increment quantity yy position  1 1 1 method constantzone history name 'yyRate'  strain-rate quantity yy position  1 1 1 method constantzone history name 'zzInc' strain-increment quantity zz position  1 1 1 method constantzone history name 'zzRate'  strain-rate quantity zz position  1 1 1 method constant
model gravity 9.81model solve 
history export 'xxInc' table '1'history export 'xxRate' table '2'history export 'yyInc' table '3'history export 'yyRate' table '4'history export 'zzInc' table '5'history export 'zzRate' table '6'
fish def _dataTrans    loop i(1,72)        j = i - 1        table.x(7,i) = i        table.x(8,i) = i        table.x(9,i) = i        if i = 1            table.y(7,i) = table.y(1,i)            table.y(8,i) = table.y(3,i)            table.y(9,i) = table.y(5,i)        endif        if i > 1            table.y(7,i) = table.y(1,i) - table.y(1,j) ;xx rate            table.y(8,i) = table.y(3,i) - table.y(3,j) ;yy rate             table.y(9,i) = table.y(5,i) - table.y(5,j) ;zz rate        endif    endloopend@_dataTrans