In the analysis of finite deformation problems, the principle equation of virtual work becomes a nonlinear equation regarding the displacement-strain relation. To solve the nonlinear equation, an iterative algorithm is generally applied. When implementing an incremental analysis for a finite deformation problem, whether to refer to the initial status as a reference layout, or refer to the starting point of the increments can be selected. The former is called the total Lagrange method, and the latter is called the updated Lagrange method. Both the total Lagrange and updated Lagrange methods are available in the program. This section discusses the various geometrically nonlinear options available, including the large strain.
Decomposition of increments of virtual work equation¶
Given the solid deformation at time t is known, the status at time t'=t+\triangle t is unknown. The equilibrium equation, dynamic boundary condition, and external boundary condition can be expressed as
where \(^{t'}\sigma\), \(^{t'}\mathbf{b}\), \(^{t'}\mathbf{n}\), \(^{t'}\mathbf{t}\), \(^{t'}\mathbf{u}\) are the Cauchy stress, body force, outward normal vector of the object's surface, fixed surface force, and fixed displacement in each time t'.
The equation \(\eqref{eq:ch5_nonlinear_gov2}\) needs to be solved referring to layout V at time 0, or layout \(^{t}v\) at time t. The following sections will introduce these two algorithms: total Lagrange method and updated Lagrange method, respectively.
where \(_{0}^{t'}\mathbf{S}\) and \(_{0}^{t'}\mathbf{E}\) are the 2nd order Piola-Kirchhoff stress tensor and the Green-Lagrange strain tensor at time t', respectively. The initial domain at time 0 is called the reference domain. The body force \(_{0}^{t'}\mathbf{b}\) and nominal surface force vector \(_{0}^{t'}\mathbf{t}\) are
Substituting equations \(\eqref{eq:ch5_nonlinear_total_lag_gov3}\), \(\eqref{eq:ch5_nonlinear_total_lag_gov4}\), \(\eqref{eq:ch5_nonlinear_total_lag_gov5}\), and \(\eqref{eq:ch5_nonlinear_total_lag_gov6}\) into governing equations \(\eqref{eq:ch5_nonlinear_total_lag_gov1}\) and \(\eqref{eq:ch5_nonlinear_total_lag_gov2}\), we have
where \(_{0}^{t}\mathbf{K}_{L}\), \(_{0}^{t}\mathbf{K}_{NL}\), \(_{0}^{t'}\mathbf{F}\), \(_{0}^{t}\mathbf{Q}\) are the initial material stiffness matrix, initial geometric stiffness (stress) matrix, external force vector, and internal force vector, respectively. The recursive algorithm to calculate the deformation status at time t' from time t is given:
In addition to the total Lagrange algorithm, the updated Lagrange algorithm is also widely applied in the nonlinear structural model computation. The principle virtual work equation at time t' uses the current domain at time t as reference domain.
The tensors \(_{t}^{t'}\mathbf{S}\), \(_{t}^{t'}\mathbf{E}\) and vectors \(_{t}^{t'}\mathbf{t}\), \(_{t}^{t'}\mathbf{b}\) are using the current time domain t as the reference domain. Therefore, the Green-Lagrange strain does not contain the initial displacement (the displacement at the time t) \(^{t}\mathbf{u}\);
Substituting equations \(\eqref{eq:ch5_nonlinear_updated_lag_gov3}\) and \(\eqref{eq:ch5_nonlinear_updated_lag_gov2}\) into governing equations \(\eqref{eq:ch5_nonlinear_updated_lag_gov1}\), we have
where \(_{t}^{t}\mathbf{K}_{L}\), \(_{t}^{t}\mathbf{K}_{NL}\), \(_{t}^{t'}\mathbf{F}\), \(_{t}^{t}\mathbf{Q}\) are the initial material stiffness matrix, initial geometric stiffness (stress) matrix, external force vector, and internal force vector, respectively. The recursive algorithm to calculate the deformation status at time t' from time t is given: