For 3D beam elements, when considering the geometry nonlinearity, the final equation becomes (Ke+Kg)*deltaU=deltaF, Kg the geometric stiffness matrix is a function of total axial force acting on the element in the reference configuration. The character of symmetric positive-definite for this matrix always changes with the element force?