Hi Bill, Two questions. I was checking the function where the forces are computed (in InsParameter.C), and I wonder why do you multiply the pressure by the fluid density (see getNormalForce)? Between the euler and leapfrog method in RigodBodyMotion, which one gives better results? regards, Joel ----- Original Message ---- From: Bill Henshaw <henshaw@xxxxxxxx> To: overture@xxxxxxxxxxxxx Sent: Thu, April 1, 2010 8:47:29 PM Subject: [overture] Re: rigid body motion Yes! Joel Guerrero wrote: > Hi Bill, > > > The forces used in RigidBodyMotion.C are the forces computed in MovingGrids.C? > > > Joel > > > > ----- Original Message ---- > From: Bill Henshaw <henshaw@xxxxxxxx> > To: overture@xxxxxxxxxxxxx > Sent: Thu, April 1, 2010 5:52:36 AM > Subject: [overture] Re: rigid body motion > > Joel, > It seems that you just need to change the Newton-Euler quations > RigidBodyMotion.C > in functions integrate and correct. > > Regards, > Bill > > Joel Guerrero wrote: >> Hi, >> >> Specifically I want to to model a spring-damper-mass system, where the >> forces (let say the vertical force) is calculated as follow: m*a+k*v+c*y=fy >> plus an additional derivate in the RHS. Currently and if I am not mistaken, >> the forces are computed as fy = m*a, hence, I want to modify the way fy is >> computed by adding the additional terms. >> >> Regards, >> >> Joel >> >> >> >> ----- Original Message ---- >> From: Bill Henshaw <henshaw@xxxxxxxx> >> To: overture@xxxxxxxxxxxxx >> Sent: Wed, March 31, 2010 12:42:50 AM >> Subject: [overture] Re: rigid body motion >> >> Hi Joel, >> I wonder if you can be more specific on what exactly you want to do? >> >> Regards, >> Bill >> >> Joel Guerrero wrote: >>> Hi, >>> >>> I have to modify the way forces are computed for rigid body motion, I >>> wonder which files should I modify (rigidbodymotion or movinggrids)? What >>> I need is to add three more derivatives to the way that forces are compute. >>> Any suggestion will be appreciated. >>> >>> Cheers, >>> >>> joel >>> >>> >>> >>> >> >> >> > > > > >