# Biotissue Macroscale Nonlinear Analysis

(Difference between revisions)
 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, $\sigma$, 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, $\tilde{C}$, which replaces $C B_L$ in the calculation of the linear strain portion of the stiffness matrix, $K_L$. The microscale RVE solver returns the six unique components of the stress tensor, $\sigma$, 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, $\tilde{C}$, 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. ==Initial linear elastic solve== ==Initial linear elastic solve==

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

Back to Biotissue Project

## 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):

(KL + KNLU = PF

where KL is the linear strain portion of the stiffness matrix, KNL 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:

$K_L = \int_V B_L^T C B_L dV$

$K_{NL} = \int_V B_{NL}^T \tilde{\sigma} B_{NL} dV$

$F = \int_V B_L^T \hat{\sigma} dV$

where BL is the linear strain-displacement matrix, C is the stress-strain matrix, BNL is the nonlinear strain-displacement matrix, $\tilde{\sigma}$ is a matrix of Cauchy stresses, and $\hat{\sigma}$ is a vector of Cauchy stresses. The matrices are given by:

$B_L = \begin{bmatrix} N_{1,1} & 0 & 0 & N_{2,1} & \cdots & 0 \\ 0 & N_{1,2} & 0 & 0 & \cdots & 0 \\ 0 & 0 & N_{1,3} & 0 & \cdots & N_{4,3} \\ N_{1,2} & N_{1,1} & 0 & N_{2,2} & \cdots & 0 \\ 0 & N_{1,3} & N_{1,2} & 0 & \cdots & N_{4,2} \\ N_{1,3} & 0 & N_{1,1} & N_{2,3} & \cdots & N_{4,1} \\ \end{bmatrix}$

$B_{NL} = \begin{bmatrix} \tilde{B}_{NL} & 0 & 0 \\ 0 & \tilde{B}_{NL} & 0 \\ 0 & 0 & \tilde{B}_{NL} \\ \end{bmatrix}$

$\tilde{B}_{NL} = \begin{bmatrix} N_{1,1} & 0 & 0 & N_{2,1} & \cdots & N_{4,1} \\ N_{1,2} & 0 & 0 & N_{2,2} & \cdots & N_{4,2} \\ N_{1,3} & 0 & 0 & N_{2,3} & \cdots & N_{4,3} \\ \end{bmatrix}$

$\tilde{\sigma} = \begin{bmatrix} \sigma & 0 & 0 \\ 0 & \sigma & 0 \\ 0 & 0 & \sigma \\ \end{bmatrix}$

$\sigma = \begin{bmatrix} \sigma_{11} & \sigma_{12} & \sigma_{13} \\ \sigma_{21} & \sigma_{22} & \sigma_{23} \\ \sigma_{31} & \sigma_{32} & \sigma_{33} \\ \end{bmatrix}$

$\hat{\sigma} = \begin{bmatrix} \sigma_{11} & \sigma_{22} & \sigma_{33} & \sigma_{12} & \sigma_{23} & \sigma_{31} \\ \end{bmatrix}^T$

where we've assumed 4 nodes and Ni,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, KNL, and the nodal point forces, F. It also returns nodal stress derivatives, $\tilde{C}$, which replaces CBL in the calculation of the linear strain portion of the stiffness matrix, KL.

## 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.