Parameter Identification of Weakly Nonlinear Vibration System in Frequency Domain

Total Page:16

File Type:pdf, Size:1020Kb

Parameter Identification of Weakly Nonlinear Vibration System in Frequency Domain Shock and Vibration 11 (2004) 685–692 685 IOS Press Parameter identification of weakly nonlinear vibration system in frequency domain Jiehua Penga,b,∗, Jiashi Tanga and Zili Chena,c aDepartment of Mechanics, Hunan University, Changsha, 410082, P.R. China bDepartment of physics, Shaoyang University, Shaoyang, 422000, P.R. China cCollege of Architectural Engineering, Central South Forestry University, Zhuzhou, 412008, P.R. China Received 1 November 2001 Revised 11 December 2003 Abstract. A new method of identifying parameters of nonlinearly vibrating system in frequency domain is presented in this paper. The problems of parameter identification of the nonlinear dynamic system with nonlinear elastic force or nonlinear damping force are discussed. In the method, the mathematic model of parameter identification is frequency response function. Firstly, by means of perturbation method the frequency response function of weakly nonlinear vibration system is derived. Next, a parameter transformation is made and the frequency response function becomes a linear function of the new parameters. Then, based on this function and with the least square method, physical parameters of the system are identified. Finally, the applicability of the proposed technique is confirmed by numerical simulation. Keywords: Parameter identification, weakly nonlinear vibration system, least square method, frequency domain 1. Introduction The parameter identification is one of the important problems of the vibration research. There are many methods of identifying parameters of the linear system. But the parameter identification of nonlinear system is much more difficult than the linear system. Hence, this problem attracts much attention. A. H. Nayfeh has presented an excellent method of parametric identification of nonlinear system that exploits nonlinear resonances and comparisons of the system to be identified with those of known system [1], but the difficult usually appears when we find the known system. Tang et al. [2–4] have suggested time domain method and phase plane method for identifying parameters of nonlinear system which are effective for integrable nonlinear system. K. Yasuda et al. [5,6] have proposed an experiment method of parameter identification which is useful for geometrically nonlinear systems. R. Bachmayer et al. [7] have established adaptive parameter identification method for marine thruster which is an on-line technique for adaptive identification. There are some other methods of parameter identification, such as orthogonal function method [8], wavelet-based technique [9], Volterraseries method [10], cross-correlation analysis [11], power spectrum technique [12], and so on. A new method of identifying parameters of the weakly nonlinear vibration system in the frequency domain is presented in this paper. We discuss two kinds of nonlinear vibration systems with nonlinear elastic force or nonlinear damping force respectively. First, by using the multiple scale method, we obtain the frequency response function of the system which is regarded as a mathematical model of the curve fitting. Next, a new set of parameters are introduced to replace the parameters to be identified and the frequency response function is transformed into the linear function of the new parameters so that the least square method can be used to identify the parameters. This ∗Corresponding author. E-mail: [email protected]. ISSN 1070-9622/04/$17.00 2004 – IOS Press and the authors. All rights reserved 686 J. Peng et al. / Parameter identification of weakly nonlinear vibration system in frequency domain is the key step of our method. Then, with the least square method the normal equations are established and they are linear equations of the new parameters. From the solution of normal equations we can get the values of the parameters to be identified. The method we present is a combination of the perturbation method and the least square method, hence, the method has the same precision with those of perturbation method and least square method. On the other hand, the perturbation method is an effective method to deal with weakly nonlinear system and the least square method is an excellent method of curve fitting such that this parameter identification method has a wide application in the study of weakly nonlinear vibration system. Although frequency response function is the first order approximation of the relationship between amplitude and frequency of nonlinear system, it is the most useful and important mathematical model to study the property of nonlinear vibration system [13]. It can be determined by experiment with good accuracy because in practice any complicated nonlinear system will finally move in a steady motion state if the excitation is harmonic and the system doesn’t come into chaotic motion. Compared with response function in time domain, frequency response function is independent of initial conditions and the test to determine this relationship possesses good repeatability and stability, such that this method of parameter identification can be easily realized by experiment. 2. Principle of the method The governing equation of the vibration system with nonlinear damping force is as follows mx¨ + cx˙ + f(˙x)+kx = G cos Ωt (1) The governing equation of the vibration system with nonlinear elastic force is as follows mx¨ + cx˙ + kx + f(x)=G cos Ωt (2) where f(˙x) is a nonlinear function of x˙, f(x) is a nonlinear function of x. Rewrite Eqs (1) and (2) into the standard form 2 x¨ + ω x = ε[−2nx˙ − f1(˙x)+g cos Ωt] (3) and 2 x¨ + ω x = ε[−2nx˙ − f1(x)+g cos Ωt] (4) k c where ε being a small positive parameter, ω = m being the linear natural frequency of the system, n = 2m , f(˙x) f(x) G f1(˙x)= m , f1(x)= m and g = m . When Ω approaches to ω, we introduce a detuning parameter σ, letting Ω=ω + εσ (5) We can find the approximate frequency response function of the nonlinear system Eqs (3) and (4) by using approximate analytical method, such as the methods of iteration, multiple scales, KBM and so on [13–15]. The general approximate frequency response function of Eqs (3) and (4) can be written in the form F (ω,n,...,a,σ)=0 (6) Where ω,n,..., are physical parameters to be identified of the system Eqs (3) and (4), a is the amplitude, σ is the detuning parameter whose dimension is the same with frequency. Equation (6) plays key role in revealing the properties of nonlinear vibration system and here we take it as a mathematical model of the curve fitting. But F (ω,n,...,a,σ) is a nonlinear function of parameters (ω,n,...,) we must replace the parameters (ω,n,...,) with a new set of parameters (A,B,...,) before employing the least square method to identify the parameters, such that F (ω,n,...,a,σ) can be transformed into a linear function of (A,B,...,). That is, if we chose (ω,n,...,)=R(A,B,...,) (7) then J. Peng et al. / Parameter identification of weakly nonlinear vibration system in frequency domain 687 F (ω,n,...,a,σ)=F1(A,B,...,a,σ)=0 (8) Where F1(A,B,...,a,σ) is a linear function of parameters (A,B,...,). Transform function R(A,B,...,) is determined by the frequency response function of the systems. We use the test data (σ i,ai)(i =1, 2,...,m) for fitting the curve. Substituting σi,ai into Eq. (8) yields small errors Ei = F1(A,B,...,ai,σi) i =1, 2,...,m The sum of the square of all errors is m m 2 2 E = Ei = [F1(A,B,...,ai,σi)] (9) i=1 i=1 It is obvious that the parameters ω,n...will fit to make the totle square errors E get the minimum. This means that we can use the method of the minimum squares to identify the parameters of the system. By determining the extreme value of E and solving a set of normal equations resulted from the extreme value of E, we can identify parameters ω,n,...of the nonlinear system. 3. The system with nonlinear damping Let us consider a weakly nonlinear vibration system with a small cubic nonlinear damping force, its governing equation is as follows. 3 mx¨ + cx˙ + c1x˙ + kx = G cos Ωt (10) Rewrite the standard equation in the form mx¨ + ω2x = ε(−2nx˙ − αx˙ 3 + g cos Ωt) (11) where k g c G ω2 = , 2n = ,α= ,g= (12) m m m m ω,n and α are the physical parameters to be identified of the system Eq. (11). We begin by assuming an expansion of the solution having the form x(t, ε)=x0(T0,T1)+εx1(T0,T1)+... (13) New independent variables are introduced according to n Tn = ε t for n =0, 1,... (14) Substituting Eq. (13) into Eq. (11) and equating the coefficients of ε i, (i =0, 1, 2,...) to zero, we obtain 2 2 D0x0 + ω x0 =0 (15) 2 2 3 D0x1 + ω x1 = −2D0D1x0 − 2nD0x0 − α(D0x0) + g cos(ωT0 + σT1) (16) The solution of the Eq. (15) is x0 = a cos(ωT0 + ϕ) (17) Substituting Eq. (17) into Eq. (16) leads to 1 2 + 2 =2 sin( + )+2 cos( + )+2 sin( + )+ 3 3[3 sin( + ) D0x1 ω x1 ωa ωT0 ϕ ωaβ ωT0 ϕ nωa ωT0 ϕ 4αa ω ωT0 ϕ − sin 3(ωT0 + ϕ)] + g cos(ωT0 + ϕ)cos(σT1 − ϕ) − g sin(ωT0 + ϕ)sin(σT1 − ϕ) To eliminate secular terms from x1, we must put 688 J. Peng et al. / Parameter identification of weakly nonlinear vibration system in frequency domain 3 3 3 2ωa +2naω + 4 αa ω −g sin(σT1 − ϕ)=0 (18) 2ωaϕ + g cos(σT1 − ϕ)=0 Letting σT1 − ϕ = γ (19) Determining the derivatives with respect to T1 yields ϕ = σ − γ . Considering the steady-state response, that is, a =0,γ =0,wehave 2naω + 3 αa3ω3 − g sin γ =0 4 (20) 2ωaσ + g cos γ =0 The result by eliminating γ in Eq.
Recommended publications
  • List of Medical Device Clinical Trial Filing Institutions
    List Of Medical Device Clinical Trial Filing Institutions Serial Record number Institution name number Beijing: 5 6 Ge Mechanical temporary 1 agency Beijing Tsinghua Chang Gung Memorial Hospital preparation 201800003 Mechanical temporary 2 agency Plastic Surgery Hospital of Chinese Academy of Medical Sciences preparation 201800008 Mechanical temporary 3 agency Beijing Youan Hospital, Capital Medical University preparation 201800019 Mechanical temporary 4 agency Peking University Shougang Hospital preparation 201800044 Mechanical temporary 5 agency Beijing Cancer Hospital preparation 201800048 Mechanical temporary 6 agency Eye Hospital of China Academy of Chinese Medical Sciences preparation 201800077 Mechanical temporary Beijing Traditional Chinese Medicine Hospital Affiliated to Capital Medical 7 agency University preparation 201800086 Mechanical temporary 8 agency Beijing Anorectal Hospital (Beijing Erlong Road Hospital) preparation 201800103 Mechanical temporary 9 agency Cancer Hospital of Chinese Academy of Medical Sciences preparation 201800108 Serial Record number Institution name number Mechanical temporary Peking Union Medical College Hospital, Chinese Academy of Medical 10 agency Sciences preparation 201800119 Mechanical temporary 11 agency Beijing Luhe Hospital, Capital Medical University preparation 201800128 Mechanical temporary 12 agency Beijing Huilongguan Hospital preparation 201800183 Mechanical temporary 13 agency Beijing Children's Hospital, Capital Medical University preparation 201800192 Mechanical temporary 14 agency
    [Show full text]
  • A Complete Collection of Chinese Institutes and Universities For
    Study in China——All China Universities All China Universities 2019.12 Please download WeChat app and follow our official account (scan QR code below or add WeChat ID: A15810086985), to start your application journey. Study in China——All China Universities Anhui 安徽 【www.studyinanhui.com】 1. Anhui University 安徽大学 http://ahu.admissions.cn 2. University of Science and Technology of China 中国科学技术大学 http://ustc.admissions.cn 3. Hefei University of Technology 合肥工业大学 http://hfut.admissions.cn 4. Anhui University of Technology 安徽工业大学 http://ahut.admissions.cn 5. Anhui University of Science and Technology 安徽理工大学 http://aust.admissions.cn 6. Anhui Engineering University 安徽工程大学 http://ahpu.admissions.cn 7. Anhui Agricultural University 安徽农业大学 http://ahau.admissions.cn 8. Anhui Medical University 安徽医科大学 http://ahmu.admissions.cn 9. Bengbu Medical College 蚌埠医学院 http://bbmc.admissions.cn 10. Wannan Medical College 皖南医学院 http://wnmc.admissions.cn 11. Anhui University of Chinese Medicine 安徽中医药大学 http://ahtcm.admissions.cn 12. Anhui Normal University 安徽师范大学 http://ahnu.admissions.cn 13. Fuyang Normal University 阜阳师范大学 http://fynu.admissions.cn 14. Anqing Teachers College 安庆师范大学 http://aqtc.admissions.cn 15. Huaibei Normal University 淮北师范大学 http://chnu.admissions.cn Please download WeChat app and follow our official account (scan QR code below or add WeChat ID: A15810086985), to start your application journey. Study in China——All China Universities 16. Huangshan University 黄山学院 http://hsu.admissions.cn 17. Western Anhui University 皖西学院 http://wxc.admissions.cn 18. Chuzhou University 滁州学院 http://chzu.admissions.cn 19. Anhui University of Finance & Economics 安徽财经大学 http://aufe.admissions.cn 20. Suzhou University 宿州学院 http://ahszu.admissions.cn 21.
    [Show full text]
  • University of Leeds Chinese Accepted Institution List 2021
    University of Leeds Chinese accepted Institution List 2021 This list applies to courses in: All Engineering and Computing courses School of Mathematics School of Education School of Politics and International Studies School of Sociology and Social Policy GPA Requirements 2:1 = 75-85% 2:2 = 70-80% Please visit https://courses.leeds.ac.uk to find out which courses require a 2:1 and a 2:2. Please note: This document is to be used as a guide only. Final decisions will be made by the University of Leeds admissions teams.
    [Show full text]
  • Analysis and Research on Defrosting Ventilation Pipe Deformation of Automobile
    Analysis and Research on Defrosting Ventilation Pipe Deformation of Automobile Guangming Li1, 2, Congrui Zuo3, a, *, Yi Ma1, 2 and Qike Huang1, 2 1Department of Mechanical and Energy Engineering, Shaoyang University, Shaoyang 422000,China 2Key Laboratory of Hunan Province for Efficient Power System and Intelligent Manufacturing, Shaoyang University,Shaoyang 422000,China 3Institute Of Thermal And Energy Metrology, Hunan Institute Of Metrology And Test ,Changsha 410000,China Keywords: Car, defrosting ventilation pipe, Deformation, The finite element. Abstract: According to the market feedback, automobile defrosting ventilation pipe does not fit in the top 10 of xiangtan geely automobile market quality problems. According to the detection of vehicles in this condition, the reason for the deformation of automobile defrosting ventilation pipe position is that the gap between instrument cover and ventilation pipe is less than 4mm. Ventilation tube ends around not level off, the present state of collapse, and the intensity of this state collapse is due to the lack of overall deformation, for deformation reasons of automobile ventilation pipe outlet increase 11 on the back of the improved design scheme of reinforcement, using the method of three-dimensional modeling, design parameters, according to auto defrost ventilation tube through the Pro/E 3 d software to complete the ventilation pipe before and after the improvement model, using ANSYS Workbench software of finite element analysis was carried out on the ventilation pipe before and after improvement, the maximum deformation analysis ventilation pipe outlet, After the improvement, the deformation amount decreased from 10.643mm to 10.283mm, and the maximum stress value of the ventilation pipe decreased from 85.783MPa to 55.781MPa, less than the ultimate strength of 70MPa, verifying the rationality of the improved design, Solve the quality problem for the enterprise, have very good reference value.
    [Show full text]
  • Comparative Clinical Characteristics Among Different
    Comparative Clinical Characteristics Among Different Age Group of Adult COVID-19 Patients yating peng second xiangya hospital of central south university https://orcid.org/0000-0001-8084-8471 Aiyuan Zhou second xiangya hospital Xin Liao aliated shaoyang central hospital of university of south china Dingding Deng the rst attached hospital of shaoyang univerisity Peng Huang zhuzhou central hospital Wenlong Liu yueyang second people's hospital Mingyan Jiang xiangtan central hospital Xudong Xiang emergency of second xiangya hospital Qingcui Shuang the rst attached hospital of shaoyang university Shan Cai respiratory department of second xiangya hospital Ping Chen ( [email protected] ) https://orcid.org/0000-0001-6707-8636 Xucai Liao the rst attached hospital of shaoyang university Research article Keywords: COVID-19, severe or critical type, age Posted Date: May 15th, 2020 DOI: https://doi.org/10.21203/rs.3.rs-26179/v1 License: This work is licensed under a Creative Commons Attribution 4.0 International License. Read Full License Page 1/11 Abstract Introduction: Coronavirus disease (COVID-19) is a global infectious disease with a large burden of illness and high healthcare costs. Objectives: This study aimed to compare clinical features among adult COVID-19 patients in different age groups. Methods: Laboratory-conrmed adult COVID-19 infection cases between Dec 31, 2019 to March 8, 2020 obtained from Neighboring Cities. Patients were divided into ve age groups: age30, 30–40, 40–50, 50–65, ≥ 65 y (elderly). Age, sex, history of chronic disease and epidemiology, symptoms, laboratory tests, and outcomes were compared among different age groups. Binary logistic regression analysis was conducted to evaluate associated factors for severe or critical type.
    [Show full text]
  • Electric Power Steering System for Agricultural Truck Matching Design
    International Conference on Power Electronics and Energy Engineering (PEEE 2015) Electric Power Steering System for Agricultural Truck Matching Design Z.G. Chen Q. Zhang Mechanical and energy engineering department Mechanical and energy engineering department Shaoyang University Shaoyang University Shaoyang, Hunan, China Shaoyang, Hunan, China W. Zhang T.M. Zhou Mechanical and energy engineering department Zhuzhou Elite Electro Mechanical Co, Ltd Shaoyang University Zhuzhou, Hunan, China Shaoyang, Hunan, China Abstract—this thesis aimed at the characteristics and the manipulate and rotate the steering wheel. At the same time, structural of a type of truck. It confirmed that its chosen torque sensor will receive signals of engine speed and vehicle power-assisted form of electric power steering system is Column speed. Then it will transform into digital signal and input EPS——that is C-EPS. On the platform of a sort of car we control unit. After that it calculate we will get and make combined with the technical requirements of a certain model car driving cycle so good at adapting assisting torque. Next, steering system performance upgrades and matched design of the controlling assisted motor, outputting the corresponding size motor, the decelerating device, and the transducer. Finally device and direction, increasing torque by the decelerating device, it of the EPS in a certain automobile model truck for the finally applied to the gear and rack so that it perform characteristics test of input and output, it have shown that the power-assisted steering[1]. EPS has a good matching design of the booster effect. Keywords-EPS; matching; design; test I. INTRODUCTION With the development of products of the automotive industry, energy efficiency and environmental protection will become a trend .automotive electric power steering system( Electric Power Steering,referred to as EPS ) is the first choice for every large automobile company both at home FIGURE I.
    [Show full text]
  • Anti‑Inflammatory Activities of Gardenia Jasminoides Extracts in Retinal Pigment Epithelial Cells and Zebrafish Embryos
    EXPERIMENTAL AND THERAPEUTIC MEDICINE 22: 700, 2021 Anti‑inflammatory activities ofGardenia jasminoides extracts in retinal pigment epithelial cells and zebrafish embryos JIANRONG CHEN1, GABRIEL MBUTA TCHIVELEKETE2, XINZHI ZHOU2, WEIZHUO TANG1, FANG LIU1, MINZHUO LIU1, CHENXI ZHAO1, XINHUA SHU2‑4 and ZHIHONG ZENG1,5 1College of Biological and Environmental Engineering, Changsha University, Changsha, Hunan 410022, P.R. China; Departments of 2Biological and Biomedical Sciences and 3Vision Science, Glasgow Caledonian University, Glasgow, G4 0BA, UK; 4School of Basic Medical Sciences, Shaoyang University, Shaoyang, Hunan 422000; 5Hunan Provincial Key Laboratory of Nutrition and Quality Control Aquatic Animals, Changsha, Hunan 410022, P.R. China Received July 10, 2020; Accepted March 19, 2021 DOI: 10.3892/etm.2021.10132 Abstract. Age‑related macular degeneration (AMD) is the most and is the third leading cause of blindness after glaucoma and common cause of visual impairment in developed countries. cataracts (1). The prevalence of AMD by 2040 is projected Inflammation serves a critical role in the pathogenesis of AMD. to significantly increase to ~288 million cases (2). Different Gardenia jasminoides is found in several regions of China factors may contribute toward AMD, including age, envi‑ and is traditionally used as an organic yellow dye but has also ronmental factors and genetic susceptibility. Ageing may been widely used as a therapeutic agent in numerous diseases, strengthen the problem of visual difficulty associated with including inflammation, depression, hepatic and vascular disor‑ AMD. Patients with AMD have extracellular deposits called ders, which may reflect the variability of functional compounds drusen, which are accumulated within Bruch's membrane that are present in Gardenia jasminoides extracts (GJE).
    [Show full text]
  • Research on Creative Education, Integration of Industry and Education In
    3rd International Conference on Management, Education, Information and Control (MEICI 2015) Research on Creative Education, Integration of Industry and Education in the Application Talents Cultivation Jianefeng HU Jiangxi University of Technology, 330098 Nanchang, China [email protected] Keyword: Creative Education; Integration of Industry and Education; Application talents cultivation; Iocal college and university Abstract. Training application talents is the primary mission and the fundamental task of the local college and university. Deepen the reform of creative education in local college and university is a breakthrough point of the comprehensive reform of higher education. Integration of industry and education in college and university is an inevitable approach to cultivate application talents. This paper puts forward some thinking about the cultivation of application talents in local college and university from two perspectives: creative education, integration of industry and education. Introduction Cultivating application talents is the primary mission and basic task of local colleges and universities. Local colleges and universities must highlight the characteristics of "local", take the initiative to integrate into the local economic and social development needs, broaden the space for colleges and universities, in order to achieve the sustainable development of local colleges and universities. The relationship between higher education and economic social development has three kind of realm. The first is the active service, the second is comprehensive support, the third is leading innovation. Deepening the reform of creative and entrepreneurship education is the urgent need to accelerate the implementation of the innovation drive development strategy, is a breakthrough point to push forward the comprehensive reform of higher education, is great significance to promote graduates more high quality employment.
    [Show full text]
  • Experimental Researches on Friction Coefficient for Disc Cutters
    GUIJU ZHANG et al: EXPERIMENTAL RESEARCHES ON FRICTION COEFFICIENT FOR DISC CUTTERS Experimental Researches on Friction Coefficient for Disc Cutters Guiju Zhang 1,2, Caiyuan Xiao 1,* 1 Department of Mechanical and Energy Engineering Shaoyang University Shaoyang Hunan, China 2 College of Mechanical and Electrical Engineering Central South University, Changsha,Hunan, China Abstract — The disc cutter is the main tool for Tunnel Boring Machine (TBM for short), both the rock breaking efficiency and service life of which are directly related to the project schedule and cost. On the basis of contact stress characteristic for disc cutter and its movement rules, the mechanical model of disc is established, together with the micro model of friction coefficient. By using the MMW-1 Vertical universal friction and wear testing machine, the dry friction test bench is designed to simulate friction counterparts between disc cutter and rock. By experiments, the changing rules of friction coefficient with time are separately analyzed between H13 steel and sandstone, H13 and marble as well as H13 and granite. Finally, the theoretical and experimental results of friction coefficient’s changing rules with 3 different rocks are obtained, which can provide a theoretical reference for the optimization design of disc cutter and improve the study of disc cutter rock breaking efficiency and service life. Keywords- TBM; Disc cutter; Friction Coefficient; Rock breaking efficiency; Service life to analyze changing rules of friction coefficient, it provides I. INTRODUCTION some theoretical basis for optimization design of TBM disc Tunnel boring machine (TBM), which use rotating tool to cutter and improves the rock breaking efficiency and service tunnel and crush cave surrounding rock and boring at the life, so as to better guide engineering practice.
    [Show full text]
  • 1 Please Read These Instructions Carefully
    PLEASE READ THESE INSTRUCTIONS CAREFULLY. MISTAKES IN YOUR CSC APPLICATION COULD LEAD TO YOUR APPLICATION BEING REJECTED. Visit http://studyinchina.csc.edu.cn/#/login to CREATE AN ACCOUNT. • The online application works best with Firefox or Internet Explorer (11.0). Menu selection functions may not work with other browsers. • The online application is only available in Chinese and English. 1 • Please read this page carefully before clicking on the “Application online” tab to start your application. 2 • The Program Category is Type B. • The Agency No. matches the university you will be attending. See Appendix A for a list of the Chinese university agency numbers. • Use the + by each section to expand on that section of the form. 3 • Fill out your personal information accurately. o Make sure to have a valid passport at the time of your application. o Use the name and date of birth that are on your passport. Use the name on your passport for all correspondences with the CLIC office or Chinese institutions. o List Canadian as your Nationality, even if you have dual citizenship. Only Canadian citizens are eligible for CLIC support. o Enter the mailing address for where you want your admission documents to be sent under Permanent Address. Leave Current Address blank. Contact your home or host university coordinator to find out when you will receive your admission documents. Contact information for you home university CLIC liaison can be found here: http://clicstudyinchina.com/contact-us/ 4 • Fill out your Education and Employment History accurately. o For Highest Education enter your current degree studies.
    [Show full text]
  • The PI Index of Polyomino Chains$
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Elsevier - Publisher Connector Applied Mathematics Letters 21 (2008) 1101–1104 www.elsevier.com/locate/aml The PI Index of polyomino chains$ Lixing Xua,∗, Shubo Chenb a Department of Science and Information Science, Shaoyang University, Shaoyang, Hunan 422000, PR China b Department of Mathematics and Computer Science, Hunan City University, Yiyang, Hunan 413000, PR China Received 5 October 2006; received in revised form 11 November 2007; accepted 11 December 2007 Abstract The PI index is a graph invariant defined as the summation of the sums of neu(e|G) and nev(e|G) over all the edges e = uv of P a connected graph G, i.e., PI(G) = e∈E(G)[neu(e|G) + nev(e|G)], where neu(e|G) is the number of edges of G lying closer to u than to v and nev(e|G) is the number of edges of G lying closer to v than to u. An efficient formula for calculating the PI index of polyomino chains is given, and the bounds for the PI index of polyomino chains are established. c 2008 Elsevier Ltd. All rights reserved. Keywords: Graph invariant; PI index; Polyomino chain; Square lattice; Segment 1. Introduction The structure of a molecule could be represented in a variety of ways. The information on the chemical constitution of molecule is conventionally represented by a molecular graph. And graph theory was successfully provided the chemist with a variety of very useful tools, namely, topological indices.
    [Show full text]
  • A Conservative Difference Scheme for the Riesz Space-Fractional Sine-Gordon Equation Zhiyong Xing1,2* and Liping Wen1
    Xing and Wen Advances in Difference Equations (2018)2018:238 https://doi.org/10.1186/s13662-018-1689-5 R E S E A R C H Open Access A conservative difference scheme for the Riesz space-fractional sine-Gordon equation Zhiyong Xing1,2* and Liping Wen1 *Correspondence: [email protected] Abstract 1School of Mathematics and Computational Science, Xiangtan In this paper, we study a conservative difference scheme for the sine-Gordon University, Xiangtan, P.R. China equation (SGE) with the Riesz space fractional derivative. We rigorously establish the 2Department of Mathematics, conservation property and solvability of the difference scheme. We discuss the Shaoyang University, Shaoyang, P.R. China stability and convergence of the difference scheme in the L∞ norm. To reduce the computational complexity, we introduce a revised Newton method for implementing the difference scheme. Finally, we provide several numerical experiments to support the theoretical results. Keywords: Space-fractional sine-Gordon equation; Riesz fractional derivative; Conservation law; Newton method; Convergence and stability 1 Introduction The nonlinear SGE ∂2u – u + sin u =0, x =(x ,...,x ) ∈ Rn, (1.1) ∂t2 1 n arises in many different areas, such as stability of fluid motions, differential geometry, Josephson junctions, models of particle physics [1], the propagation of fluxon [2], the motion of a rigid pendulum attached to a stretched wire [3], the phenomenon of supra- transmission in nonlinear media [4], and so on. Therefore the investigation for the SGE has attracted attention of some researchers, and many significant achievements have been made. Aremarkablepropertyof(1.1) is the energy conservation: 1 ∂u 2 ∂u 2 + +2P(u) dx =const, 2 R ∂t ∂x where P(u)=1–cos u; it was studied by many researchers [5–9].
    [Show full text]