Dynamic Characteristics of Functionally Graded Timoshenko Beams by Improved Differential Quadrature Method

This study proposes an effective method to enhance the accuracy of the Differential Quadrature Method(DQM)for calculating the dynamic characteristics of functionally graded beams by improving the form of discrete node distribution.Firstly,based on the first-order shear deformation theory,the governing equation of free vibration of a functionally graded beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam axial displacement,transverse displacement,and cross-sectional rotation angle by considering the effects of shear deformation and rotational inertia of the beam cross-section.Then,ignoring the shear deformation of the beam section and only considering the effect of the rotational inertia of the section,the governing equation of the beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam transverse displacement.Based on the differential quadrature method theory,the eigenvalue problem of ordinary differential equations is transformed into the eigenvalue problem of standard generalized algebraic equations.Finally,the first several natural frequencies of the beam can be calculated.The feasibility and accuracy of the improved DQM are verified using the finite element method(FEM)and combined with the results of relevant literature.

A Hybrid Level Set Optimization Design Method of Functionally Graded Cellular Structures Considering Connectivity

With the continuous advancement in topology optimization and additive manufacturing(AM)technology,the capability to fabricate functionally graded materials and intricate cellular structures with spatially varying microstructures has grown significantly.However,a critical challenge is encountered in the design of these structures–the absence of robust interface connections between adjacent microstructures,potentially resulting in diminished efficiency or macroscopic failure.A Hybrid Level Set Method(HLSM)is proposed,specifically designed to enhance connectivity among non-uniform microstructures,contributing to the design of functionally graded cellular structures.The HLSM introduces a pioneering algorithm for effectively blending heterogeneous microstructure interfaces.Initially,an interpolation algorithm is presented to construct transition microstructures seamlessly connected on both sides.Subsequently,the algorithm enables the morphing of non-uniform unit cells to seamlessly adapt to interconnected adjacent microstructures.The method,seamlessly integrated into a multi-scale topology optimization framework using the level set method,exhibits its efficacy through numerical examples,showcasing its prowess in optimizing 2D and 3D functionally graded materials(FGM)and multi-scale topology optimization.In essence,the pressing issue of interface connections in complex structure design is not only addressed but also a robust methodology is introduced,substantiated by numerical evidence,advancing optimization capabilities in the realm of functionally graded materials and cellular structures.

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.

A Radial Basis Function Method with Improved Accuracy for Fourth Order Boundary Value Problems

Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.

AHermitian C^(2) Differential Reproducing Kernel Interpolation Meshless Method for the 3D Microstructure-Dependent Static Flexural Analysis of Simply Supported and Functionally Graded Microplates

This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredependent static flexural behavior of a functionally graded(FG)microplate subjected to mechanical loads and placed under full simple supports.In the formulation,we select the transverse stress and displacement components and their first-and second-order derivatives as primary variables.Then,we set up the differential reproducing conditions(DRCs)to obtain the shape functions of the Hermitian C^(2) differential reproducing kernel(DRK)interpolant’s derivatives without using direct differentiation.The interpolant’s shape function is combined with a primitive function that possesses Kronecker delta properties and an enrichment function that constituents DRCs.As a result,the primary variables and their first-and second-order derivatives satisfy the nodal interpolation properties.Subsequently,incorporating ourHermitianC^(2)DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C^(2) DRKIM method is confirmed to be accurate and fast in its convergence rate by comparing the solutions it produces with the relevant 3D solutions available in the literature.Finally,the impact of essential factors on the transverse stresses,in-plane stresses,displacements,and couple stresses that are induced in the loaded microplate is examined.These factors include the length-to-thickness ratio,the material length-scale parameter,and the inhomogeneity index,which appear to be significant.

A comparison study on structure-function relationship of polysaccharides obtained from sea buckthorn berries using different methods:antioxidant and bile acid-binding capacity

In this study,the structural characters,antioxidant activities and bile acid-binding ability of sea buckthorn polysaccharides(HRPs)obtained by the commonly used hot water(HRP-W),pressurized hot water(HRP-H),ultrasonic(HRP-U),acid(HRP-C)and alkali(HRP-A)assisted extraction methods were investigated.The results demonstrated that extraction methods had significant effects on extraction yield,monosaccharide composition,molecular weight,particle size,triple-helical structure,and surface morphology of HRPs except for the major linkage bands.Thermogravimetric analysis showed that HRP-U with filamentous reticular microstructure exhibited better thermal stability.The HRP-A with the lowest molecular weight and highest arabinose content possessed the best antioxidant activities.Moreover,the rheological analysis indicated that HRPs with higher galacturonic acid content and molecular weight showed higher viscosity and stronger crosslinking network(HRP-C,HRP-W and HRP-U),which exhibited stronger bile acid binding capacity.The present findings provide scientific evidence in the preparation technology of sea buckthorn polysaccharides with good antioxidant and bile acid binding capacity which are related to the structure affected by the extraction methods.

Random Green's Function Method for Large-Scale Electronic Structure Calculation

We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projecting onto Krylov subspace. With the rGF method, the Fermi–Dirac operator can be obtained directly, avoiding the polynomial expansion to Fermi–Dirac function. To demonstrate the applicability, we implement the rGF method with the density-functional tight-binding method. It is shown that the Krylov subspace can maintain at small size for materials with different gaps at zero temperature, including H_(2)O and Si clusters. We find with a simple deflation technique that the rGF self-consistent calculation of H_(2)O clusters at T = 0 K can reach an error of~ 1 me V per H_(2)O molecule in total energy, compared to deterministic calculations. The rGF method provides an effective stochastic method for large-scale electronic structure simulation.

Extraction of Acoustic Normal Mode Depth Functions Using Range-Difference Method with Vertical Linear Array Data

Data-derived normal mode extraction is an effective method for extracting normal mode depth functions in the absence of marine environmental data.However,when the corresponding singular vectors become nonunique when two or more singular values obtained from the cross-spectral density matrix diagonalization are nearly equal,this results in unsatisfactory extraction outcomes for the normal mode depth functions.To address this issue,we introduced in this paper a range-difference singular value decomposition method for the extraction of normal mode depth functions.We performed the mode extraction by conducting singular value decomposition on the individual frequency components of the signal's cross-spectral density matrix.This was achieved by using pressure and its range-difference matrices constructed from vertical line array data.The proposed method was validated using simulated data.In addition,modes were successfully extracted from ambient noise.

Implementation of the Integrated Green’s Function Method for 3D Poisson’s Equation in a Large Aspect Ratio Computational Domain

The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, the integrated Green’s function method has been adopted to solve the 3D Poisson equation subject to open boundary conditions. In this paper, we report on the efficient implementation of this method, which can save more than a factor of 50 computing time compared with the direct brute force implementation and its improvement under certain extreme conditions.

A non-probabilistic reliability topology optimization method based on aggregation function and matrix multiplication considering buckling response constraints

A non-probabilistic reliability topology optimization method is proposed based on the aggregation function and matrix multiplication.The expression of the geometric stiffness matrix is derived,the finite element linear buckling analysis is conducted,and the sensitivity solution of the linear buckling factor is achieved.For a specific problem in linear buckling topology optimization,a Heaviside projection function based on the exponential smooth growth is developed to eliminate the gray cells.The aggregation function method is used to consider the high-order eigenvalues,so as to obtain continuous sensitivity information and refined structural design.With cyclic matrix programming,a fast topology optimization method that can be used to efficiently obtain the unit assembly and sensitivity solution is conducted.To maximize the buckling load,under the constraint of the given buckling load,two types of topological optimization columns are constructed.The variable density method is used to achieve the topology optimization solution along with the moving asymptote optimization algorithm.The vertex method and the matching point method are used to carry out an uncertainty propagation analysis,and the non-probability reliability topology optimization method considering buckling responses is developed based on the transformation of non-probability reliability indices based on the characteristic distance.Finally,the differences in the structural topology optimization under different reliability degrees are illustrated by examples.

Impacts of Aggregation Methods and Trophospecies Number on the Structure and Function of Marine Food Webs

Aggregation of species with similar ecological properties is one of the effective methods to simplify food web researches.However,species aggregation will affect not only the complexity of modeling process but also the accuracy of models’outputs.Selection of aggregation methods and the number of trophospecies are the keys to study the simplification of food web.In this study,three aggregation methods,including taxonomic aggregation(TA),structural equivalence aggregation(SEA),and self-organizing maps(SOM),were analyzed and compared with the linear inverse model–Markov Chain Monte Carlo(LIM-MCMC)model.Impacts of aggregation methods and trophospecies number on food webs were evaluated based on the robustness and unitless of ecological net-work indices.Results showed that aggregation method of SEA performed better than the other two methods in estimating food web structure and function indices.The effects of aggregation methods were driven by the differences in species aggregation principles,which will alter food web structure and function through the redistribution of energy flow.According to the results of mean absolute percentage error(MAPE)which can be applied to evaluate the accuracy of the model,we found that MAPE in food web indices will increase with the reducing trophospecies number,and MAPE in food web function indices were smaller and more stable than those in food web structure indices.Therefore,trade-off between simplifying food webs and reflecting the status of ecosystem should be con-sidered in food web studies.These findings highlight the importance of aggregation methods and trophospecies number in the analy-sis of food web simplification.This study provided a framework to explore the extent to which food web models are affected by dif-ferent species aggregation,and will provide scientific basis for the construction of food webs.

Continuum Skyrme Hartree-Fock-Bogoliubov theory with Green’s function method for neutron-rich Ca,Ni,Zr,and Sn isotopes

The possible exotic nuclear properties in the neutron-rich Ca,Ni,Zr,and Sn isotopes are examined with the continuum Skyrme Hartree-Fock-Bogoliubov theory in the framework of the Green’s function method.The pairing correlation,the couplings with the continuum,and the blocking effects for the unpaired nucleon in odd-A nuclei are properly treated.The Skyrme interaction SLy4 is adopted for the ph channel and the density-dependentinteraction is adopted for the pp chan-nel,which well reproduce the experimental two-neutron separation energies S_(2n)and one-neutron separation energies Sn.It is found that the criterion S_(n)>0 predicts a neutron drip line with neutron numbers much smaller than those for S_(2n)>0.Owing to the unpaired odd neutron,the neutron pairing energies−E_(pair)in odd-A nuclei are much lower than those in the neighbor-ing even-even nuclei.By investigating the single-particle structures,the possible halo structures in the neutron-rich Ca,Ni,and Sn isotopes are predicted,where sharp increases in the root-mean-square(rms)radii with significant deviations from the traditional rA^(1∕3)rule and diffuse spatial density distributions are observed.Analyzing the contributions of various partial waves to the total neutron densityρlj(r)∕ρ(r)reveals that the orbitals located around the Fermi surface-particularly those with small angular momenta-significantly affect the extended nuclear density and large rms radii.The number of neutrons Nλ(N_(0))occupying above the Fermi surfacen(continuum threshold)is discussed,whose evolution as a function of the mass number A in each isotope is consistent with that of the pairing energy,supporting the key role of the pairing correlation in halo phenomena.

An Improved Preisach Distribution Function Identification Method Considering the Reversible Magnetization

This paper presents an identification method of the scalar Preisach model to consider the effect of reversible magnetization in the process of distribution function identification.By reconsidering the identification process by stripping the influence of reversible components from the measurement data,the Preisach distribution function is identified by the pure irreversible components.In this way,the simulation accuracy of both limiting hysteresis loops and the inner internal symmetrical small hysteresis loop is ensured.Furthermore,through a discrete Preisach plane with a hybrid discretization method,the irreversible magnetic flux density components are computed more efficiently through the improved Preisach model.Finally,the proposed method results are compared with the traditional method and the traditional method considering reversible magnetization and validated by the laboratory test for the B30P105 electrical steel by Epstein frame.

Generation of endurance time excitation functions using spectral representation method

In this study,application of the spectral representation method for generation of endurance time excitation functions is introduced.Using this method,the intensifying acceleration time series is generated so that its acceleration response spectrum in any desired time duration is compatible with a time-scaled predefined acceleration response spectrum.For this purpose,simulated stationary acceleration time series is multiplied by the time dependent linear modulation function,then using a simple iterative scheme,it is forced to match a target acceleration response spectrum.It is shown that the generated samples have excellent conformity in low frequency,which is useful for nonlinear endurance time analysis.In the second part of this study,it is shown that this procedure can be extended to generate a set of spatially correlated endurance time excitation functions.This makes it possible to assess the performance of long structures under multi-support seismic excitation using endurance time analysis.

Comprehensive Evaluation of New Maize Varieties for Grain and Fodder Based on Membership Function Method

To screen new maize(Zea mays L.)varieties suitable for food and fodder dual-purpose in Du'an Yao Autonomous County of Guangxi,the agronomic characters,yield and quality indexes of 12 new maize varieties were measured,and the correlation between various indexes were analyzed,and the comprehensive performance of tested varieties was evaluated by membership function method.The results showed that Guidan 671 had the highest grain yield and whole-plant biomass at 10908 and 49965 kg/hm^(2),respectively,and the second was Zhaoyu 215 with a grain yield and whole-plant biomass of 10086 and 47175 kg/hm^(2),respectively.Grain yield was highly significantly positively correlated with ear diameter and 100-grain weight(P<0.01),and significantly correlated with whole-plant biomass,starch content,ear length and grain number per row(P<0.05);and the whole-plant biomass was highly significantly correlated with the number of grains per row(P<0.01),and significantly correlated with starch content,panicle length,plant height and panicle height(P<0.05).The comprehensive performance scores of the tested varieties from high to low were Guidan 671,Zhaoyu 215,Guidan 669,Guidan 6208,Guidan 666,Guidan 6205,Guidan 660,Guidan 6203,Guidan 6206,Guidan 162,Guidan 668 and Guidan 673.According to the values of membership function and combined with various indexes,Guidan 671 and Zhaoyu 215 had good comprehensive performance,and could be used as the first choice for food and fodder dual-purpose maize varieties in Du'an Yao Autonomous County.

Galerkin Method for Numerical Solution of Volterra Integro-Differential Equations with Certain Orthogonal Basis Function

This paper concerns the implementation of the orthogonal polynomials using the Galerkin method for solving Volterra integro-differential and Fredholm integro-differential equations. The constructed orthogonal polynomials are used as basis functions in the assumed solution employed. Numerical examples for some selected problems are provided and the results obtained show that the Galerkin method with orthogonal polynomials as basis functions performed creditably well in terms of absolute errors obtained.

"Meridian-sinew Acupuncture Method"Improves Deformity Joint Dysfunction in Rheumatoid Arthritis

Joint deformity and dysfunction are common and serious complications in the late stage of rheumatoid arthritis,which seriously affect the quality of life of patients.Traditional Chinese medicine(TCM)believes that joint deformity and dysfunction in some patients with rheumatoid arthritis are closely related to the apraxia of meridians and tendons due to enduring illness.Based on the theory of meridians and tendons circulation,using the local and nearby therapeutic effect of acupoints as the treatment method in clinical practice,we conducted penetration needling of Houxi,Baxie,Wailaogong as well as Ashi points of interphalangeal joints of both hands through bilateral Sanjian,and used the uniform reinforcing-reducing method to soothe tendons and meridians,thus effectively improving the dysfunction of deformed joints.

Exactness of penalization for exact minimax penalty function method in nonconvex programming

The exact minimax penalty function method is used to solve a noncon- vex differentiable optimization problem with both inequality and equality constraints. The conditions for exactness of the penalization for the exact minimax penalty function method are established by assuming that the functions constituting the considered con- strained optimization problem are invex with respect to the same function η （with the exception of those equality constraints for which the associated Lagrange multipliers are negative these functions should be assumed to be incave with respect to η）. Thus, a threshold of the penalty parameter is given such that, for all penalty parameters exceeding this threshold, equivalence holds between the set of optimal solutions in the considered constrained optimization problem and the set of minimizer in its associated penalized problem with an exact minimax penalty function. It is shown that coercivity is not suf- ficient to prove the results.

Crack Fault Diagnosis and Location Method for a Dual-Disk Hollow Shaft Rotor System Based on the Radial Basis Function Network and Pattern Recognition Neural Network

The crack fault is one of the most common faults in the rotor system,and researchers have paid close attention to its fault diagnosis.However,most studies focus on discussing the dynamic response characteristics caused by the crack rather than estimating the crack depth and position based on the obtained vibration signals.In this paper,a novel crack fault diagnosis and location method for a dual-disk hollow shaft rotor system based on the Radial basis function(RBF)network and Pattern recognition neural network(PRNN)is presented.Firstly,a rotor system model with a breathing crack suitable for a short-thick hollow shaft rotor is established based on the finite element method,where the crack's periodic opening and closing pattern and different degrees of crack depth are considered.Then,the dynamic response is obtained by the harmonic balance method.By adjusting the crack parameters,the dynamic characteristics related to the crack depth and position are analyzed through the amplitude-frequency responses and waterfall plots.The analysis results show that the first critical speed,first subcritical speed,first critical speed amplitude,and super-harmonic resonance peak at the first subcritical speed can be utilized for the crack fault diagnosis.Based on this,the RBF network and PRNN are adopted to determine the depth and approximate location of the crack respectively by taking the above dynamic characteristics as input.Test results show that the proposed method has high fault diagnosis accuracy.This research proposes a crack detection method adequate for the hollow shaft rotor system,where the crack depth and position are both unknown.

Properties of a family of merit functions and a merit function method for the NCP

A family of merit functions are proposed, which are the generalization of several existing merit functions. A number of favorable properties of the proposed merit functions are established. By using these properties, a merit function method for solving nonlinear complementarity problem is investigated, and the global convergence of the proposed algorithm is proved under some standard assumptions. Some preliminary numerical results are given.