asdeasoft.net

Code: Select all

proc modal { num_modes filename1 {eig_solver -genBandArpack}} {

	# begin
	puts "\nRunning modal analysis ..."

	# get all node tags
	set nodes [getNodeTags]
	if {[llength $nodes] == 0} {
		error "modal - Error: no node in model"
	}

	# check problem size (2D or 3D) from the first node, we do not support mixed dimesions!!
	set ndm [llength [nodeCoord [lindex $nodes 0]]]

	# compute total masses
	if {$ndm == 3} { 
		set ndf_max 6 
		set total_mass {0.0 0.0 0.0 0.0 0.0 0.0}
		set mass_labels {"MX" "MY" "MZ" "MRX" "MRY" "MRZ"}
		set mass_labels1 {"MODE" "MX" "MY" "MZ" "MRX" "MRY" "MRZ"}
	} else {
		set ndf_max 3
		set total_mass {0.0 0.0 0.0}
		set mass_labels {"MX" "MY" "MRZ"}
		set mass_labels1 {"MODE" "MX" "MY" "MRZ"}
	}
	foreach node $nodes {
		set indf [llength [nodeDisp $node]]
		for {set i 0} {$i < $indf} {incr i} {
			set imass [nodeMass $node [expr $i+1]]
			set imass_total [lindex $total_mass $i]
			lset total_mass $i [expr $imass_total + $imass]
		}
	}

	# some constants
	set pi [expr acos(-1.0)]
	
	#modalDamping 0.02 0.02
	
	# solve the eigenvalue problem
	set lambdas [eigen $eig_solver $num_modes]
	if {[llength $lambdas] != $num_modes} {
		error "modal - Error: something went wrong in the eigen analysis"
	}

	#open a new document to record the values of the eigenvectors
	set fp1 [open $filename1 w]

	#set nodes_x {8} where we are interested in recording the eigenvectors TO BE EDITED!!!!
	set nodes_x {12 9 7 4 6 11}

	# process each mode of vibration
	for {set imode 0} {$imode < $num_modes} {incr imode} {

		# compute i-mode eigen analysis derived data
		set lambda [lindex $lambdas $imode]
		set omega [expr {sqrt($lambda)}]
		set omega [expr {sqrt($lambda)}]
		set frequency [expr $omega / 2.0 / $pi]
		set period [expr 1.0 / $frequency]

		foreach node $nodes_x {
			# get eigenvector
			puts $node
			set V [nodeEigenvector $node [expr $imode+1]]
			puts $fp1 "$node  $V"
		}

	}
	# done
	close $fp1
	puts "\nModal Analysis done\n"
}

#ask for the number of modes you want to extract
modal 4 "_eigenvectors.txt" "_frequencies.txt"
record