Week - 2 - Explicit and Implicit Analysis

EXPLICIT AND IMPLICIT ANALYSIS

INTRODUCTION

Numerical analysis employs both explicit and implicit techniques to arrive at numerical approximations for the solutions of time-dependent ordinary and partial differential equations that call for computer simulations of physical processes. While implicit techniques arrive at a solution by solving an equation involving both the current state of the system and the later one, explicit methods determine the state of a system at a later time from the state of the system at the current time. Incremental loads are required for all dynamic and nonlinear analyses. To tackle these issues, "implicit" and "explicit" strategies are frequently employed. Mass (inertia) and damping have little impact on static analysis. Nodal forces related to mass/inertia and damping are taken into account in dynamic analysis. Both the explicit solver and the implicit solver can be used for dynamic analysis. Each step's solution in a nonlinear implicit analysis requires iterations in order to achieve equilibrium within a predetermined tolerance. Since the nodal accelerations are immediately solved in explicit analysis, iteration is not necessary.

Difference between Implicit and Explicit

Implicit and Explicit analysis dissociate in the approach to time incrementation. In Implicit analysis, each time increment has to converge, but users can set quite long time increments. Conversely, explicit analysis doesn’t have to converge each increment, but time increments must be super small for the accuracy of solution.

Implicit transient analysis has no time step limit. Implicit time steps are most often larger than explicit time steps.

Formulations:

The matrix equation is reduced by the partial differential equations in linear problems to:

[K]{x} = {f}

Where

k = stiffness matrix

x = displacement/deflection

F = force

For static non-linear situations like:

[K(x)]

{x} = {f} → [k0 + k1• x + k2 • x2 +…] {x} = {f}

The matrix equations change in dynamic issues to:

[M]{x´´} + [C]{x´} + [K]{x} = {f}

x´ = velocity

x´´ = acceleration

C = damping matrix

M = mass matrix

When to use Implicit?

When the structure is being affected by the boundary conditions slowly and the effects of strain rates are minimal, implicit methods should be employed. After it is possible to determine the increase in stress as a function of strain, geometry can be examined utilising implicit techniques. Each time increment in the model establishes the global equilibrium. This implies that every increment must converge. The solver determines all local finite element variables (stresses, strain, etc.) for each increment after the global equilibrium has converged.

When to use Explicit?

When the strain rates are equal to or greater than 10^-3 per second, explicit approaches should be utilised. Examples of these occurrences include vehicular accidents, ballistic explosions, drops, and more. In these situations, in addition to the variation in stress with strain, the rate of strain must also be taken into account. There are no iterative processes or convergence criteria to evaluate. The local finite element variables are the solver's primary emphasis. For a particular increment, the solver calculates all of the local finite element variables before moving on to the next one.

AIM

As demonstrated in the PDF, F(u) = u3+9u2+4u use this equation and solve using both Explicit and Implicit Methods ( have a tolerance of 0.001).

1. Write a detailed report about Implicit and explicit methods.
2. Do the calculations in a spreadsheet, try to plot the graph and compare them (optional).

EXPLANATION

Assume the instance of a simple bar in tension. Consider a scenario where the force in the bar is a nonlinear function of the displacement.

f(u)=u+u^3/3 (1)

As a result of (1), the stiffness of the bar is

k(u)=df/du=1+u^2 (2)

u

F

Figure 1 A simple bar in tension

Explicit Analysis

Using equation (1), it is evidently possible to plot a graph of force, f, versus displacement, u. However, if the stiffness is known but not the nonlinear force against displacement relationship, it is still possible to numerically create the force versus displacement graph. Utilizing an incremental explicit load control strategy is one approach. Next, this is demonstrated. The tension bar is loaded in this example using three equal load increments (or steps) of ∆F = 1 unit each. The variables are: k for stiffness, u for displacement, f for internal force in the bar, and f for external force acting on the bar. ∆u or ∆F are used to represent incremental displacements or incrementally applied external forces, respectively. The equation ∆F = k∆u is put to use. The analysis goes on like this:

Step 1

u₀ = 0.0, ⇒ k (u₀) = 1, ∆F = 1, ∆u₁ = ∆F/k (u₀) = 1.0, u₁ = u₀ + ∆u₁ = 1.0

Step 2

⇒ k (u₁) = 1 + 1² = 2, ∆F = 1, ∆u₂ = ∆F/k(u₁) = 1/2, u₂ = u₁ + ∆u₂ = 1.5

Step 3

