# Biotissue Macroscale Nonlinear Analysis

### From ScorecWiki

Revision as of 13:27, 6 March 2014Fovard (Talk | contribs) ← Previous diff |
Revision as of 19:04, 18 March 2014Fovard (Talk | contribs) Next diff → |
||

Line 68: |
Line 68: | ||

The microscale RVE solver returns the six unique components of the stress tensor, <math>\sigma</math>, which contributes to the nonlinear strain portion of the stiffness matrix, <math>K_{NL}</math>, and the nodal point forces, <math>F</math>. It also returns nodal stress derivatives, <math>\tilde{C}</math>, which replaces <math>C B_L</math> in the calculation of the linear strain portion of the stiffness matrix, <math>K_L</math>. | The microscale RVE solver returns the six unique components of the stress tensor, <math>\sigma</math>, which contributes to the nonlinear strain portion of the stiffness matrix, <math>K_{NL}</math>, and the nodal point forces, <math>F</math>. It also returns nodal stress derivatives, <math>\tilde{C}</math>, which replaces <math>C B_L</math> in the calculation of the linear strain portion of the stiffness matrix, <math>K_L</math>. | ||

+ | ==Fiber-Matrix Parallel== | ||

+ | |||

+ | The ''fiber-matrix parallel'' model utilizes the fiber only RVEs and simply adds the stress found from the fiber network with stress computed from a neo-Hookean material law that approximates the response of the interfibrillar matrix. | ||

==Initial linear elastic solve== | ==Initial linear elastic solve== |

## Revision as of 19:04, 18 March 2014

Back to Biotissue Project

## Contents |

## Overview

Generally, we follow the finite element nonlinear analysis of solids presented in Bathe's Finite Element Procedures. We use the static updated Lagrangian formulation (section 6.3.2):

(*K*_{L} + *K*_{NL})Δ*U* = *P* − *F*

where *K*_{L} is the linear strain portion of the stiffness matrix, *K*_{NL} is the nonlinear strain portion of the stiffness matrix (for geometric nonlinearity), Δ*U* is the incremental displacement vector, *P* is external virtual work, and *F* is nodal point forces. The portions of these matrices and vectors corresponding to certain nodes are computed at integration points in the source file macro/src/MultiscaleIntegrator.h in the atPoint function.

The system above is solved for the incremental displacements, Δ*U*. The terms are given by:

where *B*_{L} is the linear strain-displacement matrix, *C* is the stress-strain matrix, *B*_{NL} is the nonlinear strain-displacement matrix, is a matrix of Cauchy stresses, and is a vector of Cauchy stresses. The matrices are given by:

where we've assumed 4 nodes and *N*_{i,j} denotes a shape function derivative at node *i* wrt spatial coordinate *j*.

## Incorporating Microscale Data

The microscale RVE solver returns the six unique components of the stress tensor, σ, which contributes to the nonlinear strain portion of the stiffness matrix, *K*_{NL}, and the nodal point forces, *F*. It also returns nodal stress derivatives, , which replaces *C**B*_{L} in the calculation of the linear strain portion of the stiffness matrix, *K*_{L}.

## Fiber-Matrix Parallel

The *fiber-matrix parallel* model utilizes the fiber only RVEs and simply adds the stress found from the fiber network with stress computed from a neo-Hookean material law that approximates the response of the interfibrillar matrix.