FiberSection3d question

[jwaugh]

In FiberSection3d.cpp the strain for a trial state is calculated for each fiber in the line:

FiberSection3d::SetTrialSectionDeformation(const Vector &deforms)
double strain = d0 + y*d1 + z*d2;

Now I would assume this means that d0 is the axial strain, d1 is the curvature about the z axis, and d2 is hte curvature about the y axis. If one wanted to make the strains vary linear across the width of a section under bending about the y-axis then you could use the following expression:

double strain = d0 + y*d1 + z*d2 - 1/2*fabs(y/(flangeWidth/2))*z*d2

However if I do that, the strain doesn't vary across the section like it should. Instead of varying across the width it is constant like in the normal FiberSection3d. Is the strain calculated somewhere else? I am trying to make it so that on a wide section the strains at the outside edges is lower than the strain at the centroid. And it is not doing this even though it seems like it should.

THanks for all the help.