Element Size Research Articles

Permeability of porous membrane polymers modified by supercritical carbon dioxide

An approach to predict the gas permeability of membrane polymers after supercritical CO2 treatment is proposed. The approach is based on the connection of the temperatures of secondary relaxation transitions with the effective sizes of mobile free volume elements in the polymers. The correlation between permeability of nitrogen and effective sizes of mobile holes for a set of polymers is established. The effect of supercritical CO2 on the nitrogen permeability of polycarbonate, polysulfone, polyvinylbutyral is analyzed by FTIR spectroscopy of low-molecular weight conformationally-inhomogeneous compounds introduced in the polymers. Membranes were exposed at 40 MPa and 333 K for 4 h through static treatment and dynamic treatment separately. For polyvinylbutyral, the nitrogen permeability did not change after supercritical CO2 modification while for polysulphone, the effective volume of mobile holes increased, but the nitrogen permeability decreased. For polycarbonate after supercritical CO2, the effective volume of mobile holes and the nitrogen permeability increased.

Investigating the flexural behavior of ASR‐damaged RC beam through internal stress and cracking development using 3D Rigid Body Spring Model

AbstractThe alkali‐silica reaction (ASR) is one of the durability issues that affect the safety of concrete structures. Numerical simulation is useful for monitoring the internal condition of a reinforced concrete (RC) structure with ASR damage and predicting its long‐term behavior. Since discrete methods of simulation suit situations where concrete undergoes durability issues associated with expansion and cracking, a mesoscale discrete analysis method known as the three‐dimensional Rigid Body Spring Model (3D RBSM), as already used to simulate concrete ASR damage at the material scale, has been further developed to simulate ASR damage at the structural scale. In this development, the aggregate is not separately modeled and expansive elements are included to simulate ASR expansion, allowing a larger element size (1–2 cm). The number of elements and computation time are reduced to less than 1%. This proposed model is used to simulate ASR damage in a RC beam and the beam's residual flexural capacity. The simulation allows the internal stress and cracking condition to be visualized, leading to an explanation as to why the loading capacity of a RC beam is not affected by ASR damage: first, almost no cracks form at the core of the beam section due to the stirrup confinement, so the effective depth of the cross section is maintained after ASR damage; second, the tensile reinforcements have already yielded although ASR damage exists. Besides, the inhibited development of shear cracks may be one of the reasons for the slight improvement of flexural capacity of the damaged‐RC beam in experiment and simulation.

Two-step printing bionic self-lubricating composite surface fabricated by selective laser melting ink-printed metal nanoparticles

The combination of self-lubricating materials and surface texture was a promising strategy to improve tribology performance. The designed parameters and manufacturing methods of the textured composite surface were pivotal in fully harnessing this synergistic effect. Inspired by natural structures such as the scapharca subcrenata shell, trimeresurus stejnegeri snakeskin, crocodile, and tree frog toe pad surface, bionic self-lubricating composite surfaces (BSCSs) with different textures were prepared on AISI 4140 steel substrate to enhance tribological properties. These BSCSs combined self-lubricating composite material (SLCM) with copper (Cu) material in bionic patterns using the selective laser melting ink-printing metal nanoparticles manufacturing method. Tribological properties and worn surfaces of the BSCSs with various design parameters were analyzed under dry sliding conditions. Results indicated that the synergy effects of printed Cu and SLCM texture effectively reduced friction coefficient and improved wear resistance of the BSCSs. The SLCM texture density and the size of array element in the BSCSs were found to influence the types of lubricating film formed during the sliding process, thereby affecting the tribological performance. This investigation provided a potential surface design strategy and a more flexible manufacturing method for surface engineering.

Some results on conjugacy class sizes of primary elements of a finite group

Some results on conjugacy class sizes of primary elements of a finite group

Determining the Size of Representative Volume Elements for a Two-Dimensional Random Aggregate Numerical Model of Asphalt Mortar without Damage.

The size of the representative volume element (RVE) for the two-dimensional (2D) random aggregate numerical model of asphalt mortar in a non-destructive state, which directly affects the time required to simulate the linear viscoelastic behavior from asphalt mastic to asphalt mortar. However, in the existing literature, limited research has been conducted on the size determination of the numerical model RVE for asphalt mortar. To provide a recommended size for the typical 2D random aggregate numerical model RVE of asphalt mortar in a nondestructive state, this paper first applies the virtual specimen manufacturing method of asphalt concrete 2D random aggregate to asphalt mortar. Then, it generates numerical model RVEs of asphalt mortar with different maximum particle sizes, after which geometric and numerical analyses are conducted on these models. Finally, based on the geometric and numerical analysis results, the recommended minimum sizes of RVE for the 2D asphalt mortar numerical model are provided.

