s_angle = sin(angle_r) c_angle = cos(angle_r) diag = sqrt(2.) ts = TS( tse = tseval, 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 = [ # constraint for all left dofs in x-direction: BCDof(var='u', dof = domain.get_bottom_left_dofs()[0,0], value = 0.)] + # constraint for all left dofs in y-direction: [BCDof(var='u', dof = domain.get_bottom_left_dofs()[0,1], value = 0.)] + [BCDof(var='u', dof = domain.get_bottom_right_dofs()[0,0], value = 2.*c_angle - 2.)]+ [BCDof(var='u', dof = domain.get_bottom_right_dofs()[0,1], value = 2.*s_angle)]+ [BCDof(var='u', dof = domain.get_top_left_dofs()[0,0], value = -1.*s_angle)]+ [BCDof(var='u', dof = domain.get_top_left_dofs()[0,1], value = c_angle - 1.)]+ [BCDof(var='u', dof = domain.get_top_right_dofs()[0,0], value = 2.*c_angle - s_angle - 2.)]+ [BCDof(var='u', dof = domain.get_top_right_dofs()[0,1], value = 2.*s_angle + c_angle - 1.)] , 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'), # RTraceDomainField(name = 'Deformation' , # var = 'eps', idx = 0, # record_on = 'update'), RTraceDomainField(name = 'Displacement' , var = 'u', idx = 1, record_on = 'update',