求助！岩盐挤压套管

李文良

岩盐由于蠕变逐渐挤压套管，在刚接触时，我认为由于冲击效应（岩盐流变具有速度），致使套管所受法向力波动剧烈，怎么样能减小这种波动或者消除之？
config creep
;---casing geometry---
gene zone cshell group casing p0 0 0 0 p1 0.1223 0 0 p2 0 0.02 0 p3 0 0 0.1223 size 3 1 30 30 dimension 0.1085 0.1085 0.1085 0.1085
;---formation geometry---
gene zone cshell group salt p0 0 0 0 p1 0.5 0 0 p2 0 0.02 0 p3 0 0 0.5 size 15 1 30 30 dimension 0.1556 0.1556 0.1556 0.1556 ratio 1.2 1 1 1
pause
;---interface geometry---
interface 1 face range cylinder end1 0 0 0 end2 0 0.02 0 rad 0.1225 &
                       cylinder end1 0 0 0 end2 0 0.02 0 rad 0.1220 not
;---property of fomation---
model power range group salt
prop bu 1e9 shear 3e8
prop a_1=1e-25 n_1=3
;---property of casing---
model el range group casing
prop bulk 140e9 shear 80e9
;---property of interface---
interface 1 prop kn 6e14 ks 6e14 tens 1e10 SBRATIO = 1
;---boundary condition---
fix y
fix x range x 0
fix z range z 0
apply nstress -100e6 range cylinder end1 0 0 0 end2 0 0.02 0 rad 0.51 &
                           cylinder end1 0 0 0 end2 0 0.02 0 rad 0.49 not
solve

;---max-mises determination---
;call maxt-cr---
;the max_mises=1.93e8
;max_crt=(1.98e8)^(-2)/A*G
;---creep model---
ini xv 0 yv 0 zv 0
;---fish function---
def in_stress
p_in=i_node_head(i_head)
loop while p_in # null
if in_id(p_in)=1 then
    in_stress=in_nstr(p_in)
end_if
p_in = in_next(p_in)
end_loop
end
in_stress
;---hist---
hist nstep 50
hist zone sxx id 1
hist crtime
hist in_stress
set creep on
set cr dt 1.0e-4
set cr lfob 1.0e-5 lmul 1.01 ufob 5.0e-5 umul 0.9
set cr mindt 1.0e-4
set cr maxdt 0.75
set cr dt auto on
set large

solve age 500