A micro-meso coupled model for coral reef rocks based on CT Scanning

The coral reef rock has a rich porous structure, and establishing an entity model with a realistic porous structure will do great good to study the mechanical behavior and its mechanisms. However, the pore's characterized size of coral rocks spans four orders of magnitude from micro to meso scales (1 μm to 10 mm). This requires a large specimen size and a very small element size in the model, which leads to extremely high computational costs. To overcome this fundamental challenge, this paper uses CT to scan a macro coral specimen to obtain an entity model and establishes a micro-meso scale coupled model: (1) retaining pore’s equivalent diameter larger than 1mm in the overall CT entity, and regarding the rest part as a uniformly dense matrix material with unknown material’s properties, which is called the meso model; (2) extracting a small entity containing pore’s equivalent diameter less than 1mm from the CT entity, which is called the micro model also with unknown matrix material properties; (3) the unknown matrix material parameters in meso model are derived from mechanical response of the micro model, and the unknown matrix material parameters in micro model are confined by that meso model's mechanical responses should be coincident with that of a real coral specimen. Comparing with the uniaxial compression test result, this micro-meso coupled model has high precision and retains the influence of pores at various scales. The computational result shows that it can reduce the computational complexity of a single CT entity model by four orders of magnitude. It provides a new approach for studying similar materials' mechanical behavior and mechanisms under complex loading conditions.

Effects of 2D hydrodynamic model resolution on habitat estimates for rearing Coho Salmon in contrasting channel forms

AbstractEstimating the impacts of water allocation decisions on fish populations and habitat availability is an important part of environmental flow assessments, especially in locations where water resources are limited. Two‐dimensional hydrodynamic models (2DHMs) are commonly coupled with biological models to estimate fish habitat quality, area, and capacity across a range of proposed streamflows. Increasingly, resource managers are relying on landscape‐scale model domains with coarse model resolutions to maintain feasible computational loads, but this may affect habitat estimates if the mesh element size of the model exceeds the spatial scale relevant to the organism. We investigated how coarsening the resolution of a 2DHM influences the area and spatial distribution of estimated Coho Salmon (Oncorhynchus kisutch) fry habitats. We used an interpolation scheme that upscaled mesh elements from a high‐resolution (0.25 m2) 2DHM to quantify and visualize the effects of 2DHM resolution on estimates of Coho Salmon fry habitat for two contrasting channel morphologies and across a broad range of streamflows. Estimates of Coho Salmon fry habitat at increasingly coarser resolutions led to 20%–50% reductions in weighted usable habitat area (WUA) across several streamflow scenarios for a complex channel type, but did not impact estimates in a confined, flume‐like channel. Additionally, flow‐to‐habitat area relationships were not congruent at a given streamflow when resolution coarsened. Along with almost 500% more high‐quality habitat area estimated in the complex channel type over the confined, discrepancies in habitat area increased with higher flows in areas defined as optimal for rearing Coho Salmon fry. Considering that complex channel types contain critical habitat for Coho Salmon fry, this study suggests coarse 2DHM resolutions may exclude important wetted edge and off‐channel habitats from environmental flow assessments.

Discrete-continuum-discrete approach for the modeling of the dynamic behavior of 2D lattice systems

This work presents a new methodology to model elastic lattice systems through a two-step approach that permits to reliably capture its dynamic behavior with a lower computational cost than modeling the lattice explicitly. The first step consists of a non-standard continualization accounting for scale effects. Several methods are explored to derive new continuum models, whose dispersive and vibrational behaviors are compared to that of the lattice, considered as a reference. Non-classical models with micro-inertia reveal high accuracy, not presenting physical inconsistencies. The second step follows a FEM spatial discretization of the developed continua, accounting for micro inertia terms in the mass matrix. Finally, the FEM formulation allows the use of element sizes larger (up to four times) than the physical length scale of the lattice system, thus significantly reducing the computational cost while maintaining accuracy and enabling a versatile application to materials, geometries and boundary conditions. The methodology is tested here for a 2D system with displacements in the plane, but can be extended to other lattice typologies as well.

Guidelines for element size and type selection for the finite element simulation of laser-induced elastic waves in thermoelastic laser ultrasonic testing

