ANALYSIS OF NONLINEAR PIEZOELECTRIC CIRCULAR SHALLOW SPHERICAL SHELLS BY DIFFERENTIAL QUADRATURE ELEMENT METHOD

The static behavior of piezoelectric circular spherical shallow shells under both electrical and mechanical loads is studied by using the differential quadrature element method (DQEM). Geometrical nonlinearity effect is considered. Detailed formulations and procedures are given for the first time. Several examples are analyzed and accurate results are obtained by the DQEM. Based on the results in this paper, one may conclude that the DQEM is a useful tool for obtaining solutions of structural elements. It can be seen that the shell shape may be theore tically controlled and snap through may occur when the applied voltage reaches a critical value even without mechanical load for certain geometric configurations.

Cl/S与Na相互作用对Shell气化炉合成气冷却器入口积灰机制的影响

Shell干粉煤加压气化是煤炭洁净高效利用的重要技术之一,由碱金属化合物引起的合成气冷却器入口积灰结垢是导致气化炉非正常停工检修的主要原因。以添加不同含量的Na、Cl和S的Shell气化炉飞灰为原料,利用自主设计的高温竖直炉中沉积探针模拟Shell气化炉合成气冷却器入口管路,通过对积灰进行内、外分层研究,探讨内外层积灰质量的变化,并结合ICP-MS、IC、SEM-EDS和XRD等表征手段对内外层积灰的理化性质进行比较分析,获得Na、Cl、S和Fe等不同元素之间的相互作用对积灰行为的影响。结果表明,内层积灰质量随时间延长而增大,含S化合物的添加会降低内外层积灰质量,且外层积灰质量会随着时间延长而减小。Na更多以铝硅酸盐形式在外层积灰中存在,促进积灰增长;Cl通常以碱金属氯化物的形式集中在初始黏性层;S的存在会减缓管路积灰;当Cl和S共同存在时,Fe易与灰中的Si、Al和Na形成多种低温共熔物促进内、外层积灰熔融。Shell气化炉合成气冷却器入口积灰形成机制为:飞灰颗粒组分在Na、Cl、Si和Al的共同作用下,于内层形成碱金属氯化物和铝硅酸盐共晶;同时Cl、S的存在促使Fe和Na迁移到这些共晶中,形成Fe-O-Si、Fe-O-S和Fe-Na-O-Al-S共熔体。进而,铝硅酸盐与多种低温共熔体相互熔融使灰颗粒尺寸增加,促进积灰的进一步生长。

NUMERICAL ANALYSIS OF MINDLIN SHELL BY MESHLESS LOCAL PETROV-GALERKIN METHOD

The objectives of this study are to employ the meshless local Petrov-Galerkin method （MLPGM） to solve three-dimensional shell problems. The computational accuracy of MLPGM for shell problems is affected by many factors, including the dimension of compact support domain, the dimension of quadrture domain, the number of integral cells and the number of Gauss points. These factors＇ sensitivity analysis is to adopt the Taguchi experimental design technology and point out the dimension of the quadrature domain with the largest influence on the computational accuracy of the present MLPGM for shells and give out the optimum combination of these factors. A few examples are given to verify the reliability and good convergence of MLPGM for shell problems compared to the theoretical or the finite element results.

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.

Method of Semi Analytic Perturbation Weighted Residuals for Nonlinear Bending of Shallow Shells

The semi? analytic perturbation weighted residuals method was used to solve the nonlinear bending problem of shallow shells, and the fifth order B spline was taken as trial function to seek an efficient method for nonlinear bending problem of shallow shells. The results from the present method are in good agreement with those derived from other methods. The present method is of higher accuracy, lower computing time and wider adaptability. In addition, the design of computer program is simple and it is easy to be programmed.

Key Techniques and Applications of Adaptive Growth Method for Stiffener Layout Design of Plates and Shells

