摘要 |
For each (ij)-th component of a tensor, an equation that is denoted by DeltaF1 (ij) and that uses epsilon1 is computed on the basis of a function W(F) of an inputted tensor amount F and a value (F=F^) of the tensor amount F (312). For each (kl)-th component of a tensor, an equation that is denoted by ~DeltaF2 (k1) and that uses epsilon1 and epsilon2 is computed on the basis of the value (F=F^) of the tensor amount F (314). For each combination of an (ij)-th component and a (kl)-th component of a tensor, a function W(F^+DeltaF1 (ij)+~DeltaF2 (k1)) is computed by using the computed equation that is denoted by DeltaF1 (ij) and the computed equation that is denoted by ~DeltaF2 (k1) (316). For each (ij)-th component of a tensor, a coefficient of epsilon1 in the computed function W(F^+DeltaF1 (ij)+~DeltaF2 (k1)) is taken-out, and stress, that is based on a first order derivative with respect to the tensor amount F of the function W(F), is computed. For each combination of an (ij)-th component and a (kl)-th component of a tensor, a coefficient of epsilon1sepsilon2 in the function W(F^+DeltaF1 (ij)+~DeltaF2 (k1)) is taken-out, and a material Jacobian, that is based on a second order derivative with respect to the tensor amount F of the function W(F), is computed (318, 320). |