Abstract In this article, a series of nearly analytic symplectic partitioned Runge–Kutta (NSPRK) methods for the 3D seismic‐wave equation are developed. First, the Hamiltonian formulations for acoustic and elastic‐wave equations are presented, and then the spatial derivatives are discretized by a nearly analytic discrete operator to obtain a semidiscrete Hamiltonian system. The second‐, third‐, and fourth‐order symplectic partitioned Runge–Kutta schemes are then applied as the time integrator. A semianalytic procedure to facilitate the analysis of stability conditions and numerical dispersion relations is presented, and subsequent theoretical analysis shows that the NSPRK schemes preserve the wave velocity better than conventional symplectic schemes, especially on coarser grids. The numerical solutions computed by NSPRK schemes are compared with analytic solutions for the 3D acoustic and elastic cases. We implemented the NSPRK and some conventional schemes for a 3D acoustic wave propagation simulation in a parallel computer and compared their computational efficiencies. To generate comparable results, the NSPRK schemes require much less computer memory, central processing unit time, and communication time, which substantially accelerates the computation speed. The final simulation in the two‐layer acoustic model shows that the NSPRK schemes can suppress numerical dispersion and preserve the waveforms better than conventional symplectic schemes.
Background: Prostheses for the reconstruction of periacetabular bone tumors are prone to instigate stress shielding. The purpose of this study is to design 3D-printed prostheses with topology optimization (TO) for the reconstruction of periacetabular bone tumors and to add porous structures to reduce stress shielding and facilitate integration between prostheses and host bone. Methods: Utilizing patient CT data, we constructed a finite element analysis (FEA) model. Subsequent phases encompassed carrying out TO on the designated area, utilizing the solid isotropic material penalization model (SIMP), and this optimized removal area was replaced with a porous structure. Further analyses included preoperative FEA simulations to comparatively evaluate parameters, including maximum stress, stress distribution, strain energy density (SED), and the relative micromotion of prostheses before and after TO. Furthermore, FEA based on patients’ postoperative CT data was conducted again to assess the potential risk of stress shielding subsequent to implantation. Ultimately, preliminary follow-up findings from two patients were documented. Results: In both prostheses, the SED before and after TO increased by 143.61% (from 0.10322 to 0.25145 mJ/mm 3 ) and 35.050% (from 0.30964 to 0.41817 mJ/mm 3 ) respectively, showing significant differences ( p < 0.001). The peak stress in the Type II prosthesis decreased by 10.494% (from 77.227 to 69.123 MPa), while there was no significant change in peak stress for the Type I prosthesis. There were no significant changes in stress distribution or the proportion of regions with micromotion less than 28 μm before and after TO for either prosthesis. Postoperative FEA verified results showed that the stress in the pelvis and prostheses remained at relatively low levels. The results of follow-up showed that the patients had successful osseointegration and their MSTS scores at the 12th month after surgery were both 100%. Conclusion: These two types of 3D-printed porous prostheses using TO for periacetabular bone tumor reconstruction offer advantages over traditional prostheses by reducing stress shielding and promoting osseointegration, while maintaining the original stiffness of the prosthesis. Furthermore, in vivo experiments show that these prostheses meet the requirements for daily activities of patients. This study provides a valuable reference for the design of future periacetabular bone tumor reconstruction prostheses.
The Turkey–Syria earthquake on 6 February 2023 resulted in losses such as casualties, road damage, and building collapses. We mapped and quantified the areas impacted by the earthquake at different distances and directions using NOAA-20 VIIRS nighttime light (NTL) data. We then explored the relationship between the average changes in the NTL intensity, population density, and building density using the bivariate local indicators of the spatial association (LISA) method. In Turkey, Hatay, Gaziantep, and Sanliurfa experienced the largest NTL losses. Ar Raqqah was the most affected city in Syria, with the highest NTL loss rate. A correlation analysis showed that the number of injured populations in the provinces in Turkey and the number of pixels with a decreased NTL intensity exhibited a linear correlation, with an R-squared value of 0.7395. Based on the changing value of the NTL, the areas with large NTL losses were located 50 km from the earthquake epicentre in the east-by-south and north-by-west directions and 130 km from the earthquake epicentre in the southwest direction. The large NTL increase areas were distributed 130 km from the earthquake epicentre in the north-by-west and north-by-east directions and 180 km from the earthquake epicentre in the northeast direction, indicating a high resilience and effective earthquake rescue. The areas with large NTL losses had large populations and building densities, particularly in the areas approximately 130 km from the earthquake epicentre in the south-by-west direction and within 40 km of the earthquake epicentre in the north-by-west direction, which can be seen from the low–high (L-H) pattern of the LISA results. Our findings provide insights for evaluating natural disasters and can help decision makers to plan post-disaster reconstruction and determine risk levels on a national or regional scale.
Abstract We report on interfacial characteristics and chemistry of bonded Mg-Fe interfaces welded using friction stir assisted scribe technique (FaST). Two pairs of dissimilar joints: (AZ31-DP590) and (Pure Mg-DP590) were studied to shed light on joining mechanisms responsible for bonding of “immiscible” pairs of Mg and Fe. We present first direct experimental evidence of presence of oxide layer, Al segregation by atom probe tomography and nano steel grains close to interface by transmission electron microscopy study.
A large amount of traffic crash investigations have shown that rear-end collisions are the main type collisions on the freeway. The purpose of this study is to investigate the rear-end collision risk on the freeway. Firstly, a new framework was proposed to develop the rear-end collision probability (RCP) model between two vehicles based on Generalized Pareto Distribution (GPD). Secondly, the freeway rear-end collision risk (F-RCR) was defined as the sum of the rear-end collision probability of each vehicle and divided into three levels which was high, median, and low rear-end collision risk. Then, different machine learning algorithms were used to model F-RCR under the condition of an unbalanced dataset. The result of the RCP model showed continuous change and can identify the dangerous condition quickly compared to the traditional models even when the speed of the leading vehicle is faster than the following vehicle. When the vehicle distribution was unbalanced on road and the speed difference between adjacent lanes and the traffic volume was large, F-RCR will increase. Multi-Layer Perceptron (MLP) was found to be more suitable for modeling F-RCR. The framework provided in this research was transferrable and can be used in the freeway proactive traffic safety management system.
Precipitation hardening has been recently validated as a new mechanism for domain wall pinning and mechanical loss reduction in piezoelectrics. While anisometric precipitates have high pinning strengths, there is limited knowledge about the electrical anisotropy of the precipitation-hardened piezoceramics. In the present work, we successfully orient the precipitates in Li0.18Na0.82NbO3 piezoceramics by applying a uniaxial stress during the aging and studied its electrical anisotropy. Predicted by mechanical simulation and verified by transmission electron microscopy, it is demonstrated that the precipitate variant with its long axis perpendicular to the applied stress is energetically favored. The electrical anisotropy of the stress-assisted aged Li0.18Na0.82NbO3 is studied by applying electrical fields parallel or perpendicular to the stress axis. The domain wall contribution to permittivity is found to vary by more than a factor of two depending on orientation. In addition, the domain walls are more difficult to be activated by increasing the temperature when the electric field is perpendicular to the stress axis. Our work highlights the precipitate variant selection induced by stress-assisted aging and the related electrical anisotropy in piezoceramics. This technique enables the precipitate orientation in piezoceramics and the utilization of its anisotropy, providing fundamental insight into precipitate-domain-wall interactions and setting the ground for leveraging precipitation hardening effect in piezoceramics. The authors achieve precipitate alignment in piezoceramic by applying a uniaxial stress during aging process, resulting in electrical anisotropy. The optimized material has succinct advantages for high-power piezoelectric applications.
Significance Downscaling material structure into ever-decreasing levels for modern technological advancements demands small-scale mechanical characterizations of materials, such as nanoindentation. Thus, a comprehensive understanding of indentation size effect (ISE) has a far-reaching impact, but such understanding remains incomplete at the submicron scale despite a decades-long quest since the 1980s. Here, we explored the origins of this effect experimentally by linking dislocation behavior evolution with nanoindentation mechanical data at various depths. Importantly, the dislocation behavior is found to be strongly depth dependent and gives rise to subgrain formation progressively. Aided by subgrain boundaries, the ISE mechanism transitions from dislocation source starvation to dislocation interaction as indentation depth increases. The critical mechanism transition elucidates the distinctive ISE behavior at the submicron scale.