I tried to curve-fit viscoelastic material data, but the fit I get is very different than the experimental data I supply. I tried curvefitting with several terms, but without luck. I only read-in data for shear. I'm probably doing something wrong when curve-fitting, but I don't know what. Could you please have a look and try to create a good curve-fit with the data I attach?

The data
0.0001 100e6
.05 43e6
0.1 22e6
.15 15e6
0.2 11e6
.25 10e6


This an example of that old adage in computer science - garbage in = garbage out. In this case, the guess that you are providing is likely sending the algorithm off in the wrong direction.

ANSYS curve fitting is using the Levenberg-Marquette algorithm to perform a non linear regression on the data. One of the problems with Levenberg-Marquette is that it does not guarantee that the minimum energy condition is satisfied, that the minimum obtained using the algorithm may be a local minimum rather than the global minimum. Because of this difficulty, the judgment of the analyst is required to determine when "good enough" has been found.

For the data that you supplied, visual inspection implies that Ao=1e7, A1=9e7 and tau1=.05 should give a good approximation using a 1 term Prony series. The logic here is that the Prony series will look like
G = Ginfinity + G1 * exp(-time/tau1)
where Ginfinity would be the value of the modulus at infinite time, G1 would be the amount of decay (in this case the difference between the modulus at zero time and the modulus at infinite time), and tau1 would be the time constant where 5 time constants represents the required time.

To see this, consider the following two cases. Guess 1 assumes that the default values of 1 for all constants can be used for the curve fitting. Figure guess1a.jpg shows the setup screen from the curve fitting GUI and guess1b.jpg shows the resulting plot comparing the experimental data with the curve fit. Clearly these resultsare not acceptable. Guess 2 assumes that the initial values for the curve fit will be Ao=1e7, A1=9e7 and tau1=1, as shown in Figure guess2a.jpg. The resulting plot comparing the experimental data with the curve fit shows good agreement.

The conclusion is that for non linear curve fitting, the analyst must take care to insure that the progress of the curve fitting algorithm will be in the direction of the global minimum. That poor choice for the initial guesses for the series constants will lead to poor correlations with the experimental data. Unfortunately, the ability to make good choices come from experience, and a lot of that expericence will come from making poor choices. Some suggestions can be provide. Selecting Ao as the infinite modulus will be a good choice. Selecting Ai equal to the total decay divided by the number of terms in the Prony series will be a good choice. Selecting time constants that span the entire range and that are at least a decade apart will be a good choice. Fewer terms in the series, rather than more terms in the series, is probably a good choice.

One final word. The ANSYS auto time stepping algorithm does not survey the viscoelastic behavior in determining whether the time step size should be increased or decreased. At least for now, it is the analyst's responsibility to set the incremental time step size for the problem. Experience suggests that the number of substeps should be no greater than 1/10 of the currently active Prony time constant.





Show Form
No comments yet. Be the first to add a comment!