Geometric integrator

Template:Continuum mechanics In continuum mechanics, the finite strain theory—also called large strain theory, or large deformation theory—deals with deformations in which both rotations and strains are arbitrarily large, i.e. invalidates the assumptions inherent in infinitesimal strain theory. In this case, the undeformed and deformed configurations of the continuum are significantly different and a clear distinction has to be made between them. This is commonly the case with elastomers, plastically-deforming materials and other fluids and biological soft tissue.

Displacement

Figure 1. Motion of a continuum body.

The displacement of a body has two components: a rigid-body displacement and a deformation.

• A rigid-body displacement consists of a simultaneous translation and rotation of the body without changing its shape or size.
• Deformation implies the change in shape and/or size of the body from an initial or undeformed configuration ${\kappa }_0(\mathcal{ℬ})$ to a current or deformed configuration ${\kappa }_t(\mathcal{ℬ})$ (Figure 1).

A change in the configuration of a continuum body can be described by a displacement field. A displacement field is a vector field of all displacement vectors for all particles in the body, which relates the deformed configuration with the undeformed configuration. Relative displacement between particles occurs if and only if deformation has occurred. If displacement occurs without deformation, then it is deemed a rigid-body displacement.

Material coordinates (Lagrangian description)

The displacement of particles indexed by variable Buying, selling and renting HDB and personal residential properties in Singapore are simple and transparent transactions. Although you are not required to engage a real property salesperson (generally often known as a "public listed property developers In singapore agent") to complete these property transactions, chances are you'll think about partaking one if you are not accustomed to the processes concerned.

Professional agents are readily available once you need to discover an condominium for hire in singapore In some cases, landlords will take into account you more favourably in case your agent comes to them than for those who tried to method them by yourself. You need to be careful, nevertheless, as you resolve in your agent. Ensure that the agent you are contemplating working with is registered with the IEA – Institute of Estate Brokers. Whereas it might sound a hassle to you, will probably be worth it in the end. The IEA works by an ordinary algorithm and regulations, so you'll protect yourself in opposition to probably going with a rogue agent who prices you more than they should for his or her service in finding you an residence for lease in singapore.

There isn't any deal too small. Property agents who are keen to find time for any deal even if the commission is small are the ones you want on your aspect. Additionally they present humbleness and might relate with the typical Singaporean higher. Relentlessly pursuing any deal, calling prospects even without being prompted. Even if they get rejected a hundred times, they still come again for more. These are the property brokers who will find consumers what they need eventually, and who would be the most successful in what they do. 4. Honesty and Integrity

This feature is suitable for you who need to get the tax deductions out of your PIC scheme to your property agency firm. It's endorsed that you visit the correct site for filling this tax return software. This utility must be submitted at the very least yearly to report your whole tax and tax return that you're going to receive in the current accounting 12 months. There may be an official website for this tax filling procedure. Filling this tax return software shouldn't be a tough thing to do for all business homeowners in Singapore.

A wholly owned subsidiary of SLP Worldwide, SLP Realty houses 900 associates to service SLP's fast rising portfolio of residential tasks. Real estate is a human-centric trade. Apart from offering comprehensive coaching applications for our associates, SLP Realty puts equal emphasis on creating human capabilities and creating sturdy teamwork throughout all ranges of our organisational hierarchy. Worldwide Presence At SLP International, our staff of execs is pushed to make sure our shoppers meet their enterprise and investment targets. Under is an inventory of some notable shoppers from completely different industries and markets, who've entrusted their real estate must the expertise of SLP Worldwide.

If you're looking for a real estate or Singapore property agent online, you merely need to belief your instinct. It is because you don't know which agent is sweet and which agent will not be. Carry out research on a number of brokers by looking out the internet. As soon as if you find yourself certain that a selected agent is dependable and trustworthy, you'll be able to choose to utilize his partnerise find you a house in Singapore. More often than not, a property agent is considered to be good if she or he places the contact data on his web site. This is able to imply that the agent does not thoughts you calling them and asking them any questions regarding properties in Singapore. After chatting with them you too can see them of their office after taking an appointment.

Another method by way of which you could find out whether the agent is sweet is by checking the feedback, of the shoppers, on the website. There are various individuals would publish their comments on the web site of the Singapore property agent. You can take a look at these feedback and the see whether it will be clever to hire that specific Singapore property agent. You may even get in contact with the developer immediately. Many Singapore property brokers know the developers and you may confirm the goodwill of the agent by asking the developer. may be expressed as follows. The vector joining the positions of a particle in the undeformed configuration $P_i$ and deformed configuration $p_i$ is called the displacement vector. Using $\mathbf{X}$ in place of $P_i$ and $\mathbf{x}$ in place of $p_i$, both of which are vectors from the origin of the coordinate system to each respective point, we have the Lagrangian description of the displacement vector:

$\mathbf{u}(\mathbf{X},t)=u_i{\mathbf{e}}_i$

Where ${\mathbf{e}}_i$ is the unit vector that defines the basis of the material (body-frame) coordinate system.

Expressed in terms of the material coordinates, the displacement field is:

$\mathbf{u}(\mathbf{X},t)=\mathbf{b}(\mathbf{X},t)+\mathbf{x}(\mathbf{X},t)-\mathbf{X}\text{or}{u}_i={\alpha }_{iJ}{b}_J+x_i-{\alpha }_{iJ}{X}_J$

The partial derivative of the displacement vector with respect to the material coordinates yields the material displacement gradient tensor ${\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u}$. Thus we have,

$\begin{array}{rlrl}\mathbf{u}(\mathbf{X},t)& =\mathbf{x}(\mathbf{X},t)-\mathbf{X}& \text{or}& u_i=x_i-{\delta }_{iJ}{X}_J=x_i-X_i\\ {\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u}& ={\mathrm{\nabla }}_{\mathbf{X}}\mathbf{x}-\mathbf{I}=\mathbf{F}-\mathbf{I}& \text{or}& \frac{\partial u_i}{\partial X_K}=\frac{\partial x_i}{\partial X_K}-{\delta }_{iK}=F_{iK}-{\delta }_{iK}\end{array}$

where $\mathbf{F}$ is the deformation gradient tensor.

Spatial coordinates (Eulerian description)

In the Eulerian description, the vector joining the positions of a particle $P$ in the undeformed configuration and deformed configuration is called the displacement vector:

$\mathbf{U}(\mathbf{x},t)=U_i{\mathbf{E}}_i$

Where ${\mathbf{E}}_i$ is the unit vector that defines the basis of the spatial (lab-frame) coordinate system.

Expressed in terms of spatial coordinates, the displacement field is:

$\mathbf{U}(\mathbf{x},t)=\mathbf{b}(\mathbf{x},t)+\mathbf{x}-\mathbf{X}(\mathbf{x},t)\text{or}{U}_J=b_J+{\alpha }_{Ji}{x}_i-X_J$

The partial derivative of the displacement vector with respect to the spatial coordinates yields the spatial displacement gradient tensor ${\mathrm{\nabla }}_{\mathbf{x}}\mathbf{U}$. Thus we have,

$\begin{array}{rlrl}\mathbf{U}(\mathbf{x},t)& =\mathbf{x}-\mathbf{X}(\mathbf{x},t)& \text{or}& U_J={\delta }_{Ji}{x}_i-X_J=x_J-X_J\\ {\mathrm{\nabla }}_{\mathbf{x}}\mathbf{U}& =\mathbf{I}-{\mathrm{\nabla }}_{\mathbf{x}}\mathbf{X}=\mathbf{I}-{\mathbf{F}}^{-1}& \text{or}& \frac{\partial U_J}{\partial x_k}={\delta }_{Jk}-\frac{\partial X_J}{\partial x_k}={\delta }_{Jk}-F_{Jk}^{-1}.\end{array}$

Relationship between the material and spatial coordinate systems

${\alpha }_{Ji}$ are the direction cosines between the material and spatial coordinate systems with unit vectors ${\mathbf{E}}_J$ and ${\mathbf{e}}_i$, respectively. Thus

${\mathbf{E}}_J\cdot {\mathbf{e}}_i={\alpha }_{Ji}={\alpha }_{iJ}$

The relationship between $u_i$ and $U_J$ is then given by

$u_i={\alpha }_{iJ}{U}_J\text{or}{U}_J={\alpha }_{Ji}{u}_i$

Knowing that

${\mathbf{e}}_i={\alpha }_{iJ}{\mathbf{E}}_J$

then

$\mathbf{u}(\mathbf{X},t)=u_i{\mathbf{e}}_i=u_i({\alpha }_{iJ}{\mathbf{E}}_J)=U_J{\mathbf{E}}_J=\mathbf{U}(\mathbf{x},t)$

Combining the coordinate systems of deformed and undeformed configurations

It is common to superimpose the coordinate systems for the deformed and undeformed configurations, which results in $\mathbf{b}=0$, and the direction cosines become Kronecker deltas, i.e.

${\mathbf{E}}_J\cdot {\mathbf{e}}_i={\delta }_{Ji}={\delta }_{iJ}$

Thus in material (deformed) coordinates, the displacement may be expressed as:

$\mathbf{u}(\mathbf{X},t)=\mathbf{x}(\mathbf{X},t)-\mathbf{X}\text{or}{u}_i=x_i-{\delta }_{iJ}{X}_J$

