Ziegler, Moritz O.; Ziebarth, Malte; Reiter, Karsten (2019): Python Script Apple PY. V. 1.02. GFZ Data Services. https://doi.org/10.5880/wsm.2019.001
Status
I N R E V I E W : Ziegler, Moritz O.; Ziebarth, Malte; Reiter, Karsten (2019): Python Script Apple PY. V. 1.02. GFZ Data Services. https://doi.org/10.5880/wsm.2019.001
In geosciences the discretization of complex 3D model volumes into finite elements can be a time-consuming task and often needs experience with a professional software. Especially outcropping or out-pinching geological units, i.e. geological layers that are represented in the model volume, pose serious challenges. Changes in the geometry of a model may occur well into a project at a point, when re-meshing is not an option anymore or would involve a significant amount of additional time to invest.
In order to speed up and automate the process of discretization, Apple PY (Automatic Portioning Preventing Lengthy manual Element assignment for PYthon) separates the process of mesh-generation and unit assignment. It requires an existing uniform mesh together with separate information on the depths of the interfaces between geological units (herein called horizons). These two pieces of information are combined and used to assign the individual elements to different units. The uniform mesh is created with a standard meshing software and contains no or only very few and simple structures. The mesh has to be available as an Abaqus input file. The information on the horizons depths and lateral variations in the depths is provided in a text file. Apple PY compares the element location and depth with that of the horizons in order to assign each element to a corresponding geological unit below or above a certain horizon.
Additional Information
Version History:
Version 1.02 (8 September 2020):
Bug fixes - no change in functionality. Manual for Version 1.0 remains valid except for certain cases where only single units should be processed. Then the added workaround should be used.
- If you are working with only one unit/element set in the variable strata or elems_exclude use either: squared brackets, e.g. elems_exclude = ['disregarded_unit'] or append a comma after the string, e.g. elems_exclude = ('disegarded_unit',)
Version 1.01 (29 August 2019):
Bug fixes - no change in functionality. Manual for Version 1.0 remains valid
- elems_exclude works now as designed and described in the manual.
- commenting out elems_exclude does not crash the script anymore.
- create_horizon_file does not create two instances of the uppermost horizon.
Authors
Ziegler, Moritz O.;GFZ German Research Centre for Geosciences, Potsdam, Germany
Ziebarth, Malte;GFZ German Research Centre for Geosciences, Potsdam, Germany
Reiter, Karsten;Technical University Darmstadt, Darmstadt, Germany
Contact
Ziegler, Moritz O.; GFZ German Research Centre for Geosciences, Potsdam, Germany;
Funders
Bundesministerium für Bildung und Forschung:
SUBI (FONA/GEO:N) (03G0869A)
Helmholtz-Gemeinschaft:
iCross
Keywords
geomechanical-numerical model, stress, in-situ stress, modelling tool, finite-elemeent model
affiliation (affiliationIdentifier=0000-0002-2288-2820 affiliationIdentifierScheme=ORCID): GFZ German Research Centre for Geosciences, Potsdam, Germany
affiliation (affiliationIdentifier=0000-0002-5190-4478 affiliationIdentifierScheme=ORCID): GFZ German Research Centre for Geosciences, Potsdam, Germany
affiliation (affiliationIdentifier=0000-0002-2288-2820 affiliationIdentifierScheme=ORCID): GFZ German Research Centre for Geosciences, Potsdam, Germany
contributor (contributorType=Researcher)
contributorName (nameType=Personal): Ziegler, Moritz O.
affiliation (affiliationIdentifier=0000-0002-2288-2820 affiliationIdentifierScheme=ORCID): GFZ German Research Centre for Geosciences, Potsdam, Germany
affiliation (affiliationIdentifier=0000-0002-5190-4478 affiliationIdentifierScheme=ORCID): GFZ German Research Centre for Geosciences, Potsdam, Germany
CharacterString: In geosciences the discretization of complex 3D model volumes into finite elements can be a time-consuming task and often needs experience with a professional software. Especially outcropping or out-pinching geological units, i.e. geological layers that are represented in the model volume, pose serious challenges. Changes in the geometry of a model may occur well into a project at a point, when re-meshing is not an option anymore or would involve a significant amount of additional time to invest.
In order to speed up and automate the process of discretization, Apple PY (Automatic Portioning Preventing Lengthy manual Element assignment for PYthon) separates the process of mesh-generation and unit assignment. It requires an existing uniform mesh together with separate information on the depths of the interfaces between geological units (herein called horizons). These two pieces of information are combined and used to assign the individual elements to different units. The uniform mesh is created with a standard meshing software and contains no or only very few and simple structures. The mesh has to be available as an Abaqus input file. The information on the horizons depths and lateral variations in the depths is provided in a text file. Apple PY compares the element location and depth with that of the horizons in order to assign each element to a corresponding geological unit below or above a certain horizon.
CI_OnLineFunctionCode (codeList=http://www.isotc211.org/2005/resources/Codelist/gmxCodelists.xml#CI_OnLineFunctionCode codeListValue=http://www.isotc211.org/2005/resources/Codelist/gmxCodelists.xml#CI_OnLineFunctionCode_information): information