The application of the adaptive growth method is limited because several key techniques during the design process need manual intervention of designers. Key techniques of the method including the ground structure construction and seed selection are studied, so as to make it possible to improve the effectiveness and applicability of the adaptive growth method in stiffener layout design optimization of plates and shells. Three schemes of ground structures, which are comprised by different shell elements and beam elements, are proposed. It is found that the main stiffener layouts resulted from different ground structures are almost the same, but the ground structure comprised by 8-nodes shell elements and both 3-nodes and 2-nodes beam elements can result in clearest stiffener layout, and has good adaptability and low computational cost. An automatic seed selection approach is proposed, which is based on such selection rules that the seeds should be positioned on where the structural strain energy is great for the minimum compliance problem, and satisfy the dispersancy requirement. The adaptive growth method with the suggested key techniques is integrated into an ANSYS-based program, which provides a design tool for the stiffener layout design optimization of plates and shells. Typical design examples, including plate and shell structures to achieve minimum compliance and maximum bulking stability are illustrated. In addition, as a practical mechanical structural design example, the stiffener layout of an inlet structure for a large-scale electrostatic precipitator is also demonstrated. The design results show that the adaptive growth method integrated with the suggested key techniques can effectively and flexibly deal with stiffener layout design problem for plates and shells with complex geometrical shape and loading conditions to achieve various design objectives, thus it provides a new solution method for engineering structural topology design optimization.

Thermal Modal Analysis of Doubly Curved Shell Based on Rayleigh⁃Ritz Method

The doubly curved shell(DCS)is a common structure in the engineering field.In a thermal environment,the vibration characteristics of the DCS will be affected by the thermal effect.The research on the vibration characteristics of DCS in thermal environment is relatively limited.In this paper,the thermal strain and the change of Young’s modulus caused by the changing of temperature are studied,and the DCS energy equation is established systematically.The displacement tolerance function of the DCS is constructed by the spectral geometry method,and the natural frequencies and mode shapes of the DCS with different structural parameters,such as thicknesses,ratios of R_(a)/R_(b) and a/b,at different temperatures are solved by the Rayleigh-Ritz method.The results show that the natural frequency of the DCS decreases with the increasing temperature,R_(a)/R_(b) and a/b ratios,and increases with the increasing thickness.

Free vibration and critical speed of moderately thick rotating laminated composite conical shell using generalized differential quadrature method

The generalized differential quadrature method （GDQM） is employed to con- sider the free vibration and critical speed of moderately thick rotating laminated compos- ite conical shells with different boundary conditions developed from the first-order shear deformation theory （FSDT）. The equations of motion are obtained applying Hamilton＇s concept, which contain the influence of the centrifugal force, the Coriolis acceleration, and the preliminary hoop stress. In addition, the axial load is applied to the conical shell as a ratio of the global critical buckling load. The governing partial differential equations are given in the expressions of five components of displacement related to the points ly- ing on the reference surface of the shell. Afterward, the governing differential equations are converted into a group of algebraic equations by using the GDQM. The outcomes are achieved considering the effects of stacking sequences, thickness of the shell, rotating velocities, half-vertex cone angle, and boundary conditions. Furthermore, the outcomes indicate that the rate of the convergence of frequencies is swift, and the numerical tech- nique is superior stable. Three comparisons between the selected outcomes and those of other research are accomplished, and excellent agreement is achieved.

Neumann＇s method for boundary problems of thin elastic shells

The possibility of using Neumann＇s method to solve the boundary problems for thin elastic shells is studied. The variational statement of the static problems for the shells allows for a problem examination within the distribution space. The convergence of Neumann＇s method is proven for the shells with holes when the boundary of the domain is not completely fixed. The numerical implementation of Neumann＇s method normally requires significant time before any reliable results can be achieved. This paper suggests a way to improve the convergence of the process, and allows for parallel computing and evaluation during the calculations.

Natural frequency characteristics of thin-walled homogeneous and manifold layered cylindrical shells under pressure using energy method

