def run(): fets_eval = FETS2D4Q( mats_eval = MATS2DElastic() ) # Discretization fe_grid = FEGrid( coord_max = ( 10., 4., 0. ), shape = ( 10, 3 ), fets_eval = fets_eval ) bcg = BCDofGroup( var = 'u', value = 0., dims = [0], get_dof_method = fe_grid.get_left_dofs ) bcg.setup( None ) print 'labels', bcg._get_labels() print 'points', bcg._get_mvpoints() mf = MFnLineArray( #xdata = arange(10), ydata = array( [0, 1, 2, 3] ) ) right_dof = 2 tstepper = TS( sdomain = fe_grid, bcond_list = [ BCDofGroup( var = 'u', value = 0., dims = [0, 1], get_dof_method = fe_grid.get_left_dofs ), BCDofGroup( var = 'u', value = .005, dims = [1], time_function = mf.get_value, get_dof_method = fe_grid.get_right_dofs ) ], ) # Add the time-loop control tloop = TLoop( tstepper = tstepper, KMAX = 300, tolerance = 1e-4, tline = TLine( min = 0.0, step = 1.0, max = 1.0 ) ) U_k = tloop.eval() print 'dir', tloop.rtrace_mngr.dir # RTrace should not contain backward link to RTraceMngr # The definition should be forward. # rt1 = RTraceGraph( sd = tstepper.sdomain, rmgr = tstepper.rtrace_mngr, name = 'Fi,right over u_right (iteration)' , var_y = 'F_int', idx_y = right_dof, var_x = 'U_k', idx_x = right_dof, record_on = 'update' ) print 'dir', rt1.dir
def example_with_new_domain(): from ibvpy.api import \ TStepper as TS, RTraceGraph, RTraceDomainListField, \ RTraceDomainListInteg, TLoop, \ TLine, BCDof, IBVPSolve as IS, DOTSEval from ibvpy.mats.mats2D.mats2D_elastic.mats2D_elastic import MATS2DElastic from ibvpy.mats.mats2D.mats2D_sdamage.mats2D_sdamage import MATS2DScalarDamage from ibvpy.api import BCDofGroup mats_eval = MATS2DElastic() fets_eval = FETS2D4Q(mats_eval = mats_eval) #fets_eval = FETS2D4Q(mats_eval = MATS2DScalarDamage()) print fets_eval.vtk_node_cell_data from ibvpy.mesh.fe_grid import FEGrid from ibvpy.mesh.fe_refinement_grid import FERefinementGrid from ibvpy.mesh.fe_domain import FEDomain from mathkit.mfn import MFnLineArray # Discretization fe_grid = FEGrid(coord_max = (10., 4., 0.), shape = (10, 3), fets_eval = fets_eval) bcg = BCDofGroup(var = 'u', value = 0., dims = [0], get_dof_method = fe_grid.get_left_dofs) bcg.setup(None) print 'labels', bcg._get_labels() print 'points', bcg._get_mvpoints() mf = MFnLineArray(#xdata = arange(10), ydata = array([0, 1, 2, 3])) right_dof = 2 tstepper = TS(sdomain = fe_grid, bcond_list = [ BCDofGroup(var = 'u', value = 0., dims = [0, 1], get_dof_method = fe_grid.get_left_dofs), # BCDofGroup( var='u', value = 0., dims = [1], # get_dof_method = fe_grid.get_bottom_dofs ), BCDofGroup(var = 'u', value = .005, dims = [1], time_function = mf.get_value, get_dof_method = fe_grid.get_right_dofs) ], rtrace_list = [ RTraceGraph(name = 'Fi,right over u_right (iteration)' , var_y = 'F_int', idx_y = right_dof, var_x = 'U_k', idx_x = right_dof, record_on = 'update'), RTraceDomainListField(name = 'Stress' , var = 'sig_app', idx = 0, position = 'int_pnts', record_on = 'update'), # RTraceDomainListField(name = 'Damage' , # var = 'omega', idx = 0, # # record_on = 'update', # warp = True), RTraceDomainListField(name = 'Displacement' , var = 'u', idx = 0, record_on = 'update', warp = True), RTraceDomainListField(name = 'Strain energy' , var = 'strain_energy', idx = 0, record_on = 'update', warp = False), RTraceDomainListInteg(name = 'Integ strain energy' , var = 'strain_energy', idx = 0, record_on = 'update', warp = False), # RTraceDomainListField(name = 'N0' , # var = 'N_mtx', idx = 0, # record_on = 'update') ] ) # Add the time-loop control #global tloop tloop = TLoop(tstepper = tstepper, KMAX = 300, tolerance = 1e-4, tline = TLine(min = 0.0, step = 1.0, max = 1.0)) #import cProfile #cProfile.run('tloop.eval()', 'tloop_prof' ) #print tloop.eval() #import pstats #p = pstats.Stats('tloop_prof') #p.strip_dirs() #print 'cumulative' #p.sort_stats('cumulative').print_stats(20) #print 'time' #p.sort_stats('time').print_stats(20) tloop.eval() # Put the whole thing into the simulation-framework to map the # individual pieces of definition into the user interface. # from ibvpy.plugins.ibvpy_app import IBVPyApp app = IBVPyApp(ibv_resource = tloop) app.main()
def peval(self): '''Evaluation procedure. ''' #mv = MATS1DDamageView( model = mats_eval ) #mv.configure_traits() right_dof_m = self.fe_grid_m[-1, -1].dofs[0, 0, 0] right_dof_fb = self.fe_grid_fb[-1, -1, -1, -1].dofs[0, 0, 0] # Response tracers A = self.A_m + self.A_f self.sig_eps_m = RTraceGraph(name='F_u_m', var_y='F_int', idx_y=right_dof_m, var_x='U_k', idx_x=right_dof_m, transform_y='y / %g' % A) # Response tracers self.sig_eps_f = RTraceGraph(name='F_u_f', var_y='F_int', idx_y=right_dof_fb, var_x='U_k', idx_x=right_dof_fb, transform_y='y / %g' % A) self.eps_m_field = RTraceDomainListField(name='eps_m', position='int_pnts', var='eps_app', warp=False) self.eps_f_field = RTraceDomainListField(name='eps_f', position='int_pnts', var='mats_phase2_eps_app', warp=False) # Response tracers self.sig_m_field = RTraceDomainListField(name='sig_m', position='int_pnts', var='sig_app') self.sig_f_field = RTraceDomainListField(name='sig_f', position='int_pnts', var='mats_phase2_sig_app') self.omega_m_field = RTraceDomainListField(name='omega_m', position='int_pnts', var='omega', warp=False) self.shear_flow_field = RTraceDomainListField(name='shear flow', position='int_pnts', var='shear_flow', warp=False) self.slip_field = RTraceDomainListField(name='slip', position='int_pnts', var='slip', warp=False) avg_processor = None if self.avg_radius > 0.0: n_dofs = self.fe_domain.n_dofs avg_processor = RTUAvg(sd=self.fe_m_level, n_dofs=n_dofs, avg_fn=QuarticAF(radius=self.avg_radius)) ts = TStepper( u_processor=avg_processor, dof_resultants=True, sdomain=self.fe_domain, bcond_list=[ # define the left clamping BCSlice(var='u', value=0., dims=[0], slice=self.fe_grid_fb[0, 0, 0, :]), # BCSlice( var = 'u', value = 0., dims = [0], slice = self.fe_grid_m[ 0, 0 ] ), # loading at the right edge # BCSlice( var = 'f', value = 1, dims = [0], slice = domain[-1, -1, -1, 0], # time_function = ls ), BCSlice(var='u', value=self.final_displ, dims=[0], slice=self.fe_grid_fb[-1, -1, -1, :], time_function=self.time_function), # BCSlice( var = 'u', value = self.final_displ, dims = [0], slice = self.fe_grid_m[-1, -1], # time_function = self.time_function ), # fix horizontal displacement in the top layer # BCSlice( var = 'u', value = 0., dims = [0], slice = domain[:, -1, :, -1] ), # fix the vertical displacement all over the domain BCSlice(var='u', value=0., dims=[1], slice=self.fe_grid_fb[:, :, :, :]), # # Connect bond and matrix domains BCDofGroup(var='u', value=0., dims=[0], get_link_dof_method=self.fe_grid_fb.get_bottom_dofs, get_dof_method=self.fe_grid_m.get_all_dofs, link_coeffs=[1.]) ], rtrace_list=[ self.sig_eps_m, self.sig_eps_f, self.eps_m_field, self.eps_f_field, self.sig_m_field, self.sig_f_field, self.omega_m_field, self.shear_flow_field, self.slip_field, ]) # Add the time-loop control tloop = TLoop(tstepper=ts, KMAX=300, tolerance=1e-5, debug=False, verbose_iteration=True, verbose_time=False, tline=TLine(min=0.0, step=self.step_size, max=1.0)) tloop.on_accept_time_step = self.plot U = tloop.eval() self.sig_eps_f.refresh() max_sig_m = max(self.sig_eps_m.trace.ydata) return array([U[right_dof_m], max_sig_m], dtype='float_')
shape=(3, 3, 3), fets_eval=fets_eval) # Put the tseval (time-stepper) into the spatial context of the # discretization and specify the response tracers to evaluate there. # right_dof = 2 ts = TS( sdomain=domain, # conversion to list (square brackets) is only necessary for slicing of # single dofs, e.g "get_left_dofs()[0,1]" which elsewise retuns an # integer only bcond_list=[ BCDofGroup(var='u', value=0., dims=[0], get_dof_method=domain.get_left_dofs), BCDofGroup(var='u', value=0., dims=[1, 2], get_dof_method=domain.get_bottom_left_dofs), BCDofGroup(var='u', value=0.002, dims=[1], get_dof_method=domain.get_right_dofs) ], rtrace_list=[ # RTraceGraph(name = 'Fi,right over u_right (iteration)' , # var_y = 'F_int', idx_y = right_dof, # var_x = 'U_k', idx_x = right_dof, # record_on = 'update'),
def example_with_new_domain(): from ibvpy.api import \ TStepper as TS, RTDofGraph, RTraceDomainListField, \ RTraceDomainListInteg, TLoop, \ TLine, BCDof, DOTSEval from ibvpy.mats.mats2D.mats2D_elastic.mats2D_elastic import MATS2DElastic from ibvpy.mats.mats2D.mats2D_sdamage.mats2D_sdamage import MATS2DScalarDamage from ibvpy.api import BCDofGroup mats_eval = MATS2DElastic() fets_eval = FETS2D4Q(mats_eval=mats_eval) #fets_eval = FETS2D4Q(mats_eval = MATS2DScalarDamage()) print(fets_eval.vtk_node_cell_data) from ibvpy.mesh.fe_grid import FEGrid from ibvpy.mesh.fe_refinement_grid import FERefinementGrid from ibvpy.mesh.fe_domain import FEDomain from mathkit.mfn import MFnLineArray # Discretization fe_grid = FEGrid(coord_max=(10., 4., 0.), shape=(10, 3), fets_eval=fets_eval) bcg = BCDofGroup(var='u', value=0., dims=[0], get_dof_method=fe_grid.get_left_dofs) bcg.setup(None) print('labels', bcg._get_labels()) print('points', bcg._get_mvpoints()) mf = MFnLineArray( # xdata = arange(10), ydata=np.array([0, 1, 2, 3])) right_dof = 2 tstepper = TS( sdomain=fe_grid, bcond_list=[ BCDofGroup(var='u', value=0., dims=[0, 1], get_dof_method=fe_grid.get_left_dofs), # BCDofGroup( var='u', value = 0., dims = [1], # get_dof_method = fe_grid.get_bottom_dofs ), BCDofGroup(var='u', value=.005, dims=[1], time_function=mf.get_value, get_dof_method=fe_grid.get_right_dofs) ], rtrace_list=[ RTDofGraph(name='Fi,right over u_right (iteration)', var_y='F_int', idx_y=right_dof, var_x='U_k', idx_x=right_dof, record_on='update'), RTraceDomainListField(name='Stress', var='sig_app', idx=0, position='int_pnts', record_on='update'), # RTraceDomainListField(name = 'Damage' , # var = 'omega', idx = 0, # # record_on = 'update', # warp = True), RTraceDomainListField(name='Displacement', var='u', idx=0, record_on='update', warp=True), RTraceDomainListField(name='Strain energy', var='strain_energy', idx=0, record_on='update', warp=False), RTraceDomainListInteg(name='Integ strain energy', var='strain_energy', idx=0, record_on='update', warp=False), # RTraceDomainListField(name = 'N0' , # var = 'N_mtx', idx = 0, # record_on = 'update') ]) # Add the time-loop control #global tloop tloop = TLoop(tstepper=tstepper, KMAX=300, tolerance=1e-4, tline=TLine(min=0.0, step=1.0, max=1.0)) #import cProfile #cProfile.run('tloop.eval()', 'tloop_prof' ) # print tloop.eval() #import pstats #p = pstats.Stats('tloop_prof') # p.strip_dirs() # print 'cumulative' # p.sort_stats('cumulative').print_stats(20) # print 'time' # p.sort_stats('time').print_stats(20) tloop.eval() # Put the whole thing into the simulation-framework to map the # individual pieces of definition into the user interface. # from ibvpy.plugins.ibvpy_app import IBVPyApp app = IBVPyApp(ibv_resource=tloop) app.main()
def example_2d(): from ibvpy.mats.mats2D.mats2D_elastic.mats2D_elastic import MATS2DElastic from ibvpy.fets.fets2D.fets2D4q import FETS2D4Q from ibvpy.fets.fets2D.fets2D4q8u import FETS2D4Q8U from ibvpy.fets.fets2D.fets2D4q9u import FETS2D4Q9U from ibvpy.fets.fets2D.fets2D9q import FETS2D9Q fets_eval = FETS2D4Q(mats_eval=MATS2DElastic(E=1., nu=0.)) xfets_eval = FETSBimaterial(parent_fets=fets_eval, int_order=3 , mats_eval=MATS2DElastic(E=1., nu=0.), mats_eval2=MATS2DElastic(E=5., nu=0.)) # Discretization fe_domain = FEDomain() fe_level1 = FERefinementGrid(domain=fe_domain, fets_eval=fets_eval) fe_grid1 = FEGrid(coord_max=(3., 1., 0.), shape=(3, 1), fets_eval=fets_eval, level=fe_level1) fe_xdomain = XFESubDomain(domain=fe_domain, fets_eval=xfets_eval, # fe_grid_idx_slice = fe_grid1[1,0], fe_grid_slice=fe_grid1['X - 1.5']) ts = TS(dof_resultants=True, sdomain=fe_domain, bcond_list=[BCDofGroup(var='u', value=1., dims=[0], get_dof_method=fe_grid1.get_right_dofs), BCDofGroup(var='u', value=0., dims=[1], get_dof_method=fe_grid1.get_right_dofs), BCDofGroup(var='u', value=0., dims=[0, 1], get_dof_method=fe_grid1.get_left_dofs), ], rtrace_list=[ # RTDofGraph(name = 'Fi,right over u_right (iteration)' , # var_y = 'F_int', idx_y = 0, # var_x = 'U_k', idx_x = 1), # RTraceDomainListField(name = 'Stress' , # var = 'sig_app', idx = 0, warp = True ), RTraceDomainListField(name='Displacement' , var='u', idx=0, warp=True), RTraceDomainListField(name='Strain' , var='eps', idx=0, warp=True), # RTraceDomainField(name = 'N0' , # var = 'N_mtx', idx = 0, # record_on = 'update') ] ) # # # Add the time-loop control tloop = TLoop(tstepper=ts, # tolerance = 1e-4, KMAX = 4, # debug = True, RESETMAX = 2, tline=TLine(min=0.0, step=1., max=1.0)) # print "elements ",fe_xdomain.elements[0] fe_xdomain.deactivate_sliced_elems() print('parent elems ', fe_xdomain.fe_grid_slice.elems) print('parent dofs ', fe_xdomain.fe_grid_slice.dofs) print("dofmap ", fe_xdomain.elem_dof_map) print("ls_values ", fe_xdomain.dots.dof_node_ls_values) print('intersection points ', fe_xdomain.fe_grid_slice.r_i) print("triangles ", fe_xdomain.dots.rt_triangles) print("vtk points ", fe_xdomain.dots.vtk_X) print("vtk data ", fe_xdomain.dots.get_vtk_cell_data('blabla', 0, 0)) print('ip_triangles', fe_xdomain.dots.int_division) print('ip_coords', fe_xdomain.dots.ip_coords) print('ip_weigths', fe_xdomain.dots.ip_weights) print('ip_offset', fe_xdomain.dots.ip_offset) print('ip_X_coords', fe_xdomain.dots.ip_X) print('ip_ls', fe_xdomain.dots.ip_ls_values) print('vtk_ls', fe_xdomain.dots.vtk_ls_values) print('J_det ', fe_xdomain.dots.J_det_grid) print(tloop.eval()) # #ts.setup() from ibvpy.plugins.ibvpy_app import IBVPyApp ibvpy_app = IBVPyApp(ibv_resource=ts) ibvpy_app.main()
fe_grid5 = FEGrid( name = 'fe_grid5', coord_min = (0, 0, 2.), coord_max = (2., 6.,2.), shape = (1,3), fets_eval = fets_eval_4u, level = fe_rgrid3 ) fe_grid6 = FEGrid( name = 'fe_grid6', coord_min = (2., 6, 2.), coord_max = (10, 15, 2.), shape = (1,3), fets_eval = fets_eval_4u, level = fe_rgrid3 ) ts = TS( dof_resultants = True, sdomain = fe_domain, bcond_list = [BCDofGroup(var='f', value = 1., dims = [0], get_dof_method = fe_grid1.get_top_dofs ), BCDofGroup(var='u', value = 0., dims = [0,1], get_dof_method = fe_grid1.get_bottom_dofs ), ], rtrace_list = [ RTDofGraph(name = 'Fi,right over u_right (iteration)' , var_y = 'F_int', idx_y = 0, var_x = 'U_k', idx_x = 1), RTraceDomainListField(name = 'Stress', var = 'sig_app', idx = 0, warp = False ), # RTraceDomainField(name = 'Displacement' , # var = 'u', idx = 0), # RTraceDomainField(name = 'N0' , # var = 'N_mtx', idx = 0, # record_on = 'update') ] )
def _get_bc_list(self): specmn = self.specmn_fe_grid mid_specmn = self.mid_specmn_fe_grid if self.elstmr_flag: elstmr = self.elstmr_fe_grid #-------------------------------------------------------------- # boundary conditions for the symmetry #-------------------------------------------------------------- # the x-axis corresponds to the axis of symmetry along the longitudinal axis of the beam: bc_symplane_xz = BCSlice(var='u', value=0., dims=[1], slice=specmn[:, 0, :, :, 0, :]) bc_mid_symplane_xz = BCSlice(var='u', value=0., dims=[1], slice=mid_specmn[:, 0, :, :, 0, :]) bc_mid_symplane_yz = BCSlice(var='u', value=0., dims=[0], slice=mid_specmn[0, :, :, 0, :, :]) if self.elstmr_flag: bc_el_symplane_xz = BCSlice(var='u', value=0., dims=[1], slice=elstmr[:, 0, :, :, 0, :]) bc_el_symplane_yz = BCSlice(var='u', value=0., dims=[0], slice=elstmr[0, :, :, 0, :, :]) #-------------------------------------------------------------- # boundary conditions for the support #-------------------------------------------------------------- bc_support_0y0 = BCSlice(var='u', value=0., dims=[2], slice=specmn[-1, :, 0, -1, :, 0]) #-------------------------------------------------------------- # link domains #-------------------------------------------------------------- link_msp_sp = BCDofGroup(var='u', value=0., dims=[0, 1, 2], get_dof_method=mid_specmn.get_right_dofs, get_link_dof_method=specmn.get_left_dofs, link_coeffs=[1.]) # link_msp_sp_xyz = BCSlice(var = 'u', value = 0., dims = [0, 1, 2], # slice = specmn[0, :, :, 0, :, :], # link_slice = mid_specmn[-1 :, :, -1, :, :], # link_dims = [0, 1, 2], # link_coeffs = [1.]) # link_msp_sp_y = BCSlice(var = 'u', value = 0., dims = [1], # slice = specmn[0, :, :, 0, :, :], # link_slice = mid_specmn[-1 :, :, -1, :, :], # link_dims = [1], # link_coeffs = [1.]) # # link_msp_sp_z = BCSlice(var = 'u', value = 0., dims = [2], # slice = specmn[0, :, :, 0, :, :], # link_slice = mid_specmn[-1 :, :, -1, :, :], # link_dims = [2], # link_coeffs = [1.]) # link_msp_sp = [ link_msp_sp_xyz ] if self.elstmr_flag: link_el_sp = BCDofGroup( var='u', value=0., dims=[2], get_dof_method=elstmr.get_back_dofs, get_link_dof_method=mid_specmn.get_front_dofs, link_coeffs=[1.]) #-------------------------------------------------------------- # loading #-------------------------------------------------------------- w_max = self.w_max # f_max = -0.010 / 0.10 # [MN/m] if self.elstmr_flag: # apply displacement at all top node (surface load) # bc_w = BCSlice(var='u', value=w_max, dims=[2], slice=elstmr[:, :, -1, :, :, -1]) # apply a single force at the center of the beam (system origin at top of the elastomer # and us elastomer-domain as load distribution plate with a high stiffness (e.g. steel) # F_max = -0.010 #[MN] # bc_F = BCSlice(var = 'f', value = F_max, dims = [2], # slice = elstmr[0, 0, -1, 0, 0, -1]) else: # center top nodes (line load) # bc_w = BCSlice(var='u', value=w_max, dims=[2], slice=mid_specmn[0, :, -1, 0, :, -1]) # NOTE: the entire symmetry axis (yz)-plane is moved downwards # in order to avoid large indentations at the top nodes # # bc_center_w = BCSlice( var = 'w', value = w_max, dims = [2], slice = mid_specmn[0, :, :, 0, :, :] ) bc_list = [ bc_symplane_xz, bc_mid_symplane_xz, bc_mid_symplane_yz, bc_support_0y0, link_msp_sp, bc_w ] if self.elstmr_flag: bc_list_elstmr = [link_el_sp, bc_el_symplane_xz, bc_el_symplane_yz] bc_list += bc_list_elstmr return bc_list
fe_grid1 = FEGrid( coord_max = ( 1., 1., 0. ), shape = ( 1, 1 ), rt_tol = 0.1, fets_eval = fets_eval, level = fe_level1 ) # fe_grid1.deactivate( (1,0) ) # fe_grid1.deactivate( (1,1) ) fe_xdomain = XFESubDomain( domain = fe_domain, fets_eval = xfets_eval, #fe_grid_idx_slice = fe_grid1[1,0], fe_grid_slice = fe_grid1['X - 0.5 -0.1*Y'] ) ts = TS( dof_resultants = True, sdomain = fe_domain, bcond_list = [BCDofGroup( var = 'u', value = 0., dims = [0, 1], get_dof_method = fe_grid1.get_right_dofs ), BCDofGroup( var = 'u', value = 0., dims = [1], get_dof_method = fe_grid1.get_left_dofs ), BCDofGroup( var = 'u', value = -1., dims = [0], get_dof_method = fe_grid1.get_left_dofs ), ], rtrace_list = [ # RTraceGraph(name = 'Fi,right over u_right (iteration)' , # var_y = 'F_int', idx_y = 0, # var_x = 'U_k', idx_x = 1), # RTraceDomainListField(name = 'Stress' , # var = 'sig_app', idx = 0, warp = True ), RTraceDomainListField( name = 'Displacement' , var = 'u', idx = 0, warp = True ), # RTraceDomainField(name = 'N0' ,