The effect of the absolute nodal coordinate formulation (ANCF)–coupled deformation modes on the accuracy and efficiency when higher order three-dimensional beam and plate finite elements are used is investigated in this study. It is shown that while computational efficiency can be achieved in some applications by neglecting the effect of some of the ANCF-coupled deformation modes, such modes introduce geometric stiffening/softening effects that can be significant in the case of very flexible structures. As shown in previous publications, for stiff structures, the effect of the ANCF-coupled deformation modes can be neglected. For such stiff structures, the solution does not strongly depend on some of the ANCF-coupled deformation modes, and formulations that include these modes lead to numerical results that are in good agreement with formulations that exclude them. In the case of a very flexible structure, on the other hand, the inclusion of the ANCF-coupled deformation modes becomes necessary in order to obtain an accurate solution. In this case of very flexible structures, the use of the general continuum mechanics approach leads to an efficient solution algorithm and to more accurate numerical results. In order to examine the effect of the elastic force formulation on the efficiency and the coupling between different modes of deformation, three different models are used again to formulate the elastic forces in the absolute nodal coordinate formulation. These three methods are the general continuum mechanics approach, the elastic line (midsurface) approach, and the elastic line (midsurface) approach with the Hellinger–Reissner principle. Three-dimensional absolute nodal coordinate formulation beam and plate elements are used in this study. In the general continuum mechanics approach, the coupling between the cross section deformation and the beam centerline or plate midsurface displacement is considered, while in the approaches based on the elastic line and the Hellinger–Reissner principle, this coupling is neglected. In addition to the fully parametrized beam element used in this study, three different plate elements, two fully parametrized and one reduced order thin plate elements, are used. The numerical results obtained using different finite elements and elastic force formulations are compared in this study.
Computer Simulations of DNA Packing inside Bacteriophages: Elasticity, Electrostatics and Entropy
