asdeasoft.net

Code: Select all

	# ___TCL___
		# Pause time if failures occur:
		if {[llength $failure_list] > 0} {
			setTime [expr [getTime] - $timeIncrement];
		};

Code: Select all

# ___Python API___
resultsField_u_Disp = db_u.getNodalResult("Displacement", match = MpcOdb.Contains)
# Check if node still exists in model:
for Analysis_Step_u in [1,2,3,4,5,6,7,8,9,10]:
	opt_u.step = Analysis_Step_u
	field = resultsField_u_Disp.evaluate(opt_u)
	row_u_does_exist = MpcOdbResultField.node(111000)
	row_u_does_not_exist = MpcOdbResultField.node(123456)
	x_disp = field[row_u_does_exist, 0]
	y_disp = field[row_u_does_exist, 1]
	print("When Node Does Exist in model:")
	print("x("+str(row_u_does_exist.n)+"): "+str(x_disp))
	print("y("+str(row_u_does_exist.n)+"): "+str(y_disp))
	print("When Node Does NOT Exist in model:")
	try:
		y_not_exist_disp = field[row_u_does_not_exist, 1]
		x_not_exist_disp = field[row_u_does_not_exist, 0]
	except RuntimeError as ex:
		template = "An exception of type {0} occurred. Arguments:\n{1!r}"
		message = template.format(type(ex).__name__, ex.args)
		print(message)
		print(traceback.format_exc())
		print("Node Selection "+str(row_u_does_not_exist.n)+" Not in Model")
		y_not_exist_disp = None
		x_not_exist_disp = None
		pass
	print("x("+str(row_u_does_not_exist.n)+"): "+str(x_not_exist_disp))
	print("y("+str(row_u_does_not_exist.n)+"): "+str(y_not_exist_disp))

Code: Select all

from PyMpc import *
from PyMpc import MpcOdbVirtualResult as vr

App.clearTerminal()
doc = App.postDocument()
db = doc.getDatabase(1)

# the node we want to check.
# we can check the presence of a node from the mesh in
# the result field
check_node = 134

# we need a result to generate a field, and to
# get the mesh from the field
result = db.getNodalResult('Displacement', match=MpcOdb.Contains)

# create the eval options
opt = MpcOdbVirtualResultEvaluationOptions()

# keep in mind that when you remove a node or an element,
# we have a new model stage (due to the domainChange method).
# so we don't need to evaluate the result
all_stages = db.getStageIDs()
for i in range(len(all_stages)):
	# the current stage_id
	stage_id = all_stages[i]
	# get all step id and time for this stage
	steps = list(db.getStepIDs(stage_id))
	times = list(db.getStepTimes(stage_id))
	# make sure we have at least 1 step
	if len(steps) > 0:
		# take the first one
		step_id = steps[0]
		time = times[0]
		# make field
		opt.stage = stage_id
		opt.step = step_id
		field = result.evaluate(opt)
		mesh = field.mesh
		exists = check_node in mesh.nodes
		print('Stage: {}, Start Time: {}, Start Step: {}, Exists ? {}'.format(
			stage_id, time, step_id, exists))

Code: Select all

from PyMpc import *
import os
import sys
import inspect
import traceback
import pdb
import pickle
import math
from os import chdir
from inspect import getmembers, isfunction
App.clearTerminal()
doc=App.postDocument()
MPCO_filename_u = "filename_u.mpco"
MPCO_filename_d = "filename_d.mpco"
export_dir = "export_dir_full_pathname"
if not os.path.exists(export_dir):
	os.makedirs(export_dir)
App.runCommand("OpenDatabase", MPCO_filename_u)
App.setBusy(False)
App.runCommand("OpenDataBase", MPCO_filename_d)
# Set loaded databases to variable names
db_u = doc.getDatabase(1)
db_d = doc.getDatabase(2)
# Get Model Stages
stages_u = db_u.getStageIDs() # type() = int_array
stages_d = db_d.getStageIDs() # type() = int_array
# Get Analysis Step IDs
current_stage_u = stages_u[0]
current_stage_d = stages_d[0]
analysisSteps_u_curr_stage = db_u.getStepIDs(current_stage_u)
analysisSteps_d_curr_stage = db_d.getStepIDs(current_stage_d)
# Set Options
opt_u = MpcOdbVirtualResultEvaluationOptions()
opt_d = MpcOdbVirtualResultEvaluationOptions()
opt_u.stage = current_stage_u
opt_d.stage = current_stage_d
# Get timeSteps for current model stages
timeSteps_u_curr_stage = db_u.getStepTimes(current_stage_u)
timeSteps_d_curr_stage = db_d.getStepTimes(current_stage_d) 

Code: Select all

result_field_u = db_u.getNodalResult('Displacement', match=MpcOdb.Contains)
result_field_d = db_d.getNodalResult('Displacement', match=MpcOdb.Contains)

Code: Select all

result_field_u = db_u.getNodalResult('Displacement', match=MpcOdb.Contains)

Code: Select all

result_field_u = db_u.getNodalResult('Displacement', match=3)

Code: Select all

from PyMpc import *
from PyMpc import MpcOdbVirtualResult as vr
import os
import sys
import inspect
import traceback
import pdb
import pickle
import math
from os import chdir
from read_txt_to_arr import read_txt_to_arr
from inspect import getmembers, isfunction
from interr_PyMpc_obj import interr_PyMpc_obj
App.clearTerminal() # Clear the STKO Terminal

# get document
doc = App.postDocument()
# open database
input_sim_no = 1000
undamaged_filename = "MPCO_1006_UND.mpco"
App.runCommand("OpenDatabase", undamaged_filename)

# get first database
if len(doc.databases) == 0:
	raise Exception("You need a database with ID = 1 for this test")
db = doc.getDatabase(1)

# create evaluation options
# here we want to extract data for all steps of the last stage
all_stages = db.getStageIDs()  # get all model stages
last_stage = all_stages[-1]  # get the last stage (pushover)
first_stage = all_stages[0] # get the first stage
all_stage1_steps = db.getStepIDs(first_stage)
all_stage1_times = db.getStepTimes(first_stage)
all_steps = db.getStepIDs(last_stage)  # get all steps of the last stage
all_times = db.getStepTimes(last_stage)  # get all step times of the last stage

step = all_stage1_steps[-1]
print(step)
InitialReactionForce = db.getNodalResult("Displacement", match=3)

App.processEvents()

# initialize the evaluation options with the last stage
#result_field_u = db.getNodalResult("Displacement", match = MpcOdb.Contains)
result_field_u = db.getNodalResult("Displacement", match = 3)