def get_constrained_YCP(L_x, L_y, n_x, n_y):
def geo_trans(X):
x, y, z = X.T
y = y - 1.1
# y[2], y[3], y[4], y[5] = -1.428, 1.428, -1, 1
y[2], y[3], y[4], y[5] = -0.7488, 0.7488, -0.275, 0.275
# x[2], x[3], x[8] , x[9] = 3.43, 3.43, 1.93, 4.93
x[2], x[3], x[8] , x[9] = 3.6925, 3.6925, 2.045, 5.34
return np.c_[x, y, z]
cp = YoshimuraCreasePattern(L_x=L_x, L_y=L_y, n_x=n_x, n_y=n_y,
geo_transform=geo_trans)
fixed_node = fix(cp.N_h[0, 0], (0, 1, 2))
planar_front_boundary = link(cp.N_h[0, 0], 1, 1.0,
cp.N_h[1:, 0], 1, -1.0)
planar_back_boundary = link(cp.N_h[0, -1], 1, 1.0,
cp.N_h[1:, -1], 1, -1.0)
linked_left_boundary_x = link(cp.N_h[0, 0], 0, 1.0,
cp.N_h[0, 1:], 0, -1.0)
linked_left_boundary_z = link(cp.N_h[0, 0], 2, 1.0,
cp.N_h[0, 1:], 2, -1.0)
linked_left_and_right_z = link(cp.N_v[0, :], 2, 1.0,
cp.N_v[1, :], 2, -1.0)
# linked_right_boundary_x = link(ycp.N_v[-1, 0], 0, 1.0,
# ycp.N_v[-1, 1:], 0, -1.0)
cntrl_displ = [([(cp.N_v[1, 0], 0, 1.0)], -1)]
cs = fixed_node + planar_front_boundary + planar_back_boundary + linked_left_boundary_x + linked_left_boundary_z + linked_left_and_right_z + cntrl_displ
caf = CnstrTargetFace(F=[r_, s_, 4 * 0.4 * t_ * r_ * (1 - r_ / L_x) + 0.15])
n_arr = np.hstack([cp.N_h[:, :].flatten(),
cp.N_i[:, :].flatten()
])
init = Initialization(cp=cp, tf_lst=[(caf, n_arr)])
fold = Folding(source=init, n_steps=10, dof_constraints=cs)
fold.u_1
print 'nodal coordinates', fold.x_1
print 'facets', fold.F
print 'facet coordinates', fold.x_1[fold.F]
v = CreasePatternView(root=init)
v.configure_traits()
return fold
def get_constrained_YCP(L_x, L_y, n_x, n_y):
ycp = YoshimuraCreasePattern(L_x=L_x, L_y=L_y, n_x=n_x, n_y=n_y)
fixed_node = fix(ycp.N_h[0, 0], (0, 1, 2))
planar_front_boundary = link(ycp.N_h[0, 0], 1, 1.0,
ycp.N_h[1:, 0], 1, -1.0)
planar_back_boundary = link(ycp.N_h[0, -1], 1, 1.0,
ycp.N_h[1:, -1], 1, -1.0)
linked_left_boundary_x = link(ycp.N_h[0, 0], 0, 1.0,
ycp.N_h[0, 1:], 0, -1.0)
linked_left_boundary_z = link(ycp.N_h[0, 0], 2, 1.0,
ycp.N_h[0, 1:], 2, -1.0)
linked_left_and_right_z = link(ycp.N_v[0, :], 2, 1.0,
ycp.N_v[1, :], 2, -1.0)
linked_right_boundary_x = link(ycp.N_v[-1, 0], 0, 1.0,
ycp.N_v[-1, 1:], 0, -1.0)
cntrl_displ = [([(ycp.N_h[-1, 1], 0, 1.0)], -1.15)]
cs = fixed_node + planar_front_boundary + planar_back_boundary + linked_left_boundary_x + linked_left_boundary_z + linked_left_and_right_z + linked_right_boundary_x + cntrl_displ
caf = CnstrTargetFace(F=[r_, s_, 4 * 0.4 * t_ * r_ * (1 - r_ / L_x) + 0.15])
n_arr = np.hstack([ycp.N_h[:, :].flatten(),
ycp.N_i[:, :].flatten()
])
init = Initialization(cp=ycp, tf_lst=[(caf, n_arr)])
fold = Folding(source=init, n_steps=10, dof_constraints=cs)
fold.u_1
print 'nodal coordinates', fold.x_1
print 'facets', fold.F
print 'facet coordinates', fold.x_1[fold.F]
print "punkte" , fold.x_t[-1]
print "hallo", ycp.F_L
print "nodes_neighbors", ycp.N_neighbors
v = CreasePatternView(root=init)
v.configure_traits()
return fold
source=init,
n_steps=10,
MAX_ITER=500,
goal_function_type="none",
dof_constraints=fixed_node
+ planar_front_boundary
+ planar_back_boundary
+ linked_left_boundary_x
+ linked_left_boundary_z
+ linked_left_and_right_z
+ linked_right_boundary_x
+ cntrl_displ,
# tf_lst=[(caf, n_arr)]
)
print "u", lift.u_t[-1]
return init, lift
"""configure parameters:"""
init, fold = get_constrained_YCP(
L_x=6.3, L_y=4.2, n_x=4, n_y=8, d=-1.25
) # l_x length, l_y length, n_x number of elments, n_y number of Elements, d deformation of the right side
v = CreasePatternView(root=init)
v.configure_traits()
al = AbaqusLink(data=fold, n_split=10)
al.model_name = "x3y8"
al.build_inp()
n_r_h = cp.N_h[-1, (0, -1)].flatten()
n_lr_h = cp.N_h[(0, 0, -1, -1), (0, -1, 0, -1)].flatten()
n_fixed_y = cp.N_h[(0, -1), 2].flatten()
corner_slides = [([(cp.N_h[0, 0], 0, L_x), (cp.N_h[0, 0], 1, -L_y)], 0),
([(cp.N_h[-1, 0], 0, L_x), (cp.N_h[-1, 0], 1, L_y)], 0),
([(cp.N_h[-1, -1], 0, L_x), (cp.N_h[-1, -1], 1, -L_y)], 0),
([(cp.N_h[0, -1], 0, L_x), (cp.N_h[0, -1], 1, L_y)], 0),
([(cp.N_h[2, 0], 2, 1.0), (cp.N_h[2, -1], 2, -1.0)], 0),
([(cp.N_h[3, 0], 2, 1.0), (cp.N_h[3, -1], 2, -1.0)], 0),
([(cp.N_h[3, 0], 2, 1.0), (cp.N_h[2, 0], 2, -1.0)], 0),
([(cp.N_h[0, 2], 2, 1.0), (cp.N_h[-1, 2], 2, -1.0)], 0),
]
face_z_t = CnstrTargetFace(F=[r_, s_, -0.6 * t_ * (r_ * (1 - r_ / L_x) + s_ * (1 - s_ / L_y))])
init = Initialization(cp=cp, tf_lst=[(face_z_t, cp.N)], t_init=0.1)
fold = Folding(source=init,
goal_function_type='potential_energy',
n_steps=40,
MAX_ITER=1000,
#tf_lst=[(face_z_t, cp.N)],
dof_constraints=fix(n_l_h, [0], v) + fix(n_lr_h, [2]) + \
fix(n_fixed_y, [1]) + fix(n_r_h, [0], -v) + \
corner_slides)
# dof_constraints=fix([0, 1], [0], v) + fix([0, 1, 6, 7], [2]) + \
# fix([0], [1]) + fix([6, 7], [0], -v))
fold.u_t[-1]
cpw = CreasePatternView(root=init)
cpw.configure_traits()