This paper explores spatial discretization within finite element simulations of laser-induced elastic waves within the context of Laser Ultrasonic Testing (LUT). Motivated by discrepancies and oscillations detected in temperature and displacement results in the literature, we traced these issues back to spatial discretization challenges. These challenges originate from rapid localized heating and the generation and propagation of high-frequency waves across a relatively large domain. This study effectively addresses and rectifies these inaccuracies, offering guidance for selecting the appropriate element size and type. We examined two element types: four-node quadrilaterals (Q4) employing first-order Lagrange and nine-node quadrilaterals (Q9) using second-order Lagrange shape functions. Our analysis encompasses mesh refinement strategies, exploration of time and frequency domain plots for temperature and displacement, as well as an evaluation of different pulse durations. Our findings demonstrate that Q9 elements attain accuracy with grids four times larger than Q4 elements for temperature and wave propagation analyses. Furthermore, we observe that lower frequency waves exhibit reduced sensitivity to element size, emphasizing the relationship between element size and elastic wave frequency. Pulse durations in the 6 to 30 ns range affect the required element size in the heat-affected zone but exert minimal influence on wave frequency and spatial discretization in the remainder of the domain. Finally, we present a new formula for element size selection based on the dominant frequency. This study provides a comprehensive guideline for selecting element size and type, enabling the attainment of accurate results while effectively managing computational costs.

Insights into Prime Editing Technology: A Deep Dive into Fundamentals, Potentials, and Challenges.

As the most versatile and precise gene editing technology, prime editing (PE) can establish a durable cure for most human genetic disorders. Several generations of PE have been developed based on an editor machine or prime editing guide RNA (pegRNA) to achieve any kind of genetic correction. However, due to the early stage of development, PE complex elements need to be optimized for more efficient editing. Smart optimization of editor proteins as well as pegRNA has been contemplated by many researchers, but the universal PE machine's current shortcomings remain to be solved. The modification of PE elements, fine-tuning of the host genes, manipulation of epigenetics, and blockage of immune responses could be used to reach more efficient PE. Moreover, the host factors involved in the PE process, such as repair and innate immune system genes, have not been determined, and PE cell context dependency is still poorly understood. Regarding the large size of the PE elements, delivery is a significant challenge and the development of a universal viral or nonviral platform is still far from complete. PE versions with shortened variants of reverse transcriptase are still too large to fit in common viral vectors. Overall, PE faces challenges in optimization for efficiency, high context dependency during the cell cycling, and delivery due to the large size of elements. In addition, immune responses, unpredictability of outcomes, and off-target effects further limit its application, making it essential to address these issues for broader use in nonpersonalized gene editing. Besides, due to the limited number of suitable animal models and computational modeling, the prediction of the PE process remains challenging. In this review, the fundamentals of PE, including generations, potential, optimization, delivery, invivo barriers, and the future landscape of the technology are discussed.

Analysis of load distribution on the plate and lateral hinge of a valgus opening high tibial osteotomy during weight-bearing: a finite element analysis

IntroductionValgus high tibial osteotomy (HTO) is indicated for managing isolated medial knee osteoarthritis in a young patient with a metaphyseal deformity of the proximal tibia. In a medial opening HTO, maintaining the integrity of the lateral hinge is crucial for ensuring proper healing and correction retention. Using a locked plate to stabilize an HTO is common practice, allowing for earlier weight-bearing. The objective of this study was therefore to measure and track the mechanical load distribution on a locked fixation plate and the lateral hinge of an HTO using a finite element (FE) model simulating single-leg stance loading. HypothesisThe working hypothesis was that during weight-bearing, the plate and the lateral hinge absorb stress asymmetrically, predominantly on the plate. Material and methodsA numerical model of an HTO stabilized with a locked plate was developed based on the actual geometry of a healthy proximal tibia (using Autodesk Fusion 360 and Altair HyperWorks software). In this finite element simulation of loading, a mesh convergence study was conducted to optimize the accuracy of the numerical model results. The primary outcome measure was the maximum stress value in the affected areas (Von Mises stress, in MPa) of the plate and the lateral hinge. ResultsThe maximum stress intensity in the plate was approximately 20.29 MPa. The maximum stress intensity in the bony hinge was about 5.6 MPa. The results of the mesh convergence study for the hinge and the plate enabled defining the most suitable model for future FE studies: a 4 mm mesh for all model elements except for the high-stress area in the plate and the hinge, which were meshed with a 0.7 mm element size. This adaptation provided greater precision in the study. DiscussionThere is a distribution and allocation of stress both on the plate and the hinge, underlining the significance of the plate and the absolute necessity of preserving the hinge. Predictably, the plate absorbs the majority of the load, more than three times that of the hinge. ConclusionThe hypothesis is confirmed; however, additional studies would be necessary to validate these numerical results: an experimental component on instrumented cadaveric bones, as well as comparative studies of different fixation plates. Level of EvidenceV, expert opinion; controlled laboratory study.

