亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        Hamiltonian system for the inhomogeneous plane elasticity of dodecagonal quasicrystal plates and its analytical solutions

        2024-01-25 07:14:16ZhiqiangSun孫志強(qiáng)GuolinHou侯國林YanfenQiao喬艷芬andJincunLiu劉金存
        Chinese Physics B 2024年1期
        關(guān)鍵詞:侯國

        Zhiqiang Sun(孫志強(qiáng)), Guolin Hou(侯國林), Yanfen Qiao(喬艷芬), and Jincun Liu(劉金存)

        School of Mathematical Sciences,Inner Mongolia University,Hohhot 010021,China

        Keywords: Hamiltonian system,symplectic elasticity,quasicrystals,analytic solution,state function

        1.Introduction

        Quasicrystals (QCs) are a kind of solid intermediate between crystal and amorphous materials, which were first discovered by Shechtmanet al.[1]The unique quasi-periodic structure makes QCs become a new kind of functional materials and structural materials.[2,3]

        QCs have many excellent properties, such as high strength, abrasion resistance, hardness, thermal conductivity,and low adhesion.[4,5]It should be noted that the properties of QCs are affected by defects such as dislocations,cracks,holes,and inclusions.[6–12]Many methods and techniques have been developed to address the problems of elasticity and defects in QCs.Among them, symplectic elasticity approach,[13–15]the pseudo-Stroh formalism method,[16,17]the state-space method,[18,19]the integral transformation method,[20]the differential quadrature method,[21,22]the complex variable function method,[23–26]and the state-space based differential quadrature method[27]are practical to construct the analytical or numerical solutions of QCs.The integral transformation and differential quadrature methods are numerical and have been approved to be highly efficient for the rapid solution.The state-space-based differential quadrature method[27]inherits the advantages of the state-space method and the differential quadrature method.Its convergence speed is fast and can apply to different boundary conditions.However,the state-space method must be revised to solve the clamped-supported and mixed boundary conditions.The complex variable function method often leads to high-order partial differential equations,which puts forward higher requirements for the smoothness of solutions and may be difficultly solved further by separating variables.[23]The symplectic elasticity approach is the method of separating variables based on the Hamiltonian system.It can construct the general solution of the considered problem without any assumptions about the solution forms.According to the Vainberg theorem[28]in variational theory,a symmetric differential system can continually be transformed into a suitable Hamiltonian system, which provides the possibility for applications of the symplectic elasticity approach in elasticity mechanical problems because the considered mechanical equations are primarily symmetric.Recently, the symplectic elasticity approach has been applied to establish the final governing equation of Laue class 15,[14]which provides a feasible idea for solving some complex quasicrystal problems.

        This paper aims to extend the symplectic elasticity approach to solve the inhomogeneous dodecagonal QC plates.In the geometry of symplectic space, the phonon stress, phason stress, and displacements are set as unknown functions and included in the state vector, and then the formed first-order differential system,a Hamiltonian system,is introduced.Additionally, using the symplectic orthogonality relations of the eigenvectors,the general solutions of phonon and phason displacement under the given boundary conditions are obtained.Finally, the accuracy of the symplectic elasticity approach is verified by numerical examples of QC plates under uniform and concentrated loads, and the effects of different phason elastic constants on displacements and stresses are analyzed.The numerical results show that the symplectic elasticity approach is an effective tool for analyzing the mechanical behavior of QC plates.The presented method in the current paper can also be used to solve other boundary conditions.

        2.Hamiltonian system of dodecagonal QCs

        According to the linear elastic theory of QCs, the basic equations for the inhomogeneous plane elasticity problem of dodecagonal QCs are described in the following.[29]The equilibrium equations are

        and generalized Hooke’s laws have the following form:

        whereux,uyandwx,wydenote the phonon displacements and phason displacements;σxx,σxy,σyx, andσyyare the phonon stresses;Hxx,Hxy,Hyx,andHyyare the phason stresses;q1,q2,q3, andq4are the body forces;L,M,K1,K2, andK3are the elastic constants.For simplicity,K1+K2+K3is abbreviated toK4.

        To establish the Hamiltonian system,define the displacement vector as

        where the superscript T represents the transpose.The strain energy densityU(q, ˙q)can be expressed as

        where ˙()≡?/?y.Moreover,we introduce the Lagrange density function as

        Consider a dodecagonal QC rectangular domain with lengthl(0≤y ≤l)and widthh(0≤x ≤h).To build the variational operation,the equilibrium Eqs.(1a)–(1d) can be obtained by taking the variation with respect toux,uy,wx,andwy,respectively.Applying the Legendre transformation, the corresponding dual vector is given as follows:

        According to Eqs.(3)and(6),defining a state vector

        the Hamiltonian system can be expressed as

        where the Hamiltonian operator reads

        f(x,y)=?(q3,q4,q1,q2,0,0,0,0)Tis an inhomogeneous term related to body forces.Moreover, we consider Eq.(8)to satisfy the following boundary conditions shown in Fig.1,

        Remark 1 The Hamiltonian system (8) derived here is also applicable to other complex boundary conditions,such as clamped or free.The solution to those problems may require the use of the symplectic superposition methods,[30,31]which are not covered in this study.

        Fig.1.Geometry of boundary conditions(10).

        3.Eigenvalues and eigenvectors of H

        According to the theory of linear differential equations,we should first solve the following homogeneous form of Eq.(8):

        Equation (11) can be solved by the method of separating variables.[32]Assuming

        and substituting it into Eq.(11),we arrive at

        whereψ(y)is easy to get,andΦ(x)is an unknown eignvector corresponding to the unknown eigenvalueμ.The eigenvalue problem

        with the boundary condition (10) is needed to take a lot of effort to solve.According to Eq.(10),we find that the eigenvector has the following form:

        whereαandAi(i=1,2,...,8)are unknown constants which depend onμand can be calculated from Eq.(14).

        Combining Eqs.(10)and(14),the characteristic polynomial is derived as follows:

        From Eq.(15),we obtain the eigenvaluesμn=αn=±nπ/h,μ?n=?μn(n= 1,2,...), and corresponding eigenvectors with respect toμnare given by

        Similarly,the basic eigenvectors ofμ?n=?μn(n=1,2,...)are

        and first-order Jordan form eigenvectors ofμ?n=?μn(n=1,2,...)are

        Additionally,μ=0 is also a multiple eigenvalue ofH.The eigenvectors ofHwith respect to eigenvalueμ=0 include basic eigenvectors

        and first-order Jordan form eigenvectors

        4.Symplectic orthogonality and completeness

        The symplectic orthogonality of the generalized eigenvectors of the Hamiltonian operatorHis analyzed first since it plays an important role in representing the solution of the Hamiltonian system.[32]

        In addition,we have

        The following theorem is one of the main results in this paper, which characterizes the completeness of the generalized eigenvector system ofHand also ensures the feasibility of separating variables in the Hamiltonian system.It is worth noting that all of the calculations in this section can be easily achieved by symbolic software like Mathematica.

        Proof According to the symplectic orthogonal relationship(16)and(17),we take

        Thus,we have

        where the eight components in the above expression are the Fourier series expansions of functionsfi(x) (i= 1,2,...,8)under the orthogonal systemorin Hilbert spaceL2(0,h).Accordingly,

        is valid,and the proof is completed.

        5.Analytic solution of the original problem

        According to Theorem 1 and the principle of superposition of solutions, the general solution of Eq.(8)can be given in the form

        Let the Fourier expansion of the inhomogeneous termf(x,y)=(q4,q3,q1,q2,0,0,0,0)Tin Eq.(8)be

        where

        The details of(i= 1,2,3,4) are expressed in Appendix A.

        Substituting Eqs.(18)and(19)into Eq.(8),and by simple calculations,we obtain

        According to the form of the state vectorU(x,y),the analytic solutions ofHxy,Hyy,σxy,σyy,wx,wy,ux, anduycan be obtained by taking the first component through the eighth component of Eq.(20).For the sake of brevity,we do not write out their specific expressions.

        6.Numerical examples

        This section presents the behavior of solutions of the QC plates under the concentrated load,uniform load,and mechanical load based on the symplectic analytical solution obtained in Section 5.To this end,we consider the following boundary conditions in they-direction:

        For convenience of numerical simulation,we assume that the body forcesq1=q2=q3=q4.Due to the lack of comparative data, we introduce the FITM for the displacementsux,uy,wx,andwyto show the validity of the symplectic elasticity approach.The expressions of phonon and phason displacements given by the FITM can be found in Appendix B.The material constantsL=0.5741 GPa,M=0.88445 GPa,K1=1.22 GPa,K2=0.24 GPa, andK3=0.6 GPa are taken as those in Ref.[13].

        6.1.Concentrated load and uniformly distributed load

        The central concentrated loadq1=q2=q3=q4=Pδ(x ?x0,y ?y0) is examined at the (x0,y0)=(h/2,l/2)of plate,wherePis the intensity,andδ(·,·)denotes the wellknown Dirac function.As the ratio of length and width varies,the results of numerical tests under central concentrated load are shown in Tables 1 and 2 for the displacement responses at fixed horizontal coordinates (x,y)=(h/4,l/4) and (x,y)=(h/4,h/4),respectively.The uniformly distributed load of intensityqis examined, and the numerical results are given in Table 3.

        Table 1 shows that the symplectic elasticity approach performs better than the FITM since it gets a reasonable convergence rate with only 20 terms.Also,we use the first 20 terms of Eq.(20)to show the phonon and phason displacements with different length–width ratios(1≤l/h ≤2)in Fig.2.

        As shown in Fig.2, the tendency of change in displacements with a change in aspect length–width ratios of the QCs plate is different for the responses at two fixed horizontal coordinates.For the fixed valueh,the changes ofuxandwx(uyandwy)at(h/4,l/4)are increasing(decreasing)aslgoes up.This trend depends on the influences of boundary conditions (10)and (21), and the variations of distance between (h/4,l/4)and (h/2,l/2).However, the responsesux,wx,uyandwyat the fixed point (h/4,h/4) are decreasing aslbecomes larger because of the influence of boundary conditions is unchanged,and the effect of the load is gradually weakening.

        In addition,whenh=l,we have the phason displacementux=uyand the phonon displacementwx=wyfor the centralized load, which is in agreement with the case for uniformly distributed load.

        Fig.2.The phonon and phason displacements of the QC plate under a centrally concentrated load: (a) the response at (x,y)=(h/4,l/4),(b)the response at(x,y)=(h/4,h/4).

        Table 1.The calculation results under central concentrated load(x,y)=(h/4,l/4).

        Table 2.Calculation results under central concentrated load(x,y)=(h/4,h/4).

        6.2.Mechanical load

        In this subsection we analyze the effect of the phason elastic constants on the displacements and stresses.In the absence of body forces,we assume that the right side of the QC plate is subjected to the phonon and phason mechanical loads,i.e.,

        whereH0andσ0are the amplitude of the loading.Meanwhile,all other components on the left and right sides of the plate are zero.Consequently, the boundary conditions along theydirection are

        wherepandqare expressed in Eqs.(3)and(6),respectively.

        Following the previous work,[33]we takeH0=σ0=1 N/m2,α=π/h, the dimensions of the plate areh×l=0.3 m×0.3 m,and consider the responses wheny=0.75l.

        Based on the above loads, we analyze the effect of the phase elastic constantsKi(i=1,2,3)on some physical components in the QC plates.We let the phason elastic constants beKi,0.5Ki,0.1Ki(i=1,2,3),and the results of the calculations are shown in Fig.3.It can be seen that the magnitudes of the phason displacementswxandwydecrease with decreasing value of phason elastic constantsKi(i=1,2,3)[see Figs.3(c)and 3(d)], while there is no effect on the trends of the phason and phonon stresses components[see Figs.3(a)and 3(b)].Specifically, the phason displacements are more sensitive to the phason elastic constants than the phason and phonon stress.

        Fig.3.Displacements and stresses with different phason elastic constants: (a) phason stress Hyy, (b) phonon stress σxy, (c) phason displacement wx,and(d)phason displacement wy.

        Table 3.Calculation results under uniform load.

        7.Conclusions

        In summary, we have solved QC equations by employing the symplectic elasticity approach.By introducing a generalized displacement vector, the Hamiltonian system of dodecagonal QCs is derived from the strain energy density and Legendre transformation.Then, the symplectic elasticity approach is adopted to get the general solutions of the Hamiltonian system.It is worth noting that the 1st, 2nd, 3rd, and 4th components of the general solution (20) are the phonon and phason stresses,and the 5th,6th,7th,and 8th components of (20) are the phonon and phason displacements.The displacements and stresses responses of QC plates under the concentrated,uniformly distributed,and mechanical loads are discussed in detail.For the case of concentrated load,the trends of responses at the fixed point and moving point are different from the changes in length–width ratios.The numerical examples of the mechanical load show that the magnitude of the phason displacements decreases with decreasing value of phason elastic constants.The comparison data are constructed by FITM,which justifies calculating of the symplectic elasticity approach.The approach is direct and rational,providing a systematic step to solve the related problems of QCs.

        Appendix A

        Appendix B

        Substituting Eqs.(2a)–(2g)into Eqs.(1a)–(1d),we obtain

        where ?2=?2/?x2+?2/?y2.

        Based on boundary conditions(10)and(21),we can define eight double finite integral transforms with respect towx,wy,ux,uy,q1,q2,q3,andq4as follows:

        whereαm=mπ/handβn=nπ/l(m,n=1,2,3,...).

        The inverse transforms are given by

        By simple calculations, Eq.(B1) is transformed into the following form:

        Acknowledgments

        Project supported by the National Natural Science Foundation of China (Grant Nos.12261064 and 11861048), the Natural Science Foundation of Inner Mongolia,China(Grant Nos.2021MS01004 and 2022QN01008), and the High-level Talents Scientific Research Start-up Foundation of Inner Mongolia University(Grant No.10000-21311201/165).

        猜你喜歡
        侯國
        一官半職
        鴨綠江(2024年10期)2024-02-28 00:00:00
        一路沖鋒
        軍嫂(2023年9期)2023-10-12 12:52:40
        Sub-stochiometric MoOx by radio-frequency magnetron sputtering as hole-selective passivating contacts for silicon heterojunction solar cells
        西漢沛郡所轄侯國地理探賾
        三個老同學(xué)辦了個招嫖論壇
        方圓(2022年1期)2022-03-01 13:22:44
        復(fù)原西漢侯國地理全貌的扛鼎之作 《西漢侯國地理(修訂本)》評價(jià)
        南朝陳縣級封國地理與置省原因探微
        《東漢政區(qū)地理》縣級政區(qū)補(bǔ)考
        Virtual Property:the Realization of Rights and Value〔* 〕
        松花江上
        精品久久久久香蕉网| 99热婷婷一区二区三区| 国产日韩亚洲欧洲一区二区三区| 欧美巨大xxxx做受中文字幕| 日本一区二区高清精品| 日本人妻伦理在线播放| 国产精品久久人妻无码| 日韩亚洲欧美中文高清在线| 91蜜桃国产成人精品区在线| 丰满的少妇av一区二区三区| 狠狠精品久久久无码中文字幕 | 亚洲精品自产拍在线观看| 精品一级毛片| 亚洲视频精品一区二区三区| 自拍偷自拍亚洲精品第按摩| 欧美大成色www永久网站婷| 国产综合自拍| 加勒比一区二区三区av| 国产99久久久国产精品~~牛| 国产一区二区三区在线电影| 国产天堂网站麻豆| 中文字幕日本一区二区在线观看| 久久一二区女厕偷拍图| 精品精品国产自在97香蕉| 抽插丰满内射高潮视频| 91久久精品人妻一区二区| 亚洲麻豆视频免费观看| 亚洲丁香五月天缴情综合| 人妖精品视频在线观看| 日本高清一区二区不卡| 久久国产色av免费观看| 在线精品国产一区二区| 色二av手机版在线| 少妇高潮精品正在线播放| 国产免费又爽又色又粗视频| 国产人澡人澡澡澡人碰视频| 日本第一区二区三区视频| 国产一区二区视频免费在| 亚洲av永久无码精品一区二区| 中文字幕无码免费久久9一区9| 亚洲国产一区二区中文字幕|