Three-dimensional Simulation Model Research Articles

The Influence of Pre-Chamber Parameters on the Performance of a Two-Stroke Marine Dual-Fuel Low-Speed Engine

With increasing environmental pollution from ship exhaust emissions and increasingly stringent International Maritime Organization carbon regulations, there is a growing demand for cleaner and lower-carbon fuels and near-zero-emission marine engines worldwide. Liquefied natural gas is a low-carbon fuel, and when liquefied natural gas (LNG) is used on ships, dual-fuel methods are often used. The pre-chamber plays a key role in the working process of dual-fuel engines. In this paper, an effective three-dimensional simulation model based on the actual operating conditions and structural characteristics of a marine low-pressure dual-fuel engine is established. In addition, the effects of changing the Precombustion chamber (PCC) volume ratio and the PCC orifice diameter ratio on the mixture composition, engine combustion performance, and pollutant generation were thoroughly investigated. It was found that a small PPC volume ratio resulted in a higher flame jet velocity, a shorter stagnation period, and an acceleration of the combustion process in the main combustion chamber. When the PCC volume was large, the Nitrogen oxygen (NOx) ratio emission was elevated. Moreover, the angle of the PCC orifice affected the flame propagation direction of the pilot fuel. Optimizing the angle of the PCC orifice can improve combustion efficiency and reduce the generation of NOx. Furthermore, reasonable arrangement of the PCC structure can improve the stability of ignition performance and accelerate the flame jet velocity.

A Study on Three-Dimensional Multi-Cluster Fracturing Simulation under the Influence of Natural Fractures

Multi-cluster fracturing has emerged as an effective technique for enhancing the productivity of deep shale reservoirs. The presence of natural bedding planes in these reservoirs plays a significant role in shaping the evolution and development of multi-cluster hydraulic fractures. Therefore, conducting detailed research on the propagation mechanisms of multi-cluster hydraulic fractures in deep shale formations is crucial for optimizing reservoir transformation efficiency and achieving effective development outcomes. This study employs the finite discrete element method (FDEM) to construct a comprehensive three-dimensional simulation model of multi-cluster fracturing, considering the number of natural fractures present and the geo-mechanical characteristics of a target block. The propagation of hydraulic fractures is investigated in response to the number of natural fractures and the design of the multi-cluster fracturing operations. The simulation results show that, consistent with previous research on fracturing in shale oil and gas reservoirs, an increase in the number of fracturing clusters and natural fractures leads to a larger total area covered by artificial fractures and the development of more intricate fracture patterns. Furthermore, the present study highlights that an escalation in the number of fracturing clusters results in a notable reduction in the balanced expansion of the double wings of the main fracture within the reservoir. Instead, the effects of natural fractures, geo-stress, and other factors contribute to enhanced phenomena such as single-wing expansion, bifurcation, and the bending of different main fractures, facilitating the creation of complex artificial fracture networks. It is important to note that the presence of natural fractures can also significantly alter the failure mode of artificial fractures, potentially resulting in the formation of small opening shear fractures that necessitate careful evaluation of the overall renovation impact. Moreover, this study demonstrates that even in comparison to single-cluster fracturing, the presence of 40 natural main fractures in the region can lead to the development of multiple branching main fractures. This finding underscores the importance of considering natural fractures in deep reservoir fracturing operations. In conclusion, the findings of this study offer valuable insights for optimizing deep reservoir fracturing processes in scenarios where natural fractures play a vital role in shaping fracture development.

Ammonia energy fraction effect on the combustion and reduced NOX emission of ammonia/diesel dual fuel

