Quantum physics is primarily concerned with real eigenvalues,stemming from the unitarity of time evolutions.With the introduction of PT symmetry,a widely accepted consensus is that,even if the Hamiltonian of the syste...Quantum physics is primarily concerned with real eigenvalues,stemming from the unitarity of time evolutions.With the introduction of PT symmetry,a widely accepted consensus is that,even if the Hamiltonian of the system is not Hermitian,the eigenvalues can still be purely real under specific symmetry.Hence,great enthusiasm has been devoted to exploring the eigenvalue problem of non-Hermitian systems.In this work,from a distinct perspective,we demonstrate that real eigenvalues can also emerge under the appropriate recursive condition of eigenstates.Consequently,our findings provide another path to extract the real energy spectrum of non-Hermitian systems,which guarantees the conservation of probability and stimulates future experimental observations.展开更多
A numerical method is proposed to calculate the eigenvalues of the Zakharov–Shabat system based on Chebyshev polynomials. A mapping in the form of tanh(ax) is constructed according to the asymptotic of the potential ...A numerical method is proposed to calculate the eigenvalues of the Zakharov–Shabat system based on Chebyshev polynomials. A mapping in the form of tanh(ax) is constructed according to the asymptotic of the potential function for the Zakharov–Shabat eigenvalue problem. The mapping can distribute Chebyshev nodes very well considering the gradient for the potential function. Using Chebyshev polynomials, tanh(ax) mapping, and Chebyshev nodes, the Zakharov–Shabat eigenvalue problem is transformed into a matrix eigenvalue problem. This method has good convergence for the Satsuma–Yajima potential and the convergence rate is faster than the Fourier collocation method. This method is not only suitable for simple potential functions but also converges quickly for a complex Y-shape potential. It can also be further extended to other linear eigenvalue problems.展开更多
Wave propagation in horizontally layered media is a classical problem in seismic-wave theory.In semi-infinite space,a nondispersive Rayleigh wave mode exists,and the eigendisplacement decays exponentially with depth.I...Wave propagation in horizontally layered media is a classical problem in seismic-wave theory.In semi-infinite space,a nondispersive Rayleigh wave mode exists,and the eigendisplacement decays exponentially with depth.In a layered model with increasing layer velocity,the phase velocity of the Rayleigh wave varies between the S-wave velocity of the bottom half-space and that of the classical Rayleigh wave propagated in a supposed half-space formed by the parameters of the top layer.If the phase velocity is the same as the P-or S-wave velocity of the layer,which is called the critical mode or critical phase velocity of surface waves,the general solution of the wave equation is not a homogeneous(expressed by trigonometric functions)or inhomogeneous(expressed by exponential functions)plane wave,but one whose amplitude changes linearly with depth(expressed by a linear function).Theories based on a general solution containing only trigonometric or exponential functions do not apply to the critical mode,owing to the singularity at the critical phase velocity.In this study,based on the classical framework of generalized reflection and transmission coefficients,the propagation of surface waves in horizontally layered media was studied by introducing a solution for the linear function at the critical phase velocity.Therefore,the eigenvalues and eigenfunctions of the critical mode can be calculated by solving a singular problem.The eigendisplacement characteristics associated with the critical phase velocity were investigated for different layered models.In contrast to the normal mode,the eigendisplacement associated with the critical phase velocity exhibits different characteristics.If the phase velocity is equal to the S-wave velocity in the bottom half-space,the eigendisplacement remains constant with increasing depth.展开更多
The Dashuigou tellurium deposit is the world’s only known independent tellurium deposit.By restoring metamorphic rocks’protolith,we seek to understand not only the development and evolution trajectory of the region ...The Dashuigou tellurium deposit is the world’s only known independent tellurium deposit.By restoring metamorphic rocks’protolith,we seek to understand not only the development and evolution trajectory of the region but also the origin of the relevant deposits.While there are many ways to restore metamorphic rocks’protolith,we take the host metamorphic rocks of Dashuigou tellurium deposit and leverage various petrochemical eigenvalues and related diagrams previously proposed to reveal the deposit’s host metamorphic rocks’protolith.The petrochemical eigenvalues include molecular number,Niggli’s value,REE parity ratio,CaO/Al_(2)O_(3)ratio,Fe^(3+) /(Fe^(3+) -+Fe^(2+) )ratio,chondrite-normalized REE value,logarithmic REE value,various REE eigenvalues including scandium,Eu/Sm ratio,total REE amount,light and heavy REEs,δEu,Eu anomaly,Sm/Nd ratio,and silicon isotope δ^(30) SiNBS-29‰,etc.The petrochemical plots include ACMs,100 mg-c-(al+alk),SiO_(2)-(Na_(2)O+K_(2)O),(al+fm)-(c+alk)versus Si,FeO+Fe_(2)O^(3+) TiO)-Al_(2)O_(3)-MgO,c-mg,Al_(2)O_(3)-(Na_(2)O+K_(2)O),chondrite-normalized REE model,La/Yb-REE,and Sm/Nd ratio,etc.On the basis of these comprehensive analyses,the following conclusions are drawn,starting from the many mantle-derived types of basalt developed in the study area of different geological ages,combined with the previously published research results on the deposit s fluid inclusions and sulfur and lead isotopes.The deposit is formed by mantle degassing in the form of a mantle plume in the late Yanshanian orogeny.The degassed fluids are rich in nano-sc ale substances including Fe,Te,S,As,Bi,Au,Se,H_(2),CO_(2),N_(2),H_(2)O,and CH_(4),which are enriched by nano-effect,and then rise to a certain part of the crust in the form of mantle plume along the lithospheric fault to form the deposit.The ultimate power for tellurium mineralization was from H_(2)flow with high energy,which was produced through radiation from the melted iron of the Earth’s outer core.The H,flow results in the Earth’s degassing,as well as the mantle and crust’s uplift.展开更多
Letλ=(λ_(1),…,λ_(n))beβ-Jacobi ensembles with parameters p_(1),p_(2),n andβ,withβvarying with n.Set■.Suppose that■and 0≤σγ<1.We offer the large deviation for p_(1)+p_(2)/p_(1)■λ_(i)whenγ>0 via the...Letλ=(λ_(1),…,λ_(n))beβ-Jacobi ensembles with parameters p_(1),p_(2),n andβ,withβvarying with n.Set■.Suppose that■and 0≤σγ<1.We offer the large deviation for p_(1)+p_(2)/p_(1)■λ_(i)whenγ>0 via the large deviation of the corresponding empirical measure and via a direct estimate,respectively,whenγ=0.展开更多
The present paper deals with the eigenvalues of complex nonlocal Sturm-Liouville boundary value problems.The bounds of the real and imaginary parts of eigenvalues for the nonlocal Sturm-Liouville differential equation...The present paper deals with the eigenvalues of complex nonlocal Sturm-Liouville boundary value problems.The bounds of the real and imaginary parts of eigenvalues for the nonlocal Sturm-Liouville differential equation involving complex nonlocal potential terms associated with nonlocal boundary conditions are obtained in terms of the integrable conditions of coefficients and the real part of the eigenvalues.展开更多
In spite of its intrinsic complexities,the passive gait of bipedal robots on a sloping ramp is a subject of interest for numerous researchers.What distinguishes the present research from similar works is the considera...In spite of its intrinsic complexities,the passive gait of bipedal robots on a sloping ramp is a subject of interest for numerous researchers.What distinguishes the present research from similar works is the consideration of flexibility in the constituent links of this type of robotic systems.This is not a far-fetched assumption because in the transient(impact)phase,due to the impulsive forces which are applied to the system,the likelihood of exciting the vibration modes increases considerably.Moreover,the human leg bones that are involved in walking are supported by viscoelastic muscles and ligaments.Therefore,for achieving more exact results,it is essential to model the robot links with viscoelastic properties.To this end,the Gibbs-Appell formulation and Newton's kinematic impact law are used to derive the most general form of the system's dynamic equations in the swing and transient phases of motion.The most important issue in the passive walking motion of bipedal robots is the determination of the initial robot configuration with which the system could accomplish a periodic and stable gait solely under the effect of gravitational force.The extremely unstable nature of the system studied in this paper and the vibrations caused by the impulsive forces induced by the impact of robot feet with the inclined surface are some of the very serious challenges encountered for achieving the above-mentioned goal.To overcome such challenges,an innovative method that uses a combination of the linearized equations of motion in the swing phase and the algebraic motion equations in the transition phase is presented in this paper to obtain an eigenvalue problem.By solving this problem,the suitable initial conditions that are necessary for the passive gait of this bipedal robot on a sloping surface are determined.The effects of the characteristic parameters of elastic links including the modulus of elasticity and the Kelvin-Voigt coefficient on the walking stability of this type of robotic systems are also studied.The findings of this parametric study reveal that the increase in the Kelvin-Voigt coefficient enhances the stability of the robotic system,while the increase in the modulus of elasticity has an opposite effect.展开更多
This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solve...This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.展开更多
In this paper,we review computational approaches to optimization problems of inhomogeneous rods and plates.We consider both the optimization of eigenvalues and the localization of eigenfunctions.These problems are mot...In this paper,we review computational approaches to optimization problems of inhomogeneous rods and plates.We consider both the optimization of eigenvalues and the localization of eigenfunctions.These problems are motivated by physical problems including the determination of the extremum of the fundamental vibration frequency and the localization of the vibration displacement.We demonstrate how an iterative rearrangement approach and a gradient descent approach with projection can successfully solve these optimization problems under different boundary conditions with different densities given.展开更多
We consider the interior transmission eigenvalue problem corresponding to the scattering for an anisotropic medium of the scalar Helmholtz equation in the case where the boundary?Ωis split into two disjoint parts and...We consider the interior transmission eigenvalue problem corresponding to the scattering for an anisotropic medium of the scalar Helmholtz equation in the case where the boundary?Ωis split into two disjoint parts and possesses different transmission conditions.Using the variational method,we obtain the well posedness of the interior transmission problem,which plays an important role in the proof of the discreteness of eigenvalues.Then we achieve the existence of an infinite discrete set of transmission eigenvalues provided that n≡1,where a fourth order differential operator is applied.In the case of n■1,we show the discreteness of the transmission eigenvalues under restrictive assumptions by the analytic Fredholm theory and the T-coercive method.展开更多
In this paper we study the eigenvalue problem for integro-differential operators on a lasso graph.The trace formula of the operator is established by applying the residual technique in complex analysis.
This paper considers the finite difference(FD)approximations of diffusion operators and the boundary treatments for different boundary conditions.The proposed schemes have the compact form and could achieve arbitrary ...This paper considers the finite difference(FD)approximations of diffusion operators and the boundary treatments for different boundary conditions.The proposed schemes have the compact form and could achieve arbitrary even order of accuracy.The main idea is to make use of the lower order compact schemes recursively,so as to obtain the high order compact schemes formally.Moreover,the schemes can be implemented efficiently by solving a series of tridiagonal systems recursively or the fast Fourier transform(FFT).With mathematical induction,the eigenvalues of the proposed differencing operators are shown to be bounded away from zero,which indicates the positive definiteness of the operators.To obtain numerical boundary conditions for the high order schemes,the simplified inverse Lax-Wendroff(SILW)procedure is adopted and the stability analysis is performed by the Godunov-Ryabenkii method and the eigenvalue spectrum visualization method.Various numerical experiments are provided to demonstrate the effectiveness and robustness of our algorithms.展开更多
The ultimate strength of platings under compression is one of the most important factors to be addressed in the ship design.Current Rules for ship structural design generally provide explicit strength check criteria a...The ultimate strength of platings under compression is one of the most important factors to be addressed in the ship design.Current Rules for ship structural design generally provide explicit strength check criteria against buckling for simply supported and clamped platings.Nevertheless,ship platings generally exhibit an intermediate behaviour between the simple support and the clamped conditions,which implies that the torsional stiffness of supporting members should be duly considered.Hence,the main aim of this study is the development of new design formulas for the ultimate strength of platings under uniaxial compression,with short and/or long edges elastically restrained against torsion.In this respect,two benchmark studies are performed.The former is devoted to the development of new equations for the elastic buckling coefficients of platings with edges elastically restrained against torsion,based on the results of the eigenvalue buckling analysis,performed by Ansys Mechanical APDL.The latter investigates the ultimate strength of platings with elastically restrained edges,by systematically varying the plate slenderness ratio and the torsional stiffness of supporting members.Finally,the effectiveness of the new formulation is checked against a wide number of finite element(FE)simulations,to cover the entire design space of ship platings.展开更多
基金This work was supported by the National Natural Science Foundation of China(Grant No.62071248)the Natural Science Foundation of Nanjing University of Posts and Telecommunications(Grant No.NY223109)China Postdoctoral Science Foundation(Grant No.2022M721693).
文摘Quantum physics is primarily concerned with real eigenvalues,stemming from the unitarity of time evolutions.With the introduction of PT symmetry,a widely accepted consensus is that,even if the Hamiltonian of the system is not Hermitian,the eigenvalues can still be purely real under specific symmetry.Hence,great enthusiasm has been devoted to exploring the eigenvalue problem of non-Hermitian systems.In this work,from a distinct perspective,we demonstrate that real eigenvalues can also emerge under the appropriate recursive condition of eigenstates.Consequently,our findings provide another path to extract the real energy spectrum of non-Hermitian systems,which guarantees the conservation of probability and stimulates future experimental observations.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.52171251,U2106225,and 52231011)Dalian Science and Technology Innovation Fund (Grant No.2022JJ12GX036)。
文摘A numerical method is proposed to calculate the eigenvalues of the Zakharov–Shabat system based on Chebyshev polynomials. A mapping in the form of tanh(ax) is constructed according to the asymptotic of the potential function for the Zakharov–Shabat eigenvalue problem. The mapping can distribute Chebyshev nodes very well considering the gradient for the potential function. Using Chebyshev polynomials, tanh(ax) mapping, and Chebyshev nodes, the Zakharov–Shabat eigenvalue problem is transformed into a matrix eigenvalue problem. This method has good convergence for the Satsuma–Yajima potential and the convergence rate is faster than the Fourier collocation method. This method is not only suitable for simple potential functions but also converges quickly for a complex Y-shape potential. It can also be further extended to other linear eigenvalue problems.
基金supported by the National Natural Science Foundation of China(No.U1839209).
文摘Wave propagation in horizontally layered media is a classical problem in seismic-wave theory.In semi-infinite space,a nondispersive Rayleigh wave mode exists,and the eigendisplacement decays exponentially with depth.In a layered model with increasing layer velocity,the phase velocity of the Rayleigh wave varies between the S-wave velocity of the bottom half-space and that of the classical Rayleigh wave propagated in a supposed half-space formed by the parameters of the top layer.If the phase velocity is the same as the P-or S-wave velocity of the layer,which is called the critical mode or critical phase velocity of surface waves,the general solution of the wave equation is not a homogeneous(expressed by trigonometric functions)or inhomogeneous(expressed by exponential functions)plane wave,but one whose amplitude changes linearly with depth(expressed by a linear function).Theories based on a general solution containing only trigonometric or exponential functions do not apply to the critical mode,owing to the singularity at the critical phase velocity.In this study,based on the classical framework of generalized reflection and transmission coefficients,the propagation of surface waves in horizontally layered media was studied by introducing a solution for the linear function at the critical phase velocity.Therefore,the eigenvalues and eigenfunctions of the critical mode can be calculated by solving a singular problem.The eigendisplacement characteristics associated with the critical phase velocity were investigated for different layered models.In contrast to the normal mode,the eigendisplacement associated with the critical phase velocity exhibits different characteristics.If the phase velocity is equal to the S-wave velocity in the bottom half-space,the eigendisplacement remains constant with increasing depth.
基金supported by Orient Resources Ltd.College of Earth Sciences,Jilin University。
文摘The Dashuigou tellurium deposit is the world’s only known independent tellurium deposit.By restoring metamorphic rocks’protolith,we seek to understand not only the development and evolution trajectory of the region but also the origin of the relevant deposits.While there are many ways to restore metamorphic rocks’protolith,we take the host metamorphic rocks of Dashuigou tellurium deposit and leverage various petrochemical eigenvalues and related diagrams previously proposed to reveal the deposit’s host metamorphic rocks’protolith.The petrochemical eigenvalues include molecular number,Niggli’s value,REE parity ratio,CaO/Al_(2)O_(3)ratio,Fe^(3+) /(Fe^(3+) -+Fe^(2+) )ratio,chondrite-normalized REE value,logarithmic REE value,various REE eigenvalues including scandium,Eu/Sm ratio,total REE amount,light and heavy REEs,δEu,Eu anomaly,Sm/Nd ratio,and silicon isotope δ^(30) SiNBS-29‰,etc.The petrochemical plots include ACMs,100 mg-c-(al+alk),SiO_(2)-(Na_(2)O+K_(2)O),(al+fm)-(c+alk)versus Si,FeO+Fe_(2)O^(3+) TiO)-Al_(2)O_(3)-MgO,c-mg,Al_(2)O_(3)-(Na_(2)O+K_(2)O),chondrite-normalized REE model,La/Yb-REE,and Sm/Nd ratio,etc.On the basis of these comprehensive analyses,the following conclusions are drawn,starting from the many mantle-derived types of basalt developed in the study area of different geological ages,combined with the previously published research results on the deposit s fluid inclusions and sulfur and lead isotopes.The deposit is formed by mantle degassing in the form of a mantle plume in the late Yanshanian orogeny.The degassed fluids are rich in nano-sc ale substances including Fe,Te,S,As,Bi,Au,Se,H_(2),CO_(2),N_(2),H_(2)O,and CH_(4),which are enriched by nano-effect,and then rise to a certain part of the crust in the form of mantle plume along the lithospheric fault to form the deposit.The ultimate power for tellurium mineralization was from H_(2)flow with high energy,which was produced through radiation from the melted iron of the Earth’s outer core.The H,flow results in the Earth’s degassing,as well as the mantle and crust’s uplift.
基金supported by the NSFC (12171038,11871008)985 Projects.
文摘Letλ=(λ_(1),…,λ_(n))beβ-Jacobi ensembles with parameters p_(1),p_(2),n andβ,withβvarying with n.Set■.Suppose that■and 0≤σγ<1.We offer the large deviation for p_(1)+p_(2)/p_(1)■λ_(i)whenγ>0 via the large deviation of the corresponding empirical measure and via a direct estimate,respectively,whenγ=0.
基金Supported by the National Nature Science Foundation of China(12101356,12101357,12071254,11771253)the National Science Foundation of Shandong Province(ZR2021QA065,ZR2020QA009,ZR2021MA047)the China Postdoctoral Science Foundation(2019M662313)。
文摘The present paper deals with the eigenvalues of complex nonlocal Sturm-Liouville boundary value problems.The bounds of the real and imaginary parts of eigenvalues for the nonlocal Sturm-Liouville differential equation involving complex nonlocal potential terms associated with nonlocal boundary conditions are obtained in terms of the integrable conditions of coefficients and the real part of the eigenvalues.
文摘In spite of its intrinsic complexities,the passive gait of bipedal robots on a sloping ramp is a subject of interest for numerous researchers.What distinguishes the present research from similar works is the consideration of flexibility in the constituent links of this type of robotic systems.This is not a far-fetched assumption because in the transient(impact)phase,due to the impulsive forces which are applied to the system,the likelihood of exciting the vibration modes increases considerably.Moreover,the human leg bones that are involved in walking are supported by viscoelastic muscles and ligaments.Therefore,for achieving more exact results,it is essential to model the robot links with viscoelastic properties.To this end,the Gibbs-Appell formulation and Newton's kinematic impact law are used to derive the most general form of the system's dynamic equations in the swing and transient phases of motion.The most important issue in the passive walking motion of bipedal robots is the determination of the initial robot configuration with which the system could accomplish a periodic and stable gait solely under the effect of gravitational force.The extremely unstable nature of the system studied in this paper and the vibrations caused by the impulsive forces induced by the impact of robot feet with the inclined surface are some of the very serious challenges encountered for achieving the above-mentioned goal.To overcome such challenges,an innovative method that uses a combination of the linearized equations of motion in the swing phase and the algebraic motion equations in the transition phase is presented in this paper to obtain an eigenvalue problem.By solving this problem,the suitable initial conditions that are necessary for the passive gait of this bipedal robot on a sloping surface are determined.The effects of the characteristic parameters of elastic links including the modulus of elasticity and the Kelvin-Voigt coefficient on the walking stability of this type of robotic systems are also studied.The findings of this parametric study reveal that the increase in the Kelvin-Voigt coefficient enhances the stability of the robotic system,while the increase in the modulus of elasticity has an opposite effect.
基金the National Science and Tech-nology Council,Taiwan for their financial support(Grant Number NSTC 111-2221-E-019-048).
文摘This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.
基金supported by the DMS-1853701supported in part by the DMS-2208373.
文摘In this paper,we review computational approaches to optimization problems of inhomogeneous rods and plates.We consider both the optimization of eigenvalues and the localization of eigenfunctions.These problems are motivated by physical problems including the determination of the extremum of the fundamental vibration frequency and the localization of the vibration displacement.We demonstrate how an iterative rearrangement approach and a gradient descent approach with projection can successfully solve these optimization problems under different boundary conditions with different densities given.
基金supported by the National Natural Science Foundation of China(11571132,12301542)the Natural Science Foundation of Hubei(2022CFB725)the Natural Science Foundation of Yichang(A23-2-027)。
文摘We consider the interior transmission eigenvalue problem corresponding to the scattering for an anisotropic medium of the scalar Helmholtz equation in the case where the boundary?Ωis split into two disjoint parts and possesses different transmission conditions.Using the variational method,we obtain the well posedness of the interior transmission problem,which plays an important role in the proof of the discreteness of eigenvalues.Then we achieve the existence of an infinite discrete set of transmission eigenvalues provided that n≡1,where a fourth order differential operator is applied.In the case of n■1,we show the discreteness of the transmission eigenvalues under restrictive assumptions by the analytic Fredholm theory and the T-coercive method.
基金Supported by the National Natural Science Foundation of China(11871031)the National Natural Science Foundation of Jiang Su(BK20201303).
文摘In this paper we study the eigenvalue problem for integro-differential operators on a lasso graph.The trace formula of the operator is established by applying the residual technique in complex analysis.
基金supported by the NSFC grant 11801143J.Lu’s research is partially supported by the NSFC grant 11901213+3 种基金the National Key Research and Development Program of China grant 2021YFA1002900supported by the NSFC grant 11801140,12171177the Young Elite Scientists Sponsorship Program by Henan Association for Science and Technology of China grant 2022HYTP0009the Program for Young Key Teacher of Henan Province of China grant 2021GGJS067.
文摘This paper considers the finite difference(FD)approximations of diffusion operators and the boundary treatments for different boundary conditions.The proposed schemes have the compact form and could achieve arbitrary even order of accuracy.The main idea is to make use of the lower order compact schemes recursively,so as to obtain the high order compact schemes formally.Moreover,the schemes can be implemented efficiently by solving a series of tridiagonal systems recursively or the fast Fourier transform(FFT).With mathematical induction,the eigenvalues of the proposed differencing operators are shown to be bounded away from zero,which indicates the positive definiteness of the operators.To obtain numerical boundary conditions for the high order schemes,the simplified inverse Lax-Wendroff(SILW)procedure is adopted and the stability analysis is performed by the Godunov-Ryabenkii method and the eigenvalue spectrum visualization method.Various numerical experiments are provided to demonstrate the effectiveness and robustness of our algorithms.
文摘The ultimate strength of platings under compression is one of the most important factors to be addressed in the ship design.Current Rules for ship structural design generally provide explicit strength check criteria against buckling for simply supported and clamped platings.Nevertheless,ship platings generally exhibit an intermediate behaviour between the simple support and the clamped conditions,which implies that the torsional stiffness of supporting members should be duly considered.Hence,the main aim of this study is the development of new design formulas for the ultimate strength of platings under uniaxial compression,with short and/or long edges elastically restrained against torsion.In this respect,two benchmark studies are performed.The former is devoted to the development of new equations for the elastic buckling coefficients of platings with edges elastically restrained against torsion,based on the results of the eigenvalue buckling analysis,performed by Ansys Mechanical APDL.The latter investigates the ultimate strength of platings with elastically restrained edges,by systematically varying the plate slenderness ratio and the torsional stiffness of supporting members.Finally,the effectiveness of the new formulation is checked against a wide number of finite element(FE)simulations,to cover the entire design space of ship platings.