Load case sequence
The response of nonlinear materials depends not only on the magnitude and direction of the applied forces but also on the order in which the loading effects are applied. It is therefore necessary to specify the chronological sequence of loading. We start from the usual concept of a load case, by which we mean the temperature and all forces acting on the body at a given moment in time. All load cases to be used in the nonlinear computation are defined in the usual way using the RPD2/RPD3 program. From these pre-prepared loads, the HPP2/HPP3 program then selects and assembles the sequence needed to model the loading history, as shown schematically in the following figure.
It is automatically assumed that the transition from $L_i$ to $L_{i+1}$ is uniform in terms of strain. This transition between individual states is very close to a uniform change in loading, but it need not always be identical. Therefore, in cases where it is essential to strictly follow the force loading path, the relevant segment must be divided into smaller increments. Under all circumstances, the system passes through any load case that has been prescribed.
Load cases are defined as follows:
1)
In the name.i2 file, within batch AS 1 (i.e., in the first load case), the material, zero nodal displacements, any springs, the body’s initial temperature, and pairs of contact surfaces are defined in the usual way. The loads listed below (pressure, concentrated forces, etc.) are taken into account only in the linear problem. For nonlinear problems, the loading effects described in AS 1 cannot be used, and they must be specified in AS 2 and above. In the assignment /R 0 $T_o$ $T_w$/, only the initial temperature $T_o$, corresponding to the stress-free state, is significant. The value of $T_w$ has no effect.
2)
In batches AS 2, AS 3, …, all the load cases to be used in the nonlinear computation are created one after another. Their order does not matter. The loading is not specified incrementally but absolutely (with respect to zero). For example, if the assignment /R 0 $T_o$ $T_w$/ occurs, a load case arises in which the body has the temperature $T_w$. The value of $T_o$ is ignored.
3)
The input data are processed by the RPD2/RPD3 program, and the stiffness matrices are then calculated by the SRH2/SRH3 program, followed by the factorization of the system of linear equations by the FEFS program. It is also possible to calculate the elastic stresses using the STR2/STR3 program.
4)
A sequence of load cases $L_1$, $L_2$, … is assembled, where the numbers $L_i$ are the sequence numbers of the AS batches. For example, $L_2=5$ means that the second load case was defined in AS 5 (the current increment relative to the previous configuration is $L_2-L_1$). If it later turns out that a load case is missing, the entire problem must be recomputed starting from point 2.
5)
In the name.iP file, on the IP line, the following is specified:
IP KREST NLC NCYC KMOD 0 KTPR KCNT 3*0 $L_1$ $L_2$ $\dots$ $L_\mathtt{NLC}$
where
- NLC is the number of terms in the load case sequence
- NCYC is the number of cycles (repetitions of the entire sequence, the default value is $1$)
The total number of solved states is $\mathtt{NLC}\cdot\mathtt{NCYC}$.
6)
In the name.iP file, on the RP line, the times in seconds corresponding to all load cases are specified:
RP 10*0 $t_1$ $t_2$ $\dots$ $t_\mathtt{NLC}$
It is assumed that $t_\mathtt{NLC}>\ldots>t_2>t_1>0$ and that $\mathtt{NCYC}$ is always $1$ (i.e., it is also possible to write $\mathtt{NCYC}=0$).
7)
The input data are processed by the HPP2/HPP3 program, and the problem is solved by the HDYN program. If the solution was successful, we can continue at point 5 by specifying additional load cases, which, however, had to be defined in advance within the AS batches in the name.i2 file — returning to point 2 means recomputing the entire problem. In the name.iP file, it is sufficient to write $\mathtt{KREST}=2$ at the first position of the IP line:
IP 2 NLC NCYC KMOD $\dots$
For dynamic problems, it is naturally assumed that the first time $t_1$ is greater than or equal to the time at which the previous solution ended. Otherwise, the HPP2/HPP3 program reports an error. A restart can be used for solving problems with pre-stress. First, the pre-stress is calculated with $\mathtt{KREST}=1$ (e.g. residual stresses after a cooling process), and then a new series of loading (e.g. cyclic) is started with $\mathtt{KREST}=2$. In this way, it is possible to work with material properties that have been altered by their previous history.