A non-linear symmetric G1-conforming Bézier finite element formulation for the analysis of Kirchhoff beam assemblies