A stable implicit nodal integration-based particle finite element method(N-PFEM)for modelling saturated soil dynamics

In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.

Profiling Electronic and Phononic Band Structures of Semiconductors at Finite Temperatures: Methods and Applications

Semiconductor devices are often operated at elevated temperatures that are well above zero Kelvin,which is the temperature in most first-principles density functional calculations.Computational approaches to com-puting and understanding the properties of semiconductors at finite temperatures are thus in critical demand.In this review,we discuss the recent progress in computationally assessing the electronic and phononic band structures of semiconductors at finite temperatures.As an emerging semiconductor with particularly strong temperature-induced renormalization of the electronic and phononic band structures,halide perovskites are used as a representative example to demonstrate how computational advances may help to understand the band struc-tures at elevated temperatures.Finally,we briefly illustrate the remaining computational challenges and outlook promising research directions that may help to guide future research in this field.

THE SUPERCLOSENESS OF THE FINITE ELEMENT METHOD FOR A SINGULARLY PERTURBED CONVECTION-DIFFUSION PROBLEM ON A BAKHVALOV-TYPE MESH IN 2D

For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of the mesh in the layer adjacent to the transition point,resulting in a suboptimal estimate for convergence.Existing analysis techniques cannot handle these difficulties well.To fill this gap,here a novel interpolation is designed delicately for the smooth part of the solution,bringing about the optimal supercloseness result of almost order 2 under an energy norm for the finite element method.Our theoretical result is uniform in the singular perturbation parameterεand is supported by the numerical experiments.

Natural Convection and Irreversibility of Nanofluid Due to Inclined Magnetohydrodynamics(MHD)Filled in a Cavity with Y-Shape Heated Fin:FEM Computational

This study explains the entropy process of natural convective heating in the nanofluid-saturated cavity in a heated fin andmagnetic field.The temperature is constant on the Y-shaped fin,insulating the topwall while the remaining walls remain cold.All walls are subject to impermeability and non-slip conditions.The mathematical modeling of the problem is demonstrated by the continuity,momentum,and energy equations incorporating the inclined magnetic field.For elucidating the flow characteristics Finite ElementMethod(FEM)is implemented using stable FE pair.A hybrid fine mesh is used for discretizing the domain.Velocity and thermal plots concerning parameters are drawn.In addition,a detailed discussion regarding generation energy by monitoring changes in magnetic,viscous,total,and thermal irreversibility is provided.In addition,line graphs are created for the u and v components of the velocity profile to predict the flow behavior.Current simulations assume the dimensionless representative of magnetic field Hartmann number Ha between 0 and 100 and a magnetic field inclination between 0 and 90 degrees.A constant 4% volume proportion of nanoparticles is employed throughout all scenarios.

Recurrent Neural Network Inspired Finite-Time Control Design

Dear Editor,This letter is concerned with the role of recurrent neural networks(RNNs)on the controller design for a class of nonlinear systems.Inspired by the architectures of RNNs,the system states are stacked according to the dynamic along with time while the controller is represented as the neural network output.To build the bridge between RNNs and finite-time controller,a novel activation function is imposed on RNNs to drive the convergence of states at finite-time and propel the overall control process smoother.Rigorous stability proof is briefly provided for the convergence of the proposed finite-time controller.At last,a numerical simulation example is presented to illustrate the efficiency of the proposed strategy.Neural networks can be classified as static(feedforward)and dynamic(recurrent)nets[1].The former nets do not perform well in dealing with training data and using any information of the local data structure[2].In contrast to the feedforward neural networks,RNNs are constituted by high dimensional hidden states with dynamics.

On entire solutions of some Fermat type differential-difference equations

On one hand,we study the existence of transcendental entire solutions with finite order of the Fermat type difference equations.On the other hand,we also investigate the existence and growth of solutions of nonlinear differential-difference equations.These results extend and improve some previous in[5,14].

Effect of Size Change on Mechanical Properties of Monolayer Arsenene

