FALC3D单轴压缩求解轴力-轴向应变代码

n
gen zone cylind p0 0 0 0 p1 0.025 0 0 p2 0 0.1 0 p3 0 0 0.025 size 5 20 5
gen zone reflect norm 1,0,0
gen zone reflect norm 0,0,1

model mech strainsoft
pro den 2500 bulk 2e9 she 1e9 co 2e6 fric 45 ten 1e6 dil 10
pro ftab 1 ctab 2 dtab 3
table 1 0 45     0.01,46          0.02,47           0.05,48          0.2,45        0.5,42       1,40
table 2 0 2e6   0.01,3e6        0.02,3.6e6       0.05,4e6        0.2,3.2e6    0.5 1.5e6   1,1e6
table 3 0 10     0.01,11          0.02,11.6        0.05,12          0.2,10        0.5,3         1,0

;边界条件设置

fix yvelocity  range  y -0.001 0.001
fix x z range x,-0.001,0.001 y,-0.001,0.001 z,-0.001,0.001

apply yvel -1e-6 range y 0.099,0.101

;;求解轴向应力
def ax_str
   local str = 0
    local pnt = gp_head
    loop while pnt # null
       if gp_ypos(pnt) < 0.001 then
          str = str + gp_yfunbal(p