And in spatial (undeformed) coordinates, the displacement may be expressed as:

$\mathbf{U}(\mathbf{x},t)=\mathbf{x}-\mathbf{X}(\mathbf{x},t)\text{or}{U}_J={\delta }_{Ji}{x}_i-X_J$

Deformation gradient tensor

Figure 2. Deformation of a continuum body.

The deformation gradient tensor $\mathbf{F}(\mathbf{X},t)=F_{jK}{\mathbf{e}}_j\otimes {\mathbf{I}}_K$ is related to both the reference and current configuration, as seen by the unit vectors ${\mathbf{e}}_j$ and ${\mathbf{I}}_K$, therefore it is a two-point tensor.

Due to the assumption of continuity of $\chi (\mathbf{X},t)$, $\mathbf{F}$ has the inverse $\mathbf{H}={\mathbf{F}}^{-1}$, where $\mathbf{H}$ is the spatial deformation gradient tensor. Then, by the implicit function theorem, the Jacobian determinant $J(\mathbf{X},t)$ must be nonsingular, i.e. $J(\mathbf{X},t)=\det \mathbf{F}(\mathbf{X},t)\ne 0$

The material deformation gradient tensor $\mathbf{F}(\mathbf{X},t)=F_{jK}{\mathbf{e}}_j\otimes {\mathbf{I}}_K$ is a second-order tensor that represents the gradient of the mapping function or functional relation $\chi (\mathbf{X},t)$, which describes the motion of a continuum. The material deformation gradient tensor characterizes the local deformation at a material point with position vector $\mathbf{X}$, i.e. deformation at neighbouring points, by transforming (linear transformation) a material line element emanating from that point from the reference configuration to the current or deformed configuration, assuming continuity in the mapping function $\chi (\mathbf{X},t)$, i.e. differentiable function of $\mathbf{X}$ and time $t$, which implies that cracks and voids do not open or close during the deformation. Thus we have,

$\begin{array}{rlrl}d\mathbf{x}& =\frac{\partial \mathbf{x}}{\partial \mathbf{X}}d\mathbf{X}& \text{or}& d{x}_j=\frac{\partial x_j}{\partial X_K}d{X}_K\\ & =\mathrm{\nabla }\chi (\mathbf{X},t)d\mathbf{X}=\mathbf{F}(\mathbf{X},t)d\mathbf{X}& \text{or}& d{x}_j=F_{jK}d{X}_K.\end{array}$

Relative displacement vector

Consider a particle or material point $P$ with position vector $\mathbf{X}=X_I{\mathbf{I}}_I$ in the undeformed configuration (Figure 2). After a displacement of the body, the new position of the particle indicated by $p$ in the new configuration is given by the vector position $\mathbf{x}=x_i{\mathbf{e}}_i$. The coordinate systems for the undeformed and deformed configuration can be superimposed for convenience.

Consider now a material point $Q$ neighboring $P$, with position vector $\mathbf{X}+\mathrm{\Delta }\mathbf{X}=(X_I+\mathrm{\Delta }{X}_I){\mathbf{I}}_I$. In the deformed configuration this particle has a new position $q$ given by the position vector $\mathbf{x}+\mathrm{\Delta }\mathbf{x}$. Assuming that the line segments $\mathrm{\Delta }X$ and $\mathrm{\Delta }\mathbf{x}$ joining the particles $P$ and $Q$ in both the undeformed and deformed configuration, respectively, to be very small, then we can express them as $d\mathbf{X}$ and $d\mathbf{x}$. Thus from Figure 2 we have

$\begin{array}{rl}\mathbf{x}+d\mathbf{x}& =\mathbf{X}+d\mathbf{X}+\mathbf{u}(\mathbf{X}+d\mathbf{X})\\ d\mathbf{x}& =\mathbf{X}-\mathbf{x}+d\mathbf{X}+\mathbf{u}(\mathbf{X}+d\mathbf{X})\\ & =d\mathbf{X}+\mathbf{u}(\mathbf{X}+d\mathbf{X})-\mathbf{u}(\mathbf{X})\\ & =d\mathbf{X}+d\mathbf{u}\\ \end{array}$

where $\mathbf{d}\mathbf{u}$ is the relative displacement vector, which represents the relative displacement of $Q$ with respect to $P$ in the deformed configuration.

Taylor approximation

For an infinitesimal element $d\mathbf{X}$, and assuming continuity on the displacement field, it is possible to use a Taylor series expansion around point $P$, neglecting higher-order terms, to approximate the components of the relative displacement vector for the neighboring particle $Q$ as

$\begin{array}{rlrl}\mathbf{u}(\mathbf{X}+d\mathbf{X})& =\mathbf{u}(\mathbf{X})+d\mathbf{u}& \text{or}& u_i^{*}=u_i+d{u}_i\\ & \approx \mathbf{u}(\mathbf{X})+{\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u}\cdot d\mathbf{X}& \text{or}& u_i^{*}\approx u_i+\frac{\partial u_i}{\partial X_J}d{X}_J.\end{array}$

Thus, the previous equation $d\mathbf{x}=d\mathbf{X}+d\mathbf{u}$ can be written as

$\begin{array}{rl}d\mathbf{x}& =d\mathbf{X}+d\mathbf{u}\\ & =d\mathbf{X}+{\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u}\cdot d\mathbf{X}\\ & =(\mathbf{I}+{\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u})d\mathbf{X}\\ & =\mathbf{F}d\mathbf{X}\end{array}$

Time-derivative of the deformation gradient

Calculations that involve the time-dependent deformation of a body often require a time derivative of the deformation gradient to be calculated. A geometrically consistent definition of such a derivative requires an excursion into differential geometry but we avoid those issues in this article.

The time derivative of $\mathbf{F}$ is

$\dot{\mathbf{F}}=\frac{\partial \mathbf{F}}{\partial t}=\frac{\partial }{\partial t}\left[\frac{\partial \mathbf{x}(\mathbf{X},t)}{\partial \mathbf{X}}\right]=\frac{\partial }{\partial \mathbf{X}}\left[\frac{\partial \mathbf{x}(\mathbf{X},t)}{\partial t}\right]=\frac{\partial }{\partial \mathbf{X}}[\mathbf{V}(\mathbf{X},t)]$

where $\mathbf{V}$ is the velocity. The derivative on the right hand side represents a material velocity gradient. It is common to convert that into a spatial gradient, i.e.,

$\dot{\mathbf{F}}=\frac{\partial }{\partial \mathbf{X}}[\mathbf{V}(\mathbf{X},t)]=\frac{\partial }{\partial \mathbf{x}}[\mathbf{v}(\mathbf{x},t)]\cdot \frac{\partial \mathbf{x}(\mathbf{X},t)}{\partial \mathbf{X}}=\mathit{l}\cdot \mathbf{F}$

where $\mathit{l}$ is the spatial velocity gradient. If the spatial velocity gradient is constant, the above equation can be solved exactly to give

$\mathbf{F}=e^{\mathit{l}t}$

assuming $\mathbf{F}=\mathbf{1}$ at $t=0$. There are several methods of computing the exponential above.

Related quantities often used in continuum mechanics are the rate of deformation tensor and the spin tensor defined, respectively, as:

$\mathit{d}=\frac{1}{2}(\mathit{l}+{\mathit{l}}^T),\mathit{w}=\frac{1}{2}(\mathit{l}-{\mathit{l}}^T).$

The rate of deformation tensor gives the rate of stretching of line elements while the spin tensor indicates the rate of rotation or vorticity of the motion.

Transformation of a surface and volume element

To transform quantities that are defined with respect to areas in a deformed configuration to those relative to areas in a reference configuration, and vice versa, we use Nanson's relation, expressed as

$da\mathbf{n}=JdA{\mathbf{F}}^{-T}\cdot \mathbf{N}$

where $da$ is an area of a region in the deformed configuration, $dA$ is the same area in the reference configuration, and $\mathbf{n}$ is the outward normal to the area element in the current configuration while $\mathbf{N}$ is the outward normal in the reference configuration, $\mathbf{F}$ is the deformation gradient, and $J=\det \mathbf{F}$.

The corresponding formula for the transformation of the volume element is

$dv=JdV$

Derivation of Nanson's relation

To see how this formula is derived, we start with the oriented area elements in the reference and current configurations: $d\mathbf{A}=dA\mathbf{N};d\mathbf{a}=da\mathbf{n}$ The reference and current volumes of an element are $dV=d{\mathbf{A}}^T\cdot d\mathbf{L};dv=d{\mathbf{a}}^T\cdot d\mathbf{l}$ where $d\mathbf{l}=\mathbf{F}\cdot d\mathbf{L}$. Therefore, $d{\mathbf{a}}^T\cdot d\mathbf{l}=dv=JdV=Jd{\mathbf{A}}^T\cdot d\mathbf{L}$ or, $d{\mathbf{a}}^T\cdot \mathbf{F}\cdot d\mathbf{L}=dv=JdV=Jd{\mathbf{A}}^T\cdot d\mathbf{L}$ so, $d{\mathbf{a}}^T\cdot \mathbf{F}=Jd{\mathbf{A}}^T$ So we get $d\mathbf{a}=J{\mathbf{F}}^{-T}\cdot d\mathbf{A}$ or, $da\mathbf{n}=JdA{\mathbf{F}}^{-T}\cdot \mathbf{N}◻$