The Young's modulus, shear modulus and Poisson's ratio of monolayer arsenene with different sizes were calculated by finite element method, so as to explore the influence of dimension and orientation on the mechanical properties of monolayer arsenene. The calculation results show that the small size has a significant effect on the mechanical properties of the monolayer arsenene. The smaller the size, the larger the Young's modulus and Poisson's ratio of the monolayer arsenene. The size change has a great influence on the Young's modulus of the arsenene handrail direction, and the Young's modulus of the zigzag direction is not sensitive to the size change. Similarly, the size change has a significant effect on the shear modulus of arsenene in the handrail direction, while the shear modulus in the zigzag direction has no significant effect on its size change. For the Poisson's ratio, the situation is just the opposite, and the effect of the size change on the Poisson's ratio of the arsenene zigzag direction is greater than that of the handrail direction.

Tuning microstructures of TC4 ELI to improve explosion resistance

A reasonable heat treatment process for TC4 ELI titanium alloy is crucial to tune microstructures to improve its explosion resistance.However,there is limited investigation on tuning microstructures of TC4 ELI to improve explosion resistance.Moreover,the current challenge is quantifying microstructural changes'effects on explosion resistance and incorporating microstructural changes into finite element models.This work aims to tune microstructures to improve explosion resistance and elucidate their anti-explosion mechanism,and find a suitable method to incorporate microstructural changes into finite element models.In this work,we systematically study the deformation and failure characteristics of TC4 ELI plates with varying microstructures using an air explosion test and LS-DYNA finite element modeling.The Johnson-Cook(JC)constitutive parameters are used to quantify the effects of microstructural changes on explosion resistance and incorporate microstructural changes into finite element models.Because of the heat treatment,one plate has equiaxed microstructure and the other has bimodal microstructure.The convex of the plate after the explosion has a quadratic relationship with the charge mass,and the simulation results demonstrate high reliability,with the error less than 17.5%.Therefore,it is feasible to obtain corresponding JC constitutive parameters based on the differences in microstructures and mechanical properties and characterize the effects of microstructural changes on explosion resistance.The bimodal target exhibits excellent deformation resistance.The response of bimodal microstructure to the shock wave may be more intense under explosive loading.The well-coordinated structure of the bimodal target enhances its resistance to deformation.

Study on Dynamic Mechanical Behavior of Al-Mg-Si Alloy

The dynamic mechanical behavior of Al-Mg-Si alloy was investigated under different strain rates by mechanical property and microstructure characterization,constitutive behavior analysis and numerical simulation in the present study.As the strain rate increases,the yield strength,ultimate tensile strength and elongation increase first,then remain almost constant,and finally increase.The alloy always exhibits a typical ductile fracture mode,not depending on the strain rate.However,as the strain rate increases,the number of dimples gradually increases.Tensile deformation can refine grains,however,the grain structure is slightly affected by the strain rate.An optimized Johnson-Cook constitutive equation was used to describe the mechanical behavior and obtained by fitting the true stress-strain curves.The parameter C was described by a function related to the strain rate.The fitting true stress-strain curves by the JC model agree very well with the experimental true stress-strain curves.The true stress-strain curves calculated by the finite element numerical simulation agree well with the experimental true stress-strain curves.

Characteristics of the electromagnetic wave propagation in magnetized plasma sheath and practical method for blackout mitigation

“Magnetic window”is considered as an effective method to solve the communication blackout issue.COMSOL software package based on the finite element method is utilized to simulate the propagation of right-handed circularly polarized wave in the magnetized plasma sheath.We assume a double Gaussian model of electron density and an exponential attenuation model of magnetic field.The propagation characteristics of right-handed circularly polarized wave are analyzed by the observation of the reflected,transmitted and loss coefficient.The numerical results show that the propagation of right-handed circularly polarized wave in the magnetized plasma sheath varies for different incident angles,collision frequencies,non-uniform magnetic fields and non-uniform plasma densities.We notice that reducing the wave frequency can meet the propagation conditions of whistle mode in the weak magnetized plasma sheath.And the transmittance of whistle mode is less affected by the variation of the electron density and the collision frequency.It can be used as a communication window.

Modeling the Interaction between Vacancies and Grain Boundaries during Ductile Fracture

