论坛 › G0201：FLAC3D/FLAC2D › trip是什么意思？

hgshgs 发表于 2009-10-9 15:53:33

trip是什么意思？

本帖最后由 hgshgs 于 2009-10-9 15:57 编辑

; --- 采样记录 ---
his nstep 20
his unbal
his path
his sp
his lnp
his sq
his svol
his mk
his mg
his gp zdisp 0 0 1
; --- 试验 ---
Trip
; --- 结果 ---
plot his 3 vs -9
plot his 6 vs 4
plot his 7 8 vs -9
save camiso.sav

trip 试下单独运行的时候，像solve或step 1000
在手册上没找到相关的介绍,请给新手指点一下……谢谢了

hgshgs 发表于 2009-10-9 15:54:14

完整代码如下：
;-----------------------------------------------------------
;剑桥模型各向同性压缩实验例题
;-----------------------------------------------------------
new
gen zone brick p0 0 0 0 p1 1 0 0 p2 0 1 0 p3 0 0 1 size 1 1 1
title
Isotropic compression test for normally consolidated soil
; --- 材料参数 ---
model cam-clay
prop shear 250. bulk_bound 10000.
prop mm 1.02 lambda 0.2 kappa 0.05
prop mpc 5. mp1 1. mv_l 3.32
; --- 边界及初始条件 ---
fix x y z
ini sxx -5. syy -5. szz -5.
ini xvel -0.5e-4 range x 0.9 1.1
ini yvel -0.5e-4 range y 0.9 1.1
ini zvel -0.5e-4 range z 0.9 1.1
; --- 定义fish 函数 --- (numerical values for p, q, v)
def camclay_ini_p
        p_z = zone_head
        loop while p_z # null
                mean_p = -(z_sxx(p_z) + z_syy(p_z) + z_szz(p_z))/3.0 - z_pp(p_z)
                z_prop(p_z,'cam_cp') = mean_p
                p_z = z_next(p_z)
        endloop
end
camclay_ini_p
def path
        p_z = zone_head
        sp = z_prop(p_z,'cam_cp')
        sq = z_prop(p_z,'cq')
        sqcr= sp*z_prop(p_z,'mm')
        if sp = 0.0 then
                sp = 1.
        endif
        lnp = ln(sp)
        svol = z_prop(p_z,'cv')
        mk = z_prop(p_z,'bulk')
        mg = z_prop(p_z,'shear')
        cpc = z_prop(p_z,'mpc')
end
; ... 装载-卸载练习 ...
def trip
        loop i (1,5)
                command
                        ini xvel -0.5e-4 range x 0.9 1.1
                        ini yvel -0.5e-4 range y 0.9 1.1
                        ini zvel -0.5e-4 range z 0.9 1.1
                        step 300
                        ini xvel mul -.1 yvel mul -.1 zvel mul -.1
                        step 1000
                        ini xvel mul -1. yvel mul -1. zvel mul -1.
                        step 1000
                end_command
        end_loop
end
; --- 采样记录 ---
his nstep 20
his unbal
his path
his sp
his lnp
his sq
his svol
his mk
his mg
his gp zdisp 0 0 1
; --- 试验 ---
Trip
; --- 结果 ---
plot his 3 vs -9
plot his 6 vs 4
plot his 7 8 vs -9
save camiso.sav

lookcity 发表于 2009-10-9 16:53:04

本帖最后由 lookcity 于 2009-10-9 16:54 编辑

Trip是用户自定义函数
作用是分级反复加卸载

zhangruian 发表于 2009-10-9 19:30:26

学习啦
1111111111

hgshgs 发表于 2009-10-10 08:56:35

谢谢了，在积累当中…… 3# lookcity

hgshgs 发表于 2009-10-16 10:05:27

3# lookcity 这里是相当于再重新调用 前面自定义的 trip 循环式吧？相当于把def trip……再运行一此。可对？

lookcity 发表于 2009-10-16 14:04:38

函数定义：
def trip
    ....
end
函数调用：
trip

高山@流水 发表于 2010-11-12 22:09:31

调用函数。。。

wej101 发表于 2011-2-24 18:54:51

函数调用的

查看完整版本: trip是什么意思？