Dear users,
I have created a dimer particle using a raspberry model, and try to make rigid both central and surface particles (no fluctuation of distance between particles). But i am unable to do this. So I am requesting for a help.
The lines are bellow:
# First created a central particle
system.part.add(pos=pos,id=id, type=N*[0], rotation=N*[(1,1,0)],fix=N*[(0,0,1)])
# after central particle, we created surface particle and make a dimer
for i in range(0, N):
colSurfPos = np.random.randn(2)
colSurfPos = colSurfPos / np.linalg.norm(colSurfPos) * radius_col
colSurfPos = np.append(colSurfPos,0)
dip[i] = colSurfPos
system.part[i].dip = dip[i]
colSurfPos = colSurfPos + pos[i]
system.part.add(id=i+N, pos=colSurfPos, dip=dip[i], type=0)
system.part[i+N].add_bond((col_center_surface_bond, i))
# After creating a dimer we make it as a rigid using below lines. But the distance between the surface and the central particle is fluctuating.
for i in range(0,N):
m = system.part[i+N]
p = i
m.vs_auto_relate_to(p)
system.part.add(pos=pos,id=id, type=N*[0], rotation=N*[(1,1,0)],fix=N*[(0,0,1)])
# after central particle, we created surface particle and make a dimer
for i in range(0, N):
colSurfPos = np.random.randn(2)
colSurfPos = colSurfPos / np.linalg.norm(colSurfPos) * radius_col
colSurfPos = np.append(colSurfPos,0)
dip[i] = colSurfPos
system.part[i].dip = dip[i]
colSurfPos = colSurfPos + pos[i]
system.part.add(id=i+N, pos=colSurfPos, dip=dip[i], type=0)
system.part[i+N].add_bond((col_center_surface_bond, i))
# After creating a dimer we make it as a rigid using below lines. But the distance between the surface and the central particle is fluctuating.
for i in range(0,N):
m = system.part[i+N]
p = i
m.vs_auto_relate_to(p)
I am also attaching the whole lines of code.
Your help will be highly appreciated.
Thanks & Regards
Chandra Shekhar Maurya
Mechanical Eng. Dept.
Indian Institute of Technology Patna
Contact No:9793572837