The experimental results in previous studies have indicated that during the ductile fracture of pure metals,vacancies aggregate and form voids at grain boundaries.However,the physical mechanism underlying this phenomenon remains not fully understood.This study derives the equilibrium distribution of vacancies analytically by following thermodynamics and the micromechanics of crystal defects.This derivation suggests that vacancies cluster in regions under hydrostatic compression to minimize the elastic strain energy.Subsequently,a finite element model is developed for examining more general scenarios of interaction between vacancies and grain boundaries.This model is first verified and validated through comparison with some available analytical solutions,demonstrating consistency between finite element simulation results and analytical solutions within a specified numerical accuracy.A systematic numerical study is then conducted to investigate the mechanism that might govern the micromechanical interaction between grain boundaries and the profuse vacancies typically generated during plastic deformation.The simulation results indicate that the reduction in total elastic strain energy can indeed drive vacancies toward grain boundaries,potentially facilitating void nucleation in ductile fracture.

Preliminary electromagnetic analysis of the COOL blanket for CFETR

The supercritical CO_(2)cOoled Lithium-Lead(COOL)blanket has been designed as one advanced blanket candidate for the Chinese Fusion Engineering Test Reactor(CFETR).This work focuses on the electromagnetic(EM)loads(Maxwell force and Lorentz force)acting on the COOL blanket,which are important mechanical loads in further structural analysis of the COOL blanket.A 3D electromagnetic analysis is performed using the ANSYS finite element method to obtain EM loads on the COOL blanket in this study.At first,the magnetic scalar potential(MSP)method is used to obtain the magnetic field and the Maxwell force on the COOL blanket.Then,the magnetic vector potential(MVP)method is performed during a plasma disruption event to get the eddy current distribution.At last,a multi-step method is adopted for the calculation of the Lorentz force and the torque.The maximum Lorentz forces of inboard and outboard blanket structural components are 5624 kN and 2360 kN respectively.

裂纹加筋和缺口加筋对钢加筋板极限强度的影响

This paper numerically evaluates the effect of the crack position on the ultimate strength of stiffened panels.Imperfections such as notches and cracks in aged marine stiffened panels can reduce their ultimate strength.To investigate the effect of crack length and position,a series of nonlinear finite element analyses were carried out and two cases were considered,i.e.,case 1 with thin stiffeners and case 2 with thick stiffeners.In both cases,the stiffeners have the same cross-section area.To have a basis for comparison,the intact panels were modeled as well.The cracks and notches were in the longitudinal and transverse direction and were assumed to be in the middle part of the panel.The cracks and notches were assumed to be through the thickness and there is neither crack propagation nor contact between crack faces.Based on the numerical results,longitudinal cracks affect the behavior of the stiffened panels in the postbuckling region.When the stiffeners are thinner,they buckle first and provide no reserved strength after plate buckling.Thus,cracks in the stiffeners do not affect the ultimate strength in the case of the thinner stiffeners.Generally,when stiffeners are thicker,they affect the postbuckling behavior more.In that case,cracks in the stiffeners affect the buckling and failure modes of the stiffened panels.The effect of notch was also studied.In contrast to the longitudinal crack in stiffeners,a notch in the stiffeners reduces the ultimate strength of the stiffened panel for both slender and thick stiffeners.

Modelling rapid non-destructive test using light weight deflectometer on granular soils across different degrees of saturation

This study introduces an advanced finite element model for the light weight deflectometer(LWD),which integrates contact mechanics with fully coupled models.By simulating LWD tests on granular soils at various saturation levels,the model accurately reflects the dependence of the LWD modulus on dry density,water content,and effective stress.This model addresses and overcomes the limitations of previous finite element models for this specific problem.Simultaneously,this research presents the first experimentally validated fully coupled contact impact model.Furthermore,the research provides a comparative assessment of elastoplastic and nonlinear elastic models and contrasts an enriched node-tosegment method(developed in this study)with the more precise mortar technique for contact mechanics.These comparisons reveal unique advantages and challenges for each method.Moreover,the study underscores the importance of careful application of the LWD modulus,emphasising the need for sophisticated tools to interpret soil behaviour accurately.

Effects of layer interactions on instantaneous stability of finite Stokes flows

The stability analysis of a finite Stokes layer is of practical importance in flow control. In the present work, the instantaneous stability of a finite Stokes layer with layer interactions is studied via a linear stability analysis of the frozen phases of the base flow. The oscillations of two plates can have different velocity amplitudes, initial phases, and frequencies. The effects of the Stokes-layer interactions on the stability when two plates oscillate synchronously are analyzed. The growth rates of two most unstable modes when δ < 0.12 are almost equal, and δ = δ*/h*, where δ*and h*are the Stokes-layer thickness and the half height of the channel, respectively. However, their vorticities are different. The vorticity of the most unstable mode is symmetric, while the other is asymmetric. The Stokes-layer interactions have a destabilizing effect on the most unstable mode when δ < 0.68, and have a stabilizing effect when δ > 0.68. However, the interactions always have a stabilizing effect on the other unstable mode. It is explained that one of the two unstable modes has much higher dissipation than the other one when the Stokes-layer interactions are strong. We also find that the stability of the Stokes layer is closely related to the inflectional points of the base-flow velocity profile. The effects of inconsistent velocity-amplitude, initial phase, and frequency of the oscillations on the stability are analyzed. The energy of the most unstable eigenvector is mainly distributed near the plate of higher velocity amplitude or higher oscillation frequency. The effects of the initial phase difference are complicated because the base-flow velocity is extremely sensitive to the initial phase.

Application of the finite analytic numerical method to a flowdependent variational data assimilation

An anisotropic diffusion filter can be used to model a flow-dependent background error covariance matrix,which can be achieved by solving the advection-diffusion equation.Because of the directionality of the advection term,the discrete method needs to be chosen very carefully.The finite analytic method is an alternative scheme to solve the advection-diffusion equation.As a combination of analytical and numerical methods,it not only has high calculation accuracy but also holds the characteristic of the auto upwind.To demonstrate its ability,the one-dimensional steady and unsteady advection-diffusion equation numerical examples are respectively solved by the finite analytic method.The more widely used upwind difference method is used as a control approach.The result indicates that the finite analytic method has higher accuracy than the upwind difference method.For the two-dimensional case,the finite analytic method still has a better performance.In the three-dimensional variational assimilation experiment,the finite analytic method can effectively improve analysis field accuracy,and its effect is significantly better than the upwind difference and the central difference method.Moreover,it is still a more effective solution method in the strong flow region where the advective-diffusion filter performs most prominently.

A novel mitigation measure for normal fault-induced deformations on pile-raft systems

Evidence from recent earthquakes has shown destructive consequences of fault-induced permanent ground movement on structures.Such observations have increased the demand for improvements in the design of structures that are dramatically vulnerable to surface fault ruptures.In this study a novel connection between the raft and the piles is proposed to mitigate the hazards associated with a normal fault on pile-raft systems by means of 3D finite element(FE)modeling.Before embarking on the parametric study,the strain-softening constitutive law used for numerical modeling of the sand has been validated against centrifuge test results.The exact location of the fix-head and unconnected pile-raft systems relative to the outcropping fault rupture in the free-field is parametrically investigated,revealing different failure mechanisms.The performance of the proposed connection for protecting the pile-raft system against normal fault-induced deformations is assessed by comparing the geotechnical and structural responses of both types of foundation.The results indicate that the pocket connection can relatively reduce the cap rotation and horizontal and vertical displacements of the raft in most scenarios.The proposed connection decreases the bending moment response of the piles to their bending moment capacity,verging on a fault offset of 0.6 m at bedrock.

Predicting the Mechanical Behavior of a Bioinspired Nanocomposite through Machine Learning

The bioinspired nacre or bone structure represents a remarkable example of tough,strong,lightweight,and multifunctional structures in biological materials that can be an inspiration to design bioinspired high-performance materials.The bioinspired structure consists of hard grains and soft material interfaces.While the material interface has a very low volume percentage,its property has the ability to determine the bulk material response.Machine learning technology nowadays is widely used in material science.A machine learning model was utilized to predict the material response based on the material interface properties in a bioinspired nanocomposite.This model was trained on a comprehensive dataset of material response and interface properties,allowing it to make accurate predictions.The results of this study demonstrate the efficiency and high accuracy of the machine learning model.The successful application of machine learning into the material property prediction process has the potential to greatly enhance both the efficiency and accuracy of the material design process.

Numerical modeling of fracture propagation of supercritical CO_(2)compound fracturing

The exploitation of shale gas is promising due to depletion of the conventional energy and intensification of the greenhouse effect.In this paper,we proposed a heat-fluid-solid coupling damage model of supercritical CO_(2)(SC-CO_(2))compound fracturing which is expected to be an efficient and environmentally friendly way to develop shale gas.The coupling model is solved by the finite element method,and the results are in good agreement with the analytical solutions and fracturing experiments.Based on this model,the fracture propagation characteristics at the two stages of compound fracturing are studied and the influence of pressurization rate,in situ stress,bedding angle,and other factors are considered.The results show that at the SC-CO_(2)fracturing stage,a lower pressurization rate is conducive to formation of the branches around main fractures,while a higher pressurization rate inhibits formation of the branches around main fractures and promotes formation of the main fractures.Both bedding and in situ stress play a dominant role in the fracture propagation.When the in situ stress ratio(sx/sy)is 1,the presence of bedding can reduce the initiation pressure and failure pressure.Nevertheless,it will cause the fracture to propagate along the bedding direction,reducing the fracture complexity.In rocks without bedding,hydraulic fracturing has the lengthening and widening effects for SC-CO_(2)induced fracture.In shale,fractures induced at the hydraulic fracturing stage are more likely to be dominated by in situ stresses and have a shorter reorientation radius.Therefore,fracture branches propagating along the maximum principal stress direction may be generated around the main fractures induced by SC-CO_(2)at the hydraulic fracturing stage.When the branches converge with the main fractures,fracture zones are easily formed,and thus the fracture complexity and damage area can be significantly increased.The results are instructive for the design and application of SC-CO_(2)compound fracturing.

Contact Stress Reliability Analysis Model for Cylindrical Gear with Circular Arc Tooth Trace Based on an Improved Metamodel

Although there is currently no unified standard theoretical formula for calculating the contact stress of cylindrical gears with a circular arc tooth trace(referred to as CATT gear),a mathematical model for determining the contact stress of CATT gear is essential for studying how parameters affect its contact stress and building the contact stress limit state equation for contact stress reliability analysis.In this study,a mathematical relationship between design parameters and contact stress is formulated using the KrigingMetamodel.To enhance the model’s accuracy,we propose a new hybrid algorithm that merges the genetic algorithm with the Quantum Particle Swarm optimization algorithm,leveraging the strengths of each.Additionally,the“parental inheritance+self-learning”optimization model is used to fine-tune the KrigingMetamodel’s parameters.Following this,amathematicalmodel for calculating the contact stress of Variable Hyperbolic Circular-Arc-Tooth-Trace(VH-CATT)gears using the optimized Kriging model was developed.We then examined how different gear parameters affect the VH-CATT gears’contact stress.Our simulation results show:(1)Improvements in R2,RMSE,and RMAE.R2 rose from0.9852 to 0.9974(a 1.22%increase),nearing 1,suggesting the optimized Kriging Metamodel’s global error is minimized.Meanwhile,RMSE dropped from3.9210 to 1.6492,a decline of 57.94%.The global error of the GA-IQPSO-Kriging algorithm was also reduced,with RMAE decreasing by 58.69%from 0.1823 to 0.0753,showing the algorithm’s enhanced precision.In a comparison of ten experimental groups selected randomly,the GA-IQPSO-Kriging and FEM-based contact analysis methods were used to measure contact stress.Results revealed a maximum error of 12.11667 MPA,which represents 2.85%of the real value.(2)Several factors,including the pressure angle,tooth width,modulus,and tooth line radius,are inversely related to contact stress.The descending order of their impact on the contact stress is:tooth line radius>modulus>pressure angle>tooth width.(3)Complex interactions are noted among various parameters.Specifically,when the tooth line radius interacts with parameters such as pressure angle,tooth width,and modulus,the resulting stress contour is nonlinear,showcasing amultifaceted contour plane.However,when tooth width,modulus,and pressure angle interact,the stress contour is nearly linear,and the contour plane is simpler,indicating a weaker coupling among these factors.