Transfer matrix method for free and forced vibrations of multi-level functionally graded material stepped beams with different boundary conditions

Functionally graded materials(FGMs)are a novel class of composite materials that have attracted significant attention in the field of engineering due to their unique mechanical properties.This study aims to explore the dynamic behaviors of an FGM stepped beam with different boundary conditions based on an efficient solving method.Under the assumptions of the Euler-Bernoulli beam theory,the governing differential equations of an individual FGM beam are derived with Hamilton’s principle and decoupled via the separation-of-variable approach.Then,the free and forced vibrations of the FGM stepped beam are solved with the transfer matrix method(TMM).Two models,i.e.,a three-level FGM stepped beam and a five-level FGM stepped beam,are considered,and their natural frequencies and mode shapes are presented.To demonstrate the validity of the method in this paper,the simulation results by ABAQUS are also given.On this basis,the detailed parametric analyses on the frequencies and dynamic responses of the three-level FGM stepped beam are carried out.The results show the accuracy and efficiency of the TMM.

Corrigendum to“Atomic-scale electromagnetic theory bridging optics in microscopic world and macroscopic world”

The signs of the electric field markers in Figs.2 and 4 of the paper[Chin.Phys.B 32104211(2023)]have been corrected.These modifications do not affect the results derived in the paper.

Analysis profile of the fully grouted rock bolt in jointed rock using analytical and numerical methods

The purpose of this study was to investigate the effect of bolt profile on load transfer mechanism of fully grouted bolts in jointed rocks using analytical and numerical methods. Based on the analytical method with development of methods, a new model is presented. To validate the analytical model, five different profiles modeled by ANSYS software. The profile of rock bolts T3 and T4with load transfer capacity,respectively 180 and 195 kN in the jointed rocks was selected as the optimum profiles. Finally, the selected profiles were examined in Tabas Coal Mine. FLAC analysis indicates that patterns 6+7 with2 NO flexi bolt 4 m better than other patterns within the faulted zone.

Activated dissociation of H_(2) on the Cu(001)surface:The role of quantum tunneling

The activation and dissociation of hydrogen molecules(H_(2))on the Cu(001)surface are studied theoretically.Using first-principles calculations,the activation barrier for the dissociation of H_(2) on Cu(001)is determined to be~0.59 eV in height.It is found that the electron transfer from the copper substrate to H_(2) plays a key role in the activation and breaking of the H–H bond,and the formation of the Cu–H bonds.Two stationary states are identified at around the critical height of bond breaking,corresponding to the molecular and the dissociative states,respectively.Using the transfer matrix method,we also investigate the role of quantum tunneling in the dissociation process along the minimum energy pathway(MEP),which is found to be significant at or below room temperature.At a given temperature,the tunneling contributions due to the translational and the vibrational motions of H_(2) are quantified for the dissociation process.Within a wide range of temperature,the effects of quantum tunneling on the effective barriers of dissociation and the rate constants are observed.The deduced energetic parameters associated with the thermal equilibrium and non-equilibrium(molecular beam)conditions are comparable to experimental data.In the low-temperature region,the crossover from classical to quantum regime is identified.

Nonreciprocal wide-angle bidirectional absorber based on one-dimensional magnetized gyromagnetic photonic crystals

We propose magnetized gyromagnetic photonic crystals(MGPCs)composed of indium antimonide(InSb)and yttrium iron garnet ferrite(YIGF)layers,which possess the properties of nonreciprocal wide-angle bidirectional absorption.Periodical defects in the MGPCs work as filters.Absorption bands(ABs)for the positive and negative propagations arise from the optical Tamm state and resonance in cavities respectively,and they prove to share no overlaps in the studied frequency range.Givenω=2.0138 THz,for the positive propagation,the ABs in the high-frequency range are localized in the interval between 0.66ωand 0.88ω.In the angular range,the ABs for the TE and TM waves reach 60°and 51°,separately.For the negative propagation,the ABs in the low-frequency range are localized in the interval between 0.13ωand 0.3ω.The AB s extend to 60°for the TE waves and 80.4°for the TM waves.There also exists a narrow frequency band in a lower frequency range.The relevant factors,which include the external temperature,the magnetic fields applied to the YIGF,the refractive index of the impedance matching layer,and the defect thickness,are adjusted to investigate the effects on the ABs.All the numerical simulations are based on the transfer matrix method.This work provides an approach to designs of isolators and so on.

Atomic-scale electromagnetic theory bridging optics in microscopic world and macroscopic world

Atoms in the microscopic world are the basic building blocks of the macroscopic world. In this work, we construct an atomic-scale electromagnetic theory that bridges optics in the microscopic and macroscopic worlds. As the building block of the theory, we use the microscopic polarizability to describe the optical response of a single atom, solve the transport of electromagnetic wave through a single atomic layer under arbitrary incident angle and polarization of the light beam, construct the single atomic layer transfer matrix for light transport across the atomic layer. Based on this transfer matrix, we get the analytical form of the dispersion relation, refractive index, and transmission/reflection coefficient of the macroscopic medium. The developed theory can handle single-layer and few-layers of homogeneous and heterogeneous 2D materials, investigate homogeneous 2D materials with various vacancies or insertion atomic-layer defects, study compound 2D materials with a unit cell composed of several elements in both the lateral and parallel directions with respect to the light transport.

Local density of optical states calculated by the mode spectrum in stratified media

The local density of optical states(LDOS)is an important physical concept,which can characterize the spontaneous emission of microcavities.In order to calculate the LDOS,the relationship between the mode spectrum and the LDOS is established.Then,based on the transfer matrix method and the effective resonator model,the leaky loss of the leaky mode and the mode spectrum in the one-dimensional photonic bandgap crystal waveguide are calculated,results of which indicate that the mode spectrum can characterize the leaky loss of the leaky mode.At last,the density of optical states(DOS),and the LDOS in each layer are calculated.The partial DOS and the partial LDOS in the quantum well,related to the fundamental leaky mode,can be used to find out the optimal location of the quantum well in the defect layer to couple more useful photons into the lasing mode for lasers.

Antireflection Coating for Solar Cells Based on Graded-Index Materials

The conversion of sunlight into electricity via photovoltaics presents tremendous opportunities for the generation of renewable energy. However, solar cells still face several challenges and limitations to further reduce manufacturing costs and increase module efficiency. Photon management is paramount to increase the efficiency of the mainstream silicon-based cell and always includes a suitable antireflection coating (ARC) structure to decrease the reflectance (R) at the top surface. We propose a novel triple-layer anti-reflective coating (TLAR) consisting of three layers sandwiched between the upper cover (glass) and the substrate (silicon). The inner three layers are graded refractive index material (GIM) as an active layer, titanium dioxide (TiO2), and zinc sulfide (ZnS), respectively. The optical properties of the TLAR have been investigated using the transfer matrix method (TMM). The results of using GIM as the active medium lead to the reflection decaying to the minimum value, and the transmittance reaching the maximum values at a specific wavelength range. The proposed triple-layer anti-reflective coating (TLAR) structure presents a promising solution for enhancing the efficiency of solar cells. Its unique design and utilization of graded refractive index material (GIM) as the active layer make it a novel and innovative approach that holds great potential for advancing solar cell technology.

Estimation of the Forces on Rubbing Rotor-Bearing System

It*ss very important to identify the load on rotor bearing system. Based on load identification in time domain, a new method combined with transfer matrix method is applied in rubbing rotor bearing system to estimate the force on rubbing rotor for the first time. Avoiding calculating the modal parameters, the outside excitation force on multi point and the distribute of internal forces on rubbing rotor bearing system can be identified by this way at the same time. Finally, the simulation is made on partial lateral rubbing rotor bearing system. For the faulty of rotor/stator rub, the point where rubbing occurs and the change of impact forces can be detected by this method in order to improve the accuracy of rotor/stator rub faulty diagnosis.

FREE VIBRATION ANALYSIS AND PHYSICAL PARAMETER IDENTIFICATION OF NON-UNIFORM BEAM CARRYING SPRING-MASS SYSTEMS

To analyze a multibody system composed of non-uniform beam and spring-mass subsystems, the model discretization is carried on by utilizing the finite element method（FEM）, the dynamic model of non-uniform beam is developed by using the transfer matrix method of multibody system（MS-TMM）, the transfer matrix of non-u- niform beam is derived, and the natural frequencies are computed. Compared with the numerical assembly method （NAM）, the results by MS-TMM have good agreement with the results by FEM, and are better than the results by NAM. When using the high precision method, the global dynamic equations of the complex multibody system are not needed and the orders of involved system matrices are decreased greatly. For the investigation on the re- verse problem of the physical parameter identification of multibody system, MS-TMM and the optimization tech- nology based on genetic algorithms（GAs） are combined and extended. The identification problem is exchanged for an optimization problem, and it is formulated as a global minimum solution of the objective function with respect to natural frequencies of multibody system. At last, the numerical example of non-uniform beam with attach- ments is discussed, and the identification results indicate the feasibility and the effectivity of the proposed aop- proach.

Waveguide Optimization for a 9.0μm GaAs-Based Quantum Cascade Laser

Improved waveguide designs for 9.0μm GaAs-based quantum cascade laser （QCL） structures are presented. Modal losses and confinement factors are calculated for TM modes with the transfer matrix method （TMM） and effective index method （EIM）. The thicknesses of the cladding layer and waveguide layer, the ridge-width, and the cavity length are all taken into account. Appropriate thicknesses of epilayers are given with lower threshold gain and more economical material growth time.

Analysis of All-Carbon Brick Bottom and Ceramic Cup Synthetic Hearth Bottom

One of the bottlenecks of the blast furnace （BF） campaign is the life length of hearth bottom. The basic reason for the erosion of hearth bottom is its direct contact with hot metal. According to the theory of heat transfer, models of BF hearth bottom are built based on the actual examples using software and VC language, and the calculated results are in good agreement with the data of BF dissection after blowing out. The temperature distribution and the capability of the resistance to erosion for different structures of hearth bottom are analyzed, especially the two prevalent kinds of hearth bottom arrangements called ＂the method of heat transfer＂ for all-carbon brick bottom and ＂the method of heat isolation＂ for ceramic synthetic hearth bottom. Features of the two kinds of hearth bottoms are analyzed. Also the different ways of protecting the hearth bottom are clarified, according to some actual examples. After that, the same essence of prolonging life, and the fact that the existence of a ＂protective skull＂ with low thermal conductivity between the hot metal and brick layers is of utmost importance are shown.

Carbon matrix effects on the micro-structure and performance of Pt nanowire cathode prepared by decal transfer method

High performance cathode for polymer electrolyte membrane fuel cell was prepared by depositing Pt nanowires in a carbon matrix coated on a substrate, and using decal transfer method to fabricate the membrane electrode assembly. The effects of carbon and ionomer contents on the electrode micro-structure and fuel cell performance are investigated by physical characterization and single cell testing. The Pt nanowires are gradient distributed across the cathode thickness, and more Pt exists near the membrane. Both the carbon and ionomer contents can affect the Pt nanowires distribution and aggregation. In addition, the carbon loading dominates the transport distance of gas and proton, and the ionomer content affects the triple phase boundaries and porosity in the cathode. The optimal structure of Pt nanowire cathode is obtained at 0.10 mg·cm^-2 carbon loading and 10 wt% ionomer.

A SEMI-ANALYTICAL AND SEMI-NUMERICAL METHOD FOR SOLVING 2-D SOUND-STRUCTURE INTERACTION PROBLEMS

Based on the transfer matrix method and the virtual source simulation technique, this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound- structure interaction problems under a harmonic excitation.Within any integration segment, as long as its length is small enough,along the circumferential curvilinear coordinate,the non- homogeneous matrix differential equation of an elastic ring of complex geometrical shape can be rewritten in terms of the homogeneous one by the method of extended homogeneous capacity proposed in this paper.For the exterior fluid domain,the multi-circular virtual source simulation technique is adopted.The source density distributed on each virtual circular curve may be ex- panded as the Fourier's series.Combining with the inverse fast Fourier transformation,a higher accuracy and efficiency method for solving 2-D exterior Helmholtz's problems is presented in this paper.In the aspect of solution to the coupling equations,the state vectors of elastic ring induced by the given harmonic excitation and generalized forces of coefficients of the Fourier series can be obtained respectively by using a high precision integration scheme combined with the method of extended homogeneous capacity put forward in this paper.According to the superposition princi- ple and compatibility conditions at the interface between the elastic ring and fluid,the algebraic equation of system can be directly constructed by using the least square approximation.Examples of acoustic radiation from two typical fluid-loaded elastic rings under a harmonic concentrated force are presented.Numerical results show that the method proposed is more efficient than the mixed FE-BE method in common use.

Discrete time transfer matrix method for dynamics of multibody system with flexible beams moving in space

In this paper, by defining new state vectors and developing new transfer matrices of various elements mov- ing in space, the discrete time transfer matrix method of multi-rigid-flexible-body system is expanded to study the dynamics of multibody system with flexible beams moving in space. Formulations and numerical example of a rigid- flexible-body three pendulums system moving in space are given to validate the method. Using the new method to study the dynamics of multi-rigid-flexible-body system mov- ing in space, the global dynamics equations of system are not needed, the orders of involved matrices of the system are very low and the computational speed is high, irrespec- tive of the size of the system. The new method is simple, straightforward, practical, and provides a powerful tool for multi-rigid-flexible-body system dynamics.

Transfer matrix method for determination of the natural vibration characteristics of elastically coupled launch vehicle boosters

The analysis of natural vibration characteristics has become one of important steps of the manufacture and dynamic design in the aerospace industry. This paper presents a new scenario called virtual cutting in the context of the transfer matrix method of linear multibody systems closed- loop topology for computing the free vibration characteristics of elastically coupled flexible launch vehicle boosters. In this approach, the coupled system is idealized as a triple-beam system-like structure coupled by linear translational springs, where a non-uniform free-free Euler-Bemoulli beam is used. A large thrust-to-weight ratio leads to large axial accelera- tions that result in an axial inertia load distribution from nose to tail. Consequently, it causes the development of significant compressive forces along the length of the launch vehicle. Therefore, it is important to take into account this effect in the transverse vibration model. This scenario does not need the global dynamics equations of a system, and it has high computational efficiency and low memory requirements. The validity of the presented scenario is achieved through com- parison to other approaches published in the literature.

An analytical model for evaluating the dynamic response of a tunnel embedded in layered foundation soil with different saturations

This paper proposes an analytical model for evaluating the dynamic response of an underground railway tunnel in layered foundation soil with different saturations.The soil is modeled as layered media,and the circular tunnel lining is modeled as an infinite Flügge cylindrical shell.The separation of variables method is used to solve the motion equation of the shell,and the wave equation of the soil is solved using the Helmholtz decomposition theorem.A dynamic matrix reflecting the wave vectors of soil layers is established using the transfer matrix method.Based on boundary conditions,the tunnel-soil model is coupled using the transformation method of plane wave functions and cylindrical wave functions.The proposed model is validated by comparison with existing tunnel models,and the effects of saturation and the layered properties of soil on the dynamic response of a layered tunnel-soil system is demonstrated via case studies.