With stringent regulations of internal combustion engine on reducing CO2 emission, ammonia has been used as an alternative fuel. Investigating how engine-related performance is affected by partial ammonia replacement of diesel fuel is essential for understanding the combustion. Therefore, in this study, a three-dimensional numerical simulation model is developed for the burning of two fuels of diesel and ammonia based on relevant parameters (i.e., compression ratio, load, ammonia energy fraction, etc.) in a lab-made diesel engine. The consequences of load and compression proportion on combustion and pollutant emissions are investigated for ammonia energy fractions between 50% and 90%. When the ammonia portion rises, the increased ammonia equivalent ratio causes ammonia to move away from the dilute combustion boundary and accelerates the combustion rate of ammonia. An increase in compression ratio significantly increases the specified thermal performance and combustion efficacy. When the compression ratio is 16, as the ammonia energy fractions increases, due to the increase in the proportion of ammonia, that is, the proportion of nitrogen atoms increases, more NOx is generated during the combustion process. When the ammonia substitution rate is 90%, as the compression ratio increases, the cylinder pressure and temperature increase. The combustion efficiency of ammonia increases, generating more NOx and NOx emissions can reach 0.66 mg/m3. At a compression ratio of 18, the NOx emissions can reach 1.59 mg/m3. However, under medium and low load conditions, as the ammonia fraction increases, the total energy of fuel decreases, and the combustion efficiency of ammonia decreases, resulting in a decrease in the heat released during combustion and a decrease in NOx emissions. When the ammonia substitution rate is 90% and the load is 25%, NOx emissions reach 0.1 mg/m3. This research provides theoretical suggestions for the profitable and use ammonia fuel in internal combustion engines in a clean manner.

Weakly guiding approximation of a three-dimensional waveguide model for extreme ultraviolet lithography simulation

A three-dimensional (3D) waveguide model is applied in extreme ultraviolet (EUV) lithography simulations. The 3D waveguide model is equivalent to rigorous coupled-wave analysis, but fewer field components are used to solve Maxwell’s equations. The 3D waveguide model uses two components of vector potential, A x and A y , corresponding to the two polarizations. The electric field of the A x polarization is approximately parallel to the x axis, and the electric field of the A y polarization is approximately parallel to the y axis. The 3D waveguide model solves a coupled vector wave equation for two polarizations. The refractive index of conventional EUV absorbers is close to that of vacuum. The weakly guiding approximation in optical fiber theory is applied to the 3D waveguide model. The coupled vector wave equations for the two polarizations are decoupled into two independent scalar wave equations. Maxwell’s equations are simplified to a set of scalar wave equations. The weakly guiding approximation reduces the computation time to solve the equations. The computation time required to solve the weakly guiding approximation is about 1/5 of the time to solve the original 3D waveguide model.

Analysis of damage characteristics of the rectangular steel container under near-earth explosion loading

In this study, the dynamic response of a semi-buried steel vessel under a near-surface explosion shock load was investigated through full-scale model tests. A three-dimensional precise numerical simulation model of the semi-buried steel vessel structure was established, and the dynamic response and damage consequences under different load conditions were analyzed. The influence of the thickness of corrugated steel on the dynamic damage and antiknock properties of the structure was studied. The results show that the weak part of the rectangular steel vessel structure is mainly located at the connection between the blast-facing panel and the beam column under near-surface explosion conditions. The peak stress and displacement of the blasting surface decrease with the increase in the horizontal angle of the incident wave, whilst/while increase with the height of the explosion. Under the same explosive yield, the main failure mode of the structure is the tearing of the plate and beam under plastic deformation. Increasing the thickness of the corrugated steel can significantly restrain the plastic deformation and acceleration of the plate. When the thickness is increased from 2mm to 6mm, the antiknock performance is significantly improved. However, when the thickness is increased from 6mm to 8mm, the improvement in antiknock performance slows down. Increasing the height of the enclosing soil can also effectively improve the explosion-proof performance of the structure, with full burial providing the best effect. The research results provide a theoretical basis for the application of rectangular steel vessel structures in the field of protection engineering.

Interaction characteristics between the large ship-rubber fenders-wharf structure in deep-water considering the flexibility of the wharf structure

High pile-supported wharf structures in offshore deep water have significant flexibility due to long free ends and large slenderness ratios of the pile bodies, resulting in substantial deformation of the pile body under kinds of loads. In the design of this new type of wharves, it is expected to consider the characteristic of elastic deformation of pile body absorbing a portion of the ship impact energy reasonably. However, the existing calculation method still adopts the theoretical formulas of traditional wharf structures, failure to incorporate this advanced design concept fully. In order to solve this, firstly, a comprehensive three-dimensional numerical simulation model encompassing large ship, rubber fenders and wharf structure, was developed using Ansys/LS-Dyna, and the calculation model was validated against the existing research results. Subsequently, dynamic response characteristics of wharf structures during the normal berthing of ships were examined meticulously, and the variations in absorption energy and reaction forces of rubber fenders and wharf structures were analyzed systematically. And then the analysis results were compared with previous research, the limitations of the existing design code for wharf structures were discussed as well. Finally, the findings demonstrate that peak values of absorption energy and energy transfer times increase with the increase of ship loading capacities. And the elastic deformation of the wharf structure indeed absorbs considerable partial of the impact energy effectively, with about 30 % for conditions in this study, leading to lower impact forces acted on the wharf structure. Moreover, the maximum pile top displacement of the wharf structure is also smaller than that reported in the existing literature in the same load condition, indicating that considering the structural flexibility in the design of this new type of wharf structure is more in line with reality, and it is more reasonable and economical as well.

Optimization Design Method of 2D + 3D Slope Shape for Landslide Prevention in Open-Pit Coal Mine

In order to improve the stability of the slope and prevent the occurrence of landslide disasters, taking the east slope of the first mining area of Zhundong Open-pit Coal Mine as the engineering background, a combination of the two-dimensional limit equilibrium method and three-dimensional numerical simulation is used to optimize the shape of the east slope. By selecting a typical calculation profile based on the Bishop method and the residual thrust method in the two- dimensional rigid body limit equilibrium method, the stability analysis of the profile slope is carried out, and the slope shape of the profile position meeting the requirements of the safety reserve coefficient and stripping ratio is preliminarily designed. Based on the three-dimensional finite element strength reduction method, aiming at the slope shape that does not meet the requirements of safety reserve coefficient and stripping ratio, the width of the transport plate is reasonably changed, and a three-dimensional numerical simulation model of various schemes is established. The influence of different tracking distances and slope angles on slope stability is studied. The three-dimensional stability of the slope under different spatial forms is calculated, and the optimal tracking distance and optimal slope angle are determined. The slope stability and stripping ratio are further optimized. The results show that when the width of the transport plate of the DBS3 section slope is 8meters, it does not meet the requirements of safety reserve coefficient 1.2. When the width of the transport plate is set to 24 meters, it meets the requirements of a safety reserve coefficient of 1.2 and an economic stripping ratio of not more than 8.0 m3/t. The three- dimensional numerical simulation results show that the optimal tracking distance on the east side is 50 meters, and the optimal slope angle is 35°. After the optimization design of a two-dimensional and three-dimensional slope shape, 2.456 million tons of coal resources were mined, creating a profit of about 21.1268 million yuan. While preventing landslide disasters, the recovery rate of coal resources was further improved.

Optimization analysis for bolted steel I-beam: Considering various sizes of beam web with a focus on sustainability design

Abstract Steel structural construction has consistently held a pivotal role in modern construction practices. Presently, a substantial portion of research efforts has been directed towards the examination of bolted connections, while the optimization of steel beam usage has not garnered equivalent attention. This research initiative is focused on determining the ideal dimensions for the beam web through the utilization of finite element analysis (FEA). This analysis specifically investigates bolted steel I-beams featuring varying beam web sizes under concentrated loads. To achieve this, a three-dimensional (3D) simulation model of a bolted beam-to-column structure has been meticulously constructed and executed within the ABAQUS software to derive optimization insights. The beam web sizes considered span a range from 200 mm, representing a smaller size, to 400 mm, signifying a medium size, and up to 600 mm, which represents a larger beam web size. Following the simulation, the results pertaining to beam displacement and bending stress have been tabulated and compared against manually calculated values to validate the accuracy of the simulation. Remarkably, the percentage disparity between the FEA and manual calculations remained below 5 %, signifying a high level of acceptability. From the post hoc test, the efficiency rate of the 400 mm beam web size stood at 76.67 %, rendering it the most favourable choice in terms of cost-effectiveness and sustainability among the three considered sizes.

Three-dimensional numerical simulation and theoretical model of a hollow droplet impacting on a solid surface

The phenomenon of a hollow droplet impacting on a solid surface is widely found in various fields. The dynamic characteristics of hollow diesel droplets impacting on a solid surface are studied by combining numerical simulation and theoretical analysis. The dynamic contact angle model presented in this paper couples the advantages of existing dynamic contact angle models for simulating both the spreading and retracting stages. It also considers the continuous variation of the contact angle during the maximum spreading state. Compared to existing models, the maximum error has been reduced from 14.9% to 4.6%. The effects of impact velocity, impact angle, and volume ratio of a hollow droplet on the spreading and jetting characteristics are investigated by three-dimensional numerical simulations. It is found that air entrainment occurred in the counter-jet, and the presence of the impact angle increased the asymmetry of the counter-jet and spreading liquid film, promoting fingerlike splashing at the front liquid film. Based on energy conservation law, the theoretical prediction models of the maximum spreading coefficient of the hollow droplet impacting on the surface and the velocity of the counter-jet at the maximum spreading state are established using the multi-regional modeling method and the energy distribution principle. Compared with existing hollow droplet theoretical models, the proposed theoretical models exhibit a more concise expression, higher accuracy, and wider applicability range.

Static and Dynamic Analyses of Spindle Collet Made of Different Materials Using Finite Element Modeling

The spindle collet, a critical component in various machine tools, plays a pivotal role in determining the success of machining operations. This paper aims to study the static and dynamic parameters of collet structures made from three different materials using the finite element method. A three-dimensional model and computer simulation were conducted in Autodesk Inventor software. Simulations are performed using identical boundary conditions and mesh size. Static analysis is performed with varied applied forces where total deformation and Von Mises stress are measured. For the dynamic analysis, the natural frequencies and mode shapes are measured up to the first five modes. The variations in stress is minimal when the material is altered. The magnitude of deformation varies significantly with changes in material. The relative deformation values demonstrate that carbon steel deforms more than alloy steel by almost 3%, while stainless steel deforms more than alloy steel by 6%. Materials with higher Young's modulus and lower density have been found to increase the natural frequencies, reducing total deformation and von Misses stress. The use of alloy steel in the industry offers an advantage over the other two materials. The results provide improved insight into the appropriate materials for the collet.

Active Control of Longitudinal Combustion Instability in Bluff-Body Stabilized Premixed Flames with Multiple Neural Network Controller

ABSTRACT The suppression performance of multiple neural network controller on the longitudinal combustion instability is investigated in the Volvo bluff-body stabilized premixed flame. Two controlled plant models, e.g. one-dimensional (1D) reduced order acoustic network model and three-dimensional (3D) large eddy simulation model are adopted to verify this control strategy. In order to build the reduced order Volvo acoustic network model, large eddy simulation is adopted to extract the characteristics of system heat release under wide bandwidth and construct the flame transfer functions. For the 1D acoustic network model, this controller can completely eliminate the pressure oscillation using fuel valve as actuator. For the 3D large eddy simulation model that mimics the real-time oscillating characteristics, clustered prediction method and amplitude limiting method are adopted in multiple neural network control strategy to better account for the nonlinear characteristics of the 3D Volvo flame. This control strategy enables the modulation of the phase relationship between the heat release rate oscillation and pressure oscillation by introducing small perturbations to the inlet equivalence ratio, and attenuates the amplitude of the second harmonic of the unstable mode effectively.

Performance enhancements of power density and exergy efficiency for high-temperature proton exchange membrane fuel cell based on RSM-NSGA III

In this study, a three-dimensional simulation model was developed and a multi-objective optimization of a single-channel high-temperature proton exchange membrane fuel cell (PEMFC) was carried out by the response surface methodology (RSM) and the non-dominated genetic algorithm III (NSGA III). Firstly, five main influencing factors including operating temperature (T), operating pressure (p), gas diffusion layer porosity (εGDL), proton exchange membrane (PEM) thickness (Mmem) and anode stoichiometry ratio (λa) were identified. The optimization objectives include three cell performance evaluation indexes: power density (wpow), system efficiency (ηsystem), and exergy efficiency (ηexergy). Secondly, the prediction accuracies of the RSM and artificial neural network (ANN) for the three performance evaluation indexes of high-temperature PEMFC were compared. Finally, the optimal solution from the Pareto solutions derived through RSM-NSGA III. The results show that the maximum wpow, maximum ηsystem and maximum ηexergy are 0.793 W cm−2, 22.31 % and 54.69 %, respectively. The optimal result obtained by combining EWM and VIKOR is wpow = 0.6726 W cm−2, ηsystem = 19.44 % and ηexergy = 54.01 %, respectively. The corresponding condition is T = 447.18 K, p = 1.15 atm, εGDL = 0.617, Mmem = 0.0412 mm, and λa = 1.26, respectively. Moreover, the suggestions are provided for the design of more efficient high-temperature PEMFC.

Encapsulation structure and in-situ embedding test and model study of a new type of road piezoelectric energy collector

In recent years, research on road piezoelectric energy harvesters (RPEH) has made great progress in the field of road vibration energy harvesting. However, the research on the packaging structure and embedding method of RPEH is not comprehensive. In this paper, a new type of transducer unit structure is developed, which can improve the electrical energy output compared with the drum piezoelectric transducer, and its waterproof performance of the transducer unit is tested under different water depth and duration. Then, a new RPEH package structure is made, and the output performance is compared with the three commonly used package structures in the current research, and the field construction of RPEH with different embedding methods is carried out on the actual road to evaluate the electrical energy output capacity of RPEH under different embedding methods, and the best construction method of RPEH is summarized through durability test. Finally, the three-dimensional simulation model of RPEH with different construction methods is established, and the stress and voltage are analyzed by finite element method to verify the best embedding method. The test results show that the piezoelectric output efficiency of the new transducer unit is increased by 38.9 %, and it still has good electrical energy output after 24 h immersion at a depth of 1.5 m, and its waterproof level reaches IPX8. Under larger excitation frequency and load, the indoor output of RPEH in the new package structure has better output performance in indoor test. The bolt-spring package is used to apply pre-pressure to RPEH, and the output efficiency of RPEH is improved. When the excitation conditions are 10 kN, 10 Hz, the output power is increased by 11.1 %, 127.3 % and 247.4 % compared with the other three package structures. The maximum power of RPEH is 3.41 W by placing concrete on the bottom surface of RPEH sponge washer and using asphalt as the pavement layer to ensure the durability of piezoelectric pavement. The finite element analysis illustrates the strain and voltage variation of piezoelectric pavement with different embedment modes, and the simulation results are consistent with the test results. The research results provide an important basis for the design of RPEH packaging structure and embedding method in road piezoelectric engineering.

Numerical simulation and optimisation design for ventilation and heat dissipation in high-temperature and high-load indoor substations

Under high-temperature and high-load operational conditions, inadequate ventilation and suboptimal cooling arrangements within indoor substations result in high oil temperatures, posing a threat to the secure and steady operation of transformers. In this paper, the ventilation and heat dissipation effect of a 110 kV indoor substation is studied by the computational fluid dynamics method. Initially, the three-dimensional simulation model of the main transformer chamber is constructed to mirror the actual substation structure. Subsequently, the impact of the upper outlet on ventilation and heat dissipation is explored. The results show that the proximity of the upper outlet to the fan on the side wall prompts the fan to draw air from the upper outlet, diminishing the air volume entering the lower intake outlet by 38.8 %. It proves detrimental to the transformer's heat dissipation efficiency. Lastly, with the upper outlet closed, the study delves into the impact of various configurations of intake and exhaust ports on optimizing ventilation and heat dissipation. A total of six cases of two scenarios for the air inlet location and three exhaust vent positions are explored. According to the simulation results, the optimal case of substation ventilation and heat dissipation is obtained by considering four evaluation parameters.

Temperature distribution analysis and mechanical property prediction of wire loops in the Stelmor air-cooling system

During the Stelmor air-cooling process, the temperature distribution has a significant impact on the wire loops' final mechanical properties. The temperature distribution during air-cooling process is accurately solved by establishing three-dimensional model and Computational Fluid Dynamics (CFD) simulation. It is also found that the heat transfer coefficient on the cross section of the wire loop is very different, with a difference of 70–100 Wm−2K−1. The average value of the surface heat transfer coefficient of the wire loops section is used to replace the surface heat transfer coefficient of the wire loops, which improves the temperature distribution calculation accuracy, which is not considered by previous scholars. The finite difference method is used to calculate the mathematical model of the temperature distribution during the Stelmor air-cooling process. It is found that the maximum temperature difference of the wire during the cooling process is 90–110 K. The prediction model of mechanical property of wire after air-cooling was established, and distribution of mechanical property of the wire loops was obtained, and the reason for the uneven distribution of mechanical property of the wire loops was found.

A Novel Numerical Model Considering Unsaturated Soil Properties and Computational Study on Heat and Moisture Transfer Characteristics of Helix Ground Heat Exchanger

Abstract A novel three-dimensional numerical simulation model for the helix ground heat exchanger was proposed in this paper, which takes into account unsaturated soil properties. This model is more suitable for real working conditions. To validate its accuracy, a miniature model heat transfer experimental platform was constructed. Additionally, the study conducted simulation research using three types of soil with significantly different thermal and moisture characteristics. Moreover, the comprehensive thermal conductivity and water diffusion coefficient of these three soil types were determined by relevant literature and experimental tests. The aim was to comprehensively explore the impact of different soil types on the heat and mass transfer of the helix ground heat exchanger. The results indicate that the numerical model developed in this paper accurately captures the heat and mass transfer characteristics of the helix ground heat exchanger to a certain extent. Increasing the comprehensive thermal conductivity and water diffusion coefficient of the soil can significantly enhance the heat exchange capacity of the exchanger. For instance, under sandy loam conditions, the heat exchange capacity is approximately 20.73% higher compared to clay loam conditions. The study also identifies two distinct areas around the helix ground heat exchanger: the severe change region and the soft change region. In the severe change region, there is a notable decrease in soil water content near the exchanger, which inevitably weakens the thermal conductivity of the soil. It is advised to minimize this effect through measures like active water spraying.

A three-dimensional numerical model for computing the actual voltage induced in a Barkhausen surface coil

This investigation presents a novel three-dimensional simulation model designed for the analysis of Magnetic Barkhausen Noise (MBN). The model’s formulation is directly deduced from Maxwell’s equations, establishing a Boundary Value Problem (BVP) with boundary conditions tailored to a specific configuration involving a lengthy steel rod encased within a solenoid, subjected to sinusoidal current excitation, all contained within an air enclosure. This chosen configuration holds physical equivalence to the majority of setups commonly utilized in Barkhausen experiments. Utilizing the proposed three-dimensional model, it became feasible, for the first time, to computationally assess the voltage induced in the coil by the flux emanating from all spatial directions. The findings indicate that the MBN signal is primarily generated by the flux rate aligned with the direction of the applied field. Nevertheless, simulations confirm additional contributions to the signal originating from magnetic flux in directions significantly deviating from the applied field direction. It is demonstrated that these contributions from lateral flux significantly influence the angular dependence of the Barkhausen signal energy, showcasing a notable alignment between the proposed model and experimental results. The results underscore the capability of the proposed model to facilitate the simulation of Barkhausen Noise.

Numerical study on dynamic flow of flexible extension nozzle in rocket motors

The use of extension nozzles in rocket propulsion systems is one of the effective methods to increase the geometric expansion ratio and thrust of the nozzle. As a new type of extension nozzle, the flexible extension nozzle can achieve continuous changes in nozzle expansion ratio. To explore the characteristics of the flow field structure at the “bulge” structure and the step structure in the actual configuration of the flexible extension nozzle, and the impact on the performance of the nozzle, a three-dimensional numerical simulation model of the flexible extension nozzle was first established to study the flow field at the “bulge” structure, and then a two-dimensional unsteady numerical simulation model was established by using the coupling technology of dynamic grid and overlapping grid to study the flow field at the step structure. The results show that during the unfolding of the nozzle surface, a complex flow structure of “expansion fan-mixing layer-shock-vortex” were formed in the step area at the junction of the fixed section and the moving extension section. The specific impulse loss of the 3D model is larger than that of the 2D model due to the consideration of the “bulge” structure. The flow loss of the “bulge” structure and the step structure to the nozzle is 0.73% and 0.65% respectively.

Design, thermal-mechanical coupling analysis, and optimization of polymeric matrix composite sandwiches with a lattice core exposed to a high temperature

This study presents the design of a sandwich structure tailored for post-heat transfer applications subjected to out-of-plane compression. A three-dimensional finite element simulation model was developed to analyze the temperature distribution within the sandwich structure and investigate the effects of high-temperature exposure on its mechanical behaviors. The structure was subjected to a temperature of 300 °C for 400 s, and the temperature distribution at the upper connection point between the top face sheet and the struts of the core was determined. Subsequently, upon returning the sandwich to ambient temperature, a comprehensive calculation of its mechanical properties was conducted and then enhanced by applying different optimization techniques. The results demonstrate that filling the core with Saffil alumina fibers in the presence of a 1.5 mm a thermal barrier coating of Superwool 607 helps increase the mechanical properties of the sandwich structure by around 55%.

Thermal and Mechanical Stress Analysis in Aircraft Hybrid Panels with Multi-Bolt Joints.

This study investigates the thermal stress and bolt load distribution in a hybrid panel structure of an aircraft mechanical joint under extreme temperatures. The hybrid panel structure comprises two aluminum alloy splices, six T-shaped composite stringers, and two composite skins, secured together with 96 bolts. This study analyzed the strain induced by thermal stress on composite materials and metals within the structure across temperatures, employing temperature environment tests ranging from room temperature to -54 °C, alongside a carrying capacity test at -54 °C. Furthermore, a three-dimensional simulation model of the panel structure was developed, incorporating considerations for contact, metal elastoplasticity, and the progressive damage failure of composite materials. This model facilitated the determination of thermal stress and bolt load distribution patterns. The results indicate a strong consistency between the finite element analysis outcomes and the experimental data. Temperature variations exacerbate the uneven distribution of bolt loads, concentrating the load near the edges of the hybrid structure while diminishing it in the center. The bolt load distribution parallel to the mechanical load direction forms an "M" shape, with a maximum load magnitude of approximately 31 kN. Perpendicular to the mechanical load, the bolt load undergoes significant changes, especially at the edges, reaching a maximum of about 20 kN, which warrants attention. The bolt-load distribution of the structure with the increase in mechanical load at -54 °C tends to be consistent with that at room temperature.