fracture.geomp20
Syntax
- fp20 = fracture.geomp20(gset<,dfn>) (3D only)
Get the P20 of polygons. The P20 is the number of fracture intersections with the geometry polygons per unit area of the geometry polygons. If dfn is not given the calculation is undertaken over all fractures of all DFNs.
Returns: fp20 - average P20 of geometry set edges
Arguments: gset - geometry set pointer
dfn - optional DFN pointer
Usage Example
The following example illustrates how fracture.geomp20
can be
used to get the P20 on an outcrop for various DFNs, or for all fractures.
model new
model random 10000
model domain extent -5 5
;Create outcrop geometry
geometry set 'outcrop'
geometry polygon create by-position (-5,-5,0) (-5,5,0) (5,5,0) (5,-5,0)
;Generate fractures using the default template
fracture generate dfn 'p1' fracture-count 100
fracture generate dfn 'p2' fracture-count 200
fracture generate dfn 'p3' fracture-count 300
fracture generate dfn 'p4' fracture-count 400
; fish access to densities
fish define access_p20
local geom = geom.set.find('outcrop')
local dfn1 = dfn.find('p1')
local dfn2 = dfn.find('p2')
local dfn3 = dfn.find('p3')
local dfn4 = dfn.find('p4')
io.out('P20 of dfn p1 = '+string(fracture.geomp20(geom,dfn1)))
io.out('P20 of dfn p2 = '+string(fracture.geomp20(geom,dfn2)))
io.out('P20 of dfn p3 = '+string(fracture.geomp20(geom,dfn3)))
io.out('P20 of dfn p4 = '+string(fracture.geomp20(geom,dfn4)))
io.out('P20 of all dfns = '+string(fracture.geomp20(geom)))
end
@access_p20
return
Was this helpful? ... | PFC 6.0 © 2019, Itasca | Updated: Nov 19, 2021 |