def create_cp_fc_01(L_x=4, L_y=4, n_x=2, n_y=2, z0_ratio=0.1, n_steps=100): '''Create scalable rhombus crease pattern with face constraints One basic element with no general formulation ''' rcp = YoshimuraCreasePattern(n_steps=n_steps, L_x=L_x, L_y=L_y, n_x=n_x, n_y=n_y, show_iter=False, z0_ratio=z0_ratio, MAX_ITER=50) n_h = rcp.N_h n_i = rcp.N_i n_v = rcp.N_v cp = Lifting(n_steps=n_steps) cp.cp_geo(rcp) caf = CnstrTargetFace(F=[r_, s_, 4 * 0.4 * t_ * r_ * (1 - r_ / L_x) + 0.15]) n_arr = np.hstack([n_h[:, :].flatten(), #n_v[:, :].flatten(), n_i[:, :].flatten() ]) cp.init_tf_lst = [(caf, n_arr)] cp.cnstr_lhs = [[(n_h[0, 0], 1, 1.0), (n_h[1, 0], 1, -1.0)], # 1 [(n_h[0, 0], 1, 1.0), (n_h[-1, 0], 1, -1.0)], # 2 [(n_h[0, -1], 1, 1.0), (n_h[1, -1], 1, -1.0)], # 3 [(n_h[0, -1], 1, 1.0), (n_h[-1, -1], 1, -1.0)], # 4 [(n_h[0, -1], 1, 1.0)], [(n_h[1, 0], 0, 1.0)] ] print "n_h[1, 0]", n_h[1, 0] print "n_h[-1,-1]", n_h[-1, 0] print "n_h[1, -1]", n_h[1, -1] print "n_h[-1,-1]", n_h[-1, -1] print "cnstr_lhs", cp.cnstr_lhs print "cnstr_rhs", cp.cnstr_rhs A = 0.2 face_z_t = CF(Rf=z_ - 4 * A * t_ * x_ * (1 - x_ / L_x)) face_x_L2 = CF(Rf=x_ - L_x / 2) cp.cf_lst = [(face_z_t, n_h[0, :]), (face_z_t, n_h[-1, :]), (face_z_t, [n_h[1, 0]]), ] print "edge1", n_h[0, :] print "edge2", n_h[-1, :] return cp
def create_cp_fc_02(L_x=4, L_y=4, n_x=2, n_y=2, z0_ratio=0.1, n_steps=100): '''Create scalable rhombus crease pattern with face constraints One basic element with general formulation (extension in y-direction variabel) (extension in x-direction has to be adepted manually) ''' rcp = YoshimuraCreasePattern(n_steps=n_steps, L_x=L_x, L_y=L_y, n_x=n_x, n_y=n_y, show_iter=False, z0_ratio=z0_ratio, MAX_ITER=50) n_h = rcp.N_h n_i = rcp.N_i n_v = rcp.N_v cp = Lifting(n_steps=n_steps) cp.cp_geo(rcp) caf = CnstrTargetFace(F=[r_, s_, 4 * 0.4 * t_ * r_ * (1 - r_ / L_x) + 0.15]) n_arr = np.hstack([n_h[:, :].flatten(), #n_v[:, :].flatten(), n_i[:, :].flatten() ]) cp.init_tf_lst = [(caf, n_arr)] n_h_idx = n_x / 2 n_h_idx = n_x / 2 y_links = [] for n_arr in n_h[0:3, :].T: for idx, n in enumerate(n_arr[1:]): n_x = len(n_arr) y_links.append([(n_arr[0], 1, 1.0), (n, 1, -1.0)]) ''' Extension in x-direction ''' #y_links.append([(n_h[0,0], 1, 1.0), (n_h[-1,0], 1, -1.0)]) #y_links.append([(n_h[0,0], 1, 1.0), (n_h[-2,0], 1, -1.0)]) y_links.append([(n_h[0, -1], 1, 1.0)]) y_links.append([(n_h[1, 0], 0, 1.0)]) cp.cnstr_lhs = y_links print "n_h[1, 0]", n_h[1, 0] print "n_h[-1,-1]", n_h[-1, 0] print "n_h[1, -1]", n_h[1, -1] print "n_h[-1,-1]", n_h[-1, -1] print "cnstr_lhs", cp.cnstr_lhs print "cnstr_rhs", cp.cnstr_rhs A = 0.2 face_z_t = CF(Rf=z_ - 4 * A * t_ * x_ * (1 - x_ / L_x)) cp.cf_lst = [(face_z_t, n_h[1:-1, 0]), (face_z_t, n_h[0, :]), (face_z_t, n_h[-1, :]) ] print "edge1", n_h[0, :] print "edge2", n_h[-1, :] print "center", n_h[1:-1, 0] return cp
def create_cp_fc_inclined(L_x=4, L_y=4, n_x=2, n_y=4, n_steps=100): '''Create scalable rhombus crease pattern with face constraints ''' rcp = YoshimuraCreasePattern(n_steps=n_steps, L_x=L_x, L_y=L_y, n_x=n_x, n_y=n_y, show_iter=False, MAX_ITER=50) n_h = rcp.N_h n_i = rcp.N_i n_v = rcp.N_v cp = Lifting(n_steps=n_steps) cp.cp_geo(rcp) caf = CnstrTargetFace(F=[r_, s_, 4 * 0.4 * t_ * r_ * (1 - r_ / L_x) + 0.15]) n_arr = np.hstack([n_h[:, :].flatten(), #n_v[:, :].flatten(), n_i[:, :].flatten() ]) cp.init_tf_lst = [(caf, n_arr)] y_links = [] # n_h0 = n_h[(0, -1), :-1] # n_h1 = n_h[(0, -1), 1:] # for nv, nh0, nh1 in zip(n_v.T, n_h0.T, n_h1.T): # for v, h0, h1 in zip(nv, nh0, nh1): # print 'constraining', h0, h1 # y_links.append([(h0, 1, 1.0), (h1, 1, -1.0)]) n_h0 = n_h[(0, -1), :-1] n_h1 = n_h[(0, -1), 1:] for nv in n_v.T: print 'adding constraint', nv y_links.append([(nv[0], 0, 1.0), (nv[1], 0, 1.0)]) # here was a conflict @todo - resolve with Jan # for nv, nh0, nh1 in zip(n_v.T, n_h0.T, n_h1.T): # for v, h0, h1 in zip(nv, nh0, nh1): # y_links.append([(v, 1, 1.0), (h1, 1, -0.5)]) cp.cnstr_lhs = y_links print "cnstr_lhs", cp.cnstr_lhs print "cnstr_rhs", cp.cnstr_rhs # A = L_x * 0.2 A = 0.2 # face_z_t = CF(Rf = z_ - 4 * A * t_ / L_x * x_ * (1 - x_ / L_x)) face_z_t = CF(Rf=z_ - 4 * A * t_ * x_ * (1 - x_ / L_x)) face_x_L2 = CF(Rf=x_ - L_x / 2) face_y_L2 = CF(Rf=y_ - L_y / 2) # face_y_Ly = CF(Rf = y_ - L_y) #old n_h_idx = n_x / 2 z_nodes = n_h[:, :].flatten() # y_nodes = n_i[:, 0] # + list(n_v[:, :].flatten()) y_nodes = n_i[0, 0] # + list(n_v[:, :].flatten()) cp.cf_lst = [(face_y_L2, [n_i[0, 0]]), (face_z_t, z_nodes), ## (face_x_L2, n_h[2, (0, -1)].flatten()), # (face_x_L2, n_h[n_h_idx, (0, -1)].flatten()), (face_x_L2, n_h[n_h_idx, :].flatten()), ] return cp
def create_cp_fc_bow(L_x=4, L_y=4, n_x=4, n_y=2, z0_ratio=0.1, n_steps=100): '''Create scalable rhombus crease pattern with face constraints bad working ''' rcp = YoshimuraCreasePattern(n_steps=n_steps, L_x=L_x, L_y=L_y, n_x=n_x, n_y=n_y, show_iter=False, z0_ratio=z0_ratio, MAX_ITER=50) n_h = rcp.N_h n_i = rcp.N_i n_v = rcp.N_v cp = Lifting(n_steps=n_steps) cp.cp_geo(rcp) caf = CnstrTargetFace(F=[r_, s_, 4 * 0.4 * t_ * r_ * (1 - r_ / L_x) + 0.15]) n_arr = np.hstack([n_h[:, :].flatten(), #n_v[:, :].flatten(), n_i[:, :].flatten() ]) cp.init_tf_lst = [(caf, n_arr)] # y_links = [] # n_h0 = n_h[(0, -1), :-1] # n_h1 = n_h[(0, -1), 1:] # for nv, nh0, nh1 in zip(n_v.T, n_h0.T, n_h1.T): # for v, h0, h1 in zip(nv, nh0, nh1): # y_links.append([(v, 1, 1.0), (h1, 1, -0.5)]) # cp.cnstr_lhs = y_links # cp.cnstr_rhs = np.zeros((len(cp.cnstr_lhs),), dtype = float) # print "cnstr_lhs", cp.cnstr_lhs # print "cnstr_rhs", cp.cnstr_rhs # A = L_x * 0.2 A = 0.2 # face_z_t = CF(Rf = z_ - 4 * A * t_ / L_x * x_ * (1 - x_ / L_x)) face_z_t = CF(Rf=z_ - 4 * A * t_ * x_ * (1 - x_ / L_x)) face_x_L2 = CF(Rf=x_ - L_x / 2) #old # face_y_L2 = CF(Rf = y_ - L_y / 2) #new face_y_L0 = CF(Rf=y_) # face_y_Ly = CF(Rf = y_ - L_y) # n_h_idx = n_x / 2 # z_nodes = n_h[:, :].flatten() # cp.cf_lst = [(face_y_L2, [n_i[0, 0]]), # (face_z_t, z_nodes), ### (face_x_L2, n_h[2, (0, -1)].flatten()), ## (face_x_L2, n_h[n_h_idx, (0, -1)].flatten()), # (face_x_L2, n_h[n_h_idx, :].flatten()), # ] #new n_h_idx = n_x / 2 z_nodes_field1 = n_h[1:n_h_idx, 0].flatten() z_nodes_field2 = n_h[(n_h_idx + 1):-1, 0].flatten() z_nodes_field3 = n_h[1:n_h_idx, 1].flatten() z_nodes_field4 = n_h[(n_h_idx + 1):-1, 1].flatten() cp.cf_lst = [(face_y_L0, [n_h[0, 0]]), (face_y_L0, [n_h[n_h_idx, -1]]), (face_y_L0, [n_h[-1, 0]]), (face_z_t, z_nodes_field1), (face_z_t, z_nodes_field2), (face_z_t, z_nodes_field3), (face_z_t, z_nodes_field4), (face_z_t, n_h[n_h_idx, :]), (face_z_t, n_h[0, :]), (face_z_t, n_h[-1, :]), # (face_z_t, n_h[n_h_idx,1:]), # (face_x_L2, n_h[2, (0, -1)].flatten()), # (face_x_L2, n_h[n_h_idx, (0, -1)].flatten()), (face_x_L2, n_h[n_h_idx, :].flatten()), ] print "field1", z_nodes_field1 print "field2", z_nodes_field2 print "edge1", n_h[0, :] print "edge2", n_h[-1, :] print "center", n_h[n_h_idx, :] return cp