2010/4/30 Václav Šmilauer <[email protected]> > > For your code, it looks normal;
now I don't know what I must do for solving it, what is your suggestion? comment out all non-elastic things > (including plasticity!!) and run the test again, i.e. you will have just > something like > > > Real displN=geom->displacementN(); > Real Fn=phys->kn*displN; > phys->normalForce=Fn*geom->normal; > Vector3r Fs=phys->ks*geom->displacementT(); > applyForceAtContactPoint(....) > > BTW it is quite confusing to call Fs "trialFs", since trial force (or > trial stress) is coloquially one before you do stress return in case of > plastic flow (you "try" if it is outside or inside, hence "trial"): > > > applyForceAtContactPoint(phys->normalForce+trialFs,geom->contactPoint,contact->getId1(),geom->se31.position,contact->getId2(),geom->se32.position,scene); > > > BTW2, at this line: > > Real crossSection=sqrt(2)*pow(geom->refR1,2); > > order of id1 and id2 is arbitrary, and here you use just refR1 > asymetrically. > you are right, thanks. :) Nasibeh
_______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
