Python - Boreholes and soil layers set using information from a csv file

Question

Afternoon, 
 
 I have a csv file with: 
 - boreholes (1, 2, 3,..., b) characteristics ('x', head, 1st layer top level, 1st layer bottom level, 2nd layer bottom level, 3rd layer bottom level, ..., lth layer bottom level); 
 - material (1, 2, 3, ..., m) characteristics (to be assigned to each soil layer that each borehole). 
 Whilst I managed to get the Python script reading the csv files, trying to set boreholes/layers/materials in loops has been more complex (I would like the script to be run for as many boreholes/layer/materials as needed). All examples I have seen, so far, are for basic ground models with horizontal layers defined by a single borehole. Is there a way to set as many boreholes as needed, change the top and bottom layers levels, head, &lsquo;x&rsquo; coordination in a loop? 
 From what I have seen g_i.borehole(x) creates Borehole_1, and if I repeat g_i.borehole(x) Borehole_2 is created, and so on. How can I access and change the characteristics of each borehole, using a Python script, for any number of boreholes >1? 
 How can layers characteristics be set for each borehole, and a material set for each of those layer, in a loop? 
 Apologies if this question has been asked before but, I couldn&rsquo;t find it in the forum. 
 
 Cheers, 
 Joao

Micha van der Sloot · Accepted Answer

Hi Joao, 
 In order to set borehole layers from a csv we would need borehole elevation levels, in the script I see that you used thickness data? Is that correct? 
 Q1 Borehole layers 
 For the borehole data you could use this: 
 
 I used a function to read the csv file 
 I also created a small function to just retrieve the soil layer elevation from the csv file (from the excel file). If you change the boolean to False when calling the function, your attached csv file should work too 
 and I used one command to set the soillayer levels in just one line (this was added recently as a supported format in PLAXIS 2D CONNECT Edition V20 Update 3.

import csv

def read_data ( csvfilepath ):
 with open ( csvfilepath , 'r' ) as csvfile :
 csvdata = csv . reader ( csvfile )
 csvrows = [ row for row in csvdata ]
 return csvrows

def getboreholelayerinput ( cvsline , toplayerindex , isThicknessData = True ):
 """ 
 creates a list with borehole layer elevations 
 :param cvsline: the line from a csv file, a list of values 
 :param toplayerindex: index for the top layer elevation 
 :param isThicknessData: boolean to specify if data is elevation or layer thickness 
 :return: a list of elevation values top-down 
 """ 
 layerlevels = []
 # when layer levels are all defined as elevation, use: 
 if isThicknessData is False :
 layerlevels = [ float ( i ) for i in cvsline [ toplayerindex :]]
 else :
 # when layer levels are defined as y_level_top, then thinkness values, use: 
 layerlevels = [ float ( cvsline [ toplayerindex ])] # create list of elevations. top down 
 for ilvl in range ( toplayerindex + 1, len ( cvsline )):
 if 'na' in cvsline [ ilvl ] or cvsline [ ilvl ] == '' :
 break 
 else :
 layerlevels . append ( layerlevels [ - 1] - float ( cvsline [ ilvl ]))
 return layerlevels

def createboreholes ( csvfile ):
 """ 
 creates boreholes in a blank project 
 :param csvfile: csv file with borehole data per line, has header line 
 :return: a list of generated boreholes 
 """ 
 csvdata = read_data ( csvfile )

# column definition from csv file 
 indexName = 1
 indexX = 3
 indexY = None # 3D only 
 indexHead = 2
 indexLayerLevels = 4 # column for yo_top_level

# boreholes 
 boreholes = []
 for line in csvdata [1:]:
 # create borehole 
 if indexY :
 boreholes . append ( g_i . borehole ( float ( line [ indexX ]), float ( line [ indexY ])))
 else :
 boreholes . append ( g_i . borehole ( float ( line [ indexX ])))
 # set name value 
 boreholes [ - 1] . Name = line [ indexName ]

# set head value 
 boreholes [ - 1] . Head = float ( line [ indexHead ])

#set soil layer levels 
 BHlayerlevels = getboreholelayerinput ( line , indexLayerLevels , isThicknessData = False )
 # make sure we have enough soillayers: adding soillayers for the first borehole 
 if len ( boreholes ) == 1:
 for layer in range ( len ( BHlayerlevels ) - 1):
 g_i . soillayer (0) # jsut add zero thickness layers initially, correct later 
 firstborehole = False 
 #set all soil layer levels per borehole: 
 g_i . setsoillayerlevel ( boreholes [ - 1], * [( i , BHlayerlevels [ i ]) for i in range ( len ( BHlayerlevels ))])

return boreholes

# start a new project and add boreholes 
 s_i . new ()
 createboreholes ( 'Boreholes_list.csv' )

Q2 Looping over polygons 
 It seems you are looping over all Polygon names in staged construction. In that case, the Polygon created in Structures mode has been intersected and now we have a cut object. If you use the name of the polygon, we could loop over all the polygon-entities that originate from that original name. And to set materials, you would also need to tell in which phase it needs to be done. Here is a small example 
 
 polygonnames = [ 'PolygonName1' , 'PolygonName2' ]
 phase = g_i . Phase_1 
 material = g_i . MaterialB

for pgname in polygonnames :
 if hasattr ( g_i , pgname ):
 pgobject = getattr ( g_i , pgname )
 for cutobject in pgobject :
 # deactivate polygon: 
 g_i . activate ( cutobject , phase )
 # change material 
 g_i . setmaterial ( cutobject , phase , material )

I hope these examples here help to get the script working. 
 Looks to be an interesting way to automate these model generations. 
 
 Kind regards, 
 Micha

Stefanos Papavasileiou · Answer

Hello Sean, 
 Typically we use the .BoundingBox of a polygon object to identify which is which. 
 However, we are working on some further developments (still improving) on the geometry information in Staged construction I can share with you. 
 Let's say you have Polygon_1_1 in Staged construction. 
 Running the following command will give you the segments (lines) that define the polygon: 
 echo Polygon_1_1.Geometry.Segments 
 And if you want to know the points of each part, you can run this: 
 echo Polygon_1_1.Geometry.Segments[0] 
 This will give you x1, y1, x2, y2 . 
 In Python naturally, you can do a for-loop to query all the values and remove any duplicates. 
 I hope that this helps.

Python - Boreholes and soil layers set using information from a csv file

Top Replies