|
| void | readInputFile (string file_name) |
| | Read the input file.
|
| |
| const json & | getJson () |
| | Return the data file structure.
|
| |
| string | getFileName () |
| | Return the file name.
|
| |
| bool | getLoadState () |
| | Return load state activated.
|
| |
| bool | getSaveState () |
| | Return save state activated.
|
| |
| double | getSimulationTime () |
| | Return the simulation time.
|
| |
| double | getTimeStep () |
| | Return the time step.
|
| |
| double | getCriticalTimeStepMultiplier () |
| | Return critical time step multiplier.
|
| |
| Solver * | getSolver () |
| | Return the solver to be used in the model.
|
| |
| ModelSetup::InterpolationFunctionType | getInterpolationFunction () |
| | Return the interpolation functions type.
|
| |
| Vector3i | getCellsNum () |
| | Return the number of cells in each direction.
|
| |
| Vector3d | getCellDimension () |
| | Return the cell dimension.
|
| |
| Vector3d | getOrigin () |
| | Return the origin of coordinates.
|
| |
| vector< Material * > | getMaterialList () |
| | Return the material list.
|
| |
| vector< Body * > | getBodyList () |
| | Return the body list.
|
| |
| Vector3d | getGravity () |
| | Return the gravity force.
|
| |
| int | getResultNum () |
| | Return the number of results.
|
| |
| vector< string > | getResultFields () |
| | Return the fields to be written.
|
| |
| ModelSetup::DampingType | getDampingType () |
| | Read the damping definition in the input file.
|
| |
| double | getDampingValue () |
| | Read the damping value.
|
| |
| vector< Boundary::BoundaryType > | getMeshBoundaryConditions () |
| | Return the mesh boundary conditions.
|
| |
| vector< Boundary::BoundaryType > | getMeshBoundaryConditionsFluid () |
| | Return the mesh boundary conditions of fluid.
|
| |
| unsigned | getNumThreads () |
| | Return number of threads defined in the input file.
|
| |
| unsigned | getNumberPhases () |
| | Return the number of phases in the simulation.
|
| |
| vector< Loads::LoadDistributedBox > | getLoadDistributedBox () |
| | Return loads distributed in particles inside a box.
|
| |
| Loads::NodalPointLoadData | readNodalPointLoads () |
| | read nodal point loads This function reads a points from input file in the following format: "nodal_point_loads": [[p1x,p1y,p1z],[p1x,p1y,p1z]],...,[[pnx,pny,pnz],[pnx,pny,pnz]]
|
| |
| vector< Loads::PressureBox > | getPrescribedPressureBox () |
| | Return prescribed pressure in particles inside a box.
|
| |
| vector< Loads::PressureBox > | getInitialPressureBox () |
| | Return initial pore pressure in particles inside a box.
|
| |
| vector< Loads::PressureMaterial > | getInitialPressureMaterial () |
| | Return initial pore pressure in particles by material id.
|
| |
| vector< Loads::PressureBoundaryForceBox > | getPressureBoundaryForceBox () |
| | Return pore pressure force in particles inside a box.
|
| |
| Loads::SeismicData | readSeismicData (const std::string &filename, bool hasHeader) |
| |
Operations to read the input file.
In this section topics related with the input data file are considered in details.
Definition
The keywords are the way to setup the numerical model. Have keyword for each characteristic implemented in the program. Some of they are associated with numbers, for example "time":0.2, others that are associated with strings, for example "type":"elastic", or to arrays, for example "gravity":[0.0,0.0,-9.81], and with others definitions, like "mesh":{ / "cells_dimension":[1,1,1], ... Next is presented the complete list of keyword for using in the input file, organized according the characteristic to be configured.
Bodies
| Keyword | Description | Data type |
| body | define bodies | – |
| type | define a type of body | string |
| id | define a body identification | integer |
| material_id | used to setup a defined material to the body | integer |
| initial_velocity | used to setup the initial velocity to the body | integer |
Particles
| Keyword | Description | Data type |
| particles | define particles | – |
| id | identification | integer |
| position | initial position | array |
| volume | initial volume | integer |
| length | initial size | integer |
Example of Body defined by particles
The next example shows the definition of body forming by one particle
"body":
{
"soil":
{
"type":"particles",
"id":1,
"material_id":1,
"particles":
{
"id":[1],
"position":[[0.5,0.5,0.5]],
"volume":[1],
"length":[1]
}
}
}
Cuboid
| Keyword | Description | Data type |
| cuboid | cuboid body type | body type |
| point_p1 | define the a lower left point in a cuboid | array |
| point_p2 | define the a higher right point in a cuboid | array |
Cuboid definition example
The next example shows the definition of a cuboid body type witch volume extends from point P1=(2,2,2) to the point P2=(7,7,7) and witch material has id=1.
"body":
{
"cuboid":
{
"type":"cuboid",
"id":1,
"point_p1":[2,2,2],
"point_p2":[7,7,7],
"material_id":1
}
}
Polygon extrusion
| Keyword | Description | Data type |
| polygon_2d | body created by a 2D polygon | body type |
| extrude_direction | direction of the extrusion during the polygon body type | string |
| extrude_displacement | magnitude of the displacement in the extrude direction | double |
| discretization_length | cell dimension of the auxiliary mesh during the creation of the polygon body type | double |
| points | 3D coordinate points | array |
Polygon extrusion definition example
"body":
{
"soil":
{
"type":"polygon_2d",
"extrude_direction":"y",
"extrude_displacement":2,
"discretization_length":0.5,
"id":1,
"points":
[
[0, 0, 0],
[2.001, 0, 0],
[0, 0, 2.001]
],
"material_id":1
}
}
Loads
| Keyword | Description | Data type |
| gravity | define the body force intensity | array |
| load_box | distribute a load in particles inside a box | – |
Load definition examples
Gravity load Example
"gravity":[0.0,0.0,-9.81],
Distributed load in box example
"load_distributed_box":
{
"load_1":
{
"point_p1":[0,0,6],
"point_p2":[5,1,10],
"load":[0,0,-10]
}
}
Pore pressure
| Keyword | Description | Data type |
| prescribed_pressure_box | prescribe a pore pressure in particles inside a box | – |
| initial_pressure_box | configure initial pore pressure in particles inside a box | – |
Pore pressure definition examples
Prescribed pore pressure in particles inside a box example
"prescribed_pressure_box":
{
"pressure_prescibed_1":
{
"point_p1":[0,0,6],
"point_p2":[5,1,10],
"pressure":10
}
}
Initial pore pressure in particles inside a box example
"initial_pressure_box":
{
"pressure_initial_1":
{
"point_p1":[0,0,6],
"point_p2":[5,1,10],
"pressure":10
}
}
Damping
| Keyword | Description | Data type |
| damping | define damping in the model | – |
| type | define damping type | string |
| value | define the local damping value | double |
See damping type in ModelSetup::DampingType
Damping definition example
"damping":
{
"type":"local",
"value":0.0
}
Parallelization
| Keyword | Description | Data type |
| n_threads | define the number of threads in the simulation | unsigned |
Parallelization definition example
Simulation time
| Keyword | Description | Data type |
| time | define the simulation time | double |
| time_step | define the time step | double |
| critical_time_step_multiplier | define the time step by a fraction of the critical time step | double |
Simulation time definition example
Time step definition example
The time step of the simulation can be defined by the critical time multiplier:
"critical_time_step_multiplier":0.25,
So, in this case the time step will be 0.25 of the critical time step.
Two phase simulations
| Keyword | Description | Data type |
| n_phases | define number of phases in the simulations (default is 1) | int |
Two phase simulation definition example
Stress update scheme
| Keyword | Description | Data type |
| stress_scheme_update | used to define the type of stress update scheme | string |
| USL | Update Stress Last scheme type | – |
Stress update scheme definition example
"stress_scheme_update":"USL",
Shape Functions
| Keyword | Description | Data type |
| shape_function | define the type of nodal interpolation function | string |
| GIMP | generalized interpolation material point shape function type | – |
| linear | linear shape function type | – |
Shape functions definition example
Mesh
| Keyword | Description | Data type |
| mesh | used to define a mesh | – |
| cells_number | define the number of cells in each direction | array |
| cells_dimension | used to define the cell dimension in each direction | array |
| origin | used to define the origin of coordinates | array |
| boundary_conditions | used to set the mesh boundary conditions | – |
| boundary_conditions_fluid | used to set the mesh boundary conditions of fluid | – |
| sliding | boundary condition that can move perpendicular to the normal of the boundary plane | string |
| fixed | boundary condition that can not move in any direction | string |
| free | boundary condition that can move freely in all directions | string |
| plane_X0 | Plane which normal points to the negative part of the axis X | – |
| plane_Y0 | Plane which normal points to the negative part of the axis Y | – |
| plane_Z0 | Plane which normal points to the negative part of the axis Z | – |
| plane_Xn | Plane which normal points to the positive part of the axis X | – |
| plane_Yn | Plane which normal points to the positive part of the axis Y | – |
| plane_Zn | Plane which normal points to the positive part of the axis Z | – |
Mesh definition example
"mesh":
{
"cells_dimension":[1,1,1],
"cells_number":[10,10,10],
"origin":[0,0,0],
"boundary_conditions":
{
"plane_X0":"fixed",
"plane_Y0":"fixed",
"plane_Z0":"fixed",
"plane_Xn":"sliding",
"plane_Yn":"sliding",
"plane_Zn":"free"
}
},
Materials
| Keyword | Description | Data type |
| material | used to define materials | – |
| id | material identification | integer |
| density | mass density | double |
| density_fluid | mass density of fluid in mixture | double |
| porosity | porosity of mixture | double |
| bulk_fluid | volumetric modulus of fluid in mixture | double |
| hydraulic_conductivity | principal values of the hydraulic conductivity of fluid in mixture | array (1x3) |
| type | used to specify the material constitutive model | string |
Linear Elastic
| Keyword | Description | Data type |
| elastic | elastic material type | – |
| poisson | Poisson's ratio of the material | double |
| young | Young's modulus of the material | double |
Mohr-Coulomb
| Keyword | Description | Data type |
| mohr-coulomb | Mohr-Coulomb elasto plastic material type | – |
| friction | Angle of internal friction in degree | double |
| cohesion | Material Cohesion | double |
| dilation | Material dilation | double |
| tensile | Tensile cut-off | double |
Material definition example
"material":
{
"plastic":
{
"type":"mohr-coulomb",
"id":1,
"young":50e4,
"density":2000,
"poisson":0.2,
"friction":30.0,
"cohesion":1500.0
}
},
Results
| Keyword | Description | Data type |
| results | used to define the type and the number of results | – |
| print | number of results to be written | integer |
| fields | used to define the fields to be written | string array |
| id | particle identification field | – |
| displacement | particle displacement field | – |
| material | particle material field | – |
| pressure | fluid pressure | – |
| external_force | total external force in particle | – |
| plastic_strain | effective plastic strain in particle | – |
| stress | particle stress | – |
Results definition example
"results":
{
"print":50,
"fields":["id","displacement","material"]
},
Complete input file definition example
{
"stress_scheme_update":"USL",
"shape_function":"GIMP",
"n_threads":4,
"time":2,
"time_step":0.001,
"gravity":[0.0,0.0,-9.81],
"results":
{
"print":50,
"fields":["id","displacement","material"]
},
"mesh":
{
"cells_dimension":[1,1,1],
"cells_number":[10,10,10],
"origin":[0,0,0],
"boundary_conditions":
{
"plane_X0":"fixed",
"plane_Y0":"fixed",
"plane_Z0":"fixed",
"plane_Xn":"sliding",
"plane_Yn":"sliding",
"plane_Zn":"free"
}
},
"material":
{
"elastic":
{
"type":"elastic",
"id":1,
"young":50e4,
"density":2000,
"poisson":0.2
}
},
"body":
{
"cuboid":
{
"type":"cuboid",
"id":1,
"point_p1":[2,2,2],
"point_p2":[7,7,7],
"material_id":1
}
}
}
Save/Load States
These keywords allow to load and save model states, for example this can be used to load the initial stresses resulting of a previous analysis
| Keyword | Description | Data type |
| load_state | activate load model state | bool |
| save_state | activate save model state | bool |
Example To save the state of the model
And for loading the salved state
1.12.0