找回密码
 注册
Simdroid-非首页
查看: 98|回复: 8

[命令/FISH] trip是什么意思?

[复制链接]
发表于 2009-10-9 15:53:33 | 显示全部楼层 |阅读模式 来自 江苏徐州
本帖最后由 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  
在手册上没找到相关的介绍,请给新手指点一下……谢谢了
 楼主| 发表于 2009-10-9 15:54:14 | 显示全部楼层 来自 江苏徐州
Simdroid开发平台
完整代码如下:
;-----------------------------------------------------------
;  剑桥模型各向同性压缩实验例题
;-----------------------------------------------------------
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
回复 不支持

使用道具 举报

发表于 2009-10-9 16:53:04 | 显示全部楼层 来自 安徽淮南
本帖最后由 lookcity 于 2009-10-9 16:54 编辑

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

使用道具 举报

发表于 2009-10-9 19:30:26 | 显示全部楼层 来自 江苏徐州
学习啦
1111111111
回复 不支持

使用道具 举报

 楼主| 发表于 2009-10-10 08:56:35 | 显示全部楼层 来自 江苏徐州
谢谢了,在积累当中…… 3# lookcity
回复 不支持

使用道具 举报

 楼主| 发表于 2009-10-16 10:05:27 | 显示全部楼层 来自 江苏徐州
3# lookcity 这里是相当于再重新调用 前面自定义的 trip 循环式吧?相当于把def trip……再运行一此。可对?
回复 不支持

使用道具 举报

发表于 2009-10-16 14:04:38 | 显示全部楼层 来自 安徽淮南
函数定义:
def trip
    ....
end
函数调用:
trip
回复 不支持

使用道具 举报

发表于 2010-11-12 22:09:31 | 显示全部楼层 来自 四川成都
调用函数。。。
回复 不支持

使用道具 举报

发表于 2011-2-24 18:54:51 | 显示全部楼层 来自 山东青岛
函数调用的
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-9-25 06:28 , Processed in 0.039184 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表