Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F90394155
manual-cohesive_elements.tex
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Fri, Nov 1, 07:00
Size
4 KB
Mime Type
text/x-tex
Expires
Sun, Nov 3, 07:00 (2 d)
Engine
blob
Format
Raw Data
Handle
18289143
Attached To
rAKA akantu
manual-cohesive_elements.tex
View Options
\section{Cohesive Elements}
The cohesive elements that have been implemented in \akantu are based
on the work of Ortiz and Pandolfi~\cite{ortiz1999}. Their main
properties are reported in Table~\ref{tab:coh:cohesive_elements}.
\begin{table}[!htb]
\begin{center}
\begin{tabular}{l|llcc}
\toprule
Element type & Facet type & Order & \# nodes & \# quad. points \\
\midrule
\texttt{\_cohesive\_1d\_2} & \texttt{\_point\_1} & linear & 2 & 1 \\
\hline
\texttt{\_cohesive\_2d\_4} & \texttt{\_segment\_2} & linear & 4 & 1 \\
\texttt{\_cohesive\_2d\_6} & \texttt{\_segment\_3} & quadratic & 6 & 2 \\
\hline
\texttt{\_cohesive\_3d\_6} & \texttt{\_triangle\_3} & linear & 6 & 1 \\
\texttt{\_cohesive\_3d\_12} & \texttt{\_triangle\_6} & quadratic & 12 & 3 \\
\bottomrule
\end{tabular}
\end{center}
\caption{Some basic properties of the cohesive elements in \akantu.}
\label{tab:coh:cohesive_elements}
\end{table}
\begin{figure}
\centering
\includegraphics[width=.6\textwidth]{figures/cohesive2d}
\caption{Cohesive element in 2D for quadratic triangular elements
T6.}
\label{fig:smm:coh:cohesive2d}
\end{figure}
Cohesive element insertion can be either realized at the beginning of
the simulation or it can be carried out dynamically during the
simulation. The first approach is called \emph{intrinsic}, the second
one \emph{extrinsic}. When an element is present from the beginning, a
bilinear or exponential cohesive law should be used instead of a
linear one. A bilinear law works exactly like a linear one except for
an additional parameter $\delta_0$ separating an initial linear
elastic part from the linear irreversible one. For additional details
concerning cohesive laws see Section~\ref{sec:cohesive-laws}.
\begin{figure}
\centering
\includegraphics[width=.75\textwidth]{figures/insertion}
\caption{Insertion of a cohesive element.}
\label{fig:smm:coh:insertion}
\end{figure}
Extrinsic cohesive elements are dynamically inserted between two
standard elements when
\begin{equation}
\sigma_\mathrm{eff} > \sigma_\mathrm{c} \quad\text{with}\quad
\sigma_\mathrm{eff} = \sqrt{\sigma_\mathrm{n}^2 +
\frac{\tau^2}{\beta^2}}
\end{equation}
in which $\sigma_\mathrm{n}$ is the tensile normal traction and $\tau$
the resulting tangential one (Figure~\ref{fig:smm:coh:insertion}).
For the static analysis of the structures containing cohesive
elements, the stiffness of the cohesive elements should also be added
to the total stiffness of the structure. Considering a 2D quadratic
cohesive element as that in Figure~\ref{fig:smm:coh:cohesive2d}, the
opening displacement along the mid-surface can be written as:
\begin{equation}
\label{eq:opening}
\vec{\Delta}(s) = \llbracket \mat{u}\rrbracket \,\mat{N}(s) =
\begin{bmatrix}
u_3-u_0 & u_4-u_1 & u_5-u_2\\
v_3-v_0 & v_4-v_1 & v_5-v_2
\end{bmatrix}
\begin{bmatrix}
N_0(s) \\ N_1(s) \\ N_2(s)
\end{bmatrix} =
\mat{N}^\mathrm{k} \mat{A U} = \mat{PU}
\end{equation}
The \mat{U}, \mat{A} and $\mat{N}^\mathrm{k}$ are as following:
\begin{align}
\mat{U} &= \left [
\begin{array}{c c c c c c c c c c c c}
u_0 & v_0 & u_1 & v_1 & u_2 & v_2 & u_3 & v_3 & u_4 & v_4 & u_5 & v_5
\end{array}\right ]\\[1ex]
\mat{A} &= \left [\begin{array}{c c c c c c c c c c c c}
1 & 0 & 0 & 0& 0 & 0 & -1& 0 & 0 &0 &0 &0\\
0 &1& 0&0 &0 &0 &0 & -1& 0& 0 & 0 &0\\
0 &0& 1&0 &0 &0 &0 & 0& -1& 0 & 0 &0\\
0 &0& 0&1 &0 &0 &0 & 0& 0& -1 & 0 &0\\
0 &0& 0&0 &1 &0 &0 & 0& 0& 0 & -1 &0\\
0 &0& 0&0 &0 &1 &0 & 0& 0& 0 & 0 &-1
\end{array} \right ]\\[1ex]
\mat{N}^\mathrm{k} &= \begin{bmatrix}
N_0(s) & 0 & N_1(s) &0 & N_2(s) & 0\\
0 & N_0(s)& 0 &N_1(s)& 0 & N_2(s)
\end{bmatrix}
\end{align}
The consistent stiffness matrix for the element is obtained as
\begin{equation}
\label{eq:cohesive_stiffness}
\mat{K} = \int_{S_0} \mat{P}^\mathrm{T}\,
\frac{\partial{\vec{T}}} {\partial{\delta}} \mat{P} \,\mathrm{d}
S_0
\end{equation}
where $\vec{T}$ is the cohesive traction and $\delta$ the opening
displacement (for more details check
Section~\ref{tab:coh:cohesive_elements}).
%%% Local Variables:
%%% mode: latex
%%% TeX-master: "manual"
%%% End:
Event Timeline
Log In to Comment