RESISTANCE OF FERRITE-PERLITE STEEL TUBING BRANCH PIPE METAL TO DESTRUCTION IN HYDROGEN SULFIDE MEDIUM

This paper discusses the problems of resistance of the metal of the tubing branch pipe (tubing) to the process of its destruction in a medium containing hydrogen sulfide and when exposed to external forces. A 73х73 tubing nozzle made of 38G2SFA steel was selected as the material for research. The tubing fragment was destroyed during production well operation. To establish the reasons for the failure of the tubing nozzle, a visual inspection was carried out, tests were performed to assess the hardness and mechanical properties, the microstructure and microtexture of the metal were analyzed using scanning electron microscopy (SEM) methods. To analyze the nature of the distribution of introduced defects, the mutual orientation of structural components and crystallographic texture, studies were carried out by the method electron backscatter diffraction (EBSD). Studies of the suspended branch pipe, which does not have chemical protection, showed that the cause of its destruction was hydrogen sulfide stress corrosion cracking (SSCC) as a result of hydrogen embrittlement and the action of external tensile forces exceeding the yield strength of the metal. It was established that as a result of the interaction of hydrogen sulfide of the distilled product with the metal of the branch pipe, hydrogen embrittlement of the material of the branch pipe occurred. The process of hydrogen embrittlement of the fragment was proved by detecting hydrogen blisters in the metal microstructure and establishing the process of coalescence of neighboring inclusions. In addition, the results of the test to measure the hardness of the metal along the wall thickness further supported the hypothesis of hydrogen embrittlement of the tubing. EBSD analysis showed that the threaded part of the branch pipe is characterized by a small size of structural elements and an increased density of introduced defects (geometrically necessary dislocations), which are a promising site for the deposition of molecular hydrogen. It is concluded that reaching the limit concentration of molecular hydrogen reduces the binding energy of metal atoms and this process initiates the formation of pores. It has been established that the development of cracks along large-angle grain boundaries is facilitated, and small-angle and coincidence site lattice (Σ3, Σ5, Σ11 and Σ13b) grain boundaries are active barriers to crack propagation. It has been shown that by optimizing the elemental composition of steel, controlling the microstructure, crystallographic texture, as well as the type and state of grain boundaries, new pipeline steels can be produced that are more resistant to SSCC.

Impact of MoS2 Monolayers on the Thermoelastic Response of Silicon Heterostructures.

Understanding the thermoelastic response of a nanostructure is crucial for the choice of materials and interfaces in electronic devices with improved and tailored transport properties at the nanoscale. Here, we show how the deposition of a MoS2 monolayer can strongly modify the nanoscale thermoelastic dynamics of silicon substrates close to their interface. We demonstrate this by creating a transient grating with extreme ultraviolet light, using ultrashort free-electron laser pulses, whose ≈84 nm period is comparable to the size of elements typically used in nanodevices, such as electric contacts and nanowires. The thermoelastic response, featuring coherent acoustic waves and incoherent relaxation, is tangibly modified by the presence of monolayer MoS2. Namely, we observed a major reduction of the amplitude of the surface mode, which is almost suppressed, while the longitudinal mode is basically unperturbed, aside from a faster decay of the acoustic modulations. We interpret this behavior as a selective modification of the surface elasticity, and we discuss the conditions to observe such effect, which may be of immediate relevance for the design of Si-based nanoscale devices.

Multi‐Scale Seismic Imaging of the Ridgecrest, CA, Region With Waveform Inversion of Regional and Dense Array Data