Energy method for the vibration of two types of cylindrical shells,namely thin-walled homogeneous isotropic and manifold layered isotropic cylindrical shells under uniform external lateral pressure is presented.The study is carried out based on strain-displacement relationship from Love's shell theory with beam functions as axial modal function.A manifold layered cylindrical shell configuration is formed by three layers of isotropic material where the inner and outer layers are stainless steel and the middle layer is aluminum.The homogeneous cylindrical shell is made-up of isotropic one layer with stainless steel.The governing equations with uniform external lateral pressure for homogeneous isotropic and manifold layered isotropic cylindrical shells are obtained using energy functional by the Lagrangian function with Rayleigh-Ritz method.The boundary conditions that are presented at the end conditions of the cylindrical shell are simply supported-simply supported,clamped-clamped and free-free.The influences of uniform external lateral pressure and symmetrical boundary conditions on the natural frequency characteristics for both homogeneous and manifold layered isotropic cylindrical shells are examined.For all boundary conditions considered,the natural frequency of both cylindrical shells with symmetric uniform lateral pressure increases as h/R ratio increases and those considering natural frequency of the both cylindrical shells with symmetric uniform lateral pressure decrease as L/R ratio increases.

Effects of Different Drying Methods on the Functional and Structural Properties of Dietary Fiber from Peanut Shell

[Objectives] This study was conducted to investigate the drying methods,functional and structure properties of dietary fiber( DF) from peanut shells.[Methods]Peanut shells were used as a raw material to prepare peanut shell dietary fiber( DF) by hot air drying( HA) and vacuum freeze drying( VF),respectively,and their functional and structural characteristics were compared in detail. [Results]The solubility,water holding capacity,oil holding capacity and swelling capacity of HA-DF and VF-DF were 2. 15 %,7. 63 g/g,7. 73 g/g,10. 35 ml/g and 3. 85 %,14. 98 g/g,15. 25 g/g,15. 85 ml/g,respectively. The total phenol contents were 2. 623 and 5. 173 mg GAE/g,respectively. The IC(50) values of ·OH,O2^-· and DPPH free radicals were 4. 16 and 4. 09 mg/ml,7. 90 and 3. 32 mg/ml,and 3. 19 and 3. 09 mg/ml,respectively. The molecular weight of VF-DF was smaller,and it had narrow molecular weight distribution and denser particles. Electron microscopy showed that VF-DF had a porous network like honeycomb and swelled structure. [Conclusions]This study can provide a theoretical basis for the functional modification and comprehensive utilization of peanut shell dietary fiber.

h-ADAPTIVE ANALYSIS BASED ON MESHLESS LOCAL PETROV-G ALERKIN METHOD WITH B SPLINE WAVELET FOR PLATES AND SHELLS

Using the two-scale decomposition technique, the h-adaptive meshless local Petrov- Galerkin method for solving Mindlin plate and shell problems is presented. The scaling functions of B spline wavelet are employed as the basis of the moving least square method to construct the meshless interpolation function. Multi-resolution analysis is used to decompose the field variables into high and low scales and the high scale component can commonly represent the gradient of the solution according to inherent characteristics of wavelets. The high scale component in the present method can directly detect high gradient regions of the field variables. The developed adaptive refinement scheme has been applied to simulate actual examples, and the effectiveness of the present adaptive refinement scheme has been verified.

Nonlinear dynamic response and buckling of laminated cylindrical shells with axial shallow groove based on a semi-analytical method

The effect of axial shallow groove on the nonlinear dynamic response and buckling of laminated cylindrical shells subjected to radial compression loading was investigated. Based on the first-order shear deformation theory （FSDT）, the nonlinear dynamic equations involving the transverse shear deformation and initial geometric imperfections were derived with the Hamilton philosophy. The axial shallow groove of the laminated composite cylindrical shell was treated as the initial geometric imperfections in the dynamic equations. A semi-analytical method of expanding displacements and loads along the circumferential direction and employing the finite difference method along the axial direction and in the time domain is used to solve the governing equations and obtain the dynamic response of the laminated shell. The B-R criterion was employed to determine the critical loads of dynamic buckling of the shell. The effects of the parameters of the shallow groove on the dynamic response and buckling were discussed in this paper and the results show that the axial shallow grooves greatly affect the dynamic response and buckling.

FINITE ELEMENT DISPLACEMENT PERTURBATION METHOD FOR GEOMETRIC NONLINEAR BEHAVIORS OF SHELLS OF REVOLUTION OVERALL BEDING IN A MERIDIONAL PLANE AND APPLICATION TO BELLOW (Ⅱ)

The finite_element_displacement_perturbation method (FEDPM)for the geometric nonlinear behaviors of shells of revolution subjected to pure bending moments or lateral forces in one of their meridional planes (Ⅰ) was employed to calculate the stress distributions and the stiffness of the bellows. Firstly, by applying the first_order perturbation solution (the linear solution)of the FEDPM to the bellows, the obtained results were compared with those of the general solution and the initial parameter integration solution proposed by the present authors earlier, as well as of the experiments and the FEA by others.It is shown that the FEDPM is with good precision and reliability, and as it was pointed out in (Ⅰ) the abrupt changes of the meridian curvature of bellows would not affect the use of the usual straight element. Then the nonlinear behaviors of the bellows were discussed. As expected, the nonlinear effects mainly come from the bellows ring plate,and the wider the ring plate is, the stronger the nonlinear effects are. Contrarily, the vanishing of the ring plate, like the C_shaped bellows, the nonlinear effects almost vanish. In addition, when the pure bending moments act on the bellows, each convolution has the same stress distributions calculated by the linear solution and other linear theories, but by the present nonlinear solution they vary with respect to the convolutions of the bellows. Yet for most bellows, the linear solutions are valid in practice.

Sound radiation of a functionally graded material cylindrical shell in water by mobility method

Based on the fundamental dynamic equations of functionally graded material (FGM) cylindrical shell, this paper investigates the sound radiation of vibrational FGM shell in water by mobility method. This model takes into account the exterior fluid loading due to the sound press radiated by the FGM shell. The FGM cylindrical shell was excited by a harmonic line radial force uniformly distributing along the generator. The FGM shell equations of motion, the Helmholtz equation in the exterior fluid medium and the continuity equation at fluid-shell interface are used in this vibroacoustic problem. The expressions of sound radiation efficiency and sound field of the FGM shell have been derived by mobility method. Radiation efficiency, modal mobility and the directivity pattern of the sound field are solved numerically. In particular, radiation efficiency and directivity pattern with various power law index are analyzed.

Euler’s First-Order Explicit Method–Peridynamic Differential Operator for Solving Population Balance Equations of the Crystallization Process

Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.

Energy Stable Nodal DG Methods for Maxwell’s Equations of Mixed-Order Form in Nonlinear Optical Media

In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic Kerr response,and the nonlinear delayed Raman molecular vibrational response.Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al.(J Comput Phys 350:420–452,2017)and Lyu et al.(J Sci Comput 89:1–42,2021),a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part(i.e.,the auxiliary differential equations)modeling the linear and nonlinear dispersion in the material.The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization.A nodal discontinuous Galerkin(DG)method is further applied in space for efficiently handling nonlinear terms at the algebraic level,while preserving the energy stability and achieving high-order accuracy.Indeed with d_(E)as the number of the components of the electric field,only a d_(E)×d_(E)nonlinear algebraic system needs to be solved at each interpolation node,and more importantly,all these small nonlinear systems are completely decoupled over one time step,rendering very high parallel efficiency.We evaluate the proposed schemes by comparing them with the methods in Bokil et al.(2017)and Lyu et al.(2021)(implemented in nodal form)regarding the accuracy,computational efficiency,and energy stability,by a parallel scalability study,and also through the simulations of the soliton-like wave propagation in one dimension,as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric(TE)mode of the equations.

Experimental Study on Full-Surface Buckling of Variable Curvature Cylindrical Shell Using Multi-camera 3D-DIC System

