iMechanica - Inhomogeneous swelling of a gel in equilibrium with a solvent and mechanical load - Comments

I already found the answer.

Mario Juha — Thu, 10 Dec 2009 15:16:16 -0500

I already found the answer. Thanks I was measuring lambda respect to the free-swelling state and not from the dry network state. Now I have same results. But if you could send the info, I would appreciate it.

Thanks,

Mario

Material model implementation

Mario Juha — Thu, 10 Dec 2009 10:34:08 -0500

Dear Dr. Hong.

I have implemented your material model in a non-linear finite element program that I wrote. I am using a Continuum Mechanics Total Lagrangian formulation to deal with nonlinearities, and my measure for stress and strain are Second Piola-Kirchhoff stress tensor and Green-Lagrange strain tensor, respectively. In order to test the program, I would like to reproduce your results. I have followed what you nicely explained in the paper, but I feel I need an extra help.

If it is appropiate, could you give me extra information about the finite element model?, that is, what exactly boundary conditions you used and how many substeps?.

To use your model and doing the following:

1) Specify the initial chemical potential and then calculate lambda_0 using Eq.(19). This value of lambda_0 is what I am going to use in Eq.(23)

2) Increment linearly the current chemical potential using the following relation:

current_mu = mu_0 + delta * ( current_mu - mu_0 ), where delta is the increment corresponding to the load step.

current_mu is what I am going to use in Eq.(23)

3) Solve nonlinear equations using Newton-Raphson.

Thanks in advance,

Mario

Re: Future work

Wei Hong — Wed, 30 Sep 2009 09:46:33 -0400

Dear Mario,

Sure you can! Please let me know if you need more information than those on the paper.

Wei

Future work...

Mario Juha — Tue, 29 Sep 2009 16:47:13 -0400

Dear Dr. Hong.

I am really interested in this work. I have been studied it for weeks and I would like to extend your results. In section 8.3, last paragraph, you said something about instabilities that you can not solve currently with ABAQUS. I would like to know if I can collaborate with you and or your students.

thank you,

Mario

Re: Molecular incompressibility

Wei Hong — Tue, 01 Sep 2009 16:35:06 -0400

Dear Mario,

It is actually similar to the constraint you are familiar with, the only difference being the volume of solvent. When solvent migrate into a polymer network, the network expands. Here we assume the total volume to be a constant.

Your understanding of C(X) is correct.

Wei

Molecular incompressibility and the free-energy function

Mario Juha — Tue, 01 Sep 2009 15:26:00 -0400

Dear Dr. Hong. I have been busy studying your papers and I have more questions than answer (to me that is good). In: W.Hong, X. Zhao, J. Zhou, Z. Suo. A theory of coupled diffusion and large deformation in polymeric gels. Journal of the mechanics and physics 56, 1779 - 1793 (2008), section 3, you specify the molecular incompressibility as:

1 + vC = det(F)

I am familiar with the incompressibility constraint of det(F) = 1, but not with the above equation. Could you give more insight about it? Probably it is an standard constraint in gel theory, but I am not familiar with it.

Other question that I have is related with C(X), distribution of solvent molecules in the gel. Do I need to understand this notation as C(X) = C(phi^-1(x))? , where phi(x)^-1 is the inverse of the motion, lower x is the spatial coordinate and capital X is the reference coordinate. In general terms, what I have understood is that C is an inhomogeneous field, is it correct?

cordially,

Mario J. Juha

www.eng.usf.edu/~mjuha/

Re: Questions about your paper

Wei Hong — Thu, 27 Aug 2009 20:50:12 -0400

Dear Mario,

Thank you for your interest in our work.

I don't quite understand your first question though. I = F_iK F_iK is an invariant of the deformation gradient F. It is not the first invariant of F. In fact, it is the first invariant of the right Cauchy-Green deformation tensor. You can find reference on these quantities on almost any textbook on finite-deformation continuum mechanics, or even on wikipedia.

You are right that on the final published paper, the lambda_0^3 is missing in eq (23). The version here on imechanica was correct. Thank you for pointing this out, I appologize for our carelessness on proofreading the final print.

Wei

Questions about your paper

Mario Juha — Thu, 27 Aug 2009 18:01:00 -0400

Dear Dr. Hong.

I am trying to find an idea for my proposal presentation for my Ph.D and I have found interest in yor work about gels. So, I read two papers that you wrote about it.

I have two simple questions related with the paper: W. Hong, Z. Liu, and Z. Suo, Inhomonegeous swelling of a gel in equilibrium with a solvent and mechanical load. International Journal of Solids and Structures 46, 3282-3289 (2009).

1) Why do you relate I = F_ik F_ik with the first invariant of the deformation gradient (eq.(17)) ? Could you give me some reference? Now I am reading the Gerahard Holzapfel's book about continuum mechanics, and I have not found any reference to it.

2) Is there a lambda_0^3 term missing in the second log, in eq.(23) ? I have done the mathematic and I just found a discrepancy in the numerator of the second log.

I am a Ansys user and I want to reproduce your excellents results, but using Ansys.

Thank you,

Mario J. Juha

Civil and environmetal Engineering Department

University of South Florida

http://www.eng.usf.edu/~mjuha/

Dear Wei, Thanks for

marc — Thu, 29 May 2008 22:17:55 -0400