Polar decomposition of the deformation gradient tensor

Figure 3. Representation of the polar decomposition of the deformation gradient

The deformation gradient $\mathbf{F}$, like any second-order tensor, can be decomposed, using the polar decomposition theorem, into a product of two second-order tensors (Truesdell and Noll, 1965): an orthogonal tensor and a positive definite symmetric tensor, i.e.

$\mathbf{F}=\mathbf{R}\mathbf{U}=\mathbf{V}\mathbf{R}$

where the tensor $\mathbf{R}$ is a proper orthogonal tensor, i.e. ${\mathbf{R}}^{-1}={\mathbf{R}}^T$ and $\det \mathbf{R}=+1$, representing a rotation; the tensor $\mathbf{U}$ is the right stretch tensor; and $\mathbf{V}$ the left stretch tensor. The terms right and left means that they are to the right and left of the rotation tensor $\mathbf{R}$, respectively. $\mathbf{U}$ and $\mathbf{V}$ are both positive definite, i.e. $\mathbf{x}\cdot \mathbf{U}\cdot \mathbf{x}\ge 0$ and $\mathbf{x}\cdot \mathbf{V}\cdot \mathbf{x}\ge 0$, and symmetric tensors, i.e. $\mathbf{U}={\mathbf{U}}^T$ and $\mathbf{V}={\mathbf{V}}^T$, of second order.

This decomposition implies that the deformation of a line element $d\mathbf{X}$ in the undeformed configuration onto $d\mathbf{x}$ in the deformed configuration, i.e. $d\mathbf{x}=\mathbf{F}d\mathbf{X}$, may be obtained either by first stretching the element by $\mathbf{U}$, i.e. $d{\mathbf{x}}^{\prime }=\mathbf{U}d\mathbf{X}$, followed by a rotation $\mathbf{R}$, i.e. $d\mathbf{x}=\mathbf{R}d{\mathbf{x}}^{\prime }$; or equivalently, by applying a rigid rotation $\mathbf{R}$ first, i.e. $d{\mathbf{x}}^{\prime }=\mathbf{R}d\mathbf{X}$, followed later by a stretching $\mathbf{V}$, i.e. $d\mathbf{x}=\mathbf{V}d{\mathbf{x}}^{\prime }$ (See Figure 3).

It can be shown that,

$\mathbf{V}=\mathbf{R}\cdot \mathbf{U}\cdot {\mathbf{R}}^T$

so that $\mathbf{U}$ and $\mathbf{V}$ have the same eigenvalues or principal stretches, but different eigenvectors or principal directions ${\mathbf{N}}_i$ and ${\mathbf{n}}_i$, respectively. The principal directions are related by

${\mathbf{n}}_i=\mathbf{R}{\mathbf{N}}_i.$

This polar decomposition is unique as $\mathbf{F}$ is non-symmetric.

Deformation tensors

Several rotation-independent deformation tensors are used in mechanics. In solid mechanics, the most popular of these are the right and left Cauchy-Green deformation tensors.

Since a pure rotation should not induce any stresses in a deformable body, it is often convenient to use rotation-independent measures of deformation in continuum mechanics. As a rotation followed by its inverse rotation leads to no change ($\mathbf{R}{\mathbf{R}}^T={\mathbf{R}}^T\mathbf{R}=\mathbf{1}$) we can exclude the rotation by multiplying $\mathbf{F}$ by its transpose.

The Right Cauchy-Green deformation tensor

In 1839, George Green introduced a deformation tensor known as the right Cauchy-Green deformation tensor or Green's deformation tensor, defined as:

$\mathbf{C}={\mathbf{F}}^T\mathbf{F}={\mathbf{U}}^2\text{or}{C}_{IJ}=F_{kI}{F}_{kJ}=\frac{\partial x_k}{\partial X_I}\frac{\partial x_k}{\partial X_J}.$

Physically, the Cauchy-Green tensor gives us the square of local change in distances due to deformation, i.e. $d{\mathbf{x}}^2=d\mathbf{X}\cdot \mathbf{C}d\mathbf{X}$

Invariants of $\mathbf{C}$ are often used in the expressions for strain energy density functions. The most commonly used invariants are

$\begin{array}{rl}{I}_1^C& :=\text{tr}(\mathbf{C})=C_{II}={\lambda }_1^2+{\lambda }_2^2+{\lambda }_3^2\\ I_2^C& :=\frac{1}{2}[(\text{tr}\mathbf{C}{)}^2-\text{tr}({\mathbf{C}}^2)]=\frac{1}{2}[(C_{JJ}{)}^2-C_{IK}{C}_{KI}]={\lambda }_1^2{\lambda }_2^2+{\lambda }_2^2{\lambda }_3^2+{\lambda }_3^2{\lambda }_1^2\\ I_3^C& :=\det (\mathbf{C})={\lambda }_1^2{\lambda }_2^2{\lambda }_3^2.\end{array}$

where ${\lambda }_i$ are stretch ratios for the unit fibers that are initially oriented along the directions of three axis in the coordinate systems.

The Finger deformation tensor

The IUPAC recommends that the inverse of the right Cauchy-Green deformation tensor (called the Cauchy tensor in that document), i. e., ${\mathbf{C}}^{-1}$, be called the Finger tensor. However, that nomenclature is not universally accepted in applied mechanics.

$\mathbf{f}={\mathbf{C}}^{-1}={\mathbf{F}}^{-1}{\mathbf{F}}^{-T}\text{or}{f}_{IJ}=\frac{\partial X_I}{\partial x_k}\frac{\partial X_J}{\partial x_k}$

The Left Cauchy-Green or Finger deformation tensor

Reversing the order of multiplication in the formula for the right Green-Cauchy deformation tensor leads to the left Cauchy-Green deformation tensor which is defined as:

$\mathbf{B}=\mathbf{F}{\mathbf{F}}^T={\mathbf{V}}^2\text{or}{B}_{ij}=\frac{\partial x_i}{\partial X_K}\frac{\partial x_j}{\partial X_K}$

The left Cauchy-Green deformation tensor is often called the Finger deformation tensor, named after Josef Finger (1894).

Invariants of $\mathbf{B}$ are also used in the expressions for strain energy density functions. The conventional invariants are defined as

$\begin{array}{rl}{I}_1& :=\text{tr}(\mathbf{B})=B_{ii}={\lambda }_1^2+{\lambda }_2^2+{\lambda }_3^2\\ I_2& :=\frac{1}{2}[(\text{tr}\mathbf{B}{)}^2-\text{tr}({\mathbf{B}}^2)]=\frac{1}{2}(B_{ii}^2-B_{jk}{B}_{kj})={\lambda }_1^2{\lambda }_2^2+{\lambda }_2^2{\lambda }_3^2+{\lambda }_3^2{\lambda }_1^2\\ I_3& :=\det \mathbf{B}=J^2={\lambda }_1^2{\lambda }_2^2{\lambda }_3^2\end{array}$

where $J:=\det \mathbf{F}$ is the determinant of the deformation gradient.

For nearly incompressible materials, a slightly different set of invariants is used:

$({\overline{I}}_1:=J^{-2/3}{I}_1;{\overline{I}}_2:=J^{-4/3}{I}_2;J=1).$

The Cauchy deformation tensor

Earlier in 1828, Augustin Louis Cauchy introduced a deformation tensor defined as the inverse of the left Cauchy-Green deformation tensor, ${\mathbf{B}}^{-1}$. This tensor has also been called the Piola tensor and the Finger tensor in the rheology and fluid dynamics literature.

$\mathbf{c}={\mathbf{B}}^{-1}={\mathbf{F}}^{-T}{\mathbf{F}}^{-1}\text{or}{c}_{ij}=\frac{\partial X_K}{\partial x_i}\frac{\partial X_K}{\partial x_j}$

Spectral representation

If there are three distinct principal stretches ${\lambda }_i$, the spectral decompositions of $\mathbf{C}$ and $\mathbf{B}$ is given by

$\mathbf{C}=\sum _{i=1}^3{\lambda }_i^2{\mathbf{N}}_i\otimes {\mathbf{N}}_i\text{and}\mathbf{B}=\sum _{i=1}^3{\lambda }_i^2{\mathbf{n}}_i\otimes {\mathbf{n}}_i$

Furthermore,

$\mathbf{U}=\sum _{i=1}^3{\lambda }_i{\mathbf{N}}_i\otimes {\mathbf{N}}_i;\mathbf{V}=\sum _{i=1}^3{\lambda }_i{\mathbf{n}}_i\otimes {\mathbf{n}}_i$

$\mathbf{R}=\sum _{i=1}^3{\mathbf{n}}_i\otimes {\mathbf{N}}_i;\mathbf{F}=\sum _{i=1}^3{\lambda }_i{\mathbf{n}}_i\otimes {\mathbf{N}}_i$

Observe that

$\mathbf{V}=\mathbf{R}\mathbf{U}{\mathbf{R}}^T=\sum _{i=1}^3{\lambda }_i\mathbf{R}({\mathbf{N}}_i\otimes {\mathbf{N}}_i){\mathbf{R}}^T=\sum _{i=1}^3{\lambda }_i(\mathbf{R}{\mathbf{N}}_i)\otimes (\mathbf{R}{\mathbf{N}}_i)$