To achieve full-surface strain measurement of variable curvature objects,a 360°3D digital image correlation(DIC)system is proposed.The measurement system consists of four double-camera systems,which capture the object’s entire surface from multiple angles,enabling comprehensive full-surface measurement.To increase the stitching quality,a hierarchical coordinate matching method is proposed.Initially,a 3D rigid body calibration auxiliary block is employed to track motion trajectory,which enables preliminary matching of four 3D-DIC sub-systems.Subsequently,secondary precise matching is performed based on feature points on the test specimen’s surface.Through the hierarchical coordinate matching method,the local 3D coordinate systems of each double-camera system are unified into a global coordinate system,achieving 3D surface reconstruction of the variable curvature cylindrical shell,and error analysis is conducted on the results.Furthermore,axial compression buckling experiment is conducted to measure the displacement and strain fields on the cylindrical shell’s surface.The experimental results are compared with the finite element analysis,validating the accuracy and effectiveness of the proposed multi-camera 3D-DIC measuring system.

The forecasting efficiency under different selected regions by Pattern Informatics Method and seismic potential estimation in the North-South Seismic Zone

In 2022,four earthquakes with M_(S)≥6.0 including the Menyuan M_(S)6.9 and Luding M_(S)6.8 earthquakes occurred in the North-South Seismic Zone(NSSZ),which demonstrated high and strong seismicity.Pattern Informatics(PI)method,as an effective long and medium term earthquake forecasting method,has been applied to the strong earthquake forecasting in Chinese mainland and results have shown the positive performance.The earthquake catalog with magnitude above M_(S)3.0 since 1970 provided by China Earthquake Networks Center was employed in this study and the Receiver Operating Characteristic(ROC)method was applied to test the forecasting efficiency of the PI method in each selected region related to the North-South Seismic Zone systematically.Based on this,we selected the area with the best ROC testing result and analyzed the evolution process of the PI hotspot map reflecting the small seismic activity pattern prior to the Menyuan M_(S)6.9 and Luding M_(S)6.8 earthquakes.A“forward”forecast for the area was carried out to assess seismic risk.The study shows the following.1)PI forecasting has higher forecasting efficiency in the selected study region where the difference of seismicity in any place of the region is smaller.2)In areas with smaller differences of seismicity,the activity pattern of small earthquakes prior to the Menyuan M_(S)6.9 and Luding M_(S)6.8 earthquakes can be obtained by analyzing the spatio-temporal evolution process of the PI hotspot map.3)The hotspot evolution in and around the southern Tazang fault in the study area is similar to that prior to the strong earthquakes,which suggests the possible seismic hazard in the future.This study could provide some ideas to the seismic hazard assessment in other regions with high seismicity,such as Japan,Californi,Turkey,and Indonesia.

Shape Sensing of Thin Shell Structure Based on Inverse Finite Element Method

Shape sensing as a crucial component of structural health monitoring plays a vital role in real-time actuation and control of smart structures,and monitoring of structural integrity.As a model-based method,the inverse finite element method(iFEM)has been proved to be a valuable shape sensing tool that is suitable for complex structures.In this paper,we propose a novel approach for the shape sensing of thin shell structures with iFEM.Considering the structural form and stress characteristics of thin-walled structure,the error function consists of membrane and bending section strains only which is consistent with the Kirchhoff–Love shell theory.For numerical implementation,a new four-node quadrilateral inverse-shell element,iDKQ4,is developed by utilizing the kinematics of the classical shell theory.This new element includes hierarchical drilling rotation degrees-of-freedom(DOF)which enhance applicability to complex structures.Firstly,the reconstruction performance is examined numerically using a cantilever plate model.Following the validation cases,the applicability of the iDKQ4 element to more complex structures is demonstrated by the analysis of a thin wallpanel.Finally,the deformation of a typical aerospace thin-wall structure(the composite tank)is reconstructed with sparse strain data with the help of iDKQ4 element.