 鲜花( 12)  鸡蛋( 0)
|
本帖最后由 三T上人 于 2016-7-22 17:20 编辑 <br /><br />源代码如下:
title '混凝土预制管片抗压试验分析'
gen z cshell p0 0 0 0 p1 3 0 0 p2 0 1.5 0 p3 0 0 3 p4 3 1.5 0 p5 0 1.5 3 p8 2.7 0 0 p9 0 0 2.7 p10 2.7 1.5 0 p11 0 1.5 2.7 d 5 5 5 5 s 3 6 24
m m
prop b 4.2e10 s 1.668e10 coh 2.72e5 fri 44 ten 2e5
;prop ctab 1 ftab 2
;table 1 0,2.72e5 1e-4,2e5 2e-4,1.5e5 3e-4,1.03e5 1,1.03e5
;table 2 0,44 1e-4,42 2e-4,40 3e-4,38 1,38
fix x y z ran x 2.69 2.71 z -.1 .1
fix x y z ran z 2.69 2.71 x -.1 .1
;
;app xforce -215e3 zforce -215e3 ran x 1.6 1.7 y 0 1.5 z 2.45 2.55
;app xforce -215e3 zforce -215e3 ran x 2.55 2.45 y 0 1.5 z 1.6 1.7
;
def xxload
n=1
loop while n<1800
load=n*ini_load
lx_load=.0006*load
lz_load=.0008*load
rx_load=.0008*load
rz_load=.0006*load
command
app xforce lx_load zforce lz_load ran x 1.6 1.7 y 0 1.5 z 2.45 2.55
app xforce rx_load zforce rz_load ran x 2.55 2.45 y 0 1.5 z 1.6 1.7
step 1
endcommand
n=n+1
endloop
end
;
set ini_load=-1e3
his unbal
his load
;his zone szz id 189
his gp zdis 1.91 .75 1.91
his gp zdis 1.5 .75 2.243
his gp zdis 2.243 .75 1.5
;his gp xdis 1.329 .75 2.69
;his gp xdis 2.69 .75 1.329
his gp xdis .1976 .75 2.993
his gp xdis 2.993 .75 .1976
pl his -3 mark -4 mark -6 mark vs -2
SOSO.CC |
|