Loading [MathJax]/jax/output/SVG/jax.js
Research article

Reconstruction of shear force in Atomic Force Microscopy from measured displacement of the cone-shaped cantilever tip

  • Received: 13 June 2023 Revised: 01 November 2023 Accepted: 19 January 2024 Published: 25 January 2024
  • We present a new comprehensive mathematical model of the cone-shaped cantilever tip-sample interaction in Atomic Force Microscopy (AFM). The importance of such AFMs with cone-shaped cantilevers can be appreciated when its ability to provide high-resolution information at the nanoscale is recalled. It is an indispensable tool in a wide range of scientific and industrial fields. The interaction of the cone-shaped cantilever tip with the surface of the specimen (sample) is modeled by the damped Euler-Bernoulli beam equation ρA(x)utt +μ(x)ut+(r(x)uxx+κ(x)uxxt)xx=0, (x,t)(0,)×(0,T), subject to the following initial, u(x,0)=0, ut(x,0)=0 and boundary, u(0,t)=0, ux(0,t)=0, (r(x)uxx(x,t)+κ(x)uxxt)x==M(t), ((r(x)uxx+κ(x)uxxt)x)x==g(t) conditions, where M(t):=2hcosθg(t)/π is the moment generated by the transverse shear force g(t). Based on this model, we propose an inversion algorithm for the reconstruction of an unknown shear force in the AFM cantilever. The measured displacement ν(t):=u(,t) is used as additional data for the reconstruction of the shear force g(t). The least square functional J(F)=12u(,)ν2L2(0,T) is introduced and an explicit gradient formula for the Fréchet derivative of the cost functional is derived via the weak solution of the adjoint problem. Additionally, the geometric parameters of the cone-shaped tip are explicitly contained in this formula. This enables us to construct a gradient based numerical algorithm for the reconstructions of the shear force from noise free as well as from random noisy measured output ν(t). Computational experiments show that the proposed algorithm is very fast and robust. This creates the basis for developing a numerical "gadget" for computational experiments with generic AFMs.

    Citation: Alemdar Hasanov, Alexandre Kawano, Onur Baysal. Reconstruction of shear force in Atomic Force Microscopy from measured displacement of the cone-shaped cantilever tip[J]. Mathematics in Engineering, 2024, 6(1): 137-154. doi: 10.3934/mine.2024006

    Related Papers:

    [1] Mufit San, Seyma Ramazan . A study for a higher order Riemann-Liouville fractional differential equation with weakly singularity. Electronic Research Archive, 2024, 32(5): 3092-3112. doi: 10.3934/era.2024141
    [2] Eman A. A. Ziada, Hind Hashem, Asma Al-Jaser, Osama Moaaz, Monica Botros . Numerical and analytical approach to the Chandrasekhar quadratic functional integral equation using Picard and Adomian decomposition methods. Electronic Research Archive, 2024, 32(11): 5943-5965. doi: 10.3934/era.2024275
    [3] Qingcong Song, Xinan Hao . Positive solutions for fractional iterative functional differential equation with a convection term. Electronic Research Archive, 2023, 31(4): 1863-1875. doi: 10.3934/era.2023096
    [4] Jinheng Liu, Kemei Zhang, Xue-Jun Xie . The existence of solutions of Hadamard fractional differential equations with integral and discrete boundary conditions on infinite interval. Electronic Research Archive, 2024, 32(4): 2286-2309. doi: 10.3934/era.2024104
    [5] Seda IGRET ARAZ, Mehmet Akif CETIN, Abdon ATANGANA . Existence, uniqueness and numerical solution of stochastic fractional differential equations with integer and non-integer orders. Electronic Research Archive, 2024, 32(2): 733-761. doi: 10.3934/era.2024035
    [6] Meng-Xue Chang, Bang-Sheng Han, Xiao-Ming Fan . Global dynamics of the solution for a bistable reaction diffusion equation with nonlocal effect. Electronic Research Archive, 2021, 29(5): 3017-3030. doi: 10.3934/era.2021024
    [7] Vo Van Au, Jagdev Singh, Anh Tuan Nguyen . Well-posedness results and blow-up for a semi-linear time fractional diffusion equation with variable coefficients. Electronic Research Archive, 2021, 29(6): 3581-3607. doi: 10.3934/era.2021052
    [8] Mustafa Aydin, Nazim I. Mahmudov, Hüseyin Aktuğlu, Erdem Baytunç, Mehmet S. Atamert . On a study of the representation of solutions of a Ψ-Caputo fractional differential equations with a single delay. Electronic Research Archive, 2022, 30(3): 1016-1034. doi: 10.3934/era.2022053
    [9] Jun Zhou . Initial boundary value problem for a inhomogeneous pseudo-parabolic equation. Electronic Research Archive, 2020, 28(1): 67-90. doi: 10.3934/era.2020005
    [10] Yang Cao, Qiuting Zhao . Initial boundary value problem of a class of mixed pseudo-parabolic Kirchhoff equations. Electronic Research Archive, 2021, 29(6): 3833-3851. doi: 10.3934/era.2021064
  • We present a new comprehensive mathematical model of the cone-shaped cantilever tip-sample interaction in Atomic Force Microscopy (AFM). The importance of such AFMs with cone-shaped cantilevers can be appreciated when its ability to provide high-resolution information at the nanoscale is recalled. It is an indispensable tool in a wide range of scientific and industrial fields. The interaction of the cone-shaped cantilever tip with the surface of the specimen (sample) is modeled by the damped Euler-Bernoulli beam equation ρA(x)utt +μ(x)ut+(r(x)uxx+κ(x)uxxt)xx=0, (x,t)(0,)×(0,T), subject to the following initial, u(x,0)=0, ut(x,0)=0 and boundary, u(0,t)=0, ux(0,t)=0, (r(x)uxx(x,t)+κ(x)uxxt)x==M(t), ((r(x)uxx+κ(x)uxxt)x)x==g(t) conditions, where M(t):=2hcosθg(t)/π is the moment generated by the transverse shear force g(t). Based on this model, we propose an inversion algorithm for the reconstruction of an unknown shear force in the AFM cantilever. The measured displacement ν(t):=u(,t) is used as additional data for the reconstruction of the shear force g(t). The least square functional J(F)=12u(,)ν2L2(0,T) is introduced and an explicit gradient formula for the Fréchet derivative of the cost functional is derived via the weak solution of the adjoint problem. Additionally, the geometric parameters of the cone-shaped tip are explicitly contained in this formula. This enables us to construct a gradient based numerical algorithm for the reconstructions of the shear force from noise free as well as from random noisy measured output ν(t). Computational experiments show that the proposed algorithm is very fast and robust. This creates the basis for developing a numerical "gadget" for computational experiments with generic AFMs.



    As the continuous advancement of current industrial production systems in the direction of complexity, large-scale and intelligence, their functions are highly integrated, performance parameters are continuously improved, environmental conditions are more severe, and working loads are complex and changeable. Accurately assessing the reliability and resilience of complex systems is fundamental to reliable operation and safe service of the system.

    Many scholars have studied the performance change and reliability in shock environments. Ma and Trivedi [1] introduced a process to analyze the changes in performance and reliability of wireless broadcasting networks under different network parameters. Zhao et al. [2] analyzed the state changes of multi-state systems after various types of shocks, and proposed various replacement and maintenance strategies after failure. Dhulipala et al. [3] proposed a post hazard-event recovery model to forecast performance recovery trajectories from a single hazard to multi-hazard. Dehghani et al. [4] discussed the performance recovery situation of industrial systems subjected to stochastic hazards. [5] presented a novel reliability model for composite insulators, which considers the amplitude and duration of detrimental shocks, as well as the impact of self-recovery processes on the performance changes of composite insulators. Wan et al. [6] used a model to describe the dynamic performance degradation. Wang et al. [7,8] proposed more universal continuous and discrete degradation and shock models based on time and system state. Wang et al. [9] deduced the reliability of dynamic equilibrium systems and static equilibrium systems with multi-state protection equipment affected by external shock and internal degradation. Dong et al. [10] classified faults into hard fault processes and soft fault processes, and established a more generalized system reliability model under shock conditions. Zhang et al. [11] modeled the actuator degradation process through the Wiener process, and derived the corresponding actuator reliability. Zhao et al. [12] obtained indexes related to reliability indexes for multi-state equilibrium systems when subjected to external shocks. Wang et al. [13] assessed the failure probability and reliability of a k-out-n: F capacity-balancing system under shock environments with multiple sources. Ranjkesh et al. [14] assessed the system reliability considering the correlation between the appearance time of random shocks and shock strength. Anwar et al. [15] evaluated the stress-strength reliability of complex systems under the unified progressive mixed truncation scheme. Song et al. [16] analyzed the fault process and reliability of a multi-state k: F system with performance sharing, balanced requirements, and protection strategy operation mechanisms.

    Some scholars have studied the resilience and costs of complex systems. In [17], a resilience model integrating the effects of system recovery behavior, system hardening, and the severity of shock events was established to quantify the ability of microgrids to withstand shocks. Zeng et al. [18] conducted a time-static and time-dependent resilience model for energy systems under extreme shock events. Liu et al. [19] used a modified projection tracing model to assess the resilience for metro stations to flood hazards from three aspects: defense, restoration, adaptation. Somayajulu et al. [20] modeled changes in infrastructure performance under multiple disaster events. Tang et al. [21] expanded the optimization model of linear programming to model the resilience for urban railway transport systems under multiple shocks. Levitin et al. [22] evaluated and optimized the expected costs and task reliability of the warm backup system under external shocks. Dui et al. [23] discussed the relationship between resilience, importance measure and maintenance. [24] combined maintenance costs related to component and system failures with importance measure theory. Andrzejczak et al. [25] modeled corrective maintenance costs and replacement costs of vehicle and developed a simple model of vehicle fault. [26] analyzed the performance changes under the influence of maintenance costs and repair time simultaneously. Teixeira et al. [27] used generalized stochastic Petri nets to model offshore wind turbine systems and analyzed the changes in system performance and cost under various maintenance strategies. An optimal model for aircraft scheduling on air corridor ramps based on time and cost was presented [28]. Dui et al. [29] analyzed the fault and maintenance cost management based on the importance measures. Given the costs related to energy and time, Chang et al. [30] proposed effective protocols for complex networks with pinning control. [31] referred to the costs of transforming system states as response costs, and classified it into operating costs incurred by executing actions and the effect costs that can measure the consequences of shocks. Yan et al. [32] proposed three optimization models aimed at balancing the bonus of discovering defects and the costs related to ship survey.

    However, existing studies have mainly focused onthe performance degradation process and the maintenance process under one shock, or the performance degradation process after multiple shocks, and do not consider the intersecting process of multiple shocks and maintenance activities. In this paper, considering that the system performance will be affected by shocks again in the process of maintenance, the performance change processes of the system are analyzed. Then, system reliability and resilience are studied based on performance loss and recovery under the intersecting process of multiple shocks and maintenance. Finally, the direct and indirect economic losses caused by multiple shocks, as well as maintenance costs, are analyzed.

    The reminder of this paper is structured as follows. Section 2 discusses performance change and reliability of complex systems under the intersecting process of many shocks and maintenance activities. In Section 3, we analyze the performance change process of complex systems under performance degradation and recovery processes. Then, the resilience and total costs models are analyzed. Section 4 illustrates the above models with a specific example of a welding robot system and Section 5 concludes this paper.

    Consider external factors such as pressure, load changes, sudden changes in temperature and humidity, electromagnetic fluctuations, and vibration shocks as shocks that complex systems are subjected to during operation. Complex systems may be subjected to multiple shocks, and each shock will affect system performance. One reason for the performance degradation is the shocks, and the other is the natural degradation process such as fatigue and aging of the system itself. The natural degradation process is also influenced by external factors. The fault process of complex systems is shown in Figure 1.

    Figure 1.  Fault process.

    Assuming that the complex system has a total of (k + 1) states, the different states of the system are marked with integers between 0 and k. Among them, state 0 represents the optimal state, state k represents that the complex system is unable to run at all, and other states represent performance degradation states. The performance set corresponding to the states is {g0,g1,,gk}. G(t) represents the performance of the complex system at time t, G(t){g0,g1,,gk}. The probability that the complex system is in state m can be represented by pm(t). The probability set is P(t)=[p0(t),p1(t),,pk(t)].

    Even though the complex system does not receive the subsequent shock again, its performance will continue to deteriorate, and the time between state transitions during the process of performance deterioration will be shortened, owing to the continuous influence of shocks on performance. After the complex system's performance has been brought down to a fault condition, it must be repaired in order to increase the reliability. However, the system may experience shocks once more before its performance is restored to its peak level. This maintenance stops when a fresh shock occurs, and the system can only be restored after the system state drops to a failed state. The shocks and maintenance processes of the complex system intersect in this form, and the performance change of the system is shown in Figure 2.

    Figure 2.  Performance change under shock events.

    In Figure 2, g0 represents the optimal performance. gk represents the performance in the fault state. The complex system is initially performing at its optimal state. As the shocks occur, the system state will transition from state v to state z (0v<zk), and the values of system performance will also be reduced from gv to gz. Multiple shocks accelerate the transition time between states in a multi-state complex system. After multiple shocks, the performance state drops to the fault state gk. The system has to be repaired at this point. After maintenance, the system performance will be improved from gv to gz (0z<vk). Before the system performance recovers to its optimal level, it is subjected to a fourth shock, at which point maintenance is interrupted and the system performance begins to decrease again. Otherwise, the system will continue to be repaired until its performance recovers to the optimal performance g0.

    The effect of shocks and natural degradation on the performance of complex systems can be chiefly seen in the variation of the state degradation rate. Let τm,nx be the complex system's state degradation rate from state m to state n after the arrival of the x-th shock event. Considering how the effect on performance adds up, we define that the value of τm,nx is related to the initial state degradation rate τm,n1. The expression for τm,nx is

    τm,nx=(1+μmx)σxτm,n1, (1)

    where μmx is related to the system state just prior to the occurrence of the x-th shock event. μmx is not related to the number and strength of shocks, but is defined as a constant that changes with state m. σx is concerned with the strength of the x-th shock. The initial complex system is a system whose performance state is in its initial state without degradation or whose performance state has recovered to its best performance state after the shocks.

    Even if complex systems have not failed, running at very low performance can also affect normal operation. To ensure the efficiency of completing work, a performance threshold gh is defined. When the performance is less than or equal to the threshold, it is necessary to repair the complex system.

    The reliability here can be understood as the probability that the system performance is not lower than the performance threshold. So, the expression for the reliability is

    R(t)=P(G(t)gh)=hm=0pm(t). (2)

    Each shock will affect the system performance degradation rate and state transition process. Assume that both system performance degradation and maintenance processes change state by state and cannot cross intermediate states. τm,n(t) represents the performance degradation rate from state m to state n (m < n). θm,n(t) represents the performance recovery rate from state m to state n (m > n). The state transition process is shown in Figure 3.

    Figure 3.  State transition process.

    Assuming that the process of performance recovery does not depend on the quantity and strength of shocks, the performance recovery rate matrix remains the same regardless of the number of maintenances performed. The performance recovery rate matrix Wp(t) of complex systems from the beginning of the f-th maintenance to the time it is subjected to shocks again can be expressed as

    Wp(t)=[θ1,0(t)θ1,0(t)00000θk1,k2(t)θk1,k2(t)0000θk,k1(t)θk,k1(t)].
    tptmin{min(tx|txtp),min(ty|tytp)},x,p=1,2,3,

    where tp represents the starting time of the f-th maintenance, tx represents the arrival time of the x-th shock event occurs and ty represents the time when the system is repaired to its best performance state. Under various shocks, it may be repaired multiple times, so there may be multiple values for ty. The element with a value of 0 in the matrix Wp(t) represents that the system cannot change from a high state to a lower state, meaning that system performance cannot be reduced but can only be improved.

    When x=1,the initial performance degradation rate matrix B1(t) after the arrival of the first shock and before the arrival of the second shock or before repair is

    B1(t)=[τ0,11(t)τ0,11(t)0000τ1,21(t)τ1,21(t)00000τk1,k1(t)τk1,k1(t)].
    t1tmin{t2,min(tp|tpt1)},p=1,2,3,

    The element with a value of 0 in the B1(t) represents that the complex system can only undergo state degradation after being subject to shocks, meaning that the performance cannot be improved but can only be reduced.

    Under the background of multiple shocks, the performance degradation rate matrix Bx(t) from the arrival of the x-th shock to the arrival of the (x + 1)-th shock or before repair is

    Bx(t)=[τ0,1x(t)τ0,1x(t)0000τ1,2x(t)τ1,2x(t)00000τk1,kx(t)τk1,kx(t)].
    x>1,txtmin{tx+1,min(tp|tptx)},x,p=1,2,3,

    From Eq (1), the expression for Bx(t) is

    Bx(t)=(1+μmx)σxB1(t) (3)

    Under the given initial condition P(t)=[1,0,,0], the probability of the complex system being in different states could be obtained by Eq (4).

    dP(t)dt=P(t)M, (4)

    where M represents the matrix Wp(t) or matrix Bx(t),x1 of complex systems.

    At time t, the average performance ˉG(t) of complex systems is

    ˉG(t)=km=0pm(t)gm. (5)

    Set tf as the time when the performance drops to gh for the first time after multiple shocks. When time is continuous, the performance loss GL(tf) can be expressed as

    GL(tf)=tf0[g0max(gh,ˉG(t))]dt. (6)

    Set tb as the beginning time and te as the ending time for maintenance. When t[tb,te], the expression for performance recovery Gr(t) is

    Gr(t)=ttb(ˉG(t)gh)dt. (7)

    The resilience RI(t) of complex systems under continuous time can be defined as

    RI(t)=Gr(t)GL(tf)=ttb(ˉG(t)gh)dttf0[g0max(gh,ˉG(t))]dt. (8)

    When ˉG(t)=gh, RI(t)=0. This indicates that the system cannot be restored, and the system lacks any resilience. RI(t)=1 indicates that the system has returned to its optimal performance state. The closer the value of RI is to 0, the smaller the degree of performance recovery and the smaller the system resilience after maintenance.

    Complex systems incur various costs during shocks and maintenance processes. The losses directly caused by shocks are called direct economic costs. The system subjected to shocks will operate in a low-performance mode. The economic benefits generated by low-performance work will be lower than those produced by working at optimal performance. The losses of economic benefits under low performance are called indirect economic costs. For example, after a handling robot system is shocked, it may experience a decrease in speed, load capacity and accuracy, which will greatly reduce the work efficiency and thus reduce economic benefits. Repairing the complex system to higher performance also incurs economic losses, which are known as maintenance costs. The expression for the total costs generated by complex systems in (0, t) is

    C(t)=k1i=0αi,i+1Mi,i+1(t)+ki=0uiδi(t)+ki=1βi,i1Mi,i1(t), (9)

    where C(t) represents the total costs. ki=0αi,i+1Mi,i+1(t) represents the direct economic costs. αi,i+1 is the direct economic loss caused by the degradation from state i to state i + 1. Mi,i+1(t) represents how many times it has been transferred from state i to state i + 1 within time t. ki=0uiδi(t) represents the indirect economic costs. ui represents the benefit losses caused by the system running in state i per unit of time. δi(t) represents the total time the system has been in state i within time t. ki=1βi,i1Mi,i1(t) represents the maintenance costs. βi,i1 represents the costs of repairing from state i to state i-1. Mi,i1(t) represents how many times it has been transferred from state i to state i-1 within time t. The expression for δi(t) is

    δi(t)=t0Φ{G(t)=gi}dt. (10)

    If the system's performance G(t) at time t is identical to its performance in state i, then Φ{G(t)=gi} = 1, otherwise Φ{G(t)=gi} = 0.

    A welding robot, as a complex system, is mainly used in the production activity of various automotive parts and can replace manual work in welding, cutting, or thermal spraying. Positioning accuracy is commonly used to measure the working performance of a welding robot. In robot arc welding operations, due to changes in the actual welding environment, such as tooling errors and welding thermal deformation, the welding gun moving according to the original trajectory can no longer ensure accurate alignment of the weld seam, which can lead to a decrease in positioning accuracy and even the inability to maintain the normal welding process.

    Taking a six-axis welding robot as an example, suppose that the welding robot system is divided into five states: 0, 1, 2, 3, and 4. State 0 represents the optimal positioning accuracy of the system. At the initial running time, the welding robot is in the optimal performance state. State 4 represents the system has failed, indicating that the positioning accuracy error of the welding robot system is too large to finish tasks normally. States between 0 and 4 represent a gradual decrease in performance, meaning that the positioning accuracy gradually increases. To enable the welding robot to operate efficiently, state 3 is defined as the threshold state of the welding robot system. When the state is state 3 or state 4, the welding robot needs to be repaired. The positioning accuracy value corresponding to each state is j0=0.01mm, j1=0.02mm, j2=0.04mm, j3=0.06mm, j4=0.08mm. For the convenience of calculation, we simplify the positioning accuracy value jm,m=0,1,2,3,4 and obtain the performance value corresponding to positioning accuracy, which is g0=1.0,g1=0.8,g2=0.6,g3=0.3,g4=0. So, the performance thresholdgh=0.3.

    After the welding robot system is subjected to shocks or maintenance, the system state deteriorates or recovers one by one in the five states. Figure 4 shows the state degradation and recovery process.

    Figure 4.  The state degradation and recovery process of the welding robot system.

    Assuming that the performance degradation rate of the initial welding robot system is τ0,11=0.25a1,τ1,21=0.34a1,τ2,31=0.30a1,τ3,41=0.22a1, and the performance recovery rate is θ1,0=0.51a1,θ2,1=0.48a1,θ3,2=0.62a1,θ4,3=0.46a1, respectively. So, the matrix B1 and matrix Wp can be expressed as

    B1=[0.250.2500000.340.3400000.300.3000000.220.22],
    Wp=[0.510.5100000.480.4800000.620.6200000.460.46].

    μmx is only related to the state of the welding robot, and the values of μmx vary when the welding robot is in different states. Table 1 gives the values of μmx,m=0,1,2,3,4.

    Table 1.  The values of μmx.
    μmx Value
    μ0x 0.48
    μ1x 0.56
    μ2x 0.80
    μ3x 0.64
    μ4x 0

     | Show Table
    DownLoad: CSV

    Assume that the welding robot system is subjected to the first shock at the initial time, the second shock in year 4, and the third shock in year 25.

    The value of σx,x=2,3 depends on the strength of the second and third shock, with σ2=2.3 and σ3=2. According to Eq (1), it can be calculated that τ0,12=0.69a1,τ1,22=0.94a1,τ2,32=0.83a1,τ3,42=0.61a1and τ0,13=0.53a1,τ1,23=0.71a1,τ2,33=0.63a1,τ3,43=0.46a1. So, the performance degradation rate matrices after being subjected to the second and third shocks can be expressed as

    B2=(1+μm2)σ2[0.250.2500000.340.3400000.300.3000000.220.22]
    =[0.690.6900000.940.9400000.830.8300000.610.61]
    B3=(1+μm3)σ3[0.250.2500000.340.3400000.300.3000000.220.2200000]
    =[0.530.5300000.710.7100000.630.6300000.460.46]

    At first, the welding robot system is operating at its best performance, so P(0)=[1,0,0,0,0]. Then, according to Eq (4), the probability for the welding robot being in various states during three shocks and two maintenance processes can be obtained. Figure 5 shows the probability changes.

    Figure 5.  The probability changes in each state.

    In Figure 5, (a) represents the probability change process of each state of the welding robot system after the first random shock. In year 4, the system can still complete its work normally and does not require maintenance. Figure 5(b) shows that the system has degraded to the state that requires maintenance in year 18.15. Figure 5(c) shows that the system begins to be repaired in year 18.15, but has not yet returned to its optimal performance state by year 25. Figure 5(d) shows that the system is subjected to the third shock in year 25 and degrades again to the state that requires maintenance in year 42.25. Figure 5(e) shows that the system is repaired again in year 42.25, after which the performance of the welding robot system is no longer affected by the shock environments. Figure 5(f) shows the probability changes in each state during the entire process mentioned above.

    According to Eq (2), the reliability of welding robot system during three random shocks and two maintenance process can be obtained, as shown in Figure 6.

    Figure 6.  The reliability of the welding robot system.

    As shown in Figure 6, the welding robot system has the highest reliability at the initial time. After being subjected to the first and second shocks, the system reliability decreases until it reaches its lowest value in year 18.15. Then, the system is repaired, and the system performance begins to recover. The reliability gradually increases. At year 25, before the system reliability returns to its maximum value, it is subjected to the third shock, and then the system reliability gradually decreases again. After the reliability is reduced to its minimum value in year 42.25, the welding robot system is repaired again. Afterwards, the system performance is no longer affected by shocks.

    After obtaining the probability of the welding robot in various states, the average performance of the system can be obtained. According to Eq (8), the resilience curve of the welding robot system during two maintenance processes can be obtained, as shown in Figure 7.

    Figure 7.  The resilience of welding robot system.

    Before year 18.15, the system suffers two shocks without any maintenance behaviors, so resilience does not exist. Figure 7 shows that the resilience of welding robot system continuously increases after maintenance. Figure 7(a) shows the resilience change of the system during the first maintenance period. At year 25, the first maintenance is terminated and the system cannot return to its optimal state, resulting in a low resilience value. After the third shock, the system performance continues to decline. At year 42.25, ˉG(t)=gh, so RI(t)=0. Afterward, the system will no longer be shocked and its performance can continue to recover, so, as shown in Figure 7(b), the resilience continues to increase.

    Table 2 shows the parameter values related to the total costs.

    Table 2.  Parameter values related to the total costs.
    Parameter Value(¥) Parameter value(¥)
    α0,1 1.32×105 u3 1.52×105
    α1,2 2.18×105 u4 2.19×105
    α2,3 2.69×105 β1,0 1.76×105
    α3,4 2.15×105 β2,1 2.31×105
    u0 0 β3,2 2.07×105
    u1 8.20×104 β4,3 3.61×105
    u2 1.01×105

     | Show Table
    DownLoad: CSV

    According to Eq (9), the total costs of the welding robot system during three shocks and maintenance processes can be obtained, as shown in Figure 8.

    Figure 8.  The total costs of the welding robot system.

    In Figure 8, the first and third identification points represent the time points of the second and third shocks, respectively. The second and fourth identification points represent maintenance to the system, respectively. After being subjected to shocks, the system will operate in a lower state and the losses of benefits will increase. Therefore, after multiple shocks, the growth rate of the total costs will increase. After maintenance, the system performance will be improved and the losses of system benefits will be reduced. So, the growth rate of the total costs will decrease.

    We first analyze the performance change for complex systems under shocks and maintenance. Considering the effect of system state and shock strength on system performance, different state transition rates can be obtained. A resilience model is established based on the performance loss and performance recovery of the system under the cross-process of multiple shocks and maintenance. After that, the economic loss costs and maintenance costs are analyzed. Finally, combined with a case study of a welding robot system, the feasibility and effectiveness of the model can be verified, which has certain reference significance for resilience and cost analysis of complex systems.

    This paper assumes that the state transition rate is related to the strength of shocks and the system state. In the future, we can conduct sensitivity analysis on the strength parameter and state parameter, and analyze the effect of parameter values on the performance and costs of complex systems. In the future, other processes can be considered to describe the state transition process of complex systems.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The authors acknowledge the Program for young backbone teachers in Universities of Henan Province (No. 2021GGJS007).

    The authors declare there is no conflict of interest.



    [1] M. Antognozzi, Investigation of the shear force contrast mechanism in transverse dynamic force microscopy, Ph.D. Thesis, University of Bristol, UK, 2000.
    [2] M. Antognozzi, D. R. Binger, A. D. L. Humphris, P. J. James, M. J. Miles, Modeling of cylindrically tapered cantilevers for transverse dynamic force microscopy (TDFM), Ultramicroscopy, 86 (2001), 223–232. https://doi.org/10.1016/S0304-3991(00)00087-5 doi: 10.1016/S0304-3991(00)00087-5
    [3] H. T. Banks, D. J. Inman, On damping mechanisms in beams, J. Appl. Mech., 58 (1991), 716–723. https://doi.org/10.1115/1.2897253 doi: 10.1115/1.2897253
    [4] G. Binnig, C. F. Quate, C. Gerber, Atomic force microscopy, Phys. Rev. Lett., 56 (1986), 930. https://doi.org/10.1103/PhysRevLett.56.930 doi: 10.1103/PhysRevLett.56.930
    [5] W. J. Chang, T. H. Fang, C. I. Weng, Inverse determination of the cutting force on nanoscale processing using atomic force microscopy, Nanotechnology, 15 (2004), 427–430. https://doi.org/10.1088/0957-4484/15/5/004 doi: 10.1088/0957-4484/15/5/004
    [6] L. C. Evans, Partial differential equations, 2 Eds., Providence: American Mathematical Society, 2010.
    [7] B. Geist, J. R. McLaughlin, The effect of structural damping on nodes for the Euler-Bernoulli beam: a specific case study, Appl. Math. Lett., 7 (1994), 51–55. https://doi.org/10.1016/0893-9659(94)90112-0 doi: 10.1016/0893-9659(94)90112-0
    [8] A. Hasanov, O. Baysal, Identification of an unknown spatial load distribution in a vibrating cantilever beam from final overdetermination, J. Inverse Ill-Posed Probl., 23 (2015), 85–102. https://doi.org/10.1515/jiip-2014-0010 doi: 10.1515/jiip-2014-0010
    [9] A. Hasanov, O. Baysal, Identification of unknown temporal and spatial load distributions in a vibrating Euler-Bernoulli beam from Dirichlet boundary measured data, Automatica, 71 (2016), 106–117. https://doi.org/10.1016/j.automatica.2016.04.034 doi: 10.1016/j.automatica.2016.04.034
    [10] A. Hasanov, O. Baysal, C. Sebu, Identification of an unknown shear force in the Euler-Bernoulli cantilever beam from measured boundary deflection, Inverse Probl., 35 (2019), 115008. https://doi.org/10.1088/1361-6420/ab2a34 doi: 10.1088/1361-6420/ab2a34
    [11] A. Hasanov, O. Baysal, H. Itou, Identification of an unknown shear force in a cantilever Euler-Bernoulli beam from measured boundary bending moment, J. Inverse Ill-posed Probl., 27 (2019), 859–876. https://doi.org/10.1515/jiip-2019-0020 doi: 10.1515/jiip-2019-0020
    [12] A. H. Hasanoglu, A. G. Romanov, Introduction to inverse problems for differential equations, 2 Eds., New York: Springer, 2021.
    [13] G. Haugstad, Atomic force microscopy: understanding basic modes and advanced applications, Wiley, 2012.
    [14] O. Baysal, A. Hasanov, K. Sakthivel, Determination of unknown shear force in transverse dynamic force microscopy from measured final data, J. Inverse Ill-posed Probl., 2023. https://doi.org/10.1515/jiip-2023-0021 doi: 10.1515/jiip-2023-0021
    [15] S. Kumarasamy, A. Hasanov, A. Dileep, Inverse problems of identifying the unknown transverse shear force in the Euler-Bernoulli beam with Kelvin-Voigt damping, J. Inverse Ill-posed Probl., 2023. https://doi.org/10.1515/jiip-2022-0053 doi: 10.1515/jiip-2022-0053
    [16] H. L. Lee, W. J. Chang, Effects of damping on the vibration frequency of atomic force microscope cantilevers using the Timoshenko beam model, Jpn. J. Appl. Phys., 48 (2009), 065005. https://doi.org/10.1143/JJAP.48.065005 doi: 10.1143/JJAP.48.065005
    [17] J. S. Lee, J. Song, S. O. Kim, S. Kim, W. Lee, J. A. Jackman, et al., Multifunctional hydrogel nano-probes for atomic force microscopy, Nat. Commun., 7 (2016), 11566. https://doi.org/10.1038/ncomms11566 doi: 10.1038/ncomms11566
    [18] M. P. Scherer, G. Frank, A. W. Gummer, Experimental determination of the mechanical impedance of atomic force microscopy cantilevers in fluids up to 70 kHz, J. Appl. Phys., 88 (2000), 2912–2920. https://doi.org/10.1063/1.1287522 doi: 10.1063/1.1287522
    [19] K. Shen, D. C. Hurley, J. A. Turner, Dynamic behaviour of dagger-shaped cantilevers for atomic force microscopy, Nanotechnology, 15 (2004), 1582–1589. https://doi.org/10.1088/0957-4484/15/11/036 doi: 10.1088/0957-4484/15/11/036
    [20] J. A. Turner, J. S. Wiehn, Sensitivity of flexural and torsional vibration modes of atomic force microscope cantilevers to surface stiffness variations, Nanotechnology, 12 (2001), 322–330. https://doi.org/10.1088/0957-4484/12/3/321 doi: 10.1088/0957-4484/12/3/321
    [21] N. R. Wilson, J. V. Macpherson, Carbon nanotube tips for atomic force microscopy, Nature Nanotechnology, 4 (2009), 483–491. https://doi.org/10.1038/nnano.2009.154 doi: 10.1038/nnano.2009.154
    [22] K. Zhang, T. Nguyen, C. Edwards, M. Antognozzi, M. Miles, G. Herrmann, Real-time force reconstruction in a transverse dynamic force microscope, IEEE Trans. Ind. Electron., 69 (2022), 11403–11413. https://doi.org/10.1109/TIE.2021.3120487 doi: 10.1109/TIE.2021.3120487
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1154) PDF downloads(148) Cited by(0)

Figures and Tables

Figures(5)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog