
The aim of this paper is to compare the classical continuum mechanics and the peridynamic models in the structural analysis of a monolithic glass plate subjected to ball drop. Governing equations are recalled in order to highlight the differences and basic features of both approaches. In this study the behavior of glass is assumed to be linear-elastic and damage processes are ignored. The generalized Hooke's law is assumed within the classical theory, while the linear peridynamic solid constitutive model is applied within the peridynamic analysis. Mechanical models for the ball drop simulation are discussed in detail. An emphasis is placed on the discretization including finite element mesh, peridynamic node lattice and time stepping, as well as appropriate constraints and contact conditions in both finite element and non-local peridynamics models. Deflections of the plate after the ball drop are presented as functions of time and the results based on the finite element and peridynamic analysis are compared. Good agreements between the deflection values in selected points of the plate as well as deflection fields at several time points indicate, that the model assumptions for the non-local peridynamic analysis including the horizon size, the short-range force contact settings and the support conditions are well suited. The developed peridynamics models can be applied in the future to analyze damage patterns in glass plates.
Citation: Holm Altenbach, Oleksiy Larin, Konstantin Naumenko, Olha Sukhanova, Mathias Würkner. Elastic plate under low velocity impact: Classical continuum mechanics vs peridynamics analysis[J]. AIMS Materials Science, 2022, 9(5): 702-718. doi: 10.3934/matersci.2022043
[1] | Zhe Zou, Minjie Zheng . Design and stabilization analysis of luxury cruise with dynamic positioning systems based on sampled-data control. Mathematical Biosciences and Engineering, 2023, 20(8): 14026-14045. doi: 10.3934/mbe.2023626 |
[2] | Bin Hang, Beibei Su, Weiwei Deng . Adaptive sliding mode fault-tolerant attitude control for flexible satellites based on T-S fuzzy disturbance modeling. Mathematical Biosciences and Engineering, 2023, 20(7): 12700-12717. doi: 10.3934/mbe.2023566 |
[3] | Shuaiwen Sun, Zhijing Xu . Large kernel convolution YOLO for ship detection in surveillance video. Mathematical Biosciences and Engineering, 2023, 20(8): 15018-15043. doi: 10.3934/mbe.2023673 |
[4] | Xueyan Wang . A fuzzy neural network-based automatic fault diagnosis method for permanent magnet synchronous generators. Mathematical Biosciences and Engineering, 2023, 20(5): 8933-8953. doi: 10.3934/mbe.2023392 |
[5] | Wei Zou, Yanxia Shen, Lei Wang . Design of robust fuzzy iterative learning control for nonlinear batch processes. Mathematical Biosciences and Engineering, 2023, 20(11): 20274-20294. doi: 10.3934/mbe.2023897 |
[6] | Lixiang Zhang, Yian Zhu, Jie Ren, Wei Lu, Ye Yao . A method for detecting abnormal behavior of ships based on multi-dimensional density distance and an abnormal isolation mechanism. Mathematical Biosciences and Engineering, 2023, 20(8): 13921-13946. doi: 10.3934/mbe.2023620 |
[7] | Yue Song, Yi Zhang, Song Yang, Na Li . Investigation on stability and controller design for singular bio-economic systems with stochastic fluctuations. Mathematical Biosciences and Engineering, 2021, 18(3): 2991-3005. doi: 10.3934/mbe.2021150 |
[8] | Qianhong Zhang, Fubiao Lin, Xiaoying Zhong . On discrete time Beverton-Holt population model with fuzzy environment. Mathematical Biosciences and Engineering, 2019, 16(3): 1471-1488. doi: 10.3934/mbe.2019071 |
[9] | Chaoyue Wang, Zhiyao Ma, Shaocheng Tong . Adaptive fuzzy output-feedback event-triggered control for fractional-order nonlinear system. Mathematical Biosciences and Engineering, 2022, 19(12): 12334-12352. doi: 10.3934/mbe.2022575 |
[10] | Jinxin Du, Lei Liu . Adaptive fuzzy fixed time formation control of state constrained nonlinear multi-agent systems against FDI attacks. Mathematical Biosciences and Engineering, 2024, 21(3): 4724-4741. doi: 10.3934/mbe.2024207 |
The aim of this paper is to compare the classical continuum mechanics and the peridynamic models in the structural analysis of a monolithic glass plate subjected to ball drop. Governing equations are recalled in order to highlight the differences and basic features of both approaches. In this study the behavior of glass is assumed to be linear-elastic and damage processes are ignored. The generalized Hooke's law is assumed within the classical theory, while the linear peridynamic solid constitutive model is applied within the peridynamic analysis. Mechanical models for the ball drop simulation are discussed in detail. An emphasis is placed on the discretization including finite element mesh, peridynamic node lattice and time stepping, as well as appropriate constraints and contact conditions in both finite element and non-local peridynamics models. Deflections of the plate after the ball drop are presented as functions of time and the results based on the finite element and peridynamic analysis are compared. Good agreements between the deflection values in selected points of the plate as well as deflection fields at several time points indicate, that the model assumptions for the non-local peridynamic analysis including the horizon size, the short-range force contact settings and the support conditions are well suited. The developed peridynamics models can be applied in the future to analyze damage patterns in glass plates.
Ship dynamic positioning is a system used in ships and offshore vessels to automatically maintain the ship's positions and heading by thrusters and propulsion systems [1]. Eliminating the need for traditional anchoring, dynamic positioning ship (DPS) is equipped with a sophisticated system that allows it to stay in a specific location or follow a set track. And DPS has been used in various maritime applications, including offshore oil and gas exploration, subsea construction, cable laying, and so on. Recently, there is a wealth of research in the academic literature on the control for DPS ([2,3,4,5,6,7]), and the control algorithms and strategies include that proportional integral derivative (PID) control, adaptive control, intelligent control, etc. [2]. The research direction includes sensor fusion and position reference systems [3]; robust control [4]; fault tolerance control [5]; cooperative control and networked system [6]; simulation and validation [7] etc. For example, in [3], the technology of sensor fusion and the selection and optimization of position reference systems for DPS are discussed, which emphasizes the requirements of system accuracy and reliability. Reference [4] uses H∞ robust control techniques for DPS to obtain a robust controller with better performance than traditional dynamic positioning systems. Trajectory preservation is achieved through appropriate weighting functions. Reference [5] discusses the asymptotic stability conditions for DPS with partial actuator failure. The saturation proportional differential controller has been designed, and the asymptotic stability conditions are given by using the LKF approaches. In [6], a new event triggered collaborative DP control algorithm is proposed for multiple surface vessels. This algorithm adopts an undirected interactive topology to construct a communication network between surface vessels, so that all followers can obtain leader information. Reference [7] proposes a multi physics modeling method for a comprehensive electric propulsion system used for DPS, which can run in real-time and connect to control hardware.
In recent years, due to the advantages of discreteness, high efficiency, and low communication costs, the sampled-data control system has received widespread attention and it has been widely applicated in various systems, such as the hypersonic vehicle system [8,9], Lur'e system [10,11,12], neural networks system [13,14,15,16,17], Markovian jump system [18,19] and so on [20,21]. It can be applied to practical system such as fast train [22], unmanned aerial vehicle [23], large cruise liner [24,25], small unmanned ships [26], etc. Nowadays, the Takagi-Sugeno fuzzy system (TSFS) has received widespread attention due to their significant effectiveness in describing the complex nonlinear system, and considerable sampled-data control methods are reported for TSFS [27,28]. For example, in [29], the sampled-data control issue for TSFS with random time delay and actuator faults is studied using the input-delay method; In [30], the stability issue for nonlinear sampled-data TSFS is studied using the LKF method, which is validated by a single link robotic arm with an electric motor model. In [31], the stability characteristics of TSFS under sampling-data control are studied. The system is considered as a variable time delay system based on the input delay method. An appropriate LKF candidates and quadratic delay product terms are provided. Reference [32] studies the nonfragile sampling control problem for TSFS with parameter uncertainty. A new augmented LKF with sufficient sample data information has been constructed, which results in less conservative results.
Recently, the research on sampled-data control systems for DPS has attracted widespread attention. In [33], A semi-globally practically asymptotically stable controller is designed using the integration and backstepping technique of the Euler approximation models. In [34], based on the Euler approximation model and nonlinear sampling control theory, a state feedback controller and an output feedback controller based on a reduced order observer are designed for DPS to stabilize the tracking error. In [35], the sampling control system was applied to the dynamic positioning ship, and the input delay method was combined with the Wirtinger integral inequality to obtain results with less conservatism. In [36], the control problem of a sampling system for DPS with actuator failure is discussed, and a sampling controller is designed to ensure the exponential stability of the system and meet the H∞ performance index. In [37], the fuzzy robust stability criterion of sampled-control systems for DPS is studied, and a fuzzy sampling controller is designed using an input time delay method to ensure that the system remains stable under external disturbances. In [38], the reliable tracking control problem for DPS based on nonperiodic measurement information is studied and the mean square exponential stability criterion is obtained.
However, there is still much room for further study for the sampled-data DPS. Firstly, some important characteristics are still overlooked, which include the features of sampling pattern and the effective characteristics of the constructed LKF. Secondly, the positivity limit of the Lyapunov matrix is strict, which will result in the conservatism of the system.
Inspired by the above discussions, we aim to study the stability condition and controller design for sampled-data DPS with the T-S fuzzy model. In addition, we will construct an improved LKF by adding the new useful terms and introducing the reciprocally convex method. The main advantages of this article are be listed that:
(1) To capture the features of sampling pattern fully, some novel terms, like (tk+1−t)t∫tke2αs[x(s)x(tk)]T[R1R2∗R3][x(s)x(tk)]ds are introduced in the constructed LKF.
(2) The tighter technique of reciprocally convex combination is introduced during the differentiation of LKF. Hence, the conservatism of stability criteria can be greatly reduced.
(3) A second-order term (tk+1−t)(t−tk)[xT(tk)N1x(tk)] related to time t is added during the construction of V(t), thereby a less conservative result is obtained.
Notations: Rm denotes m-dimensional Euclidean space, and " T " represents the matrix transposition. "*" is a symmetric block matrix.
Consider the dynamic equations for DPS as follows:
M˙υ(t)+Dυ(t)+Gη(t)=u(t)+w(t),˙η(t)=J(ψ(t))υ(t), | (1) |
where
M=[m11000m22m230m32m33],D=[n11000n22n230n32n33],G=[g11000g22000g33],J(ψ(t))=[cos(ψ(t))−sin(ψ(t))0sin(ψ(t))cos(ψ(t))0001]. |
where η(t)=[xa(t)ya(t)ψ(t)]T is the positions and heading angle in the northeast coordinate system respectively. υ(t)=[p(t)v(t)r(t)]T represents the surge, sway, and yaw velocity vector in the attached coordinate system (see Figure 1). J(ψ) is the coordinate rotation matrix. u(t) is the control force vector; M and D denote the inertial matrix and damping matrix, respectively. G is the mooring force matrix, and w(t)=[w1(t),w2(t),w3(t)] is the environmental disturbance including wind, waves, and currents.
Remark 1. During the modelling of the ship, it can be considered that the interaction forces between the various points of action of the ship are ignored, that is, the ship is considered as a rigid body. Hence, the combination of forward, lateral, and turning motions with heave, roll, and pitch motions is neglected, and we only consider the three degrees of freedom plane motion of the DPS.
Define
x(t)=[η(t)υ(t)]T=[xp(t)yp(t)ψ(t)p(t)q(t)r(t)]T, | (2) |
Let
−M−1D=[a11000a22a230a32a33],M−1=[d11000d22d230d32d33],−M−1G=[b11000b22b230b32b33]. | (3) |
Substitute (2) and (3) into (1) such that
˙x(t)=Ax(t)+Bu(t)+Bww(t), | (4) |
where
A=[000cos(ψ(t))−sin(ψ(t))0000sin(ψ(t))cos(ψ(t))0000001b1100a11000b22b230a22a230b32b330a32a33],B=Bw=[000000000d11000d22d230d32d33], |
where a11,a22,a23,a32,a33,d11,d22,d23,d32,d33,b11,b22,b23,b32,b33 are constant.
The global fuzzy models can be described as follows:
Mode Rule i: IF z1(t) is fi1, ⋯zn(t) is fin, THEN
˙x(t)=Aix(t)+Biu(t)+Bwiw(t),i=1,2,…,n | (5) |
where z1(t),z2(t),…,zn(t) represent the premise variables, fij is the fuzzy set, and n is the rules number.
Assume the yaw angle satisfies ψ(t)=[−π/π22,π/π22], then the controller rules are considered such that:
Model Rule 1:
IF ψ(t)=0, THEN
˙x(t)=A1x(t)+B1u(t)+Bw1w(t), | (6) |
Model Rule 2:
IF ψ(t)=π2, THEN
˙x(t)=A2x(t)+B2u(t)+Bw2w(t), | (7) |
Model Rule 3:
IF ψ(t)=−π2, THEN
˙x(t)=A3x(t)+B3u(t)+Bw3w(t), | (8) |
where
A1=[0001−α0000α10000001b1100a11000b22b230a22a230b32b330a32a33],A2=[000β−100001β0000001b1100a11000b22b230a22a230b32b330a32a33],A3=[000β10000−1β0000001b1100a11000b22b230a22a230b32b330a32a33],Bi=Bwi=[000000000d11000d22d230d32d33],i=1,2,3. |
where α=sin(2∘), β=cos(88∘). The global fuzzy models are inferred such that:
˙x(t)=3∑i=1μi(z(t))[Aix(t)+Biu(t)+Bwiw(t)], | (9) |
where
μi(z(t))=ωi(z(t))3∑i=1ωi(z(t))⩾0,i=1,2,3,ωi(z(t))=n∏j=1fij(zj(t)),3∑i=1μi(z(t))=1,z(t)=[z1(t),z2(t),…,zn(t)], |
The membership function of ψ(t) is shown in Figure 2.
Assuming that the states of the DPS are obtained at each sampling instant tk, the sampling time is assumed to be
0=t0<t1<⋯<tk<⋯<limk→∞tk=+∞. |
The structure of the system is depicted in Figure 3.
The system stability process is shown in Figure 4. It can be seen that the DPS fuzzy controller calculates the ship's control force and moment based on the deviation between the desired and actual position and velocity of DPS. Then, the thrust allocation distributes control force and moment for DPS to generate the actual position and velocity under the environmental disturbances, which achieves the stability target of the system.
Let d represent the upper bound of the sampling interval, i.e.,
0<tk+1−tk⩽d,∀k⩾0,d>0, |
then we consider the following sampled-data controller
u(t)=Kx(tk),tk⩽t<tk+1, | (10) |
where K represents the controller matrix. Then the T-S fuzzy controller can be described that:
Controller Rule 1:
IF ψ(t)=0, THEN
u(t)=K1x(tk). | (11) |
Controller Rule 2:
IF ψ(t)=π2, THEN
u(t)=K2x(tk). | (12) |
Controller Rule 3:
IF ψ(t)=−π2, THEN
u(t)=K3x(tk). | (13) |
Hence, the overall fuzzy controller is inferred to be:
u(t)=3∑j=1μj(z(t))Kjx(tk),tk⩽t<tk+1,k=0,1,2,… | (14) |
Substitute (14) into (11)
˙x(t)=3∑i=13∑j=1μi(z(t))μj(z(t))[Aix(t)+BiKjx(tk)+Bwiw(t)] | (15) |
When the system has an uncertainty parameter, the system is converted into
˙x(t)=3∑i=13∑j=1μi(z(t))μj(z(t))[(Ai+ΔA)x(t)+(Bi+ΔB)Kjx(tk)+Bwiw(t)] | (16) |
where ΔA and ΔB are the uncertainty parameter matrix, and satisfy
[ΔA(t)ΔB(t)]=DF(t)[EaEb] | (17) |
D,Ea,Eb are known dimensional constant matrices, and F(t) is unknown time-varying matrices with Lebesgue measurable elements that satisfy
FT(t)F(t)⩽I,∀t | (18) |
For the further deriving process, the two lemmas are given.
Lemma 1 [39] For any matrix H>0, and a vector with appropriate dimension μ(s), the inequality satisfies
−∫t2t1μT(s)Hμ(s)ds⩽(t2−t1)ξTtFTH−1Fξt+2ξTtHT∫t2t1μ(s)ds | (19) |
where, matrix F and vector ξt are not dependent on integral variables and have any dimensions.
Lemma 2 [40] For any matrix J,R=RT>0, constant d1⩽d(t)⩽d2, and vector function ˙x:[d2,d1]→Rn, the following integral inequality holds:
(d1−d2)∫t−d1t−d2˙xT(s)R˙x(s)ds⩽−ϑT(t)[I−I00I−I]TW[I−I00I−I]ϑ(t) | (20) |
where
ϑT(t)=[xT(t−d1)xT(t−d(t))xT(t−d2)],W=[−RJ∗−R]⩽0 |
Lemma 3 [41] Given matrices ΩT=Ω,Γ,X, for any F satisfying FTF⩽L, the inequality is established
Ω+HFE+ETFTHT<0 | (21) |
and only if there is a scalar ε>0, then
Ω+εHHT+ε−1ETLE<0 | (22) |
In the section, a novel LKF is constructed, and the stabilization condition of system (15) is proposed. To begin, the notation is defined as follows:
ξT(t)=[xT(t)˙xT(t)xT(tk)xT(t−d)∫ttkxT(tk)] |
Theorem 1: The system (15) with w(t)=0 is asymptotically stable, if there exist positive-definite matrices X=[XT1XT2XT3XT4XT5]T, P=[P1P2∗P3],Q,Z,U,S,R,N1,N2,M=[M11M12∗M22], and given scale d>0, such that
[ZU∗Z]>0 | (23) |
Ψ1+dΨ3<0 | (24) |
[Ψ1+dΨ2dX∗−dM−111]<0 | (25) |
where
Ψ1=[Ξ11Ξ12Ξ13Ξ14Ξ15∗Ξ22Ξ230P2∗∗Ξ33Ξ34Ξ35∗∗∗Ξ440∗∗∗∗0]Ψ2=−IT3M22I3−IT3N1I3Ψ3=IT23MI23+IT3N1I3Ξ11=G1Ai+ATiGT1+P2+PT2+Q−Z+X1Ξ12=−G1+SRT+ATiG2+P1+XT2Ξ13=G1BiKj−MT12+Z−U+XT3−X1Ξ14=U+XT4Ξ15=PT3+XT5Ξ22=−G2−GT2+d2ZΞ23=G2BiKj−X2Ξ33=M12−2Z+UT+U−X3Ξ34=Z−U−XT4 |
Ξ35=−XT5−PT3Ξ44=−Q−ZI23=[0I000000I000],I3=[00I000],L=diag{l1,l2,⋯ln} |
Proof. Consider the following LKF:
V(t)=5∑i=1Vi(t),t∈[tk,tk+1)V1(t)=[x(t)∫ttkx(s)ds]T[P1P2PT2P3][x(t)∫ttkx(s)ds]V2(t)=∫tt−dx(s)TQx(s)dsV3(t)=(d−(t−tk))∫ttk[˙x(s)x(tk)]TM[˙x(s)x(tk)]dsV4(t)=(tk+1−t)(t−tk)[xT(tk)N1x(tk)]V5(t)=d∫0−d∫tt+θ˙xT(s)Z˙x(s)dsdθ | (26) |
Take the derivative of V(t), which yields
˙V1(t)=2[x(t)∫ttkx(s)ds]T[P1+SRTP2PT2P3][˙x(t)x(t)−x(tk)]=2xT(t)P1˙x(t)+2(∫ttkx(s)ds)TPT2˙x(t)+2xT(t)P2x(t)+2(∫ttkx(s)ds)TP3x(t)+2xT(t)SRT˙x(t)−2(∫ttkx(s)ds)TP3x(tk) | (27) |
˙V2(t)=x(t)TQx(t)−x(t−d)TQx(t−d) | (28) |
˙V3(t)=−∫ttk[˙x(s)x(tk)]TM[˙x(s)x(tk)]ds+(d−(t−tk))[˙x(t)x(tk)]TM[˙x(t)x(tk)]=−∫ttk˙x(s)M11˙x(s)ds−2x(tk)M12[x(t)−x(tk)]−(t−tk)x(tk)M22x(tk)+(d−(t−tk))ς(t)MςT(t) | (29) |
˙V4(t)=(dk−(t−tk))[xT(tk)N1x(tk)]−(t−tk)[xT(tk)N1x(tk)]⩽(d−(t−tk))[xT(tk)N1x(tk)]−(t−tk)[xT(tk)N1x(tk)] | (30) |
˙V5(t)=d2˙xT(t)Z˙x(t)−d∫tt−d˙xT(t)Z˙x(t)ds | (31) |
By Lemma 1, it can be obtained that
−∫ttk˙x(s)M11˙x(s)ds⩽(t−tk)ξT(t)XM−111XTξ(t)+2ξT(t)X[x(t)−x(tk)] | (32) |
Then
˙V3(t)⩽(t−tk)ξT(t)VM−111VTξ(t)+2ξT(t)V[x(t)−x(tk)]−2x(tk)M12[x(t)−x(tk)]−(t−tk)x(tk)M22x(tk)+(d−(t−tk))ς(t)MςT(t) |
By Lemma 2, we have
−d∫tt−d˙xT(s)Z˙x(s)ds⩽ϑT(t)[−ZZ−UU∗−2Z+UT+UZ−U∗∗−Z]ϑ(t) | (33) |
where
ϑ(t)=[xT(t)xT(tk)xT(t−d)]T |
then
˙V5(t)⩽d2˙xT(t)Z˙x(t)−ϑT(t)[Z−Z+U−U∗2Z−UT+U−Z+U∗∗Z]ϑ(t) | (34) |
For any matrixes G1,G2, we have
2[xT(t)G1+˙xT(t)G2]×[−˙x(t)+r∑i=1r∑j=1μi(θ(t))μj(θ(t))[Aix(t)+BiKjx(tk)]]=0 | (35) |
For any diagonal matrix D>0, it can be obtained that
−2[Lx(tk)]TD[−Lx(tk)]⩾0 | (36) |
Then, from (27) to (36), we have
˙V(t)⩽3∑i=13∑j=1μi(θ(t))μj(θ(t))[ξT(t)Ψ(τ(t))ξ(t)] | (37) |
where
τ(t)=t−tkξT(t)=[xT(t)˙xT(t)xT(tk)xT(t−d)∫ttkxT(tk)]Ψ(τ(t))=Ψ(t−tk)=Ψ1+(t−tk)(Ψ2+VTˉM−111V)+(d−(t−tk))Ψ3 |
Ψ1,Ψ2,Ψ3 is defined in Theorem 1.
By Lemma 2, we can derive a linear convex combination about τ(t), for all τ(t)∈[0,d], if
Ψ(0)=Ψ1+dΨ3<0 | (38) |
Ψ(d)=Ψ1+d(Ψ2+VTˉM−11V)<0 | (39) |
Then, Ψ(τ(t)) is negative definite.
Following the Schur complement, (24) and (25) imply (38) and (39), which means that ˙V(t)<0. This proof is completed.
Remark 1. In the constructed LKF, V1(t) and V2(t) are the general energy function of the system state, and it represents the energy changes of the system state and derivative from t−d to t. It is the foundation of the entire LKF. V5(t) symbolizes the relationship between state derivatives and time delay, and it is the foundation of time delay dependent stability standards. Additionally, V3(t) and V4(t) include the characteristic of time delay, and they play an important role in obtaining the final solution.
Remark 2. Compared with literature [42], two (tk,tk+1) dependent terms like (tk+1−t)t∫tke2αs[x(s)x(tk)]T[R1R2∗R3][x(s)x(tk)]ds and (tk+1−t)(t−tk)[xT(tk)N1x(tk)] are added in the LKF, which means that the available characteristic of sampling patterns are fully captured, and the LKF has a more general form.
Remark 3. In some previously published papers (such as [43,44]), they used overly relaxed delimitation techniques when estimating some integral terms such as −d∫tkt−d˙xT(s)Z˙x(s)ds. However, this paper directly uses a tighter technique by reciprocally convex combination in Lemma 2. Besides, free matrices G and G1 are proposed during the derivation procedure of LKF, which improves the design flexibility. Therefore, the conservatism of stability criteria can be greatly reduced.
To prove the robustness of the system, the following H∞ performance function is proposed
Jyw=∫t0[yT(s)y(s)−γ2wT(s)w(s)]ds,γ>0 | (40) |
where w(t)≠0 and w(t)∈L2[0,∞), and y(t) denotes the system output, and satisfies that
y(t)=3∑i=13∑j=1μi(z(t))μj(z(t))Cix(t), | (41) |
Then, based on the Theorem 2, the system (15) can be proved to satisfy the robust performance with disturbances.
Theorem 2: the system (15) is asymptotically stable, if there exist positive-definite matrices X=[XT1XT2XT3XT4XT5]T, P=[P1P2∗P3],Q,Z,U,S,R,N1,N2,M=[M11M12∗M22], and given scale d>0, such that
[ZU∗Z]>0 | (42) |
Ψ1+dΨ3<0 | (43) |
[Ψ1+dΨ2dX∗−dM−111]<0 | (44) |
where
˜Ψ1=[˜Ξ11Ξ12Ξ13Ξ14Ξ150∗Ξ22Ξ230P20∗∗Ξ33Ξ34Ξ350∗∗∗Ξ4400∗∗∗∗00∗∗∗∗∗−γ2I]˜Ξ11=G1Ai+ATiGT1+P2+PT2+Q−Z+X1+CTC |
Ξ12,Ξ13,Ξ14,Ξ15,Ξ22,Ξ23,Ξ33,Ξ34,Ξ35,Ξ44 are defined in Theorem 1.
Prove:
Jyw=∫t0[yT(s)y(s)−γ2wT(s)w(s)]ds,=∫t0[yT(s)y(s)−γ2wT(s)w(s)+˙V(s)]ds−V(t)⩽∫t0[yT(s)y(s)−γ2wT(s)w(s)+˙V(s)]ds⩽3∑i=13∑j=1μi(θ(t))μj(θ(t))[ζT(t)Θ(τ(t))ζ(t)] | (45) |
where
τ(t)=t−tkζT(t)=[xT(t)˙xT(t)xT(tk)xT(t−d)∫ttkxT(tk)wT(t)]Ψ(τ(t))=Ψ(t−tk)=˜Ψ1+(t−tk)(Ψ2+VTˉM−111V)+(d−(t−tk))Ψ3 |
The proof is finished.
Next, Theorem 3 provides the stability conditions for systems with norm bounded uncertainty.
Theorem 3: the system (15) is asymptotically stable, if there exist positive-definite matrices X=[XT1XT2XT3XT4XT5]T, P=[P1P2∗P3],Q,Z,U,S,R,N1,N2,M=[M11M12∗M22], and given scale d>0,ε>0, such that
[ZU∗Z]>0 | (46) |
Ψ1+dΨ3<0 | (47) |
[Ψ1+dΨ2dX∗−dM−111]<0 | (48) |
where
ˆΨ1=[Ξ11+εETaEaΞ12Ξ13+εETaEbΞ14Ξ15G1D∗Ξ22Ξ230P2G2D∗∗Ξ33+εETbEbΞ34Ξ350∗∗∗Ξ4400∗∗∗∗00∗∗∗∗∗−εI] |
Ξ12,Ξ13,Ξ14,Ξ15,Ξ22,Ξ23,Ξ33,Ξ34,Ξ35,Ξ44 are defined in Theorem 1.
Prove: replace A and Ad in (25) with A+DF(t)Ea and B+DF(t)Eb, respectively. It can be obtained that:
Ψ+ΓF(t)Ω+ΩTFT(t)ΓT<0 | (49) |
where
Γ=[G1DG2D000],Ω=[Ea0Eb00] |
According to Lemma 3, if there exist a scale ε>0, the following inequality holds
Ψ+εΩTΩ+ε−1ΓΓT<0 | (50) |
According to Schu's complement, Eq (48) is equivalent to Eq (50). The proof is completed.
Then, the controller method can be obtained according to Theorem 4.
Theorem 4: for a given scale d>0, if there exist positive-definite matrices X=[XT1XT2XT3XT4XT5]T, P=[P1P2∗P3],Q,Z,U,S,R,N1,N2,M=[M11M12∗M22], such that
[ZU∗Z]>0 | (51) |
ˉΨ1+dˉΨ3<0 | (52) |
[ˉΨ1+dˉΨ2dˉX∗−dˉM−111]<0 | (53) |
where
ˉΨ1=[ˉΞ11ˉΞ12ˉΞ13ˉΞ14ˉΞ15∗ˉΞ22ˉΞ230ˉP2∗∗ˉΞ33ˉΞ34ˉΞ35∗∗∗ˉΞ440∗∗∗∗0]ˉΨ2=−IT3ˉM22I3−IT3ˉN1I3ˉΨ3=IT23ˉMI23+IT3ˉN1I3ˉΞ11=AiG+GTATi+ˉP2+ˉPT2+ˉQ−Z+ˉX1ˉΞ12=−G+εGATi+SˉRT+ˉP1+ˉXT2ˉΞ13=BiˉKj−ˉMT12+ˉZ−ˉU+ˉXT3−ˉX1ˉΞ14=ˉUE+ˉXT4ˉΞ15=ˉPT3+ˉXT5ˉΞ22=−εG−εGT+d2ˉZˉΞ23=εBiˉKj−ˉX2ˉΞ33=ˉM12−2ˉZ+ˉUT+ˉU−ˉX3ˉΞ34=Z−U−ˉXT4ˉΞ35=−ˉXT5−ˉPT3ˉΞ44=−ˉQ−ˉZI23=[0I000000I000]I3=[00I000]L=diag{l1,l2,⋯ln}ξT(t)=[xT(t)˙xT(t)xT(tk)xT(t−d)∫ttkxT(tk)] |
then the system (15) with w(t)=0 is asymptotically stable. Moreover, the controller gain is
Kj=ˉKjG−1 | (54) |
Proof: Define
G1=G−1,G2=εG−1,ˉK=KG,ˉP=diag{ G,G}Pdiag{ G,G},ˉM=diag{ G,G,G}Mdiag{ G,G,G},ˉD=GDG,ˉQ=GQG,ˉN1=GN1G,ˉN2=GN2G,ˉR=GRG,ˉU=GUG,ˉZ=GZG. |
Let diag{ G,G,G,G,G,G}T and diag{ G,G,G,G,G,G}. By pre- and post-multiplying (24) respectively, (52) can be obtained. Let diag{ G,G,G,G,G,G,G}T and diag{ G,G,G,G,G,G,G}, pre- and post-multiplying (25) respectively, (53) can be obtained. This completes the proof.
Remark 4. The system (15) can obtain less conservativeness than that in [45], which is mainly due to two results. First, during the construction of V(t), a second-order term (tk+1−t)(t−tk)[xT(tk)N1x(tk)] related to time t is also added, which can reduce the conservatism of the results to some extent. Second, during the derivative of the LKF, the reciprocally convex combination method is used, which can obtain a more stringent derivative upper limit than Jensen's inequality, thereby less conservativeness will be achieved.
To demonstrate the merits of the method, a practical example is carried out. The M and D is considered as follows [46].
M=[1.08520002.0575−0.40870−0.40870.2153],D=[0.08650000.07620.151000.01510.0031],G=[0.03890000.02660000]. |
Let α=sin20 and β=cos880, then
A1=[0001.0000−0.034900000.03491.00000000001.0000−0.035800−0.0797000−0.020800−0.0818−0.12240−0.039400−0.2254−0.2468], |
A2=[0000.0349−1.000000001.00000.03490000001.0000−0.035800−0.0797000−0.020800−0.0818−0.12240−0.039400−0.2254−0.2468],A3=[0000.03491.00000000−1.00000.03490000001.0000−0.035800−0.0797000−0.020800−0.0818−0.12240−0.039400−0.2254−0.2468], |
Bi=Bwi=[0000000000.92150000.78021.481101.48117.4562],i=1,2,3. |
The allowable maximum sampling interval is an important control performance index of the sampled-data system. A longer sampling period indicates a lower communication burden and fewer time consumption of the actuator, which helps to save energy, time and so on. Table 1 lists the results which uses different methods to obtain the maximum sampling internal. From Table 1, it can be seen that the sampling internal in Theorem1 improves [36,37,38] about 189.2%, 173.86%, and 35.9%, respectively, which means that the designed fuzzy controller is beneficial to obtain a larger maximum sampling interval.
Assume the initial value of the system that xs(t)=[15m15m0.2∘0m/s0m/s0∘/s].
Then, the gain can be computed such that
K1=[−0.0281−0.0444−0.0339−0.7195−0.0002−0.00030.0708−0.00850.04570.0004−1.38770.2656−0.01420.00740.0088−0.00010.2795−0.1371]K2=[−0.01230.01660.0170−0.7473−0.00020.0000−0.02280.02810.07900.0002−1.44020.27450.00470.00020.0029−0.00010.2900−0.1421]K3=[−0.01230.01660.0170−0.7473−0.00020.0000−0.02280.02810.07900.0002−1.44020.27450.00470.00020.0029−0.00010.2900−0.1421] |
To demonstrate the control performance of the proposed method, Figures 5–10 present the comparison results with LKF in [36]. From Figures 5–10, it illustrates that the system state xs(t) reaches the expected values within approximately 10 seconds, while [36] achieves the desired target and requires a longer time. This indicates that the response speed in the proposed method is faster. Besides, the system's oscillation amplitude is smaller, which indicates that our result has good anti-interference ability under environment disturbance. In addition, Figures 5–10 further illustrate that the proposed method has more preferable control performance.
The stabilization issue for fuzzy DPS with sampled-data is studied. To reduce the conservative and improve the sampling interval, a novel LKF is constructed by adding some integral terms, which can utilize the information of actual sampling modes fully. Then a fuzzy sampling control scheme is proposed. By introducing free matrices and reciprocally convex approach, the conservatism of the system is reduced. Simulation research demonstrates that the proposed method can achieve better control performance. However, the proposed methodologies didn't consider the lower bound of the sampling interval, and some matrix in the LKF are still required to be positive definite. Additionally, how to avoid the emergence of hierarchical phenomena in the system is an urgent problem to be solved. Furthermore, we will improve the LKF and discuss the aperiodic sampled-data control problem for DPS, which both considers the lower and upper bound of the sampling interval.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
This study was funded by National Natural Science Foundation of China (52371369); the Key Projects of National Key R & D Program (No.2021YFB390150), the National Science Project of Fujian Province (No. 2020J01659, 2022J01323, 2021J01822, 2020J01660, 20230019), the Fuzhou-Xiamen-Quanzhou In-dependent Innovation Region Cooperated Special Foundation (No: 3502ZCQXT2021007); Funds of Fujian Provincial for Promoting High Quality Development of Marine and Fisheries Industry (No. FJHYF-ZH-2023-10).
The authors declare there is no conflict of interest.
[1] |
Adámek V (2018) The limits of Timoshenko beam theory applied to impact problems of layered beams. Int J Mech Sci 145: 128–137. https://doi.org/10.1016/j.ijmecsci.2018.07.001. doi: 10.1016/j.ijmecsci.2018.07.001
![]() |
[2] |
Eisenträger J, Naumenko K, Altenbach H, et al. (2015) Application of the first-order shear deformation theory to the analysis of laminated glasses and photovoltaic panels. Int J Mech Sci 96: 163–171. https://doi.org/10.1016/j.ijmecsci.2018.07.001. doi: 10.1016/j.ijmecsci.2018.07.001
![]() |
[3] |
Koutsawa Y, Daya EM (2007) Static and free vibration analysis of laminated glass beam on viscoelastic supports. Int J Solids Struct 44: 8735–8750. https://doi.org/10.1016/j.ijsolstr.2007.07.009. doi: 10.1016/j.ijsolstr.2007.07.009
![]() |
[4] |
Schulze S, Pander M, Naumenko K, et al. (2012) Analysis of laminated glass beams for photovoltaic applications. Int J Solids Struct 49(15–16): 2027–2036. https://doi.org/10.1016/j.ijsolstr.2012.03.028. doi: 10.1016/j.ijsolstr.2012.03.028
![]() |
[5] |
Naumenko K, Eremeyev VA (2014) A layer-wise theory for laminated glass and photovoltaic panels. Compos Struct 112: 283–291. https://doi.org/10.1016/j.compstruct.2014.02.009. doi: 10.1016/j.compstruct.2014.02.009
![]() |
[6] |
Naumenko K, Eremeyev VA (2017) A layer-wise theory of shallow shells with thin soft core for laminated glass and photovoltaic applications. Compos Struct 178: 434–446. https://doi.org/10.1088/1475-7516/2017/07/007. doi: 10.1088/1475-7516/2017/07/007
![]() |
[7] |
Aßmus M, Naumenko K, Altenbach H (2016) A multiscale projection approach for the coupled global–local structural analysis of photovoltaic modules. Compos Struct 158: 340–358. https://doi.org/10.1088/1475-7516/2016/09/036. doi: 10.1088/1475-7516/2016/09/036
![]() |
[8] |
Eisenträger J, Naumenko K, Altenbach H, et al. (2015) A user-defined finite element for laminated glass panels and photovoltaic modules based on a layer-wise theory. Compos Struct 133: 265–277. https://doi.org/10.1016/j.compstruct.2015.07.049. doi: 10.1016/j.compstruct.2015.07.049
![]() |
[9] |
Li G, Cinefra M, Carrera E (2020) Coupled thermo-mechanical finite element models with node-dependent kinematics for multi-layered shell structures. Int J Mech Sci 171: 105379. https://doi.org/10.1016/j.ijmecsci.2019.105379. doi: 10.1016/j.ijmecsci.2019.105379
![]() |
[10] |
Janda T, Schmidt J, Hála P, et al. (2021) Reduced order models of elastic glass plate under low velocity impact. Comput Struct 244: 106430. https://doi.org/10.1016/j.compstruc.2020.106430. doi: 10.1016/j.compstruc.2020.106430
![]() |
[11] |
Pelfrene J, Kuntsche J, Van Dam S, et al. (2016) Critical assessment of the post-breakage performance of blast loaded laminated glazing: experiments and simulations. Int J Impact Eng 88: 61–71. https://doi.org/10.1016/j.ijimpeng.2015.09.008. doi: 10.1016/j.ijimpeng.2015.09.008
![]() |
[12] |
Vedrtnam A, Pawar S (2017) Laminated plate theories and fracture of laminated glass plate–a review. Eng Fract Mech 186: 316–330. https://doi.org/10.1016/j.engfracmech.2017.10.020. doi: 10.1016/j.engfracmech.2017.10.020
![]() |
[13] | Murakami S (2012) Continuum Damage Mechanics: A Continuum Mechanics Approach to the Analysis of Damage and Fracture, Berlin: Springer. |
[14] |
Sun X, Khaleel MA (2004) Modeling of glass fracture damage using continuum damage mechanics–static spherical indentation. Int J Damage Mech 13(3): 263–285. https://doi.org/10.1177/1056789504042593. doi: 10.1177/1056789504042593
![]() |
[15] |
Wei J, Dharani L (2005) Fracture mechanics of laminated glass subjected to blast loading. Theor Appl Fract Mech 44(2): 157–167. https://doi.org/10.1016/j.tafmec.2005.06.004. doi: 10.1016/j.tafmec.2005.06.004
![]() |
[16] |
Wang Z, Fu J, Manes A (2021) Discrete fracture and size effect of aluminosilicate glass under flexural loading: Monte carlo simulations and experimental validation. Theor Appl Fract Mech 111: 102864. https://doi.org/10.1016/j.tafmec.2020.102864. doi: 10.1016/j.tafmec.2020.102864
![]() |
[17] |
Gao W, Liu X, Chen S, et al. (2020) A cohesive zone based de/fe coupling approach for interfacial debonding analysis of laminated glass. Theor Appl Fract Mech 108: 102668. https://doi.org/10.1016/j.tafmec.2020.102668. doi: 10.1016/j.tafmec.2020.102668
![]() |
[18] |
Altenbach H, Naumenko K (1997) Creep bending of thin-walled shells and plates by consideration of finite deflections. Comput Mech 19(6): 490–495. https://doi.org/10.1007/s004660050197. doi: 10.1007/s004660050197
![]() |
[19] |
Altenbach H, Kolarow G, Morachkovsky O, et al. (2000) On the accuracy of creep-damage predictions in thinwalled structures using the finite element method. Comput Mech 25(1): 87–98. https://doi.org/10.1007/s004660050018. doi: 10.1007/s004660050018
![]() |
[20] |
Peerlings RH, de Borst R, Brekelmans WM, et al. (1996) Gradient enhanced damage for quasi-brittle materials. Int J Numer Methods Eng 39: 3391–3403. https://doi.org/10.1007/978-94-011-5520-5-20. doi: 10.1007/978-94-011-5520-5-20
![]() |
[21] |
de Borst R, Verhoosel CV (2016) Gradient damage vs phase-field approaches for fracture: Similarities and differences. Comput Methods Appl Mech Eng 312: 78–94. https://doi.org/10.1016/j.cma.2016.05.015. doi: 10.1016/j.cma.2016.05.015
![]() |
[22] |
Linse T, Hennig P, Kästner M, et al. (2017) A convergence study of phase-field models for brittle fracture. Eng Fract Mech 184: 307–318. https://doi.org/10.1016/j.engfracmech.2017.09.013. doi: 10.1016/j.engfracmech.2017.09.013
![]() |
[23] |
Hansen-Dörr AC, de Borst R, Hennig P, et al. (2019) Phase-field modelling of interface failure in brittle materials. Comput Methods Appl Mech Eng 346: 25–42. https://doi.org/10.1016/j.cma.2018.11.020. doi: 10.1016/j.cma.2018.11.020
![]() |
[24] |
Rodriguez P, Ulloa J, Samaniego C, et al. (2018) A variational approach to the phase field modeling of brittle and ductile fracture. Int J Mech Sci 144: 502–517. https://doi.org/10.1016/j.ijmecsci.2018.05.009. doi: 10.1016/j.ijmecsci.2018.05.009
![]() |
[25] |
Schmidt J, Zemanová A, Zeman J, et al. (2020) Phase-field fracture modelling of thin monolithic and laminated glass plates under quasi-static bending. Materials 13(22): 5153. https://doi.org/10.3390/ma13225153. doi: 10.3390/ma13225153
![]() |
[26] |
Nase M, Rennert M, Naumenko K, et al. (2016) {Identifying traction–separation behavior of self-adhesive polymeric films from in situ digital images under T-peeling}. J Mech Phys Solids 91: 40–55. https://doi.org/10.1016/j.jmps.2016.03.001. doi: 10.1016/j.jmps.2016.03.001
![]() |
[27] |
Naumenko K, Bagheri B (2021) A direct approach to evaluate interaction forces between self-adhesive polymeric films subjected to T-peeling. Arch Appl Mech 91(2): 629–641. https://doi.org/10.1007/s00419-020-01834-9. doi: 10.1007/s00419-020-01834-9
![]() |
[28] |
Silling SA, Lehoucq RB (2010) Peridynamic theory of solid mechanics. Advs Appl Mech 44: 73–168. https://doi.org/10.1016/S0065-2156(10)44002-8. doi: 10.1016/S0065-2156(10)44002-8
![]() |
[29] |
Naumenko K, Pander M, Würkner M (2022) Damage patterns in float glass plates: Experiments and peridynamics analysis. Theor Appl Fract Mech 118: 103264. https://doi.org/10.1016/j.tafmec.2022.103264. doi: 10.1016/j.tafmec.2022.103264
![]() |
[30] |
Niazi S, Chen Z, Bobaru F (2021) Crack nucleation in brittle and quasi-brittle materials: A peridynamic analysis. Theor Appl Fract Mech 112: 102855. https://doi.org/10.1016/j.tafmec.2020.102855. doi: 10.1016/j.tafmec.2020.102855
![]() |
[31] |
Mehrmashhadi J, Bahadori M, Bobaru F (2020) On validating peridynamic models and a phase-field model for dynamic brittle fracture in glass. Eng Fract Mech 240: 107355. https://doi.org/10.1016/j.engfracmech.2020.107355. doi: 10.1016/j.engfracmech.2020.107355
![]() |
[32] |
Diana V, Ballarini R (2020) Crack kinking in isotropic and orthotropic micropolar peridynamic solids. Int J Solids Struct 196: 76–98. https://doi.org/10.1016/j.ijsolstr.2020.03.025. doi: 10.1016/j.ijsolstr.2020.03.025
![]() |
[33] |
Rahimi MN, Kefal A, Yildiz M, et al. (2020) An ordinary state-based peridynamic model for toughness enhancement of brittle materials through drilling stop-holes. Int J Mech Sci 182: 105773. https://doi.org/10.1016/j.ijmecsci.2020.105773. doi: 10.1016/j.ijmecsci.2020.105773
![]() |
[34] |
Zhang Y, Deng H, Deng J, et al. (2020) Peridynamic simulation of crack propagation of non-homogeneous brittle rock-like materials. Theor Appl Fract Mech 106: 102438. https://doi.org/10.1016/j.tafmec.2019.102438. doi: 10.1016/j.tafmec.2019.102438
![]() |
[35] |
Naumenko K, Eremeyev VA (2022) A non-linear direct peridynamics plate theory. Compos Struct 279: 114728. https://doi.org/10.1016/j.compstruct.2021.114728. doi: 10.1016/j.compstruct.2021.114728
![]() |
[36] |
Nguyen CT, Oterkus S (2021) Ordinary state-based peridynamics for geometrically nonlinear analysis of plates. Theor Appl Fract Mech 112: 102877. https://doi.org/10.1016/j.tafmec.2020.102877. doi: 10.1016/j.tafmec.2020.102877
![]() |
[37] |
Yang Z, Naumenko K, Altenbach H, et al. (2022) Some analytical solutions to peridynamic beam equations. ZAMM 2022: e202200132. https://doi.org/10.1002/zamm.202200132. doi: 10.1002/zamm.202200132
![]() |
[38] | Belytschko T, Liu WK, Moran B, et al. (2014) Nonlinear Finite Elements for Continua and Structures, New York: Wiley. |
[39] | Naumenko K, Altenbach H (2016) Modeling High Temperature Materials Behavior for Structural Analysis: Part I: Continuum Mechanics Foundations and Constitutive Models, Berlin: Springer. |
[40] |
Silling SA (2000) Reformulation of elasticity theory for discontinuities and long-range forces. J Mech Phys Solids 48(1): 175–209. https://10.1016/S0022-5096(99)00029-0. doi: 10.1016/S0022-5096(99)00029-0
![]() |
[41] | Silling SA (2016) Introduction to peridynamics, Handbook of Peridynamic Modeling, London: CRC Press, 63–98. |
[42] |
Javili A, Morasata R, Oterkus E, et al. (2019) Peridynamics review. Math Mech Solids 24(11): 3714–3739. https://doi.org/10.1177/1081286518803411. doi: 10.1177/1081286518803411
![]() |
[43] |
Silling SA, Epton M, Weckner O, et al. (2007) Peridynamic states and constitutive modeling. J Elast 88(2): 151–184. https://doi.org/10.1007/s10659-007-9125-1. doi: 10.1007/s10659-007-9125-1
![]() |
[44] |
Silling SA, Askari E (2005) A meshfree method based on the peridynamic model of solid mechanics. Comput Struct 83(17–18): 1526–1535. https://doi.org/10.1016/j.compstruc.2004.11.026. doi: 10.1016/j.compstruc.2004.11.026
![]() |
[45] | Littlewood DJ, Parks ML, Mitchell JA, et al. (2013) The peridigm framework for peridynamic simulations, 12th U.S. National Congress on Computational Mechanics, United States: Sandia National Lab. |
[46] | Littlewood D (2016) Roadmap for software implementation, Handbook of peridynamic modeling, London: CRC Press, 147–178. |