asdeasoft.net

Code: Select all

# get all nodes
set all_nodes [getNodeTags]
# find the largest dof for memory allocation
set max_dof 0
foreach node $all_nodes {
	set dofs [nodeDOFs $node]
	foreach dof $dofs {
		set max_dof [expr max($max_dof, $dof)]
	}
}
# allocate U
set U [lrepeat [expr $max_dof + 1] 0.0]
# fill U
foreach node $all_nodes {
	# get nodal dofs
	set dofs [nodeDOFs $node]
	# get nodal displacements
	set disp [nodeDisp $node]
	# process each dof
	for {set i 0} {$i < [llength $disp]} {incr i} {
		set idof [lindex $dofs $i]
		# consider only retained dofs (constrained dofs = -1)
		if {$idof >= 0} {
			# get dof displacement
			set idisp [expr [lindex $disp $i]]
			# put current dof displacement in the global vector
			lset U $idof $idisp
		}
	}
}
# print U
puts $U