All Courses
All Courses
Courses by Software
Courses by Semester
Courses by Domain
Tool-focused Courses
Machine learning
POPULAR COURSES
Success Stories
AIM: To solve the equation F(u)=u3+9u2+4uF(u)=u3+9u2+4u using both explicit and implicit methods with a residual tolerance of 0.01. In the equation, F-Force in the bar and u-displacement in the bar. The problem is solved by both explicit and implicit analysis methods. In both cases an incremental load control scheme is used.…
Amith Anoop Kumar
updated on 09 Sep 2021
AIM:
To solve the equation F(u)=u3+9u2+4uF(u)=u3+9u2+4u using both explicit and implicit methods with a residual tolerance of 0.01. In the equation, F-Force in the bar and u-displacement in the bar. The problem is solved by both explicit and implicit analysis methods. In both cases an incremental load control scheme is used.
Clearly, it is possible to plot a graph of force, f, versus displacement, u, by using the force equation. However, if the nonlinear force versus displacement relationship is not known, but the stiffness is known, it is possible to construct the graph of force versus displacement numerically. One method is to use an incremental explicit load control scheme wherein, equal load increments (or steps) of dF = 1 unit each are used to load the tension bar. Similarly, an implicit method can also be used to solve the problem by using Newton-Raphson iterations to establish convergence after each step. Solve the problem using both methods and plot the Force v/s displacement graph for both methods
THE EXPLICIT AND IMPLICIT SCHEMES
The explicit and implicit methods are simply two different approaches used to solve ODEs and PDEs used to depict physical processes as computer simulations. While the Explicit scheme uses the value from the present time step to calculate the value in the next timestep, the implicit scheme calculates the value in the next time step by using the present values and expressing the source term as a function of the next time step. This creates a need for Matrix inversion in the case of the implicit method. Having said all this, we must note that these schemes are just two different ways of solving the same problem. A fundamental difference between the two schemes is that the Explicit scheme assumes convergence for all time steps and thus doesn't check for convergence after each timestep. On the other hand, the implicit scheme checks for convergence after every timestep before proceeding to the next. This makes the implicit scheme stable but demanding in terms of computations as most of the time, convergence is reached only after several iterations. As a result, the time taken for computing for a single time step is bigger. The explicit scheme can be considered as more direct as there are no convergence checks done in between but this generates stability problems. To avoid this, the timesteps have to be kept extremely small and in most cases, the user doesn't have much control over the selection of this extremely small timestep. This is not the case for the implicit method. Since the convergence check is performed after every time step, the user has control over how big each timestep should be, in the implicit method.
Overview on Explicit and Implicit analysis: We will use Implicit and Explicit solvers to solve dynamic problems. Implicit and Explicit analysis differ in the approach to time increments.. In Implicit analysis each time increment has to converge, but we can set pretty long time increments. Explicit on the other hand doesn’t have to converge each increment, but for the solution to be accurate time increments must be super small. In general, when we want to solve a dynamic problem, we are asking the solver to calculate the series of time increments, to see what is going on in the model. In each of those increments, time slightly increases. This of course also changes the load in the model. This is why we need to have a time-dependent chart for load. So the solver knows what load should be applied at which time increment. To this extent, it works the same as load increment in static analysis. Also, since we are incrementing time, inertia effects can be taken into account when things are happening fast enough. To this extent, both implicit and explicit analysis does more or less the same thing. They divide the analysis time into small increments, and then analyze what is going on in our model one increment of time after another. Explicit analysis: Explicit analysis works in a way, that we don’t have to converge each step. There are no convergence criteria to check and no iterations. Instead of verifying the global equilibrium, solver assumes that the equilibrium simply exists.This means that the solver goes straight into calculation of local finite element variables.It calculates all of the local finite element variables for given increment, and moves to the next one. Time increment in explicit analysis is very short. It’s the time in which an elastic wave (i.e. impact wave) can travel through a single finite element of your model. Of course, this is verified for all finite elements in your model, and the smallest time is selected. Implicit analysis: Implicit analysis works in a way,that in each time increment, the global equilibrium in the model is established. This means that each increment has to converge (this is happening in iterations). Since global equilibrium is verified at each time increments, those increments can be large.Each time increment computes slowly, since iterations are needed to get to the global equilibrium. Explicit scheme: If the nonlinear force versus displacement relationship is not known, but the stiffness is known, it is possible to construct the graph of force versus displacement numerically. It can be done using an incremental explicit load control scheme. The following equation is solved using explicit method as follows, F(u)=u3+9u2+4uF(u)=u3+9u2+4u Differentiating force with respect to displacement gives the stiffness(k) of the bar (since F = kx). k(u)=dFdu=3u2+18u+4k(u)=dFdu=3u2+18u+4 Here,three equal load increments (or steps) of ΔF=1ΔF=1 unit each are used to load the tension bar. The variables used are u for displacement, f for internal force in the bar, F for external force applied to the bar and k for stiffness. Incremental displacements or incremental externally applied forces are represented as ΔuΔu or ΔFΔF, respectively. ΔF=k.ΔuΔF=k.Δu This relation is used to find the incremental displacement.
USING THE EXPLICIT SCHEME
Force function, F(u)=(u)3+9(u)2+4uF(u)=(u)3+9(u)2+4u
Stiffness, k(u)=dduF(u)=3u2+18u+4k(u)=dduF(u)=3u2+18u+4
Displacement = Force/Stiffness=Fk=Fk
Force increment=dF=1=dF=1 Displacement increment, du=dFkdu=dFk
Displacement, ui+1=ui+duiui+1=ui+dui
Step 1: dF=1
k(u0)=3(u0)2+18u0+4=3â‹…02+18â‹…0+4=4
du1=dFk(u0)=14=0.25
u1=u0+du1=0+0.25=0.25
F(u1)=(u1)3+9(u1)2+4u1=0.253+9â‹…0.252+4â‹…0.25=1.578
Fext=dF=1 uexact=0.177 { by solving the equation, Fext=(uexact)3+9(uexact)2+4uexact}
Step 2: dF=1
k(u1)=3(u1)2+18u1+4=3â‹…0.252+18â‹…0.25+4=8.688
du2=dFk(u1)=18.688=0.115
u2=u1+du2=0.25+0.115=0.365
F(u2)=(u2)3+9(u2)2+4u2=0.3653+9â‹…0.3652+4â‹…0.365=2.709
Fext=dF+dF=1+1=2
uexact=0.296 { by solving the equation, Fext=(uexact)3+9(uexact)2+4uexact}
Step 3: dF=1
k(u2)=3(u2)2+18u2+4=3â‹…0.3652+18â‹…0.365+4=10.972
du3=dFk(u2)=110.972=0.091
u3=u2+du3=0.365+0.091=0.456
F(u3)=(u3)3+9(u3)2+4u3=0.4563+9â‹…0.4562+4â‹…0.456=3.793
Fext=dF+dF+dF=1+1+1=3
uexact=0.391 { by solving the equation, Fext=(uexact)3+9(uexact)2+4uexact}
It is evident from the results that external and internal forces are NOT in equilibrium because Fexternal≠FInternal. The following table summarizes the results for explicit analysis method.
Implicit scheme: An implicit analysis is the same as the explicit analysis, except that at the end of each step Newton-Raphson iterations are used to enforce equilibrium before moving to the next step. Basically, an incremental force is applied to advance the solution forward at the beginning of a step. However, internal forces and external forces will not be in equilibrium unless the stiffness is linear for the given step. Hence, in order to achieve equilibrium, corrections must be made to the displacement. This is accomplished by using Newton-Raphson iterations to minimize the residual, R=FInternal-Fexternal
Expanding the residual as a Taylor series about the current displacement uj gives
Rj+1=Rj(uj)+dR(uj)duδuj+1+...=Rj+k(uj)δuj+1+...
By neglecting higher order terms in the series, setting Rj+1=0 and solving for δuj+1 the following correction is obtained
δuj+1=−[k(uj)]−1⋅Rj Here, the iteration variable is j and it may take several iterations for the norm of the residual, R, to be reduced below the chosen tolerance criteria of 0.01.
The implicit analysis proceeds as follows
(Step- 1)
u0=0.0,⇒k(u0)=4,ΔF=1,Δu1=ΔFk(u0)=0.25,u1=u0+Δu1=0.25
Checking the residual, Fexternal=ΔF=1
FInternal=u31+9u21+4u1=1.578
R(0)=FInternal-Fexternal=0.578
R(0)>10−2⇒Newton-Raphson iterations are necessary.
Calculate the correction to u1=u(0)1
δu(1)=−[k(u(0)1)]−1⋅R(0)=-(8.6875)-1⋅0.578=-0.0665 ⇒The updated value of u(1)1=u(0)1+δu(1)=0.25-0.0665=0.1835
Checking the residual again,
Fexternal=ΔF=1
FInternal=(u(1)1)3+9(u(1)1)2+4u(1)1=1.0432
R(1)=FInternal-Fexternal=0.0432
R(1)>10−2⇒Newton-Raphson iterations are necessary
Calculate the new total correction to u1=u(0)1
δu(2)=δu(1)−[k(u(1)1)]−1⋅R(1)=-0.0665-(7.404)-1⋅0.0432=-0.0723
The updated value of u(2)1=u(0)1+δu(2)=0.25-0.0723=0.1777
Checking the residual again,
Fexternal=ΔF=1
FInternal=(u(2)1)3+9(u(2)1)2+4u(2)1=1.0006
R(2)=FInternal-Fexternal=0.0006 R(2)<10−2⇒Newton-Raphson iterations are not necessary Therefore, the final value is u1=u(2)1=0.1777
Step-2̲
u1=0.1777,⇒k(u1)=7.2933,ΔF=1,Δu2=ΔFk(u1)=0.1371,u2=u1+Δu2=0.3148
Checking the residual, Fexternal=ΔF+ΔF=2
FInternal=u32+9u22+4u2=2.1822
R(0)=FInternal-Fexternal=0.1822
R(0)>10−2⇒Newton-Raphson iterations are necessary.
Calculate the correction to u1=u(0)1
δu(1)=−[k(u(0)2)]−1⋅R(0)=-(9.9637)-1⋅0.1822=-0.0182 ⇒The updated value of u(1)2=u(0)2+δu(1)=0.3148-0.0182=0.2966
Checking the residual again, Fexternal=ΔF+ΔF=2 FInternal=(u(1)2)3+9(u(1)2)2+4u(1)2=2.0042 R(1)=FInternal-Fexternal=0.0042 R(1)<10−2⇒Newton-Raphson iterations are not necessary
Therefore, the final value is u2=u(1)2=0.2966
Step-3̲
u2=0.2966,⇒k(u2)=9.6027,ΔF=1,Δu3=ΔFk(u2)=0.1041,u3=u2+Δu3=0.4007 Checking the residual, Fexternal=ΔF+ΔF+ΔF=3
FInternal=u33+9u23+4u3=3.1121 R(0)=FInternal-Fexternal=0.1121 R(0)>10−2⇒Newton-Raphson iterations are necessary.
Calculate the correction to u2=u(0)2
δu(2)=−[k(u(0)3)]−1⋅R(0)=-(11.6942)-1⋅0.1121=-0.00958
⇒The updated value of u(1)3=u(0)3+δu(2)=0.4007-0.00958=0.3911
Checking the residual again, Fexternal=ΔF+ΔF+ΔF=3
FInternal=(u(1)3)3+9(u(1)3)2+4u(1)3=3.0008
R(1)=FInternal-Fexternal=0.0008 R(1)<10−2⇒Newton-Raphson iterations are not necessary Therefore, the final value is u3=u(2)3=0.3911
The results of the implicit analysis are summarized in the following table.
It is evident from the table that the external forces and internal forces are coming into equilibrium.
It is evident from the plot that the explicit analysis drifts from the exact solution. The drift from the exact solution illustrates the lack of equilibrium between the internal and external forces. To correct for this problem an implicit analysis is required. The Newton-Raphson iterations correct the incremental steps so that they land on the exact solution according to the specified tolerance. Excellent agreement is achieved with just 3 increments and in each step only two Newton-Raphson iterations are required.
CONCLUSION:
As experimented we can use both implicit and explicit solvers to solve the same type of problems. The only difference is in approach to time increment. The implicit analysis allows you to select how big the time increment should be. This increment will take some time to compute since it has to iterate for global equilibrium first. But we can make it reasonably big for our analysis. Explicit time increments calculate really fast because they do not iterate for global equilibrium, but the time increment there is not the user’s choice. Solver simply assumes an acceptable time increment and goes with it. The speed of computing is exactly the difference between implicit and explicit. Of course, the size of explicit step depends on our model and computing times depend on the system hardware. The debate on what type of solver is better depends upon the type of problem and accuracy we need.
Thus, a nonlinear analysis is illustrated using an incremental explicit and implicit procedure. The results are plotted for comparison. It is evident from the results that the explicit analysis drifts from the true solution. To overcome this problem an implicit analysis is used, which includes Newton-Raphson iterations to enforce equilibrium between internal and external forces. The techniques demonstrated are for a load control scheme.Once satisfactory equilibrium is achieved, a new load increment is applied to advance the solution and the process is repeated.
Leave a comment
Thanks for choosing to leave a comment. Please keep in mind that all the comments are moderated as per our comment policy, and your email will not be published for privacy reasons. Please leave a personal & meaningful conversation.
Other comments...
Crash Box Simulation Using LS-DYNA
CRASH BOX SIMULATION AIM: To create a complete simulation file for-crash analysis of crash box having two different thickness of 1.2 mm and 1.5 mm from the given FE model of crashbox and produce the following deliverables. Input. k file and output files (D3PLOT, GLSTAT, SLEOUT, RCFORC) Animation of the final…
25 Oct 2021 06:05 AM IST
Mobile Drop test Challenge
DROP TEST CHALLENGEOBJECTIVE: To build the input deck for a drop test simulation from scratch. Note: The orientation of the cellphone is also important and there are specific guidelines on how to orient the products but for this assignment, we will ignore that. The objective of this assignment is not to get a correct simulation…
16 Sep 2021 05:33 AM IST
Week - 2 - Explicit and Implicit Analysis
AIM: To solve the equation F(u)=u3+9u2+4u using both explicit and implicit methods with a residual tolerance of 0.01. In the equation, F-Force in the bar and u-displacement in the bar. The problem is solved by both explicit and implicit analysis methods. In both cases an incremental load control scheme is used.…
09 Sep 2021 06:46 AM IST
Simulating Combustion of Natural Gas.
AIM: Perform a combustion simulation on the combustor model and plot the variation of the mass fraction of the different species’ in the simulation using line probes at different locations of the combustor as shown in Fig. You need to plot for CO2, H2O, CH4, N2, O2, NOx emissions & Soot formation. THEORY: Combustion of…
16 Aug 2021 02:06 PM IST
Related Courses
0 Hours of Content
Skill-Lync offers industry relevant advanced engineering courses for engineering students by partnering with industry experts.
© 2025 Skill-Lync Inc. All Rights Reserved.