⇒ k(u₂) = 1+1.5² = 3.25, ∆F = 1, ∆u₃ = ∆F/k(u₂) = 1/3.25, u₃ = u₂+∆u₃ ∼= 1.81

Finally, the three stages provide the results listed below.

F_ext = ∆F + ∆F + ∆F = 3.0

`f_int =u^3/3^3 =1.81+〖1.81〗^3/3 ≅ 3.79`

Because F_ext = f_int, it is clear from the results that exterior and internal forces are not in equilibrium. The key findings for the explicit analysis are presented in Table 1.

Table 1: Results of an explicit analytical summary

Step i	∆Fi	∆ui	ui	(Fext)i	(fint)i	fint − Fext = R
1	1.0	1.000	1.00	1.0	1.333	0.333
2	1.0	0.500	1.50	2.0	2.625	0.625
3	1.0	0.308	1.81	3.0	3.790	0.790

Implicit Analysis

The only difference between an implicit analysis and an explicit analysis is that before moving on to the following phase, Newton-Raphson iterations are employed to impose equilibrium at the conclusion of each stage. In essence, a small force is used to move the solution forward at the start of a step. However, unless the stiffness is linear for the specified step, internal and exterior forces won't be in equilibrium. Therefore, modifications must be made to the displacement in order to reach equilibrium. This is done by minimizing the residual using Newton-Raphson iterations, where R (u) = f_int − F_ext. The residual is expanded as a Taylor series about the present displacement u^j, which results in

`R^(j+1)=R^j (u^j )+dR(u^j )/du δu^(j+1 )+⋯= R^(j )+k(u^(j ) )δu^(j+1)+⋯`

The following adjustment is derived by ignoring higher order terms in the series, setting R^j+1 = 0 and solving for δu^j+1

^{δuj+1 = −[k(uj)]−1Rj.} (4)

The number of iterations required to lower the norm of the residual, , below the selected tolerance threshold (in the example below, a tolerance of 10²), is indicated by the iteration variable, j.

The process of the implicit analysis is as follows (results are approximations; they are more accurate when performed on a computer to machine precision):

Step 1

u₀ = 0.0, ⇒ k(u ) = 1, ∆F = 1, ∆u = ∆F/k(u ) = 1.0, u = u + ∆u = 1.0

Check the residual, F_ext = ∆F = 1, f_int = u + u³/3 = 4/3, R⁽⁰⁾ = f_int − F_ext = 1/3

R⁽⁰⁾ > 10 ² ⇒ Newton-Raphson iterations are necessary.

Calculate the correction to u₁ = u₁⁽⁰⁾

δu⁽¹⁾ = −[k(u₁⁽⁰⁾ )] ¹R⁽⁰⁾ = −2 ¹(1/3) = −0.16667.

⇒ The updated value of u₁⁽¹⁾ = u₁⁽⁰⁾ + δu⁽¹⁾ = 1 + (−0.16667) = 0.83333.

Check the residual again, F_ext = 1, f_int = u₁⁽¹⁾ + (u₁⁽¹⁾)³/3 = 1.02623

R⁽¹⁾ = f_int − F_ext = 0.02623

R⁽¹⁾ > 10² ⇒ another Newton-Raphson iteration is necessary.

Calculate the new total correction to u₁ = u₁⁽⁰⁾

δu⁽²⁾ = δu⁽¹⁾ − [k(u₁⁽¹⁾ )] ¹R⁽¹⁾ = −0.16667 − (1.6944 ¹)0.02623 = −0.1821.

The updated value of u₁⁽²⁾ = u₁⁽⁰⁾ + δu⁽²⁾ = 1 + (−0.18215) = 0.81785.

Check the residual, F_ext = 1, f_int = u₁⁽²⁾ + (u₁⁽²⁾)³/3 = 1.000197

R⁽²⁾ = f_int − F_ext = 0.000197

R⁽²⁾ < 10 ² ⇒ no further iterations are necessary.

Therefore, the final value is u₁ = u₁⁽²⁾ = 0.81785.

Step 2

k(u¹) = 1 + 0.81785² = 1.66888, ∆F = 1, ∆u² = ∆F/k(u¹) = 0.59920

u₂ = u₁ + ∆u₂ = 1.41705

Check the residual, F_ext = ∆F = 1, f_int = u₂ + u₂³/3 = 2.36554

R⁽⁰⁾ = f_int − F_ext = 0.36554

R⁽⁰⁾ > 10² ⇒ Newton-Raphson iterations are necessary.

Calculate the correction to u₂ = u₂⁽⁰⁾