A NOVEL SEMI-ANALYTICAL METHOD FOR SOLVING ACOUSTIC RADIATION FROM LONGITUDINALLY STIFFENED INFINITE NON-CIRCULAR CYLINDRICAL SHELLS IN WATER

Based on the extended homogeneous capacity high precision integration method and the spectrum method of virtual boundary with a complex radius vector, a novel semi-analytical method, which has satisfactory computation efectiveness and precision, is presented for solving the acoustic radiation from a submerged infnite non-circular cylindrical shell stifened by longitudinal ribs by means of the Fourier integral transformation and stationary phase method. In this work, besides the normal interacting force, which is commonly adopted by some researchers, the other interacting forces and moments between the longitudinal ribs and the non-circular cylindrical shell are considered at the same time. The efects of the number and the size of the cross-section of longitudinal ribs on the characteristics of acoustic radiation are investigated. Numerical results show that the method proposed is more efcient than the existing mixed FE-BE method.

Dynamics Analysis of a Parallel Mill-turn Tool Spindle Head Driven by Dual-linear Motors Using Extended Transfer Matrix Method

The hybrid dynamics of multi-rigid-body and multi-flexible-body system becomes the mainstream of multi-body dynamics.Currently there lacks a compact approach to model the hybrid dynamics,especially in modern machine tool application,due to the difficulty of solving the hybrid equations or the limitation of current software when dealing with the hybrid dynamics.The extended transfer matrix method（E-TMM）,which extends elements in three-dimensional space with higher matrixes,is proposed to simplify the modeling process of the hybrid dynamics.The E-TMM modeling approaches of 3 basic elements including 3D vibrant rigid body,joint and flexible body are studied in details.A parallel mill-turn tool spindle head unit driven by dual-linear motors is chosen as a plant to demonstrate the E-TMM modeling process.By using E-TMM,the spindle head unit is simplified as a topological network consisting of the three types of element,i.e.,3D vibrant rigid body,joint and flexible body,including 11 rigid bodies,14 joints and 1 3D-Timoshenko beam.Then the dynamic model of the system can be easily obtained by deducing the element-network by means of state vector transformation.The dynamic characteristics of the spindle head,such as natural frequencies,dynamic flexibility,etc.can be predicted by solving the obtained model.Experiment verification indicates that the E-TMM is valid with enough accuracy in the dynamic analysis of the parallel mill-turn tool spindle head.The E-TMM is capable of modeling the dynamics of machine tool structure with no requirements of deducing and solving the sophisticated differential equations.Moreover,the E-TMM provides a simple and elegant tool for hybrid dynamic analysis in future dynamic design of machine tools.

An efficient and modular modeling for launch dynamics of tubed rockets on a moving launcher

This paper develops a modular modeling and efficient formulation of launch dynamics with marching fire(LDMF)using a mixed formulation of the transfer matrix method for multibody systems(MSTMM)and Newton-Euler formulation.Taking a ground-borne multiple launch rocket systems(MLRS),the focus is on the launching subsystem comprising the rocket,flexible tube,and tube tail.The launching subsystem is treated as a coupled rigid-flexible multibody system,where the rocket and tube tail are treated as rigid bodies while the flexible tube as a beam with large motion.Firstly,the tube and tube tail can be elegantly handled by the MSTMM,a computationally efficient order-N formulation.Then,the equation of motion of the in-bore rocket with relative kinematics w.r.t.the tube using the Newton-Euler method is derived.Finally,the rocket,tube,and tube tail dynamics are coupled,yielding the equation of motion of the launching subsystem that can be regarded as a building block and further integrated with other subsystems.The deduced dynamics equation of the launching subsystem is not limited to ground-borne MLRS but also fits for tanks,self-propelled artilleries,and other air-borne and naval-borne weapons undergoing large motion.Numerical simulation results of LDMF are given and partially verified by the experiment.