Dear Wei,

Thanks for your explanation for the free-energy function. As you said,we must do some experiments on soft materials which have no corresponding free energy funtion. the acquirement of the free-energy function for a new soft material may be very difficult, and the experimental test is very importrant. For a new soft material without free-energy function,experimental tests may be the only way to describe its mechanical behavior!

THANKS

Lianhua

Limitation of the current work

Wei Hong — Wed, 28 May 2008 16:01:38 -0400

Dear Lianhua,

Thank you for your interest in our work!

As you mentioned, the current implementation does have its limitations: it is only suitable for the final equilibrium state of diffusion. Although the chemical potential is not required to be constant in the simulation, it is a predefined field. In other words, we can not solve for the chemical potential, it must be given. Therefore a steady-state calculation of a complex domain might not be possible either.

No theory could ever predict a general free-energy function, although some theoretical abstraction might give insights to some specific material behavior, Florry-Huggins, for example. The right way to investigate the mechanical behavior of a material would always be experimental. Instead of doing one experiment on one material, one should do a series of experiment on a same material, using different loading conditions, different sample shape/sizes. Instead of starting from nowhere, I think it is always better to start from a theoretical model, and see the deviation. If there is no deviation, good, we extract the material parameter; if, most likely, there is deviation, we either modify the theory to say why, or just use the test result numerically, if a result is important. Also instead of testing the static/equilibrium behavior together with the kinetic properties, I suggest to do separate tests for a same material. These are just my general thoughts. Let's keep on the discussion if you have further interest.

Thanks,

Wei

good job.

Lianhua Ma — Wed, 28 May 2008 07:12:27 -0400

Dear Wei Hong,Your investigation on mechanics of gel is very good job.

you had implemented the theory in the finite element package, ABAQUS. the simulation of large deformation of similar hyperelastic material requires the satisfaction of quasi-static mechanical behavior. so , the FEM simulation in your paper did not reflect the diffusion process of solvent in gel. In other words, you assumed that the chemical potential is constant in one simulation, In fact, the chemical potential is variable in different position of gel. the theory in your paper may be only suitable for final equilibrium state of diffusion.

another question: the key of your paper is the free energy funtion W . For other soft materials, If we have no corresponding free energy funtion put forward by predecessor, How can we investigate the mechanical behavior of soft materials? can you give me any advice? thanks for your paper! Hope to keep in touch with you.

THANKS

L.H. MA

Re: additional questions

Wei Hong — Tue, 27 May 2008 10:58:10 -0400

We don't need to specify the thermal expansion coefficient. We use the "T" just as a general field parameter, not as temperature, so it has nothing to do with thermal expansion.

Due to the definition of the chemical potential, it is always negative. It can be negative infinity to 0. (-0.05 is just an arbitrary number we picked.) As we start from -0.05 and end in 0, so it is still swelling instead of contracting. Just don't read it as temperature.

Please feel free to let me know if you have further concerns.

Wei

additional questions

Minkyoo Kang — Fri, 23 May 2008 18:31:18 -0400

Thank you for your prompt response.

I have additional questions. In defining the temperature in the predifined fields, we usually need to define thermal expansion coefficient in the property module, otherwise we don't see any deformation. I wonder in your simulation whether you also need to specify thermal expansion coefficient, and if that's the case, how does the thermal expansion coefficient relate with the chemical potential?

Your chemical potentials are in the range of -0.05 to 0. Does this means we have to use the same values for the temperture input in the predefiend field which actually means contraction not swelling in the point of temperature?

I'm bothering you with many questions and I really appreciate your help.

Thanks,

Min Kyoo

setting chemical potential

Wei Hong — Fri, 23 May 2008 14:42:40 -0400

Dear Min Kyoo,

Thank you for your interest in our work.

Yes, you should specify the initial chemical potential according to the given initial free swelling ratio (3rd material parameter)

The chemical potential is specified using pre defined fields (temperture) in abaqus input.

Let me know if you have further questions.

Wei

Question about the chemical potential

Minkyoo Kang — Thu, 22 May 2008 19:16:58 -0400

Dear Wei Hong,

Thank you for sharing your source code.

I'm trying to use it for a swelling deformation problem. I have a question about the chemical potential. Since the initial free swelling is an input parameter, shoud we specify the initial chemical potential accoroding to the free swelling equation? If so, since the chemical potential is mimicked by a temperature-like variable in ABAQUS as stated in your paper, How do we specify the increment of the chemical potential or temperature?

Thanks,

Min Kyoo Kang

Inhomogeneous swelling of a gel in equilibrium with a solvent and mechanical load

Wei Hong — Wed, 07 May 2008 16:45:32 -0400

A network of polymers can imbibe a large quantity of a solvent and swell, resulting in a gel. The swelling process can be markedly influenced by a mechanical load and geometric constraint. When the network, solvent, and mechanical load equilibrate, the gel usually swells by a field of inhomogeneous and anisotropic deformation. We show that this field in the swollen gel is equivalent to that in a hyperelastic solid. We implement this theory in the finite-element package, ABAQUS, and analyze examples of swelling-induced deformation, contact, and bifurcation. Because commercial software like ABAQUS is widely available, this work may provide a powerful tool to study complex phenomena in gels.

The source code of the UHYPER program is attached below.