I am trying to implement the constitutive model described in Storm et al to describe a collagen-based ECM. (Storm, Cornelis, et al. "Nonlinear elasticity in biological gels."

The equations for stress and elasticity tensors consider averaging over probability distribution P(r). In my case, I can model P(r) as an ellipsoidal volume matching the 3D distribution of collagen fibers. My question is how do I calculate the average over P(r) while writing a plugin? Is there a way to incorporate the ellipsoidal fiber distribution as described in the User Manual for the uncoupled continuous fiber distribution?

Thanks a lot,

Adil

image.png

*Nature*435.7039 (2005): 191-194) (See equation in attached picture).The equations for stress and elasticity tensors consider averaging over probability distribution P(r). In my case, I can model P(r) as an ellipsoidal volume matching the 3D distribution of collagen fibers. My question is how do I calculate the average over P(r) while writing a plugin? Is there a way to incorporate the ellipsoidal fiber distribution as described in the User Manual for the uncoupled continuous fiber distribution?

Thanks a lot,

Adil

image.png

## Comment