Therefore the uniqueness of the spectral decomposition also implies that ${\mathbf{n}}_i=\mathbf{R}{\mathbf{N}}_i$. The left stretch ($\mathbf{V}$) is also called the spatial stretch tensor while the right stretch ($\mathbf{U}$) is called the material stretch tensor.

The effect of $\mathbf{F}$ acting on ${\mathbf{N}}_i$ is to stretch the vector by ${\lambda }_i$ and to rotate it to the new orientation ${\mathbf{n}}_i$, i.e.,

$\mathbf{F}{\mathbf{N}}_i={\lambda }_i(\mathbf{R}{\mathbf{N}}_i)={\lambda }_i{\mathbf{n}}_i$

In a similar vein,

${\mathbf{F}}^{-T}{\mathbf{N}}_i=\frac{1}{{\lambda }_i}{\mathbf{n}}_i;{\mathbf{F}}^T{\mathbf{n}}_i={\lambda }_i{\mathbf{N}}_i;{\mathbf{F}}^{-1}{\mathbf{n}}_i=\frac{1}{{\lambda }_i}{\mathbf{N}}_i.$

Examples

Uniaxial extension of an incompressible material This is the case where a specimen is stretched in 1-direction with a stretch ratio of $\mathit{\alpha }={\mathit{\alpha }}_{\mathbf{1}}$. If the volume remains constant, the contraction in the other two directions is such that ${\mathit{\alpha }}_{\mathbf{1}}{\mathit{\alpha }}_{\mathbf{2}}{\mathit{\alpha }}_{\mathbf{3}}=\mathbf{1}$ or ${\mathit{\alpha }}_{\mathbf{2}}={\mathit{\alpha }}_{\mathbf{3}}={\mathit{\alpha }}^{\mathbf{-}\mathbf{0}.\mathbf{5}}$. Then: $\mathbf{F}=\left[\begin{array}{ccc}\alpha & 0& 0\\ 0& {\alpha }^{-0.5}& 0\\ 0& 0& {\alpha }^{-0.5}\end{array}\right]$ $\mathbf{B}=\mathbf{C}=\left[\begin{array}{ccc}{\alpha }^2& 0& 0\\ 0& {\alpha }^{-1}& 0\\ 0& 0& {\alpha }^{-1}\end{array}\right]$ Simple shear $\mathbf{F}=\left[\begin{array}{ccc}1& \gamma & 0\\ 0& 1& 0\\ 0& 0& 1\end{array}\right]$ $\mathbf{B}=\left[\begin{array}{ccc}1+{\gamma }^2& \gamma & 0\\ \gamma & 1& 0\\ 0& 0& 1\end{array}\right]$ $\mathbf{C}=\left[\begin{array}{ccc}1& \gamma & 0\\ \gamma & 1+{\gamma }^2& 0\\ 0& 0& 1\end{array}\right]$ Rigid body rotation $\mathbf{F}=\left[\begin{array}{ccc}\cos \theta & \sin \theta & 0\\ -\sin \theta & \cos \theta & 0\\ 0& 0& 1\end{array}\right]$ $\mathbf{B}=\mathbf{C}=\left[\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\end{array}\right]=\mathbf{1}$

Derivatives of stretch

Derivatives of the stretch with respect to the right Cauchy-Green deformation tensor are used to derive the stress-strain relations of many solids, particularly hyperelastic materials. These derivatives are

$\frac{\partial {\lambda }_i}{\partial \mathbf{C}}=\frac{1}{2{\lambda }_i}{\mathbf{N}}_i\otimes {\mathbf{N}}_i=\frac{1}{2{\lambda }_i}{\mathbf{R}}^T({\mathbf{n}}_i\otimes {\mathbf{n}}_i)\mathbf{R};i=1,2,3$

and follow from the observations that

$\mathbf{C}:({\mathbf{N}}_i\otimes {\mathbf{N}}_i)={\lambda }_i^2;\frac{\partial \mathbf{C}}{\partial \mathbf{C}}={\mathsf{I}}^{(s)};{\mathsf{I}}^{(s)}:({\mathbf{N}}_i\otimes {\mathbf{N}}_i)={\mathbf{N}}_i\otimes {\mathbf{N}}_i.$

Physical interpretation of deformation tensors

Let $\mathbf{X}=X^i{\mathit{E}}_i$ be a Cartesian coordinate system defined on the undeformed body and let $\mathbf{x}=x^i{\mathit{E}}_i$ be another system defined on the deformed body. Let a curve $\mathbf{X}(s)$ in the undeformed body be parametrized using $s\in [0,1]$. Its image in the deformed body is $\mathbf{x}(\mathbf{X}(s))$.

The undeformed length of the curve is given by

$l_X=\underset{0}{\overset{1}{\int }}|\frac{d\mathbf{X}}{ds}\cdot \frac{d\mathbf{X}}{ds}|ds=\underset{0}{\overset{1}{\int }}|\frac{d\mathbf{X}}{ds}\cdot \mathit{I}\cdot \frac{d\mathbf{X}}{ds}|ds$

After deformation, the length becomes

$\begin{array}{rl}{l}_x& =\underset{0}{\overset{1}{\int }}|\frac{d\mathbf{x}}{ds}\cdot \frac{d\mathbf{x}}{ds}|ds=\underset{0}{\overset{1}{\int }}|(\frac{d\mathbf{x}}{d\mathbf{X}}\cdot \frac{d\mathbf{X}}{ds})\cdot (\frac{d\mathbf{x}}{d\mathbf{X}}\cdot \frac{d\mathbf{X}}{ds})|ds\\ & =\underset{0}{\overset{1}{\int }}|\frac{d\mathbf{X}}{ds}\cdot [{\left(\frac{d\mathbf{x}}{d\mathbf{X}}\right)}^T\cdot \frac{d\mathbf{x}}{d\mathbf{X}}]\cdot \frac{d\mathbf{X}}{ds}|ds\end{array}$

Note that the right Cauchy-Green deformation tensor is defined as

$\mathit{C}:={\mathit{F}}^T\cdot \mathit{F}={\left(\frac{d\mathbf{x}}{d\mathbf{X}}\right)}^T\cdot \frac{d\mathbf{x}}{d\mathbf{X}}$

Hence,

$l_x=\underset{0}{\overset{1}{\int }}|\frac{d\mathbf{X}}{ds}\cdot \mathit{C}\cdot \frac{d\mathbf{X}}{ds}|ds$

which indicates that changes in length are characterized by $\mathit{C}$.

Finite strain tensors

The concept of strain is used to evaluate how much a given displacement differs locally from a rigid body displacement. One of such strains for large deformations is the Lagrangian finite strain tensor, also called the Green-Lagrangian strain tensor or Green – St-Venant strain tensor, defined as

$\mathbf{E}=\frac{1}{2}(\mathbf{C}-\mathbf{I})\text{or}{E}_{KL}=\frac{1}{2}(\frac{\partial x_j}{\partial X_K}\frac{\partial x_j}{\partial X_L}-{\delta }_{KL})$

or as a function of the displacement gradient tensor

$\mathbf{E}=\frac{1}{2}[({\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u}{)}^T+{\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u}+({\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u}{)}^T\cdot {\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u}]$

or

$E_{KL}=\frac{1}{2}(\frac{\partial U_K}{\partial X_L}+\frac{\partial U_L}{\partial X_K}+\frac{\partial U_M}{\partial X_K}\frac{\partial U_M}{\partial X_L})$

The Green-Lagrangian strain tensor is a measure of how much $\mathbf{C}$ differs from $\mathbf{I}$.

The Eulerian-Almansi finite strain tensor, referenced to the deformed configuration, i.e. Eulerian description, is defined as

$\mathbf{e}=\frac{1}{2}(\mathbf{I}-\mathbf{c})\text{or}{e}_{rs}=\frac{1}{2}({\delta }_{rs}-\frac{\partial X_M}{\partial x_r}\frac{\partial X_M}{\partial x_s})$

or as a function of the displacement gradients we have

$e_{ij}=\frac{1}{2}(\frac{\partial u_i}{\partial x_j}+\frac{\partial u_j}{\partial x_i}-\frac{\partial u_k}{\partial x_i}\frac{\partial u_k}{\partial x_j})$

Derivation of the Lagrangian and Eulerian finite strain tensors

