yade-users team mailing list archive

[Roxana Saghafian Larijani <question702522@xxxxxxxxxxxxxxxxxxxxx>]

New question #702522 on Yade:
https://answers.launchpad.net/yade/+question/702522

Hi,
I am using the ViscoElasticCapillar phys but in this model the capillary force is only added when the penetrationDepth is negative which means particles do not have overlap. However, the maximum capillary force which is equal to 
-(2 * M_PI * R * Gamma * cos(phys.theta));
will be exerted on the wet particles while the penetration depth is >=0. I tried to add this force but I am not sure if I am adding it in the right place as the normalForce plot does not change .
In line 205 TIMING_DELTAS_CHECKPOINT("force_calculation_penetr"); I think here the normalForce for of the particles in contact is being calculated based on the viscoelatic physics. Then here I wanted to add the maximum capillary force that I mentioned as below starting from the line " if (geom.penetrationDepth>0.0 && phys.liqBridgeCreated){"


		TIMING_DELTAS_CHECKPOINT("force_calculation_penetr");
		if (computeForceTorqueViscEl(_geom, _phys, I, force, torque1, torque2)) {
			addForce(id1, -force, scene2);
			addForce(id2, force, scene2);
			addTorque(id1, torque1, scene2);
			addTorque(id2, torque2, scene2);
			
			return true;
		} else {
			return false;
		}
		if (geom.penetrationDepth>0.0 && phys.liqBridgeCreated){

		 const auto MaxCapForce= 2 * M_PI *phys.R *phys.gamma*cos(phys.theta);

		 phys.normalForce = -(MaxCapForce) * geom.normal;

			addForce(id1, -phys.normalForce, scene2);
			addForce(id2, phys.normalForce, scene2);
        }
	}

Thanks in advance for your answer!

Regards,
Roxana Saghafian

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.