old_zu = None # The volumes of the elements in the filter mesh filtr_volumes = None # Create the array of properties rho = [2600.0] #, 1300.0, 600.0] E = [70e9] #, 35e9, 10e9] nu = [0.3] #, 0.3, 0.3] vars_per_node = 1 # Set the values of the objective array obj_array = [1.0e-4] # Create the stiffness properties object props = TMR.StiffnessProperties(rho, E, nu) for ite in range(max_iterations): # Create the TACSAssembler and TMRTopoProblem instance nlevs = mg_levels[ite] assembler, problem, filtr, varmap = createTopoProblem(props, forest, nlevels=nlevs, Xscale=Xscale) # Write out just the mesh - for visualization flag = TACS.ToFH5.NODES f5 = TACS.ToFH5(assembler, TACS.PY_SOLID, flag) f5.writeToFile(os.path.join(args.prefix, 'beam_mesh%d.f5' % (ite))) # Set the constraint type
# The volumes of the elements in the filter mesh filtr_volumes = None # Create the array of properties rho = [2600.0] #, 1300.0, 600.0] E = [70e9] #, 35e9, 10e9] nu = [0.3] #, 0.3, 0.3] vars_per_node = 1 # Set the values of the objective array obj_array = [1.0e2] # Create the stiffness properties object use_project = args.use_project props = TMR.StiffnessProperties(rho, E, nu, k0=1e-3, use_project=use_project) for ite in range(max_iterations): # Create the TACSAssembler and TMRTopoProblem instance nlevs = mg_levels[ite] assembler, problem, filtr, varmap = createTopoProblem( props, forest, nlevels=nlevs, Xscale=Xscale, ordering=TACS.PY_MULTICOLOR_ORDER) # Write out just the mesh - for visualization flag = TACS.ToFH5.NODES f5 = TACS.ToFH5(assembler, TACS.PY_SOLID, flag) f5.writeToFile(os.path.join(args.prefix, 'beam_mesh%d.f5' % (ite)))
prefix = args.prefix optimization_options['output_file'] = os.path.join(prefix, 'output_file.dat') optimization_options['tr_output_file'] = os.path.join(prefix, 'tr_output_file.dat') # Set the max nmber of outer iterations # (for continuation strategies and/or mesh refinement) mg_levels = args.mg_levels max_iterations = len(mg_levels) # Set the material properties rho = [2600.0] E = [70e9] nu = [0.3] # Create the stiffness properties object props = TMR.StiffnessProperties(rho, E, nu, k0=1e-3, eps=0.2, q=5.0) # Set the values of the objective array obj_array = [ 1.0e2 ] # Set the boundary conditions for the problem bcs = TMR.BoundaryConditions() bcs.addBoundaryCondition('fixed', [0, 1, 2], [0.0, 0.0, 0.0]) # Create the initial forest forest = create_forest(comm, args.init_depth, args.htarget) # Set the original filter to NULL orig_filter = None xopt = None