FLAC2D:
def ps3d ; put 3D principal stresses
; into ex_1,2,3
loop i (1,izones)
loop j (1,jzones)
$sdif = sxx(i,j) - syy(i,j)
$s0 = 0.5 * (sxx(i,j) + syy(i,j))
$st = 4.0 * sxy(i,j)^2
$rad = 0.5 * sqrt($sdif^2 + $st)
$si = $s0 - $rad
$sii = $s0 + $rad
section
if szz(i,j) > $sii then
; ---- szz is minor p.s. ----
ex_1(i,j) = $si
ex_2(i,j) = $sii
ex_3(i,j) = szz(i,j)
exit section
end_if
if szz(i,j) < $si then
; ---- szz is major p.s. ----
ex_1(i,j) = szz(i,j)
ex_2(i,j) = $si
ex_3(i,j) = $sii
exit section
end_if
; ---- szz is intermediate ---
ex_1(i,j) = $si
ex_2(i,j) = szz(i,j)
ex_3(i,j) = $sii
end_section
end_loop
end_loop
end
Hi. In connection with my first question, I was wondering why the following functions went haywire! Is it any problem in the loop I defined?
def install
pnt = zone_head
loop while pnt # null
_geref = 11.6e6 ; reference tangent shear mod
_keref = 25e6 ; reference tangent bulk mod
_Pa = 100e3 ; reference pressure
m = 1
z_extra (pnt, 1) = _geref*(((z_sig1(pnt)+z_sig2(pnt)+z_sig3(pnt))/3)/_Pa)^m ; elastic shear modulus
_ge = z_extra (pnt, 1)
z_extra (pnt, 2) = _keref*(((z_sig1(pnt)+z_sig2(pnt)+z_sig3(pnt))/3)/_Pa)^m ; elastic bulk modulus
_ke = z_extra (pnt, 2)
pnt = z_next(pnt)
end_loop
end
@install
pro bu @_ke sh @_ge
- Badges
- Science topic
- Similar topics
- Physical Sciences
- Waves