AbstractWe develop an inversion procedure for deriving multi‐scale velocity models with waveform inversions of earthquake and ambient noise data at multi‐frequency bands recorded by regional and dense sensor configurations. The method is applied for the area around the 2019 Ridgecrest earthquake rupture zones, utilizing data recorded by regional stations and dense 2D and 1D arrays with station spacings of ∼5 km and ∼100 m, respectively. Starting with regional Vp, Vs models and locations of Ridgecrest aftershocks, the velocity models and event locations are improved iteratively by inversions of waveforms recorded by regional stations and the 2D array, using a minimum spectral element size of ∼600 m. Waveforms from local events recorded by dense 1D arrays across the M7.1 rupture zone with frequencies of up to 10 Hz are used to resolve small‐scale features of the rupture zone and shallow crust with a local spectral element size of 80 m. The refined models provide self‐consistent descriptions of the rupture zone and the shallow crust embedded in the regional structures. The results reveal pronounced low Vs and high Vp/Vs in the M6.4 and M7.1 rupture zones coinciding with concentrations of seismicity, and also around the Garlock fault and in several local basins. We also observe clear velocity contrasts across the Garlock fault with polarity reversals along strike and with depth. The obtained multi‐scale velocity models can be used to improve derivations of earthquake source properties, simulations of dynamic ruptures and ground motions, and the understanding of fault and tectonic processes in the region.

Higher-order moments of spline chaos expansion

Spline chaos expansion, referred to as SCE, is a finite series representation of an output random variable in terms of measure-consistent orthonormal splines in input random variables and deterministic coefficients. This paper reports new results from an assessment of SCE’s approximation quality in calculating higher-order moments, if they exist, of the output random variable. A novel mathematical proof is provided to demonstrate that the moment of SCE of an arbitrary order converges to the exact moment for any degree of splines as the largest element size decreases. Complementary numerical analyses have been conducted, producing results consistent with theoretical findings. A collection of simple yet relevant examples is presented to grade the approximation quality of SCE with that of the well-known polynomial chaos expansion (PCE). The results from these examples indicate that higher-order moments calculated using SCE converge for all cases considered in this study. In contrast, the moments of PCE of an order larger than two may or may not converge, depending on the regularity of the output function or the probability measure of input random variables. Moreover, when both SCE- and PCE-generated moments converge, the convergence rate of the former is markedly faster than the latter in the presence of nonsmooth functions or unbounded domains of input random variables.

CFD study of heat transfer in power‐law fluids over a corrugated cylinder

AbstractThe computational study of power‐law fluid flow, along with heat transfer attributes over a corrugated heated cylinder, is explored using ANSYS FLUENT (Version 18.0). Fluids power‐law indices fall in the range of 0.25 ≤ n ≤ 1.5, and the Reynolds number spans in the range of 1 ≤ ReN ≤ 40. The flow is two‐dimensional, steady, and laminar. A wide range of Prandtl numbers (0.7 ≤ PrN ≤ 500) is used to cover the most industrially applied fluids. A domain height of 135Dh is used. A grid with the smallest element size of 0.04 m and 135,914 nodes was used. Flow and heat transfer attributes were studied using streamlines, isotherms, and local and average Nusselt numbers. The average Nusselt number increases with ReN and/or PrN. The heat transfer rate is significantly lower in dilatant fluids and higher in pseudoplastic fluids than in Newtonian fluids. The onset of wake formation behind the cylinder takes place at ReN = 10. The increase in Reynolds number (ReN) and power‐law index (n) causes an increase in wake size. Heat transfer increases with the Reynolds number and/or decrease in the power‐law index. The enhancement in heat transfer due to corrugation is studied in detail in terms of average Nusselt number, which has not been studied for arched corrugated cylinder, even for Newtonian fluids in low Reynolds number range. A Nusselt number correlation is also developed for the given ranges of conditions.

Design of Miniaturized and Wideband Four-Port MIMO Antenna Pair for WiFi.

A miniaturized and wideband four-port multiple-input multiple-output (MIMO) antenna pair for Wi-Fi mobile terminals application is proposed. The proposed antenna pair consists of four multi-branch antenna elements arranged orthogonally, with an overall size of 40 × 40 × 3.5 mm3 and each antenna element size of 15.2 × 3.5 mm × 0.8 mm3. The performance of the proposed antenna shows the advantages of a wide frequency band, low mutual coupling, high efficiency, and a compact structure. The wideband characteristics of the antenna elements are achieved through multi-mode resonance. The suppression of coupling is accomplished by strategically positioning the four compact antenna elements to ensure their maximum radiation directions are orthogonal, thus eliminating the need for an additional decoupling structure. In this paper, the proposed antenna is optimized in terms of the parameters then simulated and measured. The simulated results illustrate that an impedance bandwidth of the antenna is about 15% (5.06~5.88 GHz) with S11 < -10 dB, excellent port isolation exceeds 20 dB between all ports, a high radiation efficiency ranges from 51.2% to 89.9%, the maximum gain is 4.5 dBi, and the ECCs are less than 0.04. The measured results show that the -10 dB impedance bandwidth of the antenna is about 13% (5.13~5.80 GHz), the isolation between the antenna elements is better than 21 dB, the radiation efficiency ranges from 51.8% to 92.3%, the maximum gain is 5.3 dBi, and the ECCs are less than 0.05. The proposed four-port MIMO antenna works on the 5G LTE band 46 and Wi-Fi 6E operating bands. As a mobile terminal antenna, the proposed design scheme demonstrates excellent performance and applicability, fulfilling the requirements for 5G mobile terminal applications.

