Lipid vesicle

Description

Freely jointed chain model of a polymer. Attach on both sides to two big beads.

What you can do

Drag the vesicle by holding the mouse.
Press "p" key to increase the volume of the vesicle and see the resulting membrane tension variation
Press "m" key to decrease the volume of the vesicle and see the resulting membrane tension variation
Press "n" key to add small obstacles
Press "s" (left), "f" (right), "e" (up), "d" (down) to move the last adding obstacle

Note/Improvement

Code

Simulations:

function circle(x,y,N,r,k,b) { iPt=getN_Pt; iMb=getN_Mb; newPt(x+r*cos(0*2*PI/N),y+r*sin(0*2*PI/N),k,b); for (i=1;i<N;i++) { newPt(x+r*cos(i*2*PI/N),y+r*sin(i*2*PI/N),k,b); newMb(iPt+i-1,iPt+i); } newMb(iPt+N-1,iPt+0); } function line(x,y,N,r,k,b) { iPt=getN_Pt; iMb=getN_Mb; newPt(x,y,k,b); for (i=1;i<N;i++) { newPt(x+r*i,y,k,b); newMb(iPt+i-1,iPt+i); } //newMb(iPt+N-1,iPt+0); } print(commands); setMovable(0); Nv=50; R=1.5; circle(0,-0,Nv,R,200,1); //setMove(0,0);setMove(Nv-1,0); defPolyR(0,Nv,3.04*R*R); //circle(0,-3,Nv,1,100,1); //defPolyR(getN_Pt-Nv,getN_Pt,3.2); //setMove(Nv,0);setMove(2*Nv-1,0); // Pipette //newPt(-1.2,0.4,0,0);newPt(-2.2,0.4,0,0);newPt(-1.2,-0.4,0,0);newPt(-2.2,-0.4,0,0);newMb(Nv,Nv+1);newMb(Nv+2,Nv+3); //setMove(Nv/2,0);setMove(Nv-1,0); N=0;Ns=1500; // Premier graphe : setWindow(0,0); //newGraph(Ns); //colorGraph(1,0,0); //setBounds(0,0,Ns,1); //g1 = lastGraph; N_Mean_Tension=15; newGraph(Ns); colorGraph(1,1,0); setBounds(0,0,4,1); g2 = lastGraph; colorGraph(1,0,0); setBounds(0,-4,9,5.5); setAxisX(1,0,1,0,0.1); setAxisY(1,0,0.5,0,0.1); //selectGraph(g1); setT(2); setViewScale(0.1); setViewX(0.4); setViewY(-0.3); addLabel(1,1,'log(Surface tension)'); addLabel(7.5,-0.5,"Area"); colorGraph(0,0,0); setViewXYpG(0,-1); gVM=0;gTM=0; hold=0; nhold=0; hasdragged=0; distdragging=0.25; dd2=distdragging*distdragging; sync(true) { N=N+1; if (N>Ns) {N=N-Ns;} if (N%N_Mean_Tension==0) { gVM=gVM/N_Mean_Tension; gTM=gTM/N_Mean_Tension; selectGraph(g2); addPtGraph(gVM,log(gTM)/log(10)); //selectGraph(g1); //addPtGraph(N,gTM); gVM=0;gTM=0;} gVM=gVM+getV(0); gTM=gTM+getTension(0); DL=0.01; if (isMouseDragged==1) { if (hasdragged==1) { // continue to drag point kdrag=50; addFX(nhold,kdrag*(getMouseX-getPtX(nhold))); addFY(nhold,kdrag*(-getMouseY-getPtY(nhold))); //setPtY(nhold,-getMouseY); } else { // look for a point to drag found=0;i=1;gNP=getN_Pt;mx=getMouseX; my=-getMouseY; while ((found==0)&(i<gNP)) { px=getPtX(i);py=getPtY(i); if (((mx-px)*(mx-px)+(my-py)*(my-py))<dd2) { found=1; nhold=i;} i=i+1;} if (found==1) { writeInfo('Trouve'); hasdragged=1; } } } else { if (hasdragged==1) { // free dragged point hasdragged=0; addFX(nhold,0);addFY(nhold,0); } } if (kP==1) { print(getV0_Pol(0)); a=rK; if(a=='p') { setV0_Pol(0,getV0_Pol(0)*1/0.998); } if(a=='m') { setV0_Pol(0,getV0_Pol(0)*0.998); } if(a=='o') { setV0_Pol(1,getV0_Pol(1)*1/0.92); } if(a=='l') { setV0_Pol(1,getV0_Pol(1)*0.92); } if (a=='n') { print('X='+getMouseX+'Y='+getMouseY); // newPt(getMouseX,-getMouseY,10,0); circle(getMouseX,-getMouseY,4,0.1,6000,0); } if(a=='f') { setPtX(getN_Pt-1,getPtX(getN_Pt-1)+DL); setPtX(getN_Pt-2,getPtX(getN_Pt-2)+DL); setPtX(getN_Pt-3,getPtX(getN_Pt-3)+DL); setPtX(getN_Pt-4,getPtX(getN_Pt-4)+DL); } if(a=='s') { setPtX(getN_Pt-1,getPtX(getN_Pt-1)-DL); setPtX(getN_Pt-2,getPtX(getN_Pt-2)-DL); setPtX(getN_Pt-3,getPtX(getN_Pt-3)-DL); setPtX(getN_Pt-4,getPtX(getN_Pt-4)-DL); } if(a=='e') { setPtY(getN_Pt-1,getPtY(getN_Pt-1)+DL); setPtY(getN_Pt-2,getPtY(getN_Pt-2)+DL); setPtY(getN_Pt-3,getPtY(getN_Pt-3)+DL); setPtY(getN_Pt-4,getPtY(getN_Pt-4)+DL); } if(a=='d') { setPtY(getN_Pt-1,getPtY(getN_Pt-1)-DL); setPtY(getN_Pt-2,getPtY(getN_Pt-2)-DL); setPtY(getN_Pt-3,getPtY(getN_Pt-3)-DL); setPtY(getN_Pt-4,getPtY(getN_Pt-4)-DL); } } }