flac 3d中如何提取节点主应力
flac 3d中如何提取节点主应力?望高手指点:)好像只能提取单元应力吧,如果要提取节点的可能要fish程序编程,把周围单元应力相加再除以单元数。 yangxin546 发表于 2014-2-13 15:09
好像只能提取单元应力吧,如果要提取节点的可能要fish程序编程,把周围单元应力相加再除以单元数。 ...
请问如何提取单元应力? NEW
rest step18.sav
def pcyl_outp
z=0
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d1=z_sig1(p_z) ;first stress
z=-18.10
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d2=z_sig1(p_z)
z=-34.33
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d3=z_sig1(p_z)
z=-48.87
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d4=z_sig1(p_z)
z=-61.89
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d5=z_sig1(p_z)
z=-69.82
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d6=z_sig1(p_z)
z=-77.19
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d7=z_sig1(p_z)
z=-84.04
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d8=z_sig1(p_z)
z=-90.40
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d9=z_sig1(p_z)
z=-91.99
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d10=z_sig1(p_z)
z=-94.66
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d11=z_sig1(p_z)
z=-97.31
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d12=z_sig1(p_z)
z=-99.98
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d13=z_sig1(p_z)
z=-102.94
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d14=z_sig1(p_z)
z=-105.89
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d15=z_sig1(p_z)
z=-108.85
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d16=z_sig1(p_z)
z=-111.80
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d17=z_sig1(p_z)
z=-114.76
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d18=z_sig1(p_z)
z=-117.75
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d19=z_sig1(p_z)
z=-120.74
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d20=z_sig1(p_z)
z=-123.69
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d21=z_sig1(p_z)
z=-126.64
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d22=z_sig1(p_z)
z=-129.60
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d23=z_sig1(p_z)
z=-132.55
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d24=z_sig1(p_z)
z=-135.50
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d25=z_sig1(p_z)
z=-138.17
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d26=z_sig1(p_z)
z=-140.83
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d27=z_sig1(p_z)
z=-143.50
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d28=z_sig1(p_z)
z=-145.10
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d29=z_sig1(p_z)
z=-151.47
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d30=z_sig1(p_z)
z=-158.32
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d31=z_sig1(p_z)
z=-165.69
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d32=z_sig1(p_z)
z=-173.62
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d33=z_sig1(p_z)
z=-182.15
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d34=z_sig1(p_z)
z=-196.16
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d35=z_sig1(p_z)
z=-217.41
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d36=z_sig1(p_z)
z=-235.50
x=-17.75
y=0
p_z=z_nearall(x,y,z)
d37=z_sig1(p_z)
end
pcyl_outp
set logfile westwallsmax.log
set log on
set pagelength 10000
print d1 d2d3 d4 d5 d6 d7 d8 d9 d10 d11 d12 d13 d14 &
d15 d16 d17 d18 d19 d20 d21 d22 d23 d24 d25 d26 d27 d28 &
d29 d30 d31 d32 d33 d34 d35 d36 d37
set log off
你可以参考一下
页:
[1]