model.set_constr('fix', ['P4', 'P7'], 'x')
# set part material
mat = pyc.Material('steel')
mat.set_mech_props(7800, 210 * (10**9), 0.3)
model.set_matl(mat, part)
# set the element type and mesh database
model.set_eshape(eshape, 2)
model.set_etype('plstress', part, 0.1)
model.mesh(1.0, 'gmsh') # mesh 1.0 fineness, smaller is finer
model.plot_elements(proj_name + '_elem', display=show_gui)
model.plot_pressures(proj_name + '_press', display=show_gui)
model.plot_constraints(proj_name + '_constr', display=show_gui)
# make and solve the model
prob = pyc.Problem(model, 'struct')
prob.solve()
# view and query results
sx = prob.rfile.get_nmax('Sx')
print('Sx_max: %f' % sx)
# Plot results
# store the fields to plot
fields = 'Sx,Sy,S1,S2,S3,Seqv,ux,uy,utot,ex'
fields = fields.split(',')
for field in fields:
fname = proj_name + '_' + field
prob.rfile.nplot(field, fname, display=False)
def example_load_results(self, file_name):
proj_name = os.path.join('examples', file_name)
model = pyc.FeaModel(proj_name)
prob = pyc.Problem(model, 'struct')
prob.rfile.load()
