亚洲免费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孫志強GuolinHou侯國林YanfenQiao喬艷芬andJincunLiu劉金存
        Chinese Physics B 2024年1期
        關鍵詞:侯國

        Zhiqiang Sun(孫志強), 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
        西漢沛郡所轄侯國地理探賾
        三個老同學辦了個招嫖論壇
        方圓(2022年1期)2022-03-01 13:22:44
        復原西漢侯國地理全貌的扛鼎之作 《西漢侯國地理(修訂本)》評價
        藏書報(2021年19期)2021-05-30 10:48:04
        南朝陳縣級封國地理與置省原因探微
        《東漢政區(qū)地理》縣級政區(qū)補考
        南都學壇(2016年2期)2016-04-27 08:41:00
        Virtual Property:the Realization of Rights and Value〔* 〕
        學術(shù)界(2015年1期)2015-02-25 07:11:31
        松花江上
        在线人成免费视频69国产| 毛片在线视频成人亚洲| 亚洲av久播在线一区二区| 无码国产福利av私拍| 日韩亚洲欧美中文高清在线| 91久久精品一二三区蜜桃| 亚洲一区二区刺激的视频| 影音先锋色小姐| 麻豆精产国品| 日韩精品一区二区亚洲av性色| 亚洲精品98中文字幕| 久久久久亚洲精品无码蜜桃| 好男人视频在线视频| 亚洲人成网站在线播放小说| 日韩一区二区三区熟女| 青青草国产精品一区二区| 双腿张开被9个黑人调教影片| 日韩美女av二区三区四区| 久草手机视频在线观看| 人妻少妇偷人精品无码| 在线播放a欧美专区一区| 日韩国产有码精品一区二在线| 亚洲国产精品一区二区成人av| 亚洲欧洲精品无码av| 人妻无码中文专区久久综合| 久久老熟女乱色一区二区| 99国产精品久久久久久久成人热| 国产精品久久久久久久免费看| 亚洲 暴爽 AV人人爽日日碰| 在线亚洲日本一区二区| 亚洲欧美色一区二区三区| 欧美三级不卡视频| 亚洲精品中文字幕乱码人妻| 又硬又粗进去好爽免费| 精品无码国产自产野外拍在线| 精品国产1区2区3区AV| 国产综合开心激情五月| 99国产精品人妻噜啊噜| 久久人人做人人妻人人玩精| 国产激情小视频在线观看| 大学生高潮无套内谢视频|