δu⁽¹⁾ = −[k(u₂⁽⁰⁾)] ¹R⁽⁰⁾ = −3.42508¹(0.36554) = −0.10672.

⇒ The updated value of u₂⁽¹⁾ = u₂⁽⁰⁾ + δu(1) = 141705 + (−0.10672) = 1.31033.

Check the residual again, F_ext = 1, f_int = u₂⁽¹⁾ + (u₂⁽¹⁾)³/3 = 2.060260

R⁽¹⁾ = f_int − F_ext = 0.060260

R⁽¹⁾ > 10² ⇒ another Newton-Raphson iteration is necessary.

Calculate the new total correction to u₂ = u₂⁽⁰⁾

δu⁽²⁾ = δu⁽¹⁾ − [k(u₂⁽¹⁾)]¹R⁽¹⁾ = −0.10672 − (3.02729¹)0.06026 = −0.12663.

The updated value of u₂⁽²⁾ = u₂⁽⁰⁾ + δu⁽²⁾ = 1.41705 + (−0.12663) = 1.29042.

Check the residual, F_ext = 1, f_int = u₂⁽²⁾ + (u₂⁽²⁾)³/3 = 2.00663

R⁽²⁾ = f_int − F_ext = 0.00663

R⁽²⁾< 10² ⇒ no further iterations are necessary.

Therefore, the final value is u₂ = u₂⁽²⁾ = 1.29042.

Step 3 is completed similar to steps 1 and 2, the final value is u₃ = 1.6097.

Table 2 lists the outcomes of the implicit analysis.

The table clearly shows that the internal and external forces are achieving equilibrium.

Table 2: Results of an implicit analysis summary

Step i	∆Fi	ui	(Fext)i	(fint)i	fint − Fext = R
1	1.0	0.81785	1.0	1.000197	0.000197
2	1.0	1.29042	2.0	2.00663	0.00663
3	1.0	1.61	3.0	3.001	0.001

Plot of results

Figure 2 presents a visualization of the findings. The graphic clearly shows that the explicit analysis deviates from the precise solution. The explicit analysis more nearly matches the precise solution when 20 increments are employed. However, the explicit analysis still veers off the precise solution even with more increments. The deviation from the exact solution demonstrates the absence of internal and external force equilibrium. An implicit analysis is necessary to address this issue. In Figure 2, the outcomes of the implicit analysis. The incremental stages are corrected using the Newton-Raphson iterations, leading to the precise solution within the predetermined tolerance. Only three increments are needed to reach excellent agreement, and only two Newton-Raphson iterations are needed for each step. The non e.m and non i.m scripts in MATLAB are used to complete the explicit and implicit analyses.

Figure 2 Comparison of explicit and implicit results of load versus displacement

CONCLUSION

If the circumstances in your analysis occur very slowly, the implicit solver is really helpful. Ideal analyses typically take more time than one second. The ability for users to set any size time increment is an advantage.

The stiffness matrix's inversion yields the unknown x. (K).

Enforced equilibrium solutions for Newton-Raphson.

1. The Newton technique is applied.
2. More precise
3. Large models take a very long time.
4. Additional computer storage is necessary
5. Less than 10-3 strains per second are experienced.
6. A wise choice for dynamic structure types with sluggish load rates; Displacements applied gradually, Low frequency reaction, Oscillation and Vibration

For quick events, the explicit solver is excellent (mostly less than 0.1 second). The explicit solver determines the size of the time increment and sets it. The density of the substance affects the speed of sound in that medium. Mass scaling refers to this. The mesh (element size and quality), young modulus, and density all affect the explicit time step.

1. By inverting the mass matrix, an unknown x is discovered (M).
2. Direct answers.
3. Methods of central difference are employed.
4. Less precise
5. Efficient in terms of computation for big models and quick Dynamic reaction time
6. Less computer storage is needed.
7. Equal to or greater than 10-3 strains per second.
8. An excellent choice for structures subject to wave propagation, car crashes, blast effects, drop test and Stamping of metal

This question doesn't have a simple solution. The type of issues solved by implicit and explicit solutions are same. Technically, both should result in the same result in any situation. The same issue can be examined using both methods. The method of time incrementation is the key distinction between implicit and explicit.

• The time increment can be defined thanks to the implicit analysis. Users can specify the size of the time increment. For global equilibrium, this increment must be completed iteratively. Consequently, it takes some time to calculate.
• Calculation of explicit time increments happens quite quickly. Simply because global equilibrium does not require iteration. Additionally, there will be no manual adjustment of the time increment. Simple assumption is made by the solver, which results in a very minimal time increment.