The determination of bond constant in case named 'Bar Under Tension'

[pleCrazy]

The bc (bond constant) should be equal to 2.0 * emod / (3.1415*area * (delta^2)) in the case named “Bar Under Tension” .
The new results are shown in the following figure.
The original code given by Erkan is probably also incorrect.