We consider a stress field such that \(\Sigma =\Sigma _1+\Sigma _2\)  where \(\Sigma _1\in \mathcal {B}(\Omega ,\operatorname{Sym})\) , with bulk density \(\sigma \) , and \(\Sigma _2 \in \mathcal {C}(\Omega ,\operatorname{Sym})\) , with line density \(\tau \) ; a moment field such that \(M=M_1+M_2\)  where \(M_1\in \mathcal {B}(\Omega ,\operatorname{Sym})\) , with bulk density \(m\) , and \(M_2 \in \mathcal {C}(\Omega ,\operatorname{Sym})\) , with line density \(n\) ; a bending strain field such that \(\Lambda =\Lambda _1+\Lambda _2\)  where \(\Lambda _1\in \mathcal {B}(\Omega ,\operatorname{Sym})\) , with bulk density \(\lambda \) , and \(\Lambda _2 \in \mathcal {C}(\Omega ,\operatorname{Sym})\) , with line density \(\gamma \) ; and a transverse force field such that \(F=F_1+F_2\)  where \(F_1\in \mathcal {B}(\Omega )\) , with bulk density \(f_1\) , and \(F_2 \in \mathcal {C}(\Omega )\) , with line density \(f_2\) . We can then use the divergence identities from [1]} to derive the local form of (REF ):
\(\operatorname{div}\sigma = 0 ~\text{in}~\Omega -\lbrace S\cup \lbrace O\rbrace  \rbrace , \\\frac{d\tau }{ds} t - \llbracket \sigma  \rrbracket \nu =0 ~\text{on}~S-\lbrace O\rbrace , \\\tau  \nu =0 ~\text{on}~S-\lbrace O\rbrace , ~\text{and} \\\Sigma (\nabla (v^\alpha e_1))=0, ~~\Sigma (\nabla (v^\alpha e_2))=0~\text{at}~O,\) 
