Input Data

From KratosWiki
Revision as of 13:31, 1 October 2009 by JMora (Talk | contribs)
Jump to: navigation, search

The current input data consist of five archives with following extensions:

  • .node for the nodal coordinates
  • .prop which containes the elemental properties
  • .elem for elemental connectivities
  • .cond containing conditions connectivities and node's fixed dofs
  • .init for nodal initial values

All files are free format and the reading is not depend to the spaces, tabs, endlines etc.


Contents

Nodal Coordinates File

This file must have the .node extension with following format:

NODES = NodesList([
[NodeId,   X,   Y,   Z],
[NodeId,   X,   Y,   Z],
...
...
...
 [NodeId,   X,   Y,   Z],
  ])

for example:

NODES = NodesList([
[1,   0.05000,   0.00000,   0.00000],
[2,   0.05000,   0.01000,   0.00000],
...
...
...
[1000,   0.00000,   1.00000,   0.00000]
])

The NodeIds can be noconsecutive and also not in order, however the ordered input can improve the reading time.


Properties File

This file must have the .prop extension with following format:


PROPERTIES[PropertiesId][ScalarVariableName] =  Value
PROPERTIES[PropertiesId][VectorialVariableName] =  [Value1, Value2, Value3]
PROPERTIES[PropertiesId][MatricialVariableName] =  [[Value11, Value12, Value13],
                                                    [Value21, Value22, Value23],
                                                    [Value31, Value32, Value33]]

For example:

PROPERTIES[1][END_TIME] =   0.1;
PROPERTIES[1][DELTA_TIME] = 0.001;
PROPERTIES[2][DENSITY] = 2700.000000; 
PROPERTIES[2][YOUNG_MODULUS] = 7000000.000000;
PROPERTIES[2][GRAVITY] = [-39.5, 0.00, 0.00];


Elemental Connectivities File

This file must have the .elem extension and defines the elements with following format:

ELEMENTS[ElementId] = ElementName([NodeId1,NodeId2, ..,NodeIdn],PropertiesId);

For example:

ELEMENTS[1] = TotalLagrangian2D3N([1,2,3],1);
ELEMENTS[2] = TotalLagrangian2D3N([2,6,4],1);
ELEMENTS[3] = TotalLagrangian2D3N([6,11,8],1);
ELEMENTS[4] = TotalLagrangian2D3N([11,18,17],2);
ELEMENTS[5] = TotalLagrangian2D3N([18,26,24],2);
ELEMENTS[6] = TotalLagrangian2D3N([26,36,33],2);


Conditions File

This file must have the .cond extension and defines the fixed dofs in each node and boundary conditions. The fixed dofs can be define in following form:

NODES[NodeId].Fix(VariableName);

For example:

NODES[27].Fix(DISPLACEMENT_X);
NODES[27].Fix(DISPLACEMENT_T);
NODES[29].Fix(PRESSURE);

The boundary conditions can be defined similar to the elements by their connectivity as follow:

CONDITIONS[ConditionId] = ConditionName([NodeId1,NodeId2, ..,NodeIdn],PropertiesId);

For example:

CONDITIONS[1176] = Face2D([684,678],1);
CONDITIONS[1175] = Face2D([693,684],1);
CONDITIONS[1176] = Face2D([684,678],1);


Initial Values File

This file must have the .init extension and can be used to give intial values to the nodal historical data using following format:

NODES[NodeId](VariableName,0) = Value;

For example:

NODES[29](PRESSURE,0) = 100.000000;

New Input Data

A new data format for input data is under design. This format contains all ModelPart's data in a minimalistic block form and consists of 8 blocks of data in following form:

 Begin ModelPartData
 //  VARIABLE_NAME value
 End ModelPartData
 
 Begin Properties  properties_id 
 //  VARIABLE_NAME value
 End Properties
 
 Begin Nodes
 // id	  X	Y	Z
 End Nodes
 
 Begin Elements element_name
 // id prop_id	 n1	n2	n3	...
 End Elements
 
 Begin Conditions condition_name 
 // id prop_id	 n1	n2	n3	...
 End Conditions
 
 Begin NodalData VARIABLE_NAME
 //  id is_fixed value
 End NodalData
 
 Begin ElementalData VARIABLE_NAME
 //  id value
 End ElementalData
 
 Begin CondtionalData VARIABLE_NAME
 //  id value
 End CondtionalData

Each block starts with a Begin statement following by the block name and ends with the End statement again following by the block name. Some block may have some additional parameter like id or variable in their definitions.

An example is presented in the following lines:

Begin ModelPartData

//  VARIABLE_NAME value
End ModelPartData

Begin Properties 1
DENSITY 3.4E-5  //scalar
THICKNESS 19.5
ORTHOTROPIC_YOUNG_MODULUS [3] (20000,10000,8000) //vector
ORTHOTROPIC_POISSON_RATIO [3,3] ((0, 0.27,0.27),(0.087,0,0.27),(0.075,0.23,0)) // matrix
ORTHOTROPIC_SHEAR_MODULUS [3] (4225.23,4225.23,4225.23)
End Properties

Begin Nodes
1                  16                   0                   0
2                  16                 0.4                   0
3                15.6                   0                   0
...
972                   0                 7.2                   0
973                   0                 7.6                   0
974                   0                   8                   0
End Nodes

Begin Elements ASGS2D
1 0        649        661        641
2 0        765        771        785
3 0        710        700        704
...
1796 0        535        563        549
1797 0        877        857        891
1798 0        612        646        628
End Elements


Begin Conditions Condition2D
1799 0        644        650
1800 0        650        663
1801 0        663        673
...
1947 0        972        973
1948 0        973        974
End Conditions


Begin NodalData DISPLACEMENT_X
1 1 0.000000
2 1 0.000000
...
973 1 0.000000
974 1 0.000000
End NodalData

Begin NodalData DISPLACEMENT_Y
1 1 0.000000
2 1 0.000000
...
973 1 0.000000
974 1 0.000000
End NodalData

Begin NodalData DISPLACEMENT_Z
1 1 0.000000
2 1 0.000000
...
973 1 0.000000
974 1 0.000000
End NodalData

Begin NodalData VELOCITY_X
637 1 0.000000
639 1 0.000000
...
973 1 1.000000
974 1 1.000000
End NodalData

Begin NodalData VELOCITY_Y
1 1 0.000000
3 1 0.000000
...
973 1 0.000000
974 1 0.000000
End NodalData

Begin NodalData VELOCITY_Z
1 1 0.000000
3 1 0.000000
...
973 1 0.000000
974 1 0.000000
End NodalData

Begin NodalData VISCOSITY
1 0 0.010000
2 0 0.010000
...
973 0 0.010000
974 0 0.010000
End NodalData

Begin NodalData DENSITY
1 0 1.000000
2 0 1.000000
...
973 0 1.000000
974 0 1.000000
End NodalData
Personal tools
Categories