Emotional Design and Validation Study of Human–Landscape Visual Interaction

The formal beauty of “objects” is the main focus of modern rural landscapes, ignoring human interaction with the environment and the emotional reflection in this behavioral process. It is unable to satisfy the emotional needs of younger people who aspire to a high-quality life in the rural environment. The research idea of this paper is ‘first assessment—then design—then validation’. First, A 5-point Likert scale was used to investigate differences in contemporary young people’s emotional perceptions of the four rural natural landscapes in terms of instinct, behavior, and reflection. Then, using architectural design methods, a visual attraction element (viewing platform) was added by selecting samples that varied in all three dimensions (visual richness, behavioral attraction, and depth of thought). After that, a desktop eye tracker was used to record the eyeball characteristics of participants viewing the current images of natural landscapes and images of modified natural landscapes (pupil diameter, fixation duration, gaze point, etc.), and these data were combined with the subjective psychological perception scale score to determine whether or not the subjects’ positive emotions are evoked by the modified natural environment. The findings indicate that placing visually attractive elements between people and the natural world can cause subjects to feel good, think deeply, and feel more a part of the surroundings. Furthermore, we confirmed that subjects’ emotions can be evoked by 2D natural environment pictures and that the length of time subjects gaze at a picture is unaffected by the size of any individual element.

The Use of Bilinear Cohesive Zone Model to Study the Fracture of Mode I in Adhesive Joints

Adhesively bonded joints can be numericallysimulated using the Cohesive Zone Model (CZM) concepts.CZM are widely used for the strength prediction of adhesivejoints. The critical strain energy release rate and criticalinterface strength are the parameters which must be knownwhen cohesive elements in ABAQUS software are used. Theformulation of the cohesive finite elements is based on theCZM approach with the bilinear traction-separation law. Inthis work, the parameters of two industrial adhesivesHuntsman Araldite 2015 and resin LY3505/XB3405 forbonding of epoxy composites are identified. DoubleCantilever Beam (DCB) test data was used for theidentification. Finally, cohesive parameters are identifiedcomparing numerically simulated load-displacement curveswith experimental data retrieved from literature.Parametric study is performed to evaluate the variation ofinput parameters like initial stiffness, element size, peakstress and energy release rate ‘G’. From the numericalevaluation, it was noted that CZM simulation relies largelyon element size and peak cohesive strength.

Simulation of Stress Intensity Factor and J-Integral on UIC 54 Railway Failure Using A Fracture Mechanics Approach

Repeated loading on the railroad tracks will result in fatigue failure. Fatigue failure combined with a defect in the form of a crack in the railroad track will result in a decrease in strength. Defects in the form of cracks are formed due to improper manufacturing and treatment processes. One treatment that can cause the formation of defects in the form of cracks is thermite welding. Improper and non-standard thermite welding techniques can trigger the formation of defects in the form of cracks in the railroad joints. Based on these problems, this simulation aims to obtain information about the value of maximum stress, SIF, J-Intergal, number of cycles, and crack extension from variations in crack size to repeated loading. The method consists of preprocessing, processing and postprocessing. Preprocessing begins with the design of the UIC 54 railroad crack which consists of 3 variations of crack length, namely 10 mm, 15 mm, and 20 mm. The design was tested through static structural simulation using the ANSYS 2021 R2 application. Meshing is configured using an element size of 5 mm and uses curvature capture. The results of the simulation obtained maximum stress values, SIF, J-Integral, number of cycles, and crack extension. Based on the simulation of SIF 1 and J-Integral values on the specimen design with a crack length of 10 mm it shows 362.03 Mpa.mm1/2 of 0.5708 mJ/mm2, for an initial crack length of 15 mm that is equal to 482.81 Mpa.mm1/2 and 0.91738 mJ/mm2, and for an initial crack length of 20 mm, it is 600.54 Mpa.mm1/2 and 1.4465 mJ/mm2. The results show that the increase in SIF 1 and J-Integral will be proportional to the increase in the initial crack length value. .