A measure of deformation is the difference between the squares of the differential line element $d\mathbf{X}$, in the undeformed configuration, and $d\mathbf{x}$, in the deformed configuration (Figure 2). Deformation has occurred if the difference is non zero, otherwise a rigid-body displacement has occurred. Thus we have, $d{\mathbf{x}}^2-d{\mathbf{X}}^2=d\mathbf{x}\cdot d\mathbf{x}-d\mathbf{X}\cdot d\mathbf{X}\text{or}(dx{)}^2-(dX{)}^2=d{x}_{j}d{x}_j-d{X}_{M}d{X}_M$ In the Lagrangian description, using the material coordinates as the frame of reference, the linear transformation between the differential lines is $d\mathbf{x}=\frac{\partial \mathbf{x}}{\partial \mathbf{X}}d\mathbf{X}=\mathbf{F}d\mathbf{X}\text{or}d{x}_j=\frac{\partial x_j}{\partial X_M}d{X}_M$ Then we have, $\begin{array}{rl}d{\mathbf{x}}^2& =d\mathbf{x}\cdot d\mathbf{x}\\ & =\mathbf{F}\cdot d\mathbf{X}\cdot \mathbf{F}\cdot d\mathbf{X}\\ & =d\mathbf{X}\cdot {\mathbf{F}}^T\mathbf{F}\cdot d\mathbf{X}\\ & =d\mathbf{X}\cdot \mathbf{C}\cdot d\mathbf{X}\end{array}\text{or}\begin{array}{rl}(dx{)}^2& =d{x}_{j}d{x}_j\\ & =\frac{\partial x_j}{\partial X_K}\frac{\partial x_j}{\partial X_L}d{X}_{K}d{X}_L\\ & =C_{KL}d{X}_{K}d{X}_L\\ \end{array}$ where $C_{KL}$ are the components of the right Cauchy-Green deformation tensor, $\mathbf{C}={\mathbf{F}}^T\mathbf{F}$. Then, replacing this equation into the first equation we have, $\begin{array}{rl}d{\mathbf{x}}^2-d{\mathbf{X}}^2& =d\mathbf{X}\cdot \mathbf{C}\cdot d\mathbf{X}-d\mathbf{X}\cdot d\mathbf{X}\\ & =d\mathbf{X}\cdot (\mathbf{C}-\mathbf{I})\cdot d\mathbf{X}\\ & =d\mathbf{X}\cdot 2\mathbf{E}\cdot d\mathbf{X}\\ \end{array}$ or $\begin{array}{rl}(dx{)}^2-(dX{)}^2& =\frac{\partial x_j}{\partial X_K}\frac{\partial x_j}{\partial X_L}d{X}_{K}d{X}_L-d{X}_{M}d{X}_M\\ & =(\frac{\partial x_j}{\partial X_K}\frac{\partial x_j}{\partial X_L}-{\delta }_{KL})d{X}_{K}d{X}_L\\ & =2E_{KL}d{X}_{K}d{X}_L\end{array}$ where $E_{KL}$, are the components of a second-order tensor called the Green – St-Venant strain tensor or the Lagrangian finite strain tensor, $\mathbf{E}=\frac{1}{2}(\mathbf{C}-\mathbf{I})\text{or}{E}_{KL}=\frac{1}{2}(\frac{\partial x_j}{\partial X_K}\frac{\partial x_j}{\partial X_L}-{\delta }_{KL})$ In the Eulerian description, using the spatial coordinates as the frame of reference, the linear transformation between the differential lines is $d\mathbf{X}=\frac{\partial \mathbf{X}}{\partial \mathbf{x}}d\mathbf{x}={\mathbf{F}}^{-1}d\mathbf{x}=\mathbf{H}d\mathbf{x}\text{or}d{X}_M=\frac{\partial X_M}{\partial x_n}d{x}_n$ where $\frac{\partial X_M}{\partial x_n}$ are the components of the spatial deformation gradient tensor, $\mathbf{H}$. Thus we have $\begin{array}{rl}d{\mathbf{X}}^2& =d\mathbf{X}\cdot d\mathbf{X}\\ & ={\mathbf{F}}^{-1}\cdot d\mathbf{x}\cdot {\mathbf{F}}^{-1}\cdot d\mathbf{x}\\ & =d\mathbf{x}\cdot {\mathbf{F}}^{-T}{\mathbf{F}}^{-1}\cdot d\mathbf{x}\\ & =d\mathbf{x}\cdot \mathbf{c}\cdot d\mathbf{x}\end{array}\text{or}\begin{array}{rl}(dX{)}^2& =d{X}_{M}d{X}_M\\ & =\frac{\partial X_M}{\partial x_r}\frac{\partial X_M}{\partial x_s}d{x}_{r}d{x}_s\\ & =c_{rs}d{x}_{r}d{x}_s\\ \end{array}$ where the second order tensor $c_{rs}$ is called Cauchy's deformation tensor, $\mathbf{c}={\mathbf{F}}^{-T}{\mathbf{F}}^{-1}$. Then we have, $\begin{array}{rl}d{\mathbf{x}}^2-d{\mathbf{X}}^2& =d\mathbf{x}\cdot d\mathbf{x}-d\mathbf{x}\cdot \mathbf{c}\cdot d\mathbf{x}\\ & =d\mathbf{x}\cdot (\mathbf{I}-\mathbf{c})\cdot d\mathbf{x}\\ & =d\mathbf{x}\cdot 2\mathbf{e}\cdot d\mathbf{x}\\ \end{array}$ or $\begin{array}{rl}(dx{)}^2-(dX{)}^2& =d{x}_{j}d{x}_j-\frac{\partial X_M}{\partial x_r}\frac{\partial X_M}{\partial x_s}d{x}_{r}d{x}_s\\ & =({\delta }_{rs}-\frac{\partial X_M}{\partial x_r}\frac{\partial X_M}{\partial x_s})d{x}_{r}d{x}_s\\ & =2e_{rs}d{x}_{r}d{x}_s\end{array}$ where $e_{rs}$, are the components of a second-order tensor called the Eulerian-Almansi finite strain tensor, $\mathbf{e}=\frac{1}{2}(\mathbf{I}-\mathbf{c})\text{or}{e}_{rs}=\frac{1}{2}({\delta }_{rs}-\frac{\partial X_M}{\partial x_r}\frac{\partial X_M}{\partial x_s})$ Both the Lagrangian and Eulerian finite strain tensors can be conveniently expressed in terms of the displacement gradient tensor. For the Lagrangian strain tensor, first we differentiate the displacement vector $\mathbf{u}(\mathbf{X},t)$ with respect to the material coordinates $X_M$ to obtain the material displacement gradient tensor, ${\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u}$ $\begin{array}{rl}\mathbf{u}(\mathbf{X},t)& =\mathbf{x}(\mathbf{X},t)-\mathbf{X}\\ {\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u}& =\mathbf{F}-\mathbf{I}\\ \mathbf{F}& ={\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u}+\mathbf{I}\\ \end{array}\text{or}\begin{array}{rl}{u}_i& =x_i-{\delta }_{iJ}{X}_J\\ {\delta }_{iJ}{U}_J& =x_i-{\delta }_{iJ}{X}_J\\ x_i& ={\delta }_{iJ}(U_J+X_J)\\ \frac{\partial x_i}{\partial X_K}& ={\delta }_{iJ}(\frac{\partial U_J}{\partial X_K}+{\delta }_{JK})\\ \end{array}$ Replacing this equation into the expression for the Lagrangian finite strain tensor we have $\begin{array}{rl}\mathbf{E}& =\frac{1}{2}({\mathbf{F}}^T\mathbf{F}-\mathbf{I})\\ & =\frac{1}{2}[\{({\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u}{)}^T+\mathbf{I}\}({\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u}+\mathbf{I})-\mathbf{I}]\\ & =\frac{1}{2}[({\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u}{)}^T+{\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u}+({\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u}{)}^T\cdot {\mathrm{\nabla }}_{\mathbf{X}}\mathbf{u}]\\ \end{array}$ or $\begin{array}{rl}{E}_{KL}& =\frac{1}{2}(\frac{\partial x_j}{\partial X_K}\frac{\partial x_j}{\partial X_L}-{\delta }_{KL})\\ & =\frac{1}{2}[{\delta }_{jM}(\frac{\partial U_M}{\partial X_K}+{\delta }_{MK}){\delta }_{jN}(\frac{\partial U_N}{\partial X_L}+{\delta }_{NL})-{\delta }_{KL}]\\ & =\frac{1}{2}[{\delta }_{MN}(\frac{\partial U_M}{\partial X_K}+{\delta }_{MK})(\frac{\partial U_N}{\partial X_L}+{\delta }_{NL})-{\delta }_{KL}]\\ & =\frac{1}{2}[(\frac{\partial U_M}{\partial X_K}+{\delta }_{MK})(\frac{\partial U_M}{\partial X_L}+{\delta }_{ML})-{\delta }_{KL}]\\ & =\frac{1}{2}(\frac{\partial U_K}{\partial X_L}+\frac{\partial U_L}{\partial X_K}+\frac{\partial U_M}{\partial X_K}\frac{\partial U_M}{\partial X_L})\end{array}$ Similarly, the Eulerian-Almansi finite strain tensor can be expressed as $e_{ij}=\frac{1}{2}(\frac{\partial u_i}{\partial x_j}+\frac{\partial u_j}{\partial x_i}-\frac{\partial u_k}{\partial x_i}\frac{\partial u_k}{\partial x_j})$

Seth-Hill family of generalized strain tensors

B. R. Seth from the Indian Institute of Technology, Kharagpur was the first to show that the Green and Almansi strain tensors are special cases of a more general strain measure. The idea was further expanded upon by Rodney Hill in 1968. The Seth-Hill family of strain measures can be expressed as

${\mathbf{E}}_{(m)}=\frac{1}{2m}({\mathbf{U}}^{2m}-\mathbf{I})=\frac{1}{2m}[{\mathbf{C}}^m-\mathbf{I}]$

