key: cord-275124-7l53bvp1 authors: Yao, Minghui; Wang, Hongbo title: A potential treatment for COVID-19 based on modal characteristics and dynamic responses analysis of 2019-nCoV date: 2020-10-21 journal: Nonlinear Dyn DOI: 10.1007/s11071-020-06019-1 sha: doc_id: 275124 cord_uid: 7l53bvp1 The 2019-nCoV is ravaging the world, taking lots of lives, and it is emergent to find a solution to deal with this novel pneumonia. This paper provides a potential treatment for COVID-19 utilizing resonance to destroy the infection ability of 2019-nCoV. Firstly, the geometry size of 2019-nCoV is scaled up by 10,000 times. The additional mass is used to represent the effect of the fluid around a spike protein. The finite element analysis (FEA) is used to study the modal characteristics of the tuned 2019-nCoV model and mistuned 2019-nCoV model in blood, respectively. Based on FEA, the lumped parameter mechanical model of 2019-nCoV is established. Then, the dynamic responses of mistuned 2019-nCoV are investigated through harmonic response and dynamical analysis. Finally, a potential method utilizing 360° sweep excitation to cure COVID-19 is put forward. A novel coronavirus (2019-nCoV) was identified as the causative virus linked with a cluster of cases of pneumonia detected in Wuhan city by Chinese authorities on 7 January [1] . Due to the discovery of 2019-nCoV, the virus has been diagnosed quickly [2, 3] . From December 2019, coronavirus disease 2019 (COVID-19) has been ravaging the world and killed more than 181 thousand people [4] . Medical scientists around the world are exerting themselves to develop vaccines, but vaccine researches are timeconsuming, and vaccines can lose efficiency when ribonucleic acid (RNA) for this virus mutates. This paper tries to provide a potential method based on vibration to damage the infectivity of 2019-nCoV and cure COVID-19, and it will not fail if the virus mutates. 2019-nCoV is a kind of coronavirus, which consists of virus sphere and spike proteins on the sphere surface. When 2019-nCoV infect humans, they enter cells of humans by binding their spike proteins to receptors on the cell membrane, and then copy viruses in the cell. It is theoretically feasible to disable the 2019-nCoV by destroying the spike protein structures through resonance, and then it is easy for the immune system to eliminate these viruses and recover. This paper combines finite element analysis (FEA) and theoretical methods to study the modal characteristics and dynamic responses of 2019-nCoV in the blood, The geometry size of 2019-nCoV is about 200 nm, and for the convenience of FEA, the geometry size is scaled up by 10,000 times to 2 mm. All other parameters of length dimension are magnified 10,000 times, so that we can get the accurate result by scaling down length dimension of results by 10,000 times. Generally, it is necessary to utilize the nonlocal elasticity theory to consider the size effect of nanomaterials. However, based on studies given by Peddieson et al. [5] , for a cantilever under concentrated force, its bending behavior utilizing the nonlocal elasticity theory is identical to the local result. In this paper, spike proteins structure on the novel coronavirus sphere surface is a cantilever structure. Therefore, it is reasonable to apply traditional mechanical methods to conduct the analysis. The 3D 2019-nCoV model is built in Solidworks software as shown in Fig. 1 with protein material, which is input to ANSYS Workbench to do modal analysis. In the process of analyzing the modal characteristic of structures in the fluid, the additional mass is used to represent the effect of the fluid around the structure. According to fluid mechanics by Mao and Wu [6, 7] , the additional mass of the blood per unit length around a spike protein can be calculated as follows where q b is the density of the blood, R is the equivalent radius of the environment, and r is the radius of the spike protein. Since R is infinite compared to r, m bp is simplified as m bp ¼ q b pr 2 . The parameters are set as The additional mass of the blood around a spike protein is m b ¼ m bp l 1 ¼ q b pr 2 l 1 ¼ 2:6 Â 10 À22 kg in the real situation, while in the modal analysis it is added to the spike protein by an equivalent point mass 2:6 Â 10 À10 kg. Because of the relatively heavy mass of the sphere body compared to spike proteins, its modes are rigid body motion, which the sphere of COVID-19 model is fixed in modal analysis. A tuned 2019-nCoV model indicates that every spike protein is identical. The frequencies of the tuned 2019-nCoV model are shown in Table 1 . The first 200 frequencies are very close to each other, and they are all first-order bending vibrations. Figure 2 is one of the corresponding mode shapes. These frequencies are supposed to be identical according to the vibration theory, and the discrepancies are caused by the intrinsic error of finite element analysis (FEA). In the realistic situation under the electron microscope, the number of spike proteins varies from different novel coronavirus. Even though the number of spike proteins in the same kind of novel coronavirus varies, it is hard to detect the exact number currently. In this paper, the 2019-nCoV model is built based on estimation. In order to achieve mistuning, the circle of spike proteins picked to lengthen could be chosen randomly in principle. But for the sake of analysis, the circle of 30 spike proteins in our model is selected. A mistuned 2019-nCoV model is obtained by lengthening one chosen circle, which has 30 spike proteins. At the same time, keeps the other circles unchanged. The frequencies of the mistuned 2019-nCoV model are shown in Table 2 . The first 60 frequencies are very close to each other, and they are all first-order bending vibration of the lengthened spike proteins. From 61st frequency, they are the first-order bending vibration of unchanged spike proteins. Figure 3a , b is one of the corresponding mode shapes, respectively. The discrepancies of the first 60 frequencies are caused by the intrinsic error of FEA. The spike proteins in realistic 2019-nCoV are different from each other in geometry, so the frequencies of the realistic 2019-nCoV model are a range of frequencies close to each other, and every frequency corresponds to the first-order bending vibration of one spike protein. 3 Lumped parameter mechanical model of 2019-nCoV According to the mode shape of single spike protein as shown in Fig. 4 , the first-order bending mode shape is employed as the vibration form of spike proteins. Based on studies given by Weaver and Timoshenko et al. [8] , the Rayleigh method is utilized to calculate the first-order frequency to check this choice where p is circular frequency, q is the density, E is the Young's Modulus, I is the bending moment, Y is the first-order mode function, A is the area of cross section, and l is the length. Compared with the FEA result, which is 198,340,000 Hz, the relative error is 0.22%. This indicates a reasonable choice. Based on the first-order bending vibration form, the lumped mechanical model is established as illustrated in Fig. 5 . According to the equivalent frequency principle, spike proteins are simplified as mass-spring structures. Each spike protein is identical in geometry except for its location. Based on studies given by Beer et al. [9] , the equivalent bending stiffness is calculated as k eq ¼ 3EI l 3 ¼ 1:41 Â 10 À3 m/N, and the equivalent mass m eq is calculated as m eq ¼ 9:14 Â 10 À22 kg by the following Eq. (3) where T represents the kinetic energy, A 2 ðxÞ ¼ p Á ð5 Â 10 À9 Á x þ 1 Â 10 À9 Þ 2 is the area of cross section of the end part in a spike protein, v is the velocity at the indicates the location of the additional mass of the blood around a spike protein. According to studies given by Weaver et al. [8] , the forced vibration theory of multiple degrees of freedom system is used to yield the natural frequency, which the natural frequency is calculated as f ¼ 1 2p ffiffiffiffiffi k eq m eq q ¼ 1:977Â 10 8 Hz. Comparing with Eq. (2), the relative error is 0.33%, which validates the accuracy of the model. The lumped parameter mechanical model of mistuned 2019-nCoV is roughly the same with tuned 2019-nCoV. But, there are small differences on equivalent stiffness and mass. 4 Dynamic response analysis of mistuned 2019-nCoV The harmonic response is conducted after the modal analysis of the mistuned 2019-nCoV model, and acceleration excitation is set as 1 mm/s 2 , as shown in Fig. 6 . The result of the harmonic response is shown in Fig. 7 , where one of the lengthened spike proteins has the largest amplitude. The FEA amplitude-frequency curves of a lengthened spike protein and an unchanged spike protein, whose h is p=8, are illustrated in Figs. 8 and 9 , respectively. These results accord with the modal characteristics. The excitation F sin xt is imposed on the end of spike proteins of the mistuned 2019-nCoV model. There are totally 32 circles, and each circle consists of spike proteins, which have the same angle h between their displacements and excitation. The viscosity of blood is considered as viscous damping. According to studies given by Weaver et al. [8] , the vibration theory of multitude degrees of freedom model and the forced vibration equation of lumped parameter model are utilized to obtain the following equation and the steady-state response as follows Fig. 6 Acceleration excitation in the process of harmonic response analysis Fig. 7 Harmonic response analysis result where h i ¼ ðiÀ1Þp 16 , u i ¼ tan À1 r eqi x=k eq ð1 À k 2 i Þ, h i is the angle between excitation and displacement, . It is easy to find that the vibration amplitude increases with cosh i j j. The spike proteins, whose displacements parallel with the excitation direction, have the biggest vibration amplitude. When 5p=16 h i 11p=16 and 21p=16 h i 27p=16, cos h i j j 0:56, and these vibration amplitudes are relatively small, which it can somehow achieve response concentration. According to studies given by Weaver et al. [8] , the relationship between amplitude and frequency of the lumped parameter model is obtained as follows Fig. 8 The FEA amplitude-frequency curve of a lengthened spike protein Fig. 9 The FEA amplitude-frequency curve of an unchanged spike protein The amplitude-frequency curves of the spike proteins are analyzed in harmonic response, as shown in Fig. 10 , where B 1max ¼ 3:56 Â 10 À12 m, B 2max ¼ 7:01 Â 10 À13 m. Compared with the FEA results in Figs. 8 and 9, relative errors are 5.10% and 7.83%, respectively. The spike proteins of a realistic 2019-nCoV are usually different from each other in geometry, and every spike protein has a different frequency, which there is small discrepancies. Therefore, it would be possible to utilize an ultrasonic vibration exciter to 360°sweep, whose frequency range is set as the firstorder bending vibration frequencies 1.9 9 10 8 Hz to 2.0 9 10 8 Hz in order to excite every spike protein of 2019-nCoV to resonate. 360°rotating sweep excitation could evoke resonance of spike proteins in turn and vibration concentration. Given a proper excitation amplitude, it could destroy the spike protein structures and disable the infectivity of 2019-nCoV. Effects of ultrasonic resonance destruction of viral proteins are discussed further. Since the bending strength of the spike protein is unknown, for the reliability of the analysis, it is set as the same as the bending strength of common steels, which is 300 MPa. If the stress in the spike protein is bigger than 300 MPa, it would crack at the base and lose its ability to invade human cells. Then, the required destruction excitation amplitude is analyzed by FEA, as shown in Fig. 11 , which is calculated as 1.041 9 10 -5 mm and the equivalent acceleration is 4.164 m/s 2 . Based on studies given by Glaister [10] , this acceleration is safe for human body. Based on studies given by Randall et al. [11] , it is known that the range of resonant frequencies of standing human bodies was found to be from 9 to 16 Hz. Since the resonant frequencies of human bodies are below 20 Hz, the vibration energy would hardly be absorbed by human body and it would not cause much harm and discomfort to patients. Therefore, it is efficient and safe to utilize an ultrasonic vibration exciter to conduct 360 o rotating sweep excitation to human body if the vibration frequencies are set as 1.9 9 10 8 Hz to 2.0 9 10 8 Hz and the amplitude is set as larger than 1.041 9 10 -5 mm. In this paper, the modal characteristics and dynamic responses of 2019-nCoV in the blood are investigated by FEA and theoretical analysis. A possible treatment for COVID-19 based on vibration is put forward. Several findings are listed below. (1) The frequencies of the tuned 2019-nCoV model are very close to each other, and they are all first-order bending vibrations. The frequencies of the realistic mistuned 2019-nCoV model are a range of frequencies close to each other, and every frequency corresponds to the first-order bending vibration of one spike protein. (2) The first-order bending vibration is adopted as the vibration form, and the established lumped parameter mechanical model of 2019-nCoV is in good agreement with the FEA model. When applying an excitation on the 2019-nCoV model, the vibration amplitude increases with cos h i j j, and it can somehow achieve response concentration. (4) In order to help cure the COVID-19, a potential method utilizing 360 o sweep excitation to destroy the spike protein structures is put forward. When an ultrasonic vibration exciter is applied 360°rotating sweep excitation to human body, the vibration frequencies is set as 1.9 9 10 8 Hz to 2.0 9 10 8 Hz and the amplitude is selected as larger than 1.041 9 10 -5 mm, it is efficient and safe to disable the 2019-nCoV from infecting human cells, which maybe help cure COVID-19. World Health Organization. WHO statement regarding cluster of pneumonia cases in Wuhan National Health Commission of The People's Republic of China. Interim protocol of diagnosis and treatment of 2019 novel coronavirus-associated pneumonia (the second version) A novel coronavirus from patients with pneumonia in China World Health Organization. Novel Coronavirus (2019-nCoV) Situation Report-22 Application of nonlocal continuum models to nanotechnology Fluid Mechanics Fluid Mechanics Vibration Problems in Engineering Mechanics of Materials Human tolerance to impact acceleration Resonant frequencies of standing humans Conflict of interest The authors declare that there is no conflict of interest regarding the publication of this paper.