General variational many-body theory with complete self-consistency for trapped bosonic systems