flac 3d中如何提取节点主应力

hxm625344485

flac 3d中如何提取节点主应力？望高手指点

yangxin546

好像只能提取单元应力吧，如果要提取节点的可能要fish程序编程，把周围单元应力相加再除以单元数。

Donahue

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 d2  d3 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         

你可以参考一下