Search

Many plant materials in nature have the ability to change their shape to respond to external stimuli, such as humidity or moisture, to ensure their survival or safe seed release. A well-known example for this phenomenon is the pinecone, which is able to open its scales at low humidity due to the specific bilayer structures of the scale. Inspired by this, we developed a novel humidity-driven actuator based on paper. This was realized by the lamination of untreated paper made from eucalyptus fibers to a paper–carboxymethyl cellulose (CMC) composite. As observed, the hygroexpansion of the composite can be easily controlled by the amount of CMC in the impregnated paper sheet, which, thus, controls the morphologic deformation of the paper bilayer. For a more detailed understanding of these novel paper soft robots, we also studied the dynamic water vapor adsorption, polymer distribution and hygroexpansion of the paper–polymer composites. Finally, we applied a geometrically nonlinear finite element model to predict the bending behavior of paper bilayers and compared the results to experimental data. From this, we conclude that due to the complexity of structure of the paper composite, a universal prediction of the hygromorphic behavior is not a trivial matter. [ABSTRACT FROM AUTHOR]

The dynamic mechanical behavior and cushioning performance of honeycomb sandwich panels, which are extensively employed in product cushioning packaging due to their exceptional energy absorption capabilities, were examined using a combination of experimental and numerical methods. Several factors, such as maximum acceleration–static stress, cushioning coefficient–static stress, and other curves, were analyzed under various impact conditions. The simulated stress–strain, deformation modes, cushioning coefficients, and other parameters demonstrate consistency with the experimental results. The acceleration, maximum compression, and cushioning coefficient obtained from the experiment and simulation calculation were 30.68 g, 15.44 mm, and 2.65, and 31.96 g, 14.91 mm, and 2.79, respectively. The results indicate that all error values were less than 5%, confirming the precision and reliability of the model. Furthermore, the model was utilized to simulate and predict the cushioning performance of honeycomb sandwich panels with different cell structures and paper thicknesses. These results provide a solid basis for enhancing the design of subsequent honeycomb element structures. [ABSTRACT FROM AUTHOR]

The stability, robustness, accuracy, and efficiency of space-time finite element methods crucially depend on the choice of approximation spaces for test and trial functions. This is especially true for high-order, mixed finite element methods which often must satisfy an inf-sup condition in order to ensure stability. With this in mind, the primary objective of this paper and a companion paper is to provide a wide range of explicitly stated, conforming, finite element spaces in four dimensions. In this paper, we construct explicit high-order conforming finite elements on 4-cubes (tesseracts); our construction uses tools from the recently developed 'Finite Element Exterior Calculus'. With a focus on practical implementation, we provide details including Piola-type transformations, and explicit expressions for the volumetric, facet, face, edge, and vertex degrees of freedom. In addition, we establish important theoretical properties, such as the exactness of the finite element sequences, and the unisolvence of the degrees of freedom. [ABSTRACT FROM AUTHOR]

This paper introduces a topology optimization approach based on the virtual element method (VEM), incorporating uncertainties. The objective of this optimization process is to design an optimal material layout for problems governed by linear elasticity equations to minimize the volume while satisfying probabilistic compliance constraints. The VEM is used to solve the boundary value problem in reliability-based topology optimization (RBTO). In the comparison between the VEM and the standard finite element method (FEM), a key difference emerges in the absence of explicitly defined shape functions tied to discrete degrees of freedom in VEM. Unlike FEM, VEM directly constructs the discrete bilinear form and load linear form without the need for computing shape function derivatives within the elements. This flexibility accommodates meshes with intricate geometries and arbitrarily shaped elements. The paper discusses the computational efficiency of VEM RBTO and explores the geometric impact of tessellations on converged topologies, demonstrating reduced susceptibility to checkerboard patterns compared to conventional quadrilateral elements. Additionally, the single-loop approach is examined, showcasing comparable accuracy to the first-order/second-order reliability methods (FORM/SORM) of RBTO using VEM. Numerical results for several problems that demonstrate the feasibility of the proposed method are presented. [ABSTRACT FROM AUTHOR]

This article presents the characteristics of four straw-filled beetle elytron plate (BEPsc) nonload-bearing walls, with their thermal and mechanical properties ascertained by the finite element method, focusing on structural parameters reported in previous studies and further optimized in this paper. The results are as follows: 1) The sequence of the four models' mechanical and thermal insulation properties with the reported structural parameters is given. The mechanical properties are mainly affected by the core layer's out-of-plane bending moment of inertia Ix. The insulation capability is mainly governed by the minimum cross section of the core concrete. Based on this result, a BEP with nonuniform cross-sectional honeycomb walls is proposed. 2) Further optimization is conducted on the two models with the weakest mechanical and thermal properties. The I-beam beetle elytron plate (IBEPsc) is proven most suitable for nonload-bearing walls. The potential for end-trabecular beetle elytron plate (EBEPsc) in hollow floors is expected. 3) As biomimetic models emerge constantly, it is of great significance to conduct parallel comparisons between new and old models to determine their performance levels and advantages. This paper provides an example for selecting suitable BEPs for applications and an inspiration for more work on screening outstanding biomimetic models. [ABSTRACT FROM AUTHOR]

Twinning is an important deformation mechanism for magnesium alloys, which has a significant impact on the texture evolution and mechanical properties. Pre-twinning deformation has been considered as an effective method for texture regulation. For each twinning mode, there are six crystallographically equivalent variants. Properly identifying the activated twin and understanding its variant selection criteria are essential for the development of high-performance magnesium alloys. As summarized in this paper, the observed twin variant by the commonly employed electron backscatter diffraction technique can be identified by several approaches, including misorientation analysis, trace analysis and matrix method. Schmid factor analysis was commonly performed to explain the selection of twin variants. To broaden the application scope under complex stress state, a generalized Schmid factor was derived by introducing a stress tensor. The efficiency of Schmid criteria to assess twin variant selection was confirmed to be influenced by the type of the applied stress. To consider the local effect on twin activation, in particular twin-twin transfer and slip-induced twinning, displacement gradient tensor calculation and geometrical compatibility factor analysis have been employed. It was demonstrated that local strain accommodation played a critical role in selecting the variants of cross-boundary twins in magnesium alloys. Assisted with crystal plasticity finite element modeling, the resolved shear stress on twinning and a composite Schmid factor combining the global Schmid factor and geometrical compatibility factor were obtained to better explain the activated twin variants in magnesium alloys. All the above-mentioned Schmid law based criteria and some energy based criteria as well are summarized in this paper. Their applications in evaluating twin variant selection in magnesium alloys are critically reviewed and discussed. [ABSTRACT FROM AUTHOR]

The main objective of this paper is to reduce the dimensionality of unknown coefficient vectors of finite-element (FE) solutions in two-grid (CN) FE (TGCNFE) format for the nonlinear unsaturated soil water flow problem by using a proper orthogonal decomposition (POD) and to design a new reduced-dimension iteration TGCNFE (RDITGCNFE). For this objective, a new time semi-discrete CN (TSDCN) scheme for the nonlinear unsaturated soil water flow problem is first designed and the existence, stability, and error estimates of TSDCN solutions are demonstrated. Subsequently, a new TGCNFE format for the nonlinear unsaturated soil water flow problem is designed and the existence, unconditional stability, and error estimates of TGCNFE solutions are demonstrated. Next, a new RDITGCNFE format with the same FE basis functions as the TGCNFE format is built by the POD method and the existence, unconditional stability, and error estimates of RDITGCNFE solutions are discussed. Ultimately, the rightness of theory results and the superiority of the RDITGCNFE format are verified by two sets of numerical tests. It is worth noting that the RDITGCNFE format differs completely from all previous reduced-dimension methods, including the authors' previous works. Therefore, the study of this paper can not only provide a new theoretical method for the dimensionality reduction of numerical models for nonlinear problems but also provide an algorithm implementation technology for the numerical simulation of practical engineering problems. [ABSTRACT FROM AUTHOR]

The space debris in Earth's orbit has increased drastically due to the failure of spacecraft, rocket bodies, and mission-related objects. These objects in orbit increase the space waste and challenge other flying objects such as spacecraft. A hypervelocity impact of space debris on spacecraft structures can have a range of effects (mechanical damage and functional failure), raising significant concerns about spacecraft safety. This paper reviews the different studies on the performance and development of the Whipple shield against the hypervelocity impact of space debris. The study focuses on the impact mechanism, dynamic Fragmentation of materials, strength models, Equation of state, characteristics, and model of the debris cloud. The strength models (Steinberg–Guinan and Johnson–Cook) and Mie–Gruneisen equation of state, primarily used for hypervelocity impact applications, are thoroughly covered in this study. The study also reported the various experimental and numerical techniques for high and hypervelocity impact. The study concluded that mesh-based, mesh-free, and hybrid finite element methods are reliable for analyzing Whipple shield targets to resist hypervelocity impact. The study also observed that the two-stage light gas gun technique investigates most experimental analyses of hypervelocity impact on the Whipple shield. Alongside reviewing the abovementioned aspects, this paper also underlines the future scope of study in this paradigm. The authors strongly believe that this study provides more insights into the fundamentals and perceptions of the Whipple shield to protect the spacecraft against the hypervelocity impact of space debris. [ABSTRACT FROM AUTHOR]

A novel higher-order multi-scale computational approach for efficient nonlocal gradient elasticity simulation of composite materials is proposed in this paper. Based on the use of decoupling technique, raw fourth-order nonlocal gradient equations with periodically oscillatory coefficients are decoupled as the new second-order multi-scale equations for reducing the continuity requirements when implementing multi-scale simulation. Next, a novel macro-micro coupled multi-scale computational model with higher-order corrected terms is rigorously devised for high-resolution multi-scale and nonlocal analysis of composite materials via multi-scale asymptotic analysis. Then, the local error analysis of higher-order multi-scale solutions in the point-wise sense illustrates that establishing higher-order asymptotic corrected terms plays a vital role in investigating the micromechanical mechanisms. Moreover, a rigorous global error estimation with an explicit rate of higher-order multi-scale solutions is first derived in the energy norm sense. In addition, an efficient multi-scale numerical algorithm is presented to effectively simulate nonlocal gradient elasticity problems of composite materials based on finite element method (FEM). And we also derive the convergence estimation of the proposed numerical algorithm. Finally, numerical experiments are conducted to gauge the efficiency and accuracy of the proposed higher-order multi-scale method. This paper offers a unified higher-order two-scale computational framework for enabling nonlocal gradient elasticity behavior simulation of composite materials. [ABSTRACT FROM AUTHOR]

Due to the influence of many factors, distribution cables are often densely placed at the bottom of the cable trench. As a result, it is easy for distribution cables to become the thermal bottleneck of the whole transmission line. To address this dilemma, this paper establishes a finite element simulation model of a cable trench to analyze the hot spots of cables with different arrangements in the cable trench. Then, the model's accuracy is verified based on real temperature rise experiments. For an arrangement with overheating risk, the ampacity improvement method of filling the cable trench with high-thermal-conductivity material was proposed, and the ampacity improvement effect under different filling ratios was assessed. Finally, combined with the analysis of economic benefit and cost, the method of determining the optimal filling ratio was used, and the impact resistance of the cables under the impact of new energy load was analyzed. The results indicate that, for the case of the optimal filling ratio, the cables in the dense cable trench showed superior impact resistance. The investigations in this paper make significant contributions to the promotion of the maximum utilization of cables. [ABSTRACT FROM AUTHOR]

Many countries have a variety of offshore and onshore wind turbines that face extreme aging challenges. Issues with harmful vibrations that must be minimized are addressed in this paper. A new method of wind turbine tower vibration control using piezoelectricity and shunt circuits is proposed in this paper. The passive vibration control method is shown to improve the tower's structural performance under various environmental loads, like wind and seismic excitations. To examine the effectiveness of the suggested shunted piezoelectric system, a simple surrogate finite element model of a wind turbine tower is considered, and various investigations at the second eigenfrequency are carried out. An alternative way of modeling the studied structure is considered and the results demonstrate better performance. The advantages of setting up structural damping systems for decreasing tower vibrational loads and boosting their structural stability and resilience against extreme events are highlighted throughout this work. [ABSTRACT FROM AUTHOR]

The bottom platform is an important underwater sensor that can be used in communications, early warning, monitoring, and other fields. It may be affected by earthquakes, winds, waves, and other loads in the working environment, causing changes in posture and affecting its sensing function. Therefore, it is of practical engineering significance to analyze the force conditions and posture changes in the bottom platform. In order to solve the problem of postural stability of the underwater bottom platform, this paper establishes a fluid and structural simulation model of the underwater bottom platform. First, computational fluid dynamics (CFD) technology is used to solve the velocity distribution and forces in the watershed around the bottom platform under a 3 kn ocean current, where the finite element method (FEM) numerical calculation method is used to solve the initial equilibrium state of the bottom platform after it is buried. On this basis, this paper calculates the forces on the bottom platform and the posture of the bottom platform at different burial depths under the action of ocean currents. Additionally, the effects of different burial depths on the maximum displacement, deflection angle, and postural stability of the bottom platform are studied. The calculation results show that when the burial depth is greater than 0.6 m, and the deflection angle of the bottom platform under the action of the 3 kn sea current is less than 5°, the bottom platform can maintain a stable posture. This paper could be used to characterize the postural stability of underwater bottom platforms at different burial depths for the application of underwater sensors in ocean engineering. [ABSTRACT FROM AUTHOR]

To improve the utilization rate of cable lines while retaining sufficient security, the accurate thermal assessment of cable is significant for cable operation condition evaluation. The thermal rating for a cable joint, which is regarded as a hot spot of cable lines, is not covered by the scope of IEC 60287. While the existing publications for cable joint thermal evaluation also have some limitations. In this paper, the quasi-three-dimensional thermal model of the cable joint was established and the iterative solution method for the model is presented. Based on the model, an improved thermal rating method for the cable joint was proposed, which was implemented with monitored surface temperature and load data. The improved method was verified by the finite element method and the results showed an error of less than 5%. The superiority of the improved method was conducted by the comparison between the previously published method and the improved method. The improved method showed a better accuracy than the previously published method. The proposed method in this paper can be complementary to the IEC method, and is easy to use for the operating evaluation of cable joints in the field with the on-line condition monitoring technology. [ABSTRACT FROM AUTHOR]

The wear of backup rolls will have a great impact on the quality of the shape of hot rolled strip sheet. In order to overcome the limitations of the finite element method (FEM) in calculating backup roll wear in terms of efficiency and accuracy, this paper proposes a tandem FEM + ML hybrid model to optimise the predictive effect of the finite element method (FEM) on backup roll wear. Firstly, a backup roll wear model based on FEM is established. Secondly, in order to select the optimal machine learning (ML) algorithm as the finite element error compensation model, three types of finite element error compensation models were established based on the random forest (RF) algorithm, the radial basis function (RBF) neural network algorithm, and the particle swarm optimisation support vector machine (PSO-SVM) algorithm. Finally, the three types of finite element error compensation models were connected in series with the FEM model to compare the prediction performance of the three types of FEM + ML models on backup roll wear. The numerical experimental results show that the FEM + PSO-SVM model can better predict the wear of the backup roll, and the PSO-SVM algorithm is the most suitable for building the finite element error compensation model. It is proved that the FEM + ML model proposed in this paper can effectively improve the accuracy and computational efficiency of the FEM model for predicting backup roll wear without adding microelements. In addition, among the hot rolling parameters, the rolling force has the greatest influence on the backup roll wear, and excessive rolling force for a single pass should be avoided to slow down the backup roll wear. [ABSTRACT FROM AUTHOR]

In this paper, a ground-based experimental system for solar array active vibration suppression has been established. Firstly, in order to establish an accurate model of the solar array, the solar array is regarded as a flexible cantilevered thin plate, and the corresponding dynamical equations are derived using the absolute nodal coordinate method. In addition, in this paper, the more advanced MFC piezoelectric patch is used instead of the traditional PZT piezoelectric ceramic patch. The electromechanical coupling finite element model of the P1-type MFC patch is established and substituted into the kinetic equation of the solar array. Finally, the accuracy of the electromechanical coupling modeling and the control effect of active vibration suppression were verified using the PID control. A set of experimental frameworks for evaluating the active vibration suppression effect, including the free vibration test, sinusoidal perturbation test, and white noise perturbation test, as well as the analysis strategy of the test data, are established. [ABSTRACT FROM AUTHOR]

This paper has proposed an efficient analysis method to calculate interface stress intensity factors (SIFs) based on a proportional stress field of a reference problem whose exact solution is available. In the previous proportional methods, the same crack length and the same element size were applied to both reference and unknown problems so that the same FEM error can be produced. Therefore, when analyzing many unknown problems, the conventional method needs to analyze many reference problems at the same time. Since this approach is time-consuming, this paper considers how to calculate many crack lengths efficiently by using only one single reference solution modeling. For this purpose, several general relations of SIFs are derived for the unknown and the reference problems when both crack length and element size are different. To analyze many unknown problems accurately by using a single reference solution modeling, how to choose the most suitable element dimension of the reference model is clarified. The proposed method is especially useful for crack propagation analysis. [ABSTRACT FROM AUTHOR]

Kalman filter has emerged as a powerful tool for unbalance identification in rotating machinery. Recently, the augmented Kalman filter combined with the finite element model has grown up and projects its potential for complex rotor systems. This paper investigates the application of the augmented Kalman filter (AKF) to a practical turbofan engine. The current study reveals that using steady-state responses as measurements can cause fluctuation in the estimated results, even divergence for some cases, while the available signals in practice are steady-state responses generally. To the authors' knowledge, this practical problem is revealed for the first time. To address the problem, the convergence criterion is employed to improve the AKF and formulates the adaptive fading augmented Kalman filter (AFAKF) proposed in this paper. Results indicate that the increase of the amplification factor, the insufficient measurement points, and the complexity of the dynamic model can all lead to the deterioration of the estimated unbalance. The proposed AFAKF method shows favorable convergence and can achieve accurate estimation with less than 5% relative errors, and the superiority over AKF in computation cost is also observed. [ABSTRACT FROM AUTHOR]

This paper presents a forced vibration analysis of open doubly curved sandwich panels subjected to a moving constant force. In this paper, the effect of softness of the core is considered by implementing a semi-layerwise theory. To this aim, the first-order shear deformation theory is adopted for the face sheets and a higher-order theory which was obtained based on 3D elasticity theory is considered for the core. The presented formulation is general and as the deepness parameter is accounted in the strain–displacement relations, the formulation can be used for a wide range of deep as well as shallow doubly curved shells. To obtain the dynamic response of the system, the finite element method (FEM) along with the Newmark method is used. The proposed element is a higher-order one with nine nodes and each node has fifteen degrees of freedom. The effect of various parameters such as length-to-thickness ratio, in-plane aspect ratio, boundary conditions, lamination scheme, and fiber orientation angles on the dynamic response of the structure is examined. Additionally, the critical velocity of the force at which the structure experiences maximum dynamic deflection is obtained in each case. The results show that as the length-to-thickness ratio of the structure increases, the dynamic magnification factor curve increases with respect to non-dimensional velocity. This study provides insights into the dynamic behavior of doubly curved sandwich panels with soft cores and can aid in the design of such structures for specific applications. The results of this study can also serve as a benchmark for future studies on the forced vibration behavior of doubly curved sandwich panels. [ABSTRACT FROM AUTHOR]

This paper introduces a novel index for static eccentricity (SE) fault diagnosis in synchronous reluctance motors (SynRMs). Although SynRMs with rotor barriers under SE have been modeled in a few papers, any indices for the fault detection have not yet been reported. The proposed index is a specific frequency pattern in the motor current, which can also determine the fault severity. A novel analytical magnetic field investigation is applied to ascertain the proposed index. An accurate nonlinear numerical method is proposed for the motor inductances calculation. The model considers the rotor flux barriers, magnetic saturation, and the stator slots effect. The air gap flux density and the motor current are then achieved. The fast Fourier transform is exploited as a signal-processing tool to calculate motor current spectra. Then, a two-dimension time-stepping finite element method is used to attest of the proposed numerical model. Effectiveness of the proposed index is verified by simulation and experimental tests. [ABSTRACT FROM AUTHOR]

Determining the mechanical properties of a composite material represents an important stage in its design and is generally a complicated operation. These values are influenced by the topology and geometry of the resulting composite and the values of the elastic constants of the components. Due to the importance of this subject and the increasing use of composite materials, different calculation methods have been developed over the last fifty years. Some of the methods are theoretical, with results that are difficult to apply in practice due to difficulties related to numerical calculation. In the current paper, using theoretical results offered by the homogenization theory, values of engineering elastic constants are obtained. The finite element method (FEM) is used to determine the stress and strain field required in these calculations; this is an extremely powerful and verified calculation tool for the case of a material with any type of structure and geometry. In order to minimize errors, the paper proposes the method of least squares, a mathematical method that provides the best estimate for the set of values obtained by calculating FEM. It is useful to consider as many load cases as possible to obtain the best estimates. The elastic constants for a transversely isotropic material (composite reinforced with cylindrical fibers) are thus determined for a real case. [ABSTRACT FROM AUTHOR]

In order to further study the cushioning performance of concave hexagonal cores (CHCs) and expand their application range, the in-plane finite element model of CHCs is established in this paper. A dynamic cushioning coefficient method was proposed to characterize the cushioning performance of CHCs. The dynamic cushioning coefficient curve and minimum dynamic cushioning coefficient (MDCC) of CHCs with different impact velocities and structural parameters are obtained. The influence rules of structural parameters and impact velocities on the MDCC are analyzed; the deformation mode and transformation empirical formula are also obtained. The results show that when other parameters are constant, the MDCC of CHCs decreases with the increase of impact velocity, increases with the increase of wall thickness and side length ratio, and decreases with the increase of expansion angle. The theoretical analysis is consistent with the finite element results, which further verifies the reliability of the model. This paper provides a solid theoretical basis for the industrial application of the cushioning performance of CHCs and forms a key technical support. [ABSTRACT FROM AUTHOR]

For insulated gate bipolar transistor (IGBT) modules using wire bonding as the interconnection method, the main failure mechanism is cracking of the bonded interface. Studying the mechanical properties of the bonded interface is crucial for assessing the reliability of IGBT modules. In this paper, first, shear tests are conducted on the bonded interface to test the bonded interface's strength. Then, finite element–cohesive zone modeling (FE-CZM) is established to describe the mechanical behavior of the bonded interface. A novel machine learning (ML) architecture integrating a convolutional neural network (CNN) and a long short-term memory (LSTM) network is used to identify the shape and parameters of the traction separation law (TSL) of the FE-CZM model accurately and efficiently. The CNN-LSTM architecture not only has excellent feature extraction and sequence-data-processing abilities but can also effectively address the long-term dependency problem. A total of 1800 sets of datasets are obtained based on numerical computations, and the CNN-LSTM architecture is trained with load–displacement (F–δ) curves as input parameters and TSL shapes and parameters as output parameters. The results show that the error rate of the model for TSL shape prediction is only 0.186%. The performance metric's mean absolute percentage error (MAPE) is less than 3.5044% for all the predictions of the TSL parameters. Compared with separate CNN and LSTM architectures, the proposed CNN-LSTM-architecture approach exhibits obvious advantages in recognizing TSL shapes and parameters. A combination of the FE-CZM and ML methods in this paper provides a promising and effective solution for identifying the mechanical parameters of the bonded interfaces of IGBT modules. [ABSTRACT FROM AUTHOR]

The reliability of high-voltage circuit breakers (HVCBs) depends critically on the dynamic characteristics of their hydraulic operating mechanisms (OMs). However, previous analyses have been limited to discrete components due to the lack of models capturing multi-physics couplings. This paper proposes a novel co-simulation approach for hydraulic OMs that enables system-level analysis. Specifically, the hydraulic system is modelled using the lumped parameter method, and the transmission mechanism is modelled via finite element analysis. The subsystem models exchange inputs and outputs through shared memory, realizing coupling without compromising simulation performance. Experiments on a 550 kV HVCB validate that the model accurately captures OM transient responses. The co-simulation framework has diverse capabilities, including predicting dynamics under varied operating parameters and quantifying stress distributions and response evolution patterns. Overall, this paper proposes an advanced multi-physics analysis capability for hydraulic OMs to support reliabilityoriented design and optimization. [ABSTRACT FROM AUTHOR]

The injection of the Direct Current (DC) bias component will lead to oversaturation of the transformer core and affect its safe operation. In this paper, the vector magnetic potential A and winding current I are used as the variables solved in the time-step finite element method to simulate the nonlinear transient magnetization characteristics of the transformer under DC bias, and the fixed-point reluctivity is introduced to deal with the nonlinear iterative problem. The field-circuit coupling relationship is established through the electromagnetic induction law, and the two-dimensional time-step finite element model of the transformer core is found. The dynamic J–A hysteresis model is used to deal with the hysteresis effect of iron cores. The fixed-point reluctivity of difference form is introduced to solve the problem of uneven transition and multi-value of magnetization curve with hysteresis in the iterative process. The determination scheme for the equivalent thickness of the core model and the initial iterative value are discussed. A laminated core transformer is made for experiments, and the accuracy and applicability of the model proposed in this paper are verified by comparing the calculation and experimental results. The influence of different voltages and different DC bias components on the electromagnetic characteristics of the transformer is analyzed. [ABSTRACT FROM AUTHOR]

In order to solve the deficiencies of traditional board-level assembly CubeSat, such as low modularity, complex assembly process and poor scalability, a sandwich module assembly structure and method with no drive components is proposed, which can achieve quick connection and flexible assembly between modules. First, the proposed CubeSat standard module can be compatible with the standard interface of traditional CubeSat. The 4-in-1 interface of Mechanical-Electrical-Thermal-Data does not occupy the internal space of the standard module, and the electrical and thermal interfaces meet the requirements of CubeSat in-orbit work. Secondly, in order to improve the reliability of the connection between the modules, a plug-in connection scheme of the mechanical interface was designed. Finally, a genetic algorithm was used to optimize the interface size, so that it has better plug-in performance. The simulation results show that the performance of the male interface of titanium alloy material is the best, and the error between the finite element analysis and genetic algorithm results is less than 4.8 %, which is consistent with the actual tensile test results of the optimized prototype, and meets the mechanical requirements of the launch process and space orbit operation for CubeSat. The test results can guide production and provide a reference for subsequent engineering applications of CubeSat modular assembly. • In this paper, a standardized modular assembly method of CubeSat sandwich structure is designed. The assembly structure does not include driving parts, is simple in form, can be quickly connected, is flexible in assembly mode, and does not occupy the internal space of the standard module. • A standard 4-in-1 connection device was designed in this paper, including four interfaces: mechanical, electrical, thermal and data, which can achieve "plug and play". • This paper uses genetic algorithm to optimize the structure of the mechanical connection male joint, and obtains the optimal configuration of the male joint. This ensures the reliability, firmness and stability of mechanical connections, and meets the mechanical requirements of CubeSats during launch and in orbit. [ABSTRACT FROM AUTHOR]

In this paper, a software framework is presented through an application that is able to jointly optimize an electric motor and a gearbox for the design of a drive system for electric vehicles. The framework employs a global optimization method and uses both analytical and finite element method (FEM) models to evaluate the objective functions. The optimization process is supported by a statistical surrogate model, which allows a large reduction of runtime. An earlier version of this framework was only suitable for electric motor optimization. In the application presented in a previous paper, the motor of a belt-driven electric drive system was optimized. In this paper, the optimization of the same drive system is shown, but now with a combined optimization of a gear drive and motor. The objective functions of optimization are minimizing the total loss energy and the weight of the drive system. The optimization results are compared with previous results to demonstrate the further potential of joint optimization. [ABSTRACT FROM AUTHOR]

In electron optics, calculation of the electric field plays a major role in all computations and simulations. Accurate field calculation methods such as the finite element method (FEM), boundary element method and finite difference method, have been used for years. However, such methods are computationally very expensive and make the computer simulation challenging or even infeasible when trying to apply automated design of electrostatic lens systems with many free parameters. Hence, for years, electron optics scientists have been searching for a fast and accurate method of field calculation to tackle the aforementioned problem in the design and optimization of electrostatic electron lens systems. This paper presents a novel method for fast electric field calculation in electrostatic electron lens systems with reasonably high accuracy to enable the electron-optical designers to design and optimize an electrostatic lens system with many free parameters in a reasonably short time. The essence of the method is to express the off-axis potential in an axially symmetrical coordinate system in terms of derivatives of the axial potential up to the fourth order, and equate this to the potential of the electrode at that axial position. Doing this for a limited number of axial positions, we get a set of equations that can be solved to obtain the axial potential, necessary for calculating the lens properties. We name this method the fourth-order electrode method because we take the axial derivatives up to the fourth order. To solve the equations, a quintic spline approximation of the axial potential is calculated by solving three sets of linear equations simultaneously. The sets of equations are extracted from the Laplace equation and the fundamental equations that describe a quintic spline. The accuracy and speed of this method is compared with other field calculation methods, such as the FEM and second order electrode method (SOEM). The new field calculation method is implemented in design/optimization of electrostatic lens systems by using a genetic algorithm based optimization program for electrostatic lens systems developed by the authors. The effectiveness of this new field calculation method in optimizing optical parameters of electrostatic lens systems is compared with FEM and SOEM and the results are presented. It should be noted that the formulation is derived for general axis symmetrical electrostatic electron lens systems, however the examples shown in this paper are with cylindrical electrodes due to the simplicity of the implementation in the software. [ABSTRACT FROM AUTHOR]

This paper presents a computer-based methodology to support the design for additive manufacturing of metal components. Metal additive manufacturing, and in particular powder bed fusion systems, are playing a prominent role in the industry 4.0 scenario. The state of the art concerning design methods and tools to support design for additive manufacturing is reviewed by the authors. The key phases of product design and process design to achieve lightweight functional designs and reliable processes are deepened, and the computer-aided technologies to support the approaches implementation are described. Indeed, the state of the art design for additive manufacturing general workflow can be enriched by holistic approaches, use of numerical simulation, and integration and automation between the required tasks. The paper provides a methodology based on the systematic use of numerical simulation to achieve the optimization of both products and associated processes. To take advantage of the holistic perspective, the approach relies on the use of integrated product-process design platforms, allowing to streamline the digital process chain. Product design is based on the systematic integration of topology optimization and automatized tools for concept development and selection and subsequent product simulation driven design refinement. Process design is based on a systematic use of process simulation to prevent manufacturing flaws related to the high thermal gradients of metal processes and minimize residual stress and deformations. This is achieved by working on both the build cycles layouts and the 3D models' distortion compensation. An automotive use case of product and process design performed through the proposed simulation-driven integrated approach is provided to assess the actual method suitability for effective re-designs of additive manufacturing high-performance metal products. The bridged gaps are systematically outlined, and further developments are discussed. [ABSTRACT FROM AUTHOR]

Based on the classical grid theory and related regulations, a structure model of a fiber-wound composite gas cylinder was designed in this paper. Based on the design results, a finite element model of a fully wound composite cylinder of an aluminum alloy inner liner with a working pressure of 35 MPa was established based on the ABAQUS software, and its stress distribution under working pressure and minimum burst pressure was analyzed. According to engineering experience, the pressure tolerance of composite cylinders can be improved by proper autofrettage pressure before working pressure, so the influence of autofrettage pressure was analyzed in this paper. The optimum autofrettage pressure was selected by setting the autofrettage gradient, and damage analysis was carried out on the cylinder with nominal working pressure of 35 MPa based on the Hashin failure criterion. The results show the initial damage sequence: matrix stretching occurs before the fiber stretching, and the damage generally starts from the spiral-wound layer. The tensile damage first appears in the transition section between the head and the barrel body, and the damage of the spiral-wound layer develops from the inner layer of the wound layer to the outer layer, while the damage of the circumferentially wound layer develops from the outer layer to the inner layer. [ABSTRACT FROM AUTHOR]

Strip shape control is a hotspot and challenge in strip rolling, where the development trend of rolling technology is towards high strength, high toughness, and a large width-to-thickness ratio. The influence of material microstructure evolution on strip shape control is being increasingly emphasized. In this paper, a Nb-Ti microalloyed steel is taken as the research object. Thermodynamic and kinetic models focusing on the precipitation of the austenite phase are established to quantify the precipitation process. A coupled model of rolls and strips is built using ABAQUS 2022 software, where the precipitation strengthening model and high-temperature constitutive model are embedded into the finite element model (FEM) through subroutines. A two-dimensional alternating differential model is employed to acquire real-time temperature differences in the width direction of the strip. The effects of precipitation inclusion and exclusion on the strip crown under different operating conditions are compared and analyzed. The results indicate that as the temperature decreases, the strengthening effect increases, reaching around 40 MPa at temperatures above 1000 °C and 96.6 MPa at 800 °C. Furthermore, the inclusion of crown in the precipitation consideration is more sensitive to overall temperature changes, but as the strip width decreases, the sensitivity of crown to temperature decreases. The research findings of this paper provide guidance for improving strip shape control and reducing abnormalities during the rolling process. [ABSTRACT FROM AUTHOR]

Determining the properties of composite materials (knowing the properties of the component phases) is a primary objective in the design phase. Numerous methods have been developed to determine the elastic constants of a composite material. All these methods are laborious and require significant computing time. It is possible to make experimental measurements, but these too are expensive and time-consuming. In order to have a quick estimate of the value of the engineering constants of a new composite material (in our study a polymeric matrix reinforced with carbon fibers), this paper proposes a quick method for determining the homogenized material constants, using the finite element method (FEM). For this, the eigenfrequencies of a beam specimen manufactured by the studied composite material will be computed using FEM. The model will consider both phases of the composite, with the geometry and real size. The mechanical properties of the constituent's material phases are known. With the help of this model, the torsional, longitudinal and transverse vibrations of the beam are studied. Based on the eigenvalues obtained by this calculation, it now is possible to quickly estimate the values of homogenized material constants required in the design. An example for a fiber-reinforced polymer composite material is provided in the paper. [ABSTRACT FROM AUTHOR]

This paper is devoted to a study of a history-dependent mixed hemivariational-variational inequality arising in contact mechanics. The contact problem concerns the deformation of a viscoelastic body with long memory, subject to a general friction law on one part of the boundary and a frictionless Signorini condition on another part of the boundary. The solution existence and uniqueness of the history-dependent mixed hemivariational-variational inequality based on a Lagrange multiplier approach are proved. Then, a fully discrete scheme is introduced and studied. The trapezoidal rule is used to approximate the integral in the history-dependent operator. For the spatial discretization, the linear finite elements are used to discretize the displacement field, and dual basis functions are used in the approximation of the Lagrange multiplier. Optimal order error estimates are derived for the displacement and the Lagrange multiplier under appropriate solution regularity assumptions. Numerical results are presented to illustrate the theoretical prediction of the convergence orders. [ABSTRACT FROM AUTHOR]

The study investigated the fretting fatigue behavior of a bridge‐flat component, utilizing a flat specimen with two bridge‐type pads tested under cyclic loading. Numerical simulations were also carried out using the finite element method. Firstly, precise calculations of contact stress were performed with the material's elastoplastic model. Secondly, the stress distribution on the maximum stress cross‐section was extracted, and the stress gradient was computed. Finally, a stress fatigue life prediction model considering the stress gradient was proposed, and the predicted life was compared with the test life. The results of the study revealed that the specimens exhibited high stress gradients at the edges of the contact region during the bridge‐flat tests. In comparison with plain fatigue, fretting fatigue life was significantly reduced. The life prediction model proposed in this paper demonstrated good accuracy for predicting fretting fatigue. [ABSTRACT FROM AUTHOR]

In this paper, we prove the existence of weak solution and the uniqueness of strong solution to a Voigt-regularization of the three-dimensional thermally coupled inviscid, resistive MHD equations. We also propose a fully discrete scheme for the considered problem, which is proven to be stable and convergent. All computational results support the theoretical analysis and demonstrate the effectiveness of the presented scheme. [ABSTRACT FROM AUTHOR]

• An exact matrix-form solution is proposed for analyzing multiple non-circular tunnels in layered media. • Generalized complex variable theory, conformal mapping and fast Fourier transform are used in derivation. • The proposed solution avoids complicated coefficient determinations and iterative algorithm. • Effects of elastic properties, tunnels' position, equivalent layer and the layer thickness surrounding tunnels are discussed. This paper presents a matrix-form complex variable method that can predict the stress and displacement field around multiple tunnels in layered media. The proposed matrix-form method avoids complicated coefficient determinations and iterative algorithm, which effectively improves the computation efficiency. Firstly, the expressions for the stress and displacement field are given under the generalized complex variable theory. The conformal mapping technique and fast Fourier transform algorithm are used to formulate the matrix-form expressions of boundary conditions. Considering the boundary conditions of the tunnels and soil interfaces, the global matrix solution for multiple tunnels in layered media are obtained by assembling the matrix solutions of all single layers. The present theory is verified by comparison with the existing literature, finite element analysis and field data in engineering practices, and several numerical examples are given to reveal the effects of elastic properties, tunnels' position, equivalent layer and the layer thickness surrounding tunnels. [ABSTRACT FROM AUTHOR]

• A hybrid method by integrating the advantages of frequency-wavenumber method and finite element method. • The entire process simulation from the source to the local terrain has been achieved. • Seismic response of mountains due to point dislocation source. • The crustal velocity structure and source parameters influence the seismic response. This paper proposes a hybrid method to analyze seismic responses from the source to complex terrain. The hybrid method is based on the domain reduction idea, dividing the overall model into global domain (including the source and propagation path) and local domain (including complex terrain). For the propagation of seismic wave fields in the global domain, a frequency-wavenumber method based on the exact dynamic stiffness matrix is adopted. For the local domain, the finite element method combined with the artificial boundary substructure method is used to analyze the seismic response of complex terrain. The advantage of this hybrid method is that the frequency-wavenumber method can accurately and efficiently calculate the three-component ground motion of the global large-scale model, while the finite element method has a variety of elements and constitutive models, which can finely simulate the local complex terrain. We take the Zigong Xishan Park site as the research object, and verify the correctness of the hybrid method. Based on this, we use the hybrid method to simulate the seismic response of the site, focusing on the propagation process of seismic waves on the mountain surface, the spatial distribution of peak ground displacement, and the effects of crustal velocity structure and source parameters on the seismic response of the mountain. The results show that the propagation of seismic waves and the distribution of peak ground displacement are closely related to changes in terrain, and the crustal velocity structure and source parameters have significant effects on the seismic response of the mountain. [ABSTRACT FROM AUTHOR]

• This manuscript firstly establishes a half vehicle-pavement-continuous beam bridge (VPCB) coupled system analytical model. • It is innovatively that the bridge and vehicle initial conditions are considered into analytical model. • The pavement road surface local deterioration is simulated by increasing the local roughness level of the road surface. At present, the highway bridge field is in a period of vigorous development, and the application of continuous beam bridges is becoming more extensive, so the related safety research is particularly important. However, there are few studies on the vehicle-bridge coupled interaction of continuous beam bridges, and the influence of pavement is not considered. In order to ensure that the theoretical research can better serve the engineering practice, it is urgent to establish a more refined and more complex vehicle-bridge coupled model. In this paper, a vehicle-pavement-continuous beam bridge (VPCB) coupled system model is established. A half-vehicle model is adopted, the pavement is simulated by a continuous and uniform spring damper, and the bridge is a three-span continuous beam bridge. The motion equation of VPCB system is established and deduced using D 'Alembert principle and modal superposition method. The dynamic responses of VPCB system are calculated by Newmark- β method. Here, the continuous beam bridge mode functions are fitted based on the modal vector of the continuous beam calculated by the finite element model. The correctness of the method is verified by existing literatures. The effects of bridge pavement, vehicle speed, bridge second span length, the pavement equivalent spring stiffness coefficient k p and damping coefficient c p , road roughness on VPCB system are studied. The results show that the pavement can reduce the bridge dynamic amplification factor (DAF) and the vehicle acceleration to improve the ride comfort. The dynamic responses of the vehicle and the bridge contain each other's natural frequencies. The responses of VPCB system increase with the vehicle velocity, the vehicle responses are more affected than the bridge responses. The pavement equivalent stiffness and damping coefficient has little effect on the system. The responses of vehicle and bridge increase with the increase of the bridge second span length. When the vehicle travels each span of the three-span continuous beam bridge, the time history curve of the pavement deformation under the tire resembles a sine wave. When the local roughness level of the road surface increases, the system responses at the deteriorated position suddenly increases, and the responses of the subsequent position will also be affected. The larger the roughness grade, the more obvious the responses increase. The effect on the system responses of the whole pavement damage is more significant than the local pavement damage. [ABSTRACT FROM AUTHOR]

This paper presents an analysis of stress developed in flex-gear (FG) cup with conventional strain wave generating (SWG) cam and split-cam SWG in strain wave gear/harmonic drive (HD), using initial load distribution by tooth pairs as a boundary condition in finite element method (FEM). The initial load distribution by tooth pairs is achieved with the assembly of ring gear (RG) and deformed FG before the application of any external load. A closed-form solution to the analytical method describes the exact initial load distribution by tooth pairs in the presence of secondary contacts. However, with spit-cam SWG the phenomenon of secondary contact is completely eradicated. Therefore, the load is distributed by contacting tooth pairs that lie only in the line of action and the tooth pair that lies along the major axis of split-cam SWG experiences maximum load compared to conventional SWG. The strain developed in FG cup based on the load distribution with both SWG cams obtained from FEM is verified with the experiment. The comparison shows a rise of 44 microns in circumferential strain around the major axis of deformed FG, when RG is assembled. This testifies the presence of initial load, which is distributed among all interfering tooth pairs before applying external load. [ABSTRACT FROM AUTHOR]

This article introduces a unified mixed finite element framework based on a saddle-point formulation that applies to time-dependent fourth order linear and nonlinear problems with clamped, simply supported, and Cahn-Hilliard type boundary conditions. The classical mixed formulations lead to large matrix systems that demand huge storage and computational time making the schemes expensive, especially for the time-dependent problems. The proposed scheme circumvents this by employing biorthogonal basis functions that lead to sparse and positive-definite systems. The article discusses a mixed finite element method for the biharmonic problem and the time-dependent linear and nonlinear versions of the extended Fisher-Kolmogorov equations equipped with the aforementioned boundary conditions. The wellposedness of the scheme is discussed and a priori error estimates are presented for the semi-discrete and fully discrete finite element schemes. The numerical experiments validate the theoretical estimates derived in the paper. [ABSTRACT FROM AUTHOR]

The structural dynamic response reconstruction of large-scale civil engineering structures using traditional methods results in heavy computational burden with low efficiency. To improve efficiency, a novel time-domain response reconstruction method based on model condensation and modal decomposition is proposed in this paper. The finite element model of the target structure is transformed and condensed. The system matrices of the condensed model have lower dimensionality than those of the original finite element model. A checking procedure based on the modal assurance criterion is developed to ensure the accuracy of the condensed system. Measured dynamic responses are decomposed to time-domain modal responses at measured degrees of freedom (DOFs) using band-pass filters. The mode shape matrix of the condensed model is used to calculate the modal responses at the unmeasured DOFs. Following modal superposition and coordinate transformation, the dynamic responses at unmeasured DOFs are reconstructed without knowledge of the excitation. The accuracy and efficiency of the proposed methodology are verified in a numerical case study of a three-dimensional building frame model and an experimental case study of a two-story steel frame under laboratory conditions. In contrast with the numerical case study, the effects of modeling error and measurement noise are considered in the experiment verification. The results show the outstanding performance of the proposed methodology when compared with a traditional response reconstruction method. [ABSTRACT FROM AUTHOR]

With wind farms expanding into seismically active areas, wind turbines face the combined hazard of wind and earthquakes. Strong winds often accompany strong earthquakes, posing a significant threat to wind turbine safety. When subjected to multihazard loads, ignoring the correlation between loads could underestimate the risk to the structure. However, there is limited research on wind turbines under the combined effect of wind and earthquakes. In this paper, based on historical data considering the correlation between wind and earthquakes, a failure probability analysis method for wind turbine towers under the joint action of wind and earthquakes is presented. To accurately simulate the mechanical behavior of lattice-type wind turbine towers, the bolt slippage effect is taken into account. A finite element model of a lattice wind turbine tower is established, and the failure probability of the tower under the joint action of wind and earthquakes is calculated. The study also analyzes the influence of the bolt slippage effect and load correlation on the reliability of lattice-type wind turbine towers. The results indicate that ignoring load correlation can underestimate the failure probability of the structure, thus decreasing its safety. Similarly, ignoring the bolt slippage effect can also underestimate the failure probability of the structure and seriously threaten the safety of wind turbines. [ABSTRACT FROM AUTHOR]

Architectural heritage health assessment is the basis of scientific repair and maintenance. However, existing methods do not adequately take into account the fuzziness, randomness and uncertainties unique to architectural heritage assessment. In this paper, a new evaluation model of VM-NCM is constructed by combining variable weight theory and normal cloud model theory. The model enables the combination of qualitative ratings and quantitative calculation, deals with the fuzziness in the assessment process, and resolves the randomness and reflects the uncertainty to a certain extent. Based on constructing the index system combining qualitative and quantitative indexes, the structural index values are acquired by the synergistic coupling of the fine laser point cloud model and finite element structural analysis model. The acquisition of surface index values is completed by the hyperspectral intelligent detection technology of surface materials and diseases. These reduce the generation of ambiguous information in the index detection process. An evaluation study is conducted using the Yingxian wooden pagoda in China as an example. The results show that this method takes into account the fuzziness and randomness in the evaluation process, and obtains more scientific and reliable evaluation results, which provides a research paradigm for assessing the architectural heritage health status. [ABSTRACT FROM AUTHOR]

Elastic deformation in thin-walled parts during machining is affected by the coupling of force and flexibility. Obtaining flexibility information along machining tool paths is crucial for online monitoring of this deformation. However, the current finite element method (FEM) is limited by mesh nodes, hampering its ability to accurately determine flexibility along tool paths. To overcome this limitation, this paper proposes a method that combines FEM with the surrogate model to predict flexibility accurately at any position on thin-walled parts' surfaces. The surrogate model is the hybrid model K-K-CNN based on two K-nearest neighbor (KNN-KNN) algorithms and a convolutional neural network (CNN) model. Initially, an initial dataset containing positions and flexibility of mesh nodes is generated automatically through secondary development of ABAQUS. Then, the K-K-CNN hybrid model is introduced and trained on this dataset to calculate flexibility accurately at any position on the surface of thin-walled parts. The hybrid model employs a CNN to address the nonlinear spatial correlation issue in flexibility prediction. Moreover, the hybrid model incorporates two KNN algorithms to alleviate the overfitting challenge stemming from the straightforward input features and extensive dataset size. In comparison to traditional deep learning models, the K-K-CNN hybrid model presents notable benefits in predicting flexibility for complex thin-walled parts at any given position, which affirms its robustness and accuracy. The proposed prediction method for flexibility can provide high-quality data-driven information for monitoring the elastic deformation of thin-walled parts. [ABSTRACT FROM AUTHOR]

The paper presents a static analysis of the buckling and post-buckling state of thin-walled cold-formed steel (TWCFS) lipped channel section beam-columns subjected to eccentric compression. Eccentricity is taken into consideration relative to both major and minor principal axes. An analytical–numerical solution to the buckling and post-buckling problems is described. The solution is based on the theory of thin plates. Equations of equilibrium of section walls are derived from the principle of stationary energy. Then, to solve the problem, the finite difference (FDM) and Newton–Raphson methods are applied. Linear (buckling) and nonlinear (post-buckling) analyses are performed. As a result, pre- and post-buckling equilibrium paths are determined. Comparisons of the obtained numerical results, FE simulation results, and experimental test results are carried out and presented in comparative load-shortening diagrams. Additionally, a comparison of the buckling forces and buckling modes obtained from theoretical analysis and experiments is presented. [ABSTRACT FROM AUTHOR]

The accurate perception of external environment information through the robot foot is crucial for the mobile robot to evaluate its ability to traverse terrain. Adequate foot-end contact signals can provide robust support for robot motion control and decision-making processes. The shape and uncertain rotation of the wheel-legged robot foot end represent a significant challenge to sensing the robot foot-end contact state, which current foot-end sensing schemes cannot solve. This paper presents a sensing method for the tire stress field of wheel-legged robots. A finite element analysis was conducted to study the deformation characteristics of the foot-end tire under force. Based on this analysis, a heuristic contact position estimator was designed that utilizes symmetrical deformation characteristics. Strain sensors, arranged in an array, extract the deformation information on the inner surface of the tire at a frequency of 200 Hz. The contact position estimator reduces the dimensionality of the data and fits the eigenvalues to the estimated contact position. Using support vector regression, the force estimator utilizes the estimated contact position and sensor signal to estimate the normal reaction force, designated as FZ. The sensing system is capable of detecting the contact position on the wheel circumference (with a root mean square error of 1.150°), as well as the normal force of 160 N on the Z axis (with a root mean square error of 6.04%). To validate the efficacy of the sensor detection method, a series of randomized and repeated experiments were conducted on a self-constructed test platform. This novel approach offers a promising avenue for perceiving contact states in wheel-legged robots. [ABSTRACT FROM AUTHOR]

As urban population increases, integrated underground–aboveground complexes are being constructed at growing paces in major cities. The seismic analysis of such complexes is crucial for the safety and functionality in the threat of potential earthquake disasters. However, fine‐grained numerical modeling and analysis of such large and complex structures are still inefficient due to the consideration of the soil–structure interaction (SSI). To address this challenge, an efficient approach for numerical modeling of large‐scale seismic SSI analysis is presented in this paper to overcome the limitations of existing finite element analysis (FEA) software. Moreover, a multi‐graphic processing unit (GPU) parallel explicit algorithm is implemented for the nonlinear dynamic SSI problems to further increase the computational efficiency. A large underground–aboveground complex project in China is used as an example to demonstrate the capability of the integrated method. The accuracy and reliability of the multi‐GPU parallel explicit finite element algorithm for SSI analysis (GFEA‐SSIA) are verified through a comparative analysis of the linear‐elastic and nonlinear dynamic response of the building calculated by GFEA‐SSIA and common FEA software. Finally, the structural response and structural damage of the underground–aboveground complex are analyzed under multidirectional seismic motions, and the damage distributions of the structures are provided. [ABSTRACT FROM AUTHOR]

In this paper, we consider a finite element approximation of an optimal control problem governed by a time fractional wave equation. We first prove the well-posedness and regularity of the optimal control problem. Furthermore, we discuss the fully discrete scheme and analyze its stability and a priori error estimate. Numerical experiments are given to illustrate the theoretical findings. • The weak solution to the time fractional wave equation is studied under low regular data. • The well-posedness and regularity of the optimal control problem are proven. • Numerical schemes exist the first-order and second-order accuracy in time and space respectively. [ABSTRACT FROM AUTHOR]

This work studied in detail for the first time the bending of laminated composite plates subjected t mechanical variations by new theory Trigonometric and Trigonometric-Hyperbolic functions of shear deformation. From the Euler-Lagrange hypothesis and the equations of the shear deformation theory, we will develop a present method. One of the most important problems of composite plates is the analysis of their bending behavior. The correct approach used to study their bending behavior includes two trigonometric and trigonometric-hyperbolic functions satisfying the null shear stress condition at the free edges. In this paper the bending problem is solved analytically by developing a computational code and numerically solved by Finite Element Method. In order to simplify the study of the bending behavior, an approach taking into consideration the effect of the transverse shear deformation without the shear coefficient of correction with only four unknowns has been developed while requiring five or more unknowns for other theories. Convergence analysis has been carried and the results are compared to open literature available for plate bending analysis. The approach proves to be simple and useful in analyzing the bending behavior of composite layered plates. [ABSTRACT FROM AUTHOR]

Tire acoustic cavity resonance (TACR) noise is one kind of low-frequency and narrowband noise that particularly annoys passengers, especially becomes more prominent in electric vehicles. This paper demonstrates a numerical investigation on the TACR noise via the acoustic-structural coupling finite element model. The accuracy of this coupling finite element model is validated both by the analytic method based on the superposition theory of traveling waves and experimental modal tests of tire cavities. According to the simulated models, the influences of external factors including the inflation pressure and road load on the TACR noise are studied. Meanwhile, as the main constituent component in the tire model, the effect of belt cord is also discussed. To design the low-TACR noise tire, various design parameters of the inner contour in tire cavity are analyzed. By the orthogonal test and range analysis, it is found that the contour parameters can slightly affect the modal frequency and sound pressure of TACR noise. Among these, the curvature radius of the tire shoulder ρ 4 , tire sidewall ρ 5 , and half-section height H have a more obvious influence on the TACR noise than the curvature radius of the tire crown region ρ 2 and transition region ρ 3 . Meanwhile, it also illustrates that increasing the curvature radius at the tire shoulder ρ 4 and reducing the half-section height H has a reduction effect on the TACR noise. This study can contribute to the design of low road noise tires. [ABSTRACT FROM AUTHOR]