{"type":"doc","content":[{"type":"paragraph","content":[{"text":"This functions are based on the implementation of a problemtype for the ","type":"text"},{"text":"Frame3DD","type":"text","marks":[{"type":"textColor","attrs":{"color":"#0000ff"}},{"type":"link","attrs":{"href":"http://frame3dd.sourceforge.net"}}]},{"text":" ","type":"text","marks":[{"type":"link","attrs":{"href":"http://frame3dd.sourceforge.net"}}]},{"text":"solver, but can be useful to anyone.","type":"text"}]},{"type":"paragraph","content":[{"type":"hardBreak"},{"text":"Basic and advanced tcl commands for writing the input file - example","type":"text"}]},{"type":"bulletList","content":[{"type":"listItem","content":[{"type":"paragraph","content":[{"text":"How to ","type":"text"},{"text":"write a string","type":"text","marks":[{"type":"strong"}]},{"text":".","type":"text"}]}]}]},{"type":"paragraph","content":[{"text":"See how we print the units, we are asking GiD for the current active units for the different magnitudes: F for Force, L for Length, M for Mass","type":"text"}]},{"type":"codeBlock","attrs":{"language":"tcl"},"content":[{"text":"customlib::WriteString \"Input Data file for Frame3DD - 3D structural frame analysis ([gid_groups_conds::give_active_unit F],[gid_groups_conds::give_active_unit L],[gid_groups_conds::give_active_unit M])\"","type":"text"}]},{"type":"bulletList","content":[{"type":"listItem","content":[{"type":"paragraph","content":[{"text":"Nodes in the mesh","type":"text","marks":[{"type":"strong"}]}]}]}]},{"type":"paragraph","content":[{"text":"See how we get the number of nodes of the whole mesh. Also notice the way we are passing the writing format for the coordinates.","type":"text"}]},{"type":"codeBlock","attrs":{"language":"tcl"},"content":[{"text":"customlib::WriteString \"\"\ncustomlib::WriteString \"[GiD_Info Mesh NumNodes] # number of nodes\"\ncustomlib::WriteString \"#.node x y z r\"\ncustomlib::WriteString \"# [gid_groups_conds::give_active_unit L] [gid_groups_conds::give_active_unit L] [gid_groups_conds::give_active_unit L] [gid_groups_conds::give_active_unit L]\"\ncustomlib::WriteString \"\"\ncustomlib::WriteCoordinates \"%5d %14.5e %14.5e %14.5e 0.0\\n\"","type":"text"}]},{"type":"bulletList","content":[{"type":"listItem","content":[{"type":"paragraph","content":[{"text":"Elements","type":"text","marks":[{"type":"strong"}]},{"text":" assigned to the ","type":"text"},{"text":"condition","type":"text","marks":[{"type":"strong"}]},{"text":" that assigns a ","type":"text"},{"text":"material","type":"text","marks":[{"type":"strong"}]}]}]}]},{"type":"codeBlock","content":[{"text":"set condition_name \"frameData\"\n# If we want this output format \"element_id connectivity_x connectivity_y connectivity_z\"\nset condition_formats [list {\"%5d\" \"element\" \"id\"} {\"%5d\" \"element\" \"connectivities\"}]\n# If we want this output format with element id and properties\nset condition_formats [list {\"%1d\" \"element\" \"id\"} {\"%13.5e\" \"property\" \"Ax\"} {\"%13.5e\" \"property\" \"Asy\"} ... {\"%13.5e\" \"material\" \"Density\"}]\n# To get the number of elements\nset formats [customlib::GetElementsFormats $condition_name $condition_formats]\nset number_of_elements [GiD_WriteCalculationFile elements -count -elemtype Linear $formats]\n# To write the connectivities\ncustomlib::WriteConnectivities $condition_name $condition_formats \"\" active","type":"text"}]},{"type":"bulletList","content":[{"type":"listItem","content":[{"type":"paragraph","content":[{"text":"Constraints","type":"text","marks":[{"type":"strong"}]}]}]}]},{"type":"paragraph","content":[{"text":"This is the way we get the number of nodes assigned to a condition called 'constraints' in the spd.","type":"text"}]},{"type":"codeBlock","attrs":{"language":"tcl"},"content":[{"text":"set constraints_list [list \"constraints\"]\nset number_of_constraints [customlib::GetNumberOfNodes $constraints_list]","type":"text"}]},{"type":"bulletList","content":[{"type":"listItem","content":[{"type":"paragraph","content":[{"text":"Basic data","type":"text","marks":[{"type":"strong"}]}]}]}]},{"type":"paragraph","content":[{"text":"This is the way we access to the data in the tree. Notice that we are building what we call xpath, that is a kind of path in the spd, from the root to the item we want.","type":"text"},{"type":"hardBreak"},{"text":"Then we get the xml node from the document using selectNodes, and finally we get the current value with get_domnode_attribute","type":"text"}]},{"type":"codeBlock","attrs":{"language":"tcl"},"content":[{"text":"set document [$::gid_groups_conds::doc documentElement] \n set xpath \"/frame3dd_data/container\\[@n = 'extraOptions' \\]/value\\[@n = 'shear_deformation' \\]\"\n set xml_node [$document selectNodes $xpath]\n set shear_deformation [get_domnode_attribute $xml_node v]\n customlib::WriteString \"$shear_deformation # 1=Do, 0=Don't include shear deformation effects\"","type":"text"}]},{"type":"bulletList","content":[{"type":"listItem","content":[{"type":"paragraph","content":[{"text":"Single load case","type":"text","marks":[{"type":"strong"}]}]}]}]},{"type":"paragraph","content":[{"text":"This is the way we get the number of elements assigned to a condition, and then print the element id and the properties assigned to it. ","type":"text"},{"type":"hardBreak"},{"text":"set condition_name \"line_Uniform_load\"","type":"text"}]},{"type":"codeBlock","attrs":{"language":"tcl"},"content":[{"text":"set condition_name \"line_Uniform_load\"\n set condition_formats [list {\"%1d\" \"element\" \"id\"} {\"%13.5e\" \"property\" \"Load_x\"} {\"%13.5e\" \"property\" \"Load_y\"} {\"%13.5e\" \"property\" \"Load_z\"}]\n set formats [customlib::GetElementsFormats $condition_name $condition_formats]\n set number_of_elements [GiD_WriteCalculationFile elements -count -elemtype Linear $formats]\n customlib::WriteConnectivities $condition_name $formats \"\" active","type":"text"}]},{"type":"bulletList","content":[{"type":"listItem","content":[{"type":"paragraph","content":[{"text":"Several load cases","type":"text","marks":[{"type":"strong"}]}]}]}]},{"type":"paragraph","content":[{"text":"For the several load case, we will need to iterate over the xml nodes corresponding to the blockdatas that define each load case.","type":"text"},{"type":"hardBreak"},{"text":"For each load case, we need to print all the loads information, so we need to know if a load has any group assigned. ","type":"text"},{"type":"hardBreak"},{"text":"If so, we'll need to know how many elements are involved, and print them with the properties of the load.","type":"text"}]},{"type":"codeBlock","attrs":{"language":"tcl"},"content":[{"text":"set xpath \"/frame3dd_data/container\\[@n = 'staticCases' \\]/blockdata\"\n set xml_nodes [$document selectNodes $xpath]\n foreach load_case $xml_nodes {\n set xpath \"./condition\\[@n = 'uniformLoad' \\]/group\"\n set groups [$load_case selectNodes $xpath]\n set number_of_elements 0\n set formats_dict [dict create ]\n foreach group $groups {\n set group_name [get_domnode_attribute $group n]\n set Ux_node [$group selectNodes \"./value\\[@n = 'Ux'\\]\"]\n set Ux_value [get_domnode_attribute $Ux_node v]\n set Uy_node [$group selectNodes \"./value\\[@n = 'Uy'\\]\"]\n set Uy_value [get_domnode_attribute $Uy_node v]\n set Uz_node [$group selectNodes \"./value\\[@n = 'Uz'\\]\"]\n set Uz_value [get_domnode_attribute $Uz_node v]\n set format \"%5d $Ux_value $Uy_value $Uz_value\"\n set formats_dict [dict merge $formats_dict [dict create $group_name $format]]\n }\n set number_of_elements [GiD_WriteCalculationFile elements -count -elemtype Linear $formats_dict]\n if {$number_of_elements > 0} {\n GiD_WriteCalculationFile elements -elemtype Linear $formats_dict\n }\n }","type":"text"}]},{"type":"heading","attrs":{"level":2},"content":[{"text":"Alternatives","type":"text"}]},{"type":"paragraph","content":[{"text":"In this section, we will write the coordinates, connectivites, and conditions assuming that the input can be different, such as the one in ","type":"text"},{"text":"Matfem","type":"text","marks":[{"type":"link","attrs":{"href":"http://www.cimne.com/mat-fem/"}}]},{"text":" . In some codes, like matfem, the input file into the solver is not written in a text file, but in a \"code like file\". The coordinates section in the matfem solver is a matlab array, with a single comma splitting x and y (x, y) and a semicolon splitting a node from another (x1, y1; x2, y2; ...)","type":"text"}]},{"type":"bulletList","content":[{"type":"listItem","content":[{"type":"paragraph","content":[{"text":"Write coordinates:","type":"text"},{"type":"hardBreak"},{"text":"Example:","type":"text"}]},{"type":"codeBlock","attrs":{"language":"matlab"},"content":[{"text":"%\n% Coordinates\n%\nglobal coordinates\ncoordinates = [\n 0.00 , 0.00;\n 0.50 , 0.00;\n 2.50 , 1.00 ];","type":"text"}]},{"type":"paragraph","content":[{"type":"hardBreak"},{"text":"Code:","type":"text"}]},{"type":"codeBlock","attrs":{"language":"tcl"},"content":[{"text":" customlib::WriteString \"%\"\n customlib::WriteString \"% Coordinates\"\n customlib::WriteString \"%\"\n customlib::WriteString \"global coordinates\"\n set nodes [GiD_Mesh list node]\n customlib::WriteString \"coordinates = \\[\"\n for {set i 0} {$i < [llength $nodes]} {incr i} {\n set node [lindex $nodes $i]\n lassign [GiD_Mesh get node $node coordinates] x y z\n if {$i < [expr [llength $nodes] -1] } {set end \";\"} {set end \"\"}\n customlib::WriteString \"$x , $y $end\"\n\n }\n customlib::WriteString \"\\] ; \"","type":"text"}]}]}]},{"type":"paragraph"},{"type":"bulletList","content":[{"type":"listItem","content":[{"type":"paragraph","content":[{"text":"Write elements:","type":"text"},{"type":"hardBreak"},{"text":"Example:","type":"text"}]},{"type":"codeBlock","attrs":{"language":"matlab"},"content":[{"text":"%\n% Elements\n%\nglobal elements\nelements = [\n 1, 2, 7 ;\n 2, 3, 8 ;\n 18, 17, 12 ];","type":"text"}]},{"type":"paragraph","content":[{"type":"hardBreak"},{"text":"Code:","type":"text"}]},{"type":"codeBlock","attrs":{"language":"tcl"},"content":[{"text":" customlib::WriteString \"%\"\n customlib::WriteString \"% Elements\"\n customlib::WriteString \"%\"\n customlib::WriteString \"global elements\"\n set elements_t [GiD_Mesh list -element_type Triangle element]\n set elements_q [GiD_Mesh list -element_type Quadrilateral element]\n customlib::WriteString \"elements = \\[\"\n for {set i 0} {$i < [llength $elements_t] } {incr i} {\n set element [lindex $elements_t $i]\n lassign [GiD_Mesh get element $element connectivities] c1 c2 c3\n if {$i < [expr [llength $elements_t] -1] } {set end \";\"} {set end \"\"}\n customlib::WriteString \"$c1 , $c2 , $c3 $end\"\n }\n for {set i 0} {$i < [llength $elements_q] } {incr i} {\n set element [lindex $elements_q $i]\n lassign [GiD_Mesh get element $element connectivities] c1 c2 c3 c4\n if {$i < [expr [llength $elements_q] -1] } {set end \";\"} {set end \"\"}\n customlib::WriteString \"$c1 , $c2 , $c3 , $c4 $end\"\n }\n customlib::WriteString \"\\] ; \"","type":"text"}]}]}]},{"type":"paragraph"},{"type":"bulletList","content":[{"type":"listItem","content":[{"type":"paragraph","content":[{"text":"Conditions","type":"text"},{"type":"hardBreak"},{"text":"Example:","type":"text"}]},{"type":"codeBlock","attrs":{"language":"matlab"},"content":[{"text":"%\n% Point loads\n%\npointload = [\n6, 2, -1.0 ;\n12, 2, -1.0 ;\n18, 2, -1.0 ];","type":"text"}]},{"type":"paragraph","content":[{"type":"hardBreak"},{"text":"Code:","type":"text"},{"type":"hardBreak"},{"text":"spd","type":"text"}]},{"type":"codeBlock","attrs":{"language":"xml"},"content":[{"text":"\n \n \n","type":"text"}]}]}]},{"type":"paragraph","marks":[{"type":"indentation","attrs":{"level":1}}],"content":[{"text":"tcl","type":"text"}]},{"type":"codeBlock","attrs":{"language":"tcl"},"content":[{"text":"# Point load\n set root [customlib::GetBaseRoot]\n customlib::WriteString \"\"\n customlib::WriteString \"%\"\n customlib::WriteString \"% Point loads\"\n customlib::WriteString \"%\"\n customlib::WriteString \"pointload = \\[\"\n set displacement_fix_nodes [$root selectNodes \"*/condition\\[@n='PuntualLoads'\\]/group\"]\n foreach node $displacement_fix_nodes {\n set group [$node @n]\n set val_x [get_domnode_attribute [$node selectNodes \"./value\\[@n='x-force'\\]\"] v]\n set val_y [get_domnode_attribute [$node selectNodes \"./value\\[@n='y-force'\\]\"] v]\n set fix_x [expr $val_x == 0.0 ? \"false\" : \"true\"]\n set fix_y [expr $val_y == 0.0 ? \"false\" : \"true\"]\n set nodes [GiD_EntitiesGroups get $group nodes]\n set num_nodes [objarray length $nodes]\n\n for {set i 0} {$i < $num_nodes} {incr i} {\n set node_id [objarray get $nodes $i]\n if {$i < [expr $num_nodes -1] } {set end \";\"} {set end \"\"}\n if {$fix_x eq \"true\" && $fix_y eq \"true\"} {set end \";\" }\n if {$fix_x eq \"true\"} { customlib::WriteString \"$node_id , 1 , $val_x $end\" }\n if {$i < [expr $num_nodes -1] } {set end \";\"} {set end \"\"}\n if {$fix_y eq \"true\"} { customlib::WriteString \"$node_id , 2 , $val_y $end\" }\n }\n }\n customlib::WriteString \"\\] ; \"","type":"text"}]},{"type":"paragraph"}],"version":1}

Browser not supported