Input Data
(→New Input Data) |
|||
| Line 247: | Line 247: | ||
974 0 1.000000 | 974 0 1.000000 | ||
End NodalData | End NodalData | ||
| + | |||
| + | |||
| + | |||
| + | [[Category:Kernel]] | ||
Revision as of 13:35, 1 October 2009
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