For different values of $m$ we have:

$\begin{array}{rlr}{\mathbf{E}}_{(1)}& =\frac{1}{2}({\mathbf{U}}^2-\mathbf{I})=\frac{1}{2}(\mathbf{C}-\mathbf{I})& \text{Green-Lagrangian strain tensor}\\ {\mathbf{E}}_{(1/2)}& =(\mathbf{U}-\mathbf{I})={\mathbf{C}}^{1/2}-\mathbf{I}& \text{Biot strain tensor}\\ {\mathbf{E}}_{(0)}& =\ln \mathbf{U}=\frac{1}{2}\ln \mathbf{C}& \text{Logarithmic strain, Natural strain, True strain, or Hencky strain}\\ {\mathbf{E}}_{(-1)}& =\frac{1}{2}[\mathbf{I}-{\mathbf{U}}^{-2}]& \text{Almansi strain}\end{array}$

Stretch ratio

The stretch ratio is a measure of the extensional or normal strain of a differential line element, which can be defined at either the undeformed configuration or the deformed configuration.

The stretch ratio for the differential element $d\mathbf{X}=dX\mathbf{N}$ (Figure) in the direction of the unit vector $\mathbf{N}$ at the material point $P$, in the undeformed configuration, is defined as

${\mathrm{\Lambda }}_{(\mathbf{N})}=\frac{dx}{dX}$

where $dx$ is the deformed magnitude of the differential element $d\mathbf{X}$.

Similarly, the stretch ratio for the differential element $d\mathbf{x}=dx\mathbf{n}$ (Figure), in the direction of the unit vector $\mathbf{n}$ at the material point $p$, in the deformed configuration, is defined as

$\frac{1}{{\mathrm{\Lambda }}_{(\mathbf{n})}}=\frac{dX}{dx}.$

The normal strain $e_{\mathbf{N}}$ in any direction $\mathbf{N}$ can be expressed as a function of the stretch ratio,

$e_{(\mathbf{N})}=\frac{dx-dX}{dX}={\mathrm{\Lambda }}_{(\mathbf{N})}-1.$

This equation implies that the normal strain is zero, i.e. no deformation, when the stretch is equal to unity. Some materials, such as elastometers can sustain stretch ratios of 3 or 4 before they fail, whereas traditional engineering materials, such as concrete or steel, fail at much lower stretch ratios, perhaps of the order of 1.001 (reference?)

Physical interpretation of the finite strain tensor

The diagonal components $E_{KL}$ of the Lagrangian finite strain tensor are related to the normal strain, e.g.

$E_{11}=e_{({\mathbf{I}}_1)}+\frac{1}{2}{e}_{({\mathbf{I}}_1)}^2$

where $e_{({\mathbf{I}}_1)}$ is the normal strain or engineering strain in the direction ${\mathbf{I}}_1$.

The off-diagonal components $E_{KL}$ of the Lagrangian finite strain tensor are related to shear strain, e.g.

$E_{12}=\frac{1}{2}\sqrt{2E_{11}+1}\sqrt{2E_{22}+1}\sin {\varphi }_{12}$

where ${\varphi }_{12}$ is the change in the angle between two line elements that were originally perpendicular with directions ${\mathbf{I}}_1$ and ${\mathbf{I}}_2$, respectively.

Under certain circumstances, i.e. small displacements and small displacement rates, the components of the Lagrangian finite strain tensor may be approximated by the components of the infinitesimal strain tensor

Derivation of the physical interpretation of the Lagrangian and Eulerian finite strain tensors

The stretch ratio for the differential element $d\mathbf{X}=dX\mathbf{N}$ (Figure) in the direction of the unit vector $\mathbf{N}$ at the material point $P$, in the undeformed configuration, is defined as ${\mathrm{\Lambda }}_{(\mathbf{N})}=\frac{dx}{dX}$ where $dx$ is the deformed magnitude of the differential element $d\mathbf{X}$. Similarly, the stretch ratio for the differential element $d\mathbf{x}=dx\mathbf{n}$ (Figure), in the direction of the unit vector $\mathbf{n}$ at the material point $p$, in the deformed configuration, is defined as $\frac{1}{{\mathrm{\Lambda }}_{(\mathbf{n})}}=\frac{dX}{dx}$ The square of the stretch ratio is defined as ${\mathrm{\Lambda }}_{(\mathbf{N})}^2={\left(\frac{dx}{dX}\right)}^2$ Knowing that $(dx{)}^2=C_{KL}d{X}_{K}d{X}_L$ we have ${\mathrm{\Lambda }}_{(\mathbf{N})}^2=C_{KL}{N}_K{N}_L$ where $N_K$ and $N_L$ are unit vectors. The normal strain or engineering strain $e_{\mathbf{N}}$ in any direction $\mathbf{N}$ can be expressed as a function of the stretch ratio, $e_{(\mathbf{N})}=\frac{dx-dX}{dX}={\mathrm{\Lambda }}_{(\mathbf{N})}-1$ Thus, the normal strain in the direction ${\mathbf{I}}_1$ at the material point $P$ may be expressed in terms of the stretch ratio as $\begin{array}{rl}{e}_{({\mathbf{I}}_1)}=\frac{d{x}_1-d{X}_1}{d{X}_1}& ={\mathrm{\Lambda }}_{({\mathbf{I}}_1)}-1\\ & =\sqrt{C_{11}}-1=\sqrt{{\delta }_{11}+2E_{11}}-1\\ & =\sqrt{1+2E_{11}}-1\end{array}$ solving for $E_{11}$ we have $\begin{array}{rl}2E_{11}& =\frac{(d{x}_1{)}^2-(d{X}_1{)}^2}{(d{X}_1{)}^2}\\ E_{11}& =\left(\frac{d{x}_1-d{X}_1}{d{X}_1}\right)+\frac{1}{2}{\left(\frac{d{x}_1-d{X}_1}{d{X}_1}\right)}^2\\ & =e_{({\mathbf{I}}_1)}+\frac{1}{2}{e}_{({\mathbf{I}}_1)}^2\end{array}$ The shear strain, or change in angle between two line elements $d{\mathbf{X}}_1$ and $d{\mathbf{X}}_2$ initially perpendicular, and oriented in the principal directions ${\mathbf{I}}_1$ and ${\mathbf{I}}_2$, respectivelly, can also be expressed as a function of the stretch ratio. From the dot product between the deformed lines $d{\mathbf{x}}_1$ and $d{\mathbf{x}}_2$ we have $\begin{array}{rl}d{\mathbf{x}}_1\cdot d{\mathbf{x}}_2& =d{x}_{1}d{x}_2\cos {\theta }_{12}\\ \mathbf{F}\cdot d{\mathbf{X}}_1\cdot \mathbf{F}\cdot d{\mathbf{X}}_2& =\sqrt{d{\mathbf{X}}_1\cdot {\mathbf{F}}^T\cdot \mathbf{F}\cdot d{\mathbf{X}}_1}\cdot \sqrt{d{\mathbf{X}}_2\cdot {\mathbf{F}}^T\cdot \mathbf{F}\cdot d{\mathbf{X}}_2}\cos {\theta }_{12}\\ \frac{d{\mathbf{X}}_1\cdot {\mathbf{F}}^T\cdot \mathbf{F}\cdot d{\mathbf{X}}_2}{d{X}_{1}d{X}_2}& =\frac{\sqrt{d{\mathbf{X}}_1\cdot {\mathbf{F}}^T\cdot \mathbf{F}\cdot d{\mathbf{X}}_1}\cdot \sqrt{d{\mathbf{X}}_2\cdot {\mathbf{F}}^T\cdot \mathbf{F}\cdot d{\mathbf{X}}_2}}{d{X}_{1}d{X}_2}\cos {\theta }_{12}\\ {\mathbf{I}}_1\cdot \mathbf{C}\cdot {\mathbf{I}}_2& ={\mathrm{\Lambda }}_{{\mathbf{I}}_1}{\mathrm{\Lambda }}_{{\mathbf{I}}_2}\cos {\theta }_{12}\\ \end{array}$ where ${\theta }_{12}$ is the angle between the lines $d{\mathbf{x}}_1$ and $d{\mathbf{x}}_2$ in the deformed configuration. Defining ${\varphi }_{12}$ as the shear strain or reduction in the angle between two line elements that were originally perpendicular, we have ${\varphi }_{12}=\frac{\pi }{2}-{\theta }_{12}$ thus, $\cos {\theta }_{12}=\sin {\varphi }_{12}$ then ${\mathbf{I}}_1\cdot \mathbf{C}\cdot {\mathbf{I}}_2={\mathrm{\Lambda }}_{{\mathbf{I}}_1}{\mathrm{\Lambda }}_{{\mathbf{I}}_2}\sin {\varphi }_{12}$ or $\begin{array}{rl}{C}_{12}& =\sqrt{C_{11}}\sqrt{C_{22}}\sin {\varphi }_{12}\\ 2E_{12}+{\delta }_{12}& =\sqrt{2E_{11}+1}\sqrt{2E_{22}+1}\sin {\varphi }_{12}\\ E_{12}& =\frac{1}{2}\sqrt{2E_{11}+1}\sqrt{2E_{22}+1}\sin {\varphi }_{12}\end{array}$

