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',
warp = True),
RTraceDomainField(name = 'Damage' ,