Deformation tensors in curvilinear coordinates

A representation of deformation tensors in curvilinear coordinates is useful for many problems in continuum mechanics such as nonlinear shell theories and large plastic deformations. Let $\mathbf{x}=\mathbf{x}({\xi }^1,{\xi }^2,{\xi }^3)$ be a given deformation where the space is characterized by the coordinates $({\xi }^1,{\xi }^2,{\xi }^3)$. The tangent vector to the coordinate curve ${\xi }^i$ at $\mathbf{x}$ is given by

${\mathbf{g}}_i=\frac{\partial \mathbf{x}}{\partial {\xi }^i}$

The three tangent vectors at $\mathbf{x}$ form a basis. These vectors are related the reciprocal basis vectors by

${\mathbf{g}}_i\cdot {\mathbf{g}}^j={\delta }_i^j$

Let us define a second-order tensor field $\mathit{g}$ (also called the metric tensor) with components

$g_{ij}:=\frac{\partial \mathbf{x}}{\partial {\xi }^i}\cdot \frac{\partial \mathbf{x}}{\partial {\xi }^j}={\mathbf{g}}_i\cdot {\mathbf{g}}_j$

The Christoffel symbols of the first kind can be expressed as

${\mathrm{\Gamma }}_{ijk}=\frac{1}{2}[({\mathbf{g}}_i\cdot {\mathbf{g}}_k{)}_{,j}+({\mathbf{g}}_j\cdot {\mathbf{g}}_k{)}_{,i}-({\mathbf{g}}_i\cdot {\mathbf{g}}_j{)}_{,k}]$

To see how the Christoffel symbols are related to the Right Cauchy-Green deformation tensor let us define two sets of bases

${\mathbf{G}}_i:=\frac{\partial \mathbf{X}}{\partial {\xi }^i};{\mathbf{G}}_i\cdot {\mathbf{G}}^j={\delta }_i^j;{\mathbf{g}}_i:=\frac{\partial \mathbf{x}}{\partial {\xi }^i};{\mathbf{g}}_i\cdot {\mathbf{g}}^j={\delta }_i^j$

The deformation gradient in curvilinear coordinates

Using the definition of the gradient of a vector field in curvilinear coordinates, the deformation gradient can be written as

$\mathit{F}={\mathrm{\nabla }}_{\mathbf{X}}\mathbf{x}=\frac{\partial \mathbf{x}}{\partial {\xi }^i}\otimes {\mathbf{G}}^i={\mathbf{g}}_i\otimes {\mathbf{G}}^i$

The right Cauchy-Green tensor in curvilinear coordinates

The right Cauchy-Green deformation tensor is given by

$\mathit{C}={\mathit{F}}^T\cdot \mathit{F}=({\mathbf{G}}^i\otimes {\mathbf{g}}_i)\cdot ({\mathbf{g}}_i\otimes {\mathbf{G}}^i)=({\mathbf{g}}_i\cdot {\mathbf{g}}_j)({\mathbf{G}}^i\otimes {\mathbf{G}}^j)$

If we express $\mathit{C}$ in terms of components with respect to the basis {${\mathbf{G}}^i$} we have

$\mathit{C}=C_{ij}{\mathbf{G}}^i\otimes {\mathbf{G}}^j$

Therefore

$C_{ij}={\mathbf{g}}_i\cdot {\mathbf{g}}_j=g_{ij}$

and the Christoffel symbol of the first kind may be written in the following form.

${\mathrm{\Gamma }}_{ijk}=\frac{1}{2}[C_{ik,j}+C_{jk,i}-C_{ij,k}]=\frac{1}{2}[({\mathbf{G}}_i\cdot \mathit{C}\cdot {\mathbf{G}}_k{)}_{,j}+({\mathbf{G}}_j\cdot \mathit{C}\cdot {\mathbf{G}}_k{)}_{,i}-({\mathbf{G}}_i\cdot \mathit{C}\cdot {\mathbf{G}}_j{)}_{,k}]$

Some relations between deformation measures and Christoffel symbols

Let us consider a one-to-one mapping from $\mathbf{X}=\{X^1,X^2,X^3\}$ to $\mathbf{x}=\{x^1,x^2,x^3\}$ and let us assume that there exist two positive definite, symmetric second-order tensor fields $\mathit{G}$ and $\mathit{g}$ that satisfy

$G_{ij}=\frac{\partial X^{\alpha }}{\partial x^i}\frac{\partial X^{\beta }}{\partial x^j}{g}_{\alpha \beta }$

Then,

$\frac{\partial G_{ij}}{\partial x^k}=(\frac{{\partial }^2{X}^{\alpha }}{\partial x^i\partial x^k}\frac{\partial X^{\beta }}{\partial x^j}+\frac{\partial X^{\alpha }}{\partial x^i}\frac{{\partial }^2{X}^{\beta }}{\partial x^j\partial x^k})g_{\alpha \beta }+\frac{\partial X^{\alpha }}{\partial x^i}\frac{\partial X^{\beta }}{\partial x^j}\frac{\partial g_{\alpha \beta }}{\partial x^k}$

Noting that

$\frac{\partial g_{\alpha \beta }}{\partial x^k}=\frac{\partial X^{\gamma }}{\partial x^k}\frac{\partial g_{\alpha \beta }}{\partial X^{\gamma }}$

and $g_{\alpha \beta }=g_{\beta \alpha }$ we have

$\begin{array}{rl}\frac{\partial G_{ij}}{\partial x^k}& =(\frac{{\partial }^2{X}^{\alpha }}{\partial x^i\partial x^k}\frac{\partial X^{\beta }}{\partial x^j}+\frac{{\partial }^2{X}^{\alpha }}{\partial x^j\partial x^k}\frac{\partial X^{\beta }}{\partial x^i})g_{\alpha \beta }+\frac{\partial X^{\alpha }}{\partial x^i}\frac{\partial X^{\beta }}{\partial x^j}\frac{\partial X^{\gamma }}{\partial x^k}\frac{\partial g_{\alpha \beta }}{\partial X^{\gamma }}\\ \frac{\partial G_{ik}}{\partial x^j}& =(\frac{{\partial }^2{X}^{\alpha }}{\partial x^i\partial x^j}\frac{\partial X^{\beta }}{\partial x^k}+\frac{{\partial }^2{X}^{\alpha }}{\partial x^j\partial x^k}\frac{\partial X^{\beta }}{\partial x^i})g_{\alpha \beta }+\frac{\partial X^{\alpha }}{\partial x^i}\frac{\partial X^{\beta }}{\partial x^k}\frac{\partial X^{\gamma }}{\partial x^j}\frac{\partial g_{\alpha \beta }}{\partial X^{\gamma }}\\ \frac{\partial G_{jk}}{\partial x^i}& =(\frac{{\partial }^2{X}^{\alpha }}{\partial x^i\partial x^j}\frac{\partial X^{\beta }}{\partial x^k}+\frac{{\partial }^2{X}^{\alpha }}{\partial x^i\partial x^k}\frac{\partial X^{\beta }}{\partial x^j})g_{\alpha \beta }+\frac{\partial X^{\alpha }}{\partial x^j}\frac{\partial X^{\beta }}{\partial x^k}\frac{\partial X^{\gamma }}{\partial x^i}\frac{\partial g_{\alpha \beta }}{\partial X^{\gamma }}\end{array}$

Define

$\begin{array}{rl}{(x)}_{}{\mathrm{\Gamma }}_{ijk}& :=\frac{1}{2}(\frac{\partial G_{ik}}{\partial x^j}+\frac{\partial G_{jk}}{\partial x^i}-\frac{\partial G_{ij}}{\partial x^k})\\ {(X)}_{}{\mathrm{\Gamma }}_{\alpha \beta \gamma }& :=\frac{1}{2}(\frac{\partial g_{\alpha \gamma }}{\partial X^{\beta }}+\frac{\partial g_{\beta \gamma }}{\partial X^{\alpha }}-\frac{\partial g_{\alpha \beta }}{\partial X^{\gamma }})\\ \end{array}$

Hence

${(x)}_{}{\mathrm{\Gamma }}_{ijk}=\frac{\partial X^{\alpha }}{\partial x^i}\frac{\partial X^{\beta }}{\partial x^j}\frac{\partial X^{\gamma }}{\partial x^k}{}_{(X)}{\mathrm{\Gamma }}_{\alpha \beta \gamma }+\frac{{\partial }^2{X}^{\alpha }}{\partial x^i\partial x^j}\frac{\partial X^{\beta }}{\partial x^k}{g}_{\alpha \beta }$

Define

$[G^{ij}]=[G_{ij}{]}^{-1};[g^{\alpha \beta }]=[g_{\alpha \beta }{]}^{-1}$

Then

$G^{ij}=\frac{\partial x^i}{\partial X^{\alpha }}\frac{\partial x^j}{\partial X^{\beta }}{g}^{\alpha \beta }$

Define the Christoffel symbols of the second kind as

${(x)}_{}{\mathrm{\Gamma }}_{ij}^m:=G^{mk}{}_{(x)}{\mathrm{\Gamma }}_{ijk};{}_{(X)}{\mathrm{\Gamma }}_{\alpha \beta }^{\nu }:=g^{\nu \gamma }{}_{(X)}{\mathrm{\Gamma }}_{\alpha \beta \gamma }$

Then

$\begin{array}{rl}{(x)}_{}{\mathrm{\Gamma }}_{ij}^m& =G^{mk}\frac{\partial X^{\alpha }}{\partial x^i}\frac{\partial X^{\beta }}{\partial x^j}\frac{\partial X^{\gamma }}{\partial x^k}{}_{(X)}{\mathrm{\Gamma }}_{\alpha \beta \gamma }+G^{mk}\frac{{\partial }^2{X}^{\alpha }}{\partial x^i\partial x^j}\frac{\partial X^{\beta }}{\partial x^k}{g}_{\alpha \beta }\\ & =\frac{\partial x^m}{\partial X^{\nu }}\frac{\partial x^k}{\partial X^{\rho }}{g}^{\nu \rho }\frac{\partial X^{\alpha }}{\partial x^i}\frac{\partial X^{\beta }}{\partial x^j}\frac{\partial X^{\gamma }}{\partial x^k}{}_{(X)}{\mathrm{\Gamma }}_{\alpha \beta \gamma }+\frac{\partial x^m}{\partial X^{\nu }}\frac{\partial x^k}{\partial X^{\rho }}{g}^{\nu \rho }\frac{{\partial }^2{X}^{\alpha }}{\partial x^i\partial x^j}\frac{\partial X^{\beta }}{\partial x^k}{g}_{\alpha \beta }\\ & =\frac{\partial x^m}{\partial X^{\nu }}{\delta }_{\rho }^{\gamma }{g}^{\nu \rho }\frac{\partial X^{\alpha }}{\partial x^i}\frac{\partial X^{\beta }}{\partial x^j}{}_{(X)}{\mathrm{\Gamma }}_{\alpha \beta \gamma }+\frac{\partial x^m}{\partial X^{\nu }}{\delta }_{\rho }^{\beta }{g}^{\nu \rho }\frac{{\partial }^2{X}^{\alpha }}{\partial x^i\partial x^j}{g}_{\alpha \beta }\\ & =\frac{\partial x^m}{\partial X^{\nu }}{g}^{\nu \gamma }\frac{\partial X^{\alpha }}{\partial x^i}\frac{\partial X^{\beta }}{\partial x^j}{}_{(X)}{\mathrm{\Gamma }}_{\alpha \beta \gamma }+\frac{\partial x^m}{\partial X^{\nu }}{g}^{\nu \beta }\frac{{\partial }^2{X}^{\alpha }}{\partial x^i\partial x^j}{g}_{\alpha \beta }\\ & =\frac{\partial x^m}{\partial X^{\nu }}\frac{\partial X^{\alpha }}{\partial x^i}\frac{\partial X^{\beta }}{\partial x^j}{}_{(X)}{\mathrm{\Gamma }}_{\alpha \beta }^{\nu }+\frac{\partial x^m}{\partial X^{\nu }}{\delta }_{\alpha }^{\nu }\frac{{\partial }^2{X}^{\alpha }}{\partial x^i\partial x^j}\end{array}$

Therefore

${(x)}_{}{\mathrm{\Gamma }}_{ij}^m=\frac{\partial x^m}{\partial X^{\nu }}\frac{\partial X^{\alpha }}{\partial x^i}\frac{\partial X^{\beta }}{\partial x^j}{}_{(X)}{\mathrm{\Gamma }}_{\alpha \beta }^{\nu }+\frac{\partial x^m}{\partial X^{\alpha }}\frac{{\partial }^2{X}^{\alpha }}{\partial x^i\partial x^j}$

The invertibility of the mapping implies that

$\begin{array}{rl}\frac{\partial X^{\mu }}{\partial x^m}{}_{(x)}{\mathrm{\Gamma }}_{ij}^m& =\frac{\partial X^{\mu }}{\partial x^m}\frac{\partial x^m}{\partial X^{\nu }}\frac{\partial X^{\alpha }}{\partial x^i}\frac{\partial X^{\beta }}{\partial x^j}{}_{(X)}{\mathrm{\Gamma }}_{\alpha \beta }^{\nu }+\frac{\partial X^{\mu }}{\partial x^m}\frac{\partial x^m}{\partial X^{\alpha }}\frac{{\partial }^2{X}^{\alpha }}{\partial x^i\partial x^j}\\ & ={\delta }_{\nu }^{\mu }\frac{\partial X^{\alpha }}{\partial x^i}\frac{\partial X^{\beta }}{\partial x^j}{}_{(X)}{\mathrm{\Gamma }}_{\alpha \beta }^{\nu }+{\delta }_{\alpha }^{\mu }\frac{{\partial }^2{X}^{\alpha }}{\partial x^i\partial x^j}\\ & =\frac{\partial X^{\alpha }}{\partial x^i}\frac{\partial X^{\beta }}{\partial x^j}{}_{(X)}{\mathrm{\Gamma }}_{\alpha \beta }^{\mu }+\frac{{\partial }^2{X}^{\mu }}{\partial x^i\partial x^j}\end{array}$

We can also formulate a similar result in terms of derivatives with respect to $x$. Therefore

$\begin{array}{rl}\frac{{\partial }^2{X}^{\mu }}{\partial x^i\partial x^j}& =\frac{\partial X^{\mu }}{\partial x^m}{}_{(x)}{\mathrm{\Gamma }}_{ij}^m-\frac{\partial X^{\alpha }}{\partial x^i}\frac{\partial X^{\beta }}{\partial x^j}{}_{(X)}{\mathrm{\Gamma }}_{\alpha \beta }^{\mu }\\ \frac{{\partial }^2{x}^m}{\partial X^{\alpha }\partial X^{\beta }}& =\frac{\partial x^m}{\partial X^{\mu }}{}_{(X)}{\mathrm{\Gamma }}_{\alpha \beta }^{\mu }-\frac{\partial x^i}{\partial X^{\alpha }}\frac{\partial x^j}{\partial X^{\beta }}{}_{(x)}{\mathrm{\Gamma }}_{ij}^m\end{array}$

Compatibility conditions

Mining Engineer (Excluding Oil ) Truman from Alma, loves to spend time knotting, largest property developers in singapore developers in singapore and stamp collecting. Recently had a family visit to Urnes Stave Church. The problem of compatibility in continuum mechanics involves the determination of allowable single-valued continuous fields on bodies. These allowable conditions leave the body without unphysical gaps or overlaps after a deformation. Most such conditions apply to simply-connected bodies. Additional conditions are required for the internal boundaries of multiply connected bodies.

Compatibility of the deformation gradient

The necessary and sufficient conditions for the existence of a compatible $\mathit{F}$ field over a simply connected body are

$\mathrm{\nabla }\times \mathit{F}=\mathit{0}$

Compatibility of the right Cauchy-Green deformation tensor

The necessary and sufficient conditions for the existence of a compatible $\mathit{C}$ field over a simply connected body are

$R_{\alpha \beta \rho }^{\gamma }:=\frac{\partial }{\partial X^{\rho }}[{}_{(X)}{\mathrm{\Gamma }}_{\alpha \beta }^{\gamma }]-\frac{\partial }{\partial X^{\beta }}[{}_{(X)}{\mathrm{\Gamma }}_{\alpha \rho }^{\gamma }]+{}_{(X)}{\mathrm{\Gamma }}_{\mu \rho }^{\gamma }{}_{(X)}{\mathrm{\Gamma }}_{\alpha \beta }^{\mu }-{}_{(X)}{\mathrm{\Gamma }}_{\mu \beta }^{\gamma }{}_{(X)}{\mathrm{\Gamma }}_{\alpha \rho }^{\mu }=0$

We can show these are the mixed components of the Riemann-Christoffel curvature tensor. Therefore the necessary conditions for $\mathit{C}$-compatibility are that the Riemann-Christoffel curvature of the deformation is zero.

Compatibility of the left Cauchy-Green deformation tensor

No general sufficiency conditions are known for the left Cauchy-Green deformation tensor in three-dimensions. Compatibility conditions for two-dimensional $\mathit{B}$ fields have been found by Janet Blume.

See also

• Infinitesimal strain
• Compatibility (mechanics)
• Curvilinear coordinates
• Piola-Kirchhoff stress tensor, the stress tensor for finite deformations.
• Stress measures

References

43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.

Further reading

• 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
  
  My blog: http://www.primaboinca.com/view_profile.php?userid=5889534

• 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
  
  My blog: http://www.primaboinca.com/view_profile.php?userid=5889534

• 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
  
  My blog: http://www.primaboinca.com/view_profile.php?userid=5889534

• 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
  
  My blog: http://www.primaboinca.com/view_profile.php?userid=5889534

• 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
  
  My blog: http://www.primaboinca.com/view_profile.php?userid=5889534

• 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
  
  My blog: http://www.primaboinca.com/view_profile.php?userid=5889534

• 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
  
  My blog: http://www.primaboinca.com/view_profile.php?userid=5889534

• 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
  
  My blog: http://www.primaboinca.com/view_profile.php?userid=5889534

• 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
  
  My blog: http://www.primaboinca.com/view_profile.php?userid=5889534

External links

• Prof. Amit Acharya's notes on compatibility on iMechanica