Search results

This study aims to investigate the possible defects and their root causes in a soft-termination multilayered ceramic capacitor (MLCC) when subjected to a thermal reflow process.

Specimens of the capacitor assembly were subjected to JEDEC level 1 preconditioning (85 °C/85%RH/168 h) with 5× reflow at 270°C peak temperature. Then, they were inspected using a 2 µm scanning electron microscope to investigate the evidence of defects. The reliability test was also numerically simulated and analyzed using the extended finite element method implemented in ABAQUS.

Excellent agreements were observed between the SEM inspections and the simulation results. The findings showed evidence of discontinuities along the Cu and the Cu-epoxy layers and interfacial delamination crack at the Cu/Cu-epoxy interface. The possible root causes are thermal mismatch between the Cu and Cu-epoxy layers, moisture contamination and weak Cu/Cu-epoxy interface. The maximum crack length observed in the experimentally reflowed capacitor was measured as 75 µm, a 2.59% difference compared to the numerical prediction of 77.2 µm.

This work's contribution is expected to reduce the additional manufacturing cost and lead time in investigating reliability issues in MLCCs.

Despite the significant number of works on the reliability assessment of surface mount capacitors, work on crack growth in soft-termination MLCC is limited. Also, the combined experimental and numerical investigation of reflow-induced reliability issues in soft-termination MLCC is limited. These cited gaps are the novelties of this study.

1) The OE-MLPG penalty meshfree method is developed to solve cracked structure.(2) Smartening the numerical meshfree method by combining the particle swarm optimization (PSO) optimization algorithms and Voronoi computational geometric algorithm. (3). Selection of base functions, finding optimal penalty factor and distribution of appropriate nodal points to the accuracy of calculation in the meshless local Petrov–Galekrin (MLPG) meshless method.

Using appropriate shape functions and distribution of nodal points in local domains and sub-domains and choosing an approximation or interpolation method has an effective role in the application of meshless methods for the analysis of computational fracture mechanics problems, especially problems with geometric discontinuity and cracks. In this research, computational geometry technique, based on the Voronoi diagram (VD) and Delaunay triangulation and PSO algorithm, are used to distribute nodal points in the sub-domain of analysis (crack line and around it on the crack plane).

By doing this process, the problems caused by too closeness of nodal points in computationally sensitive areas that exist in general methods of nodal point distribution are also solved. Comparing the effect of the number of sentences of basic functions and their order in the definition of shape functions, performing the mono-objective PSO algorithm to find the penalty factor, the coefficient, convergence, arrangement of nodal points during the three stages of VD implementation and the accuracy of the answers found indicates, the efficiency of V-E-MLPG method with Ns = 7 and ß = 0.0037–0.0075 to estimation of 3D-stress intensity factors (3D-SIFs) in computational fracture mechanics.

The present manuscript is a continuation of the studies (Ref. [33]) carried out by the authors, about; feasibility assessment, improvement and solution of challenges, introduction of more capacities and capabilities of the numerical MLPG method have been used. In order to validate the modeling and accuracy of calculations, the results have been compared with the findings of reference article [34] and [35].

This paper aims to provide SIF calculation method for engineering application.

In this paper, the stress intensity factors (SIFs) calculation method is applied to the anisotropic Ni-based single crystal film cooling holes (FCHs) structure.

Based on contour integral, the anisotropic SIFs analysis finite element method (FEM) in Ni-based single crystal is proposed. The applicability and mesh independence of the method is assessed by comparing the calculated SIFs using mode of plate with an edge crack. Anisotropic SIFs can be calculated with excellent accuracy using the finite element contour integral approach. Then, the effect of crystal orientation and FCHs interference on the anisotropic SIFs is clarified. The SIFs of FCH edge crack in the [011] orientated Ni-based single crystal increases faster than the other two orientations. And the SIF of horizontal interference FCHs edge crack is also larger than that of the inclined interference one.

The SIFs of the FCH edge crack in the turbine air-cooled blade are innovatively computed using the sub-model method. Both the Mode I and II SIFs of FCHs edge crack in blade increase with crack growing.

In this paper, the standard Peridynamic Timoshenko beam model accounting for the shear deformation is chosen to describe the thick beam kinematics. Unfortunately, when applied to very thin beam structures, the standard Peridynamics (PD) encounters the shear locking phenomenon, leading to incorrect solutions.

PD differs from classical continuum mechanics and other nonlocal theories that do not involve spatial derivatives of the displacement field. PD is based on the integral equation instead of differential equations to handle discontinuities and other singularities.

The shear locking can be successfully alleviated using the developed selective integration method. In particular, this technique has been implemented in the standard PD, which allows an accurate result for a wide range of slenderness from very thin to thick (10 < L/t < 10³) structures. It can also accelerate the computational time for particular dynamic problems using fewer neighboring integration particles. Several numerical examples are solved to demonstrate the effectiveness of the proposed method for modeling beam structures.

The paper highlights the severe shear locking phenomenon in the Peridynamic Timoshenko beam available in the literature, especially for very thin structures. A new alternative for the alleviation of shear locking in the Peridynamic Timoshenko beam, using selective integration. Hence the developed Peridynamic Timoshenko beam model is effective for thin and thick structures. A new peridynamic formulation for the low-velocity impact beam models is presented and validated.

1. The paper highlights the severe shear locking phenomenon in the Peridynamic Timoshenko beam proposed in the literature, especially for very thin structures.

2. The developed Peridynamic Timoshenko beam model based on selective integration is effective for thin and thick structures.

3. A new peridynamic formulation for the low-velocity impact beam models is presented and validated.

The paper highlights the severe shear locking phenomenon in the Peridynamic Timoshenko beam proposed in the literature, especially for very thin structures.

The developed Peridynamic Timoshenko beam model based on selective integration is effective for thin and thick structures.

A new peridynamic formulation for the low-velocity impact beam models is presented and validated.

This paper aims to employ a multivariate nonlinear regression analysis to establish a predictive model for the final fracture area, while accounting for the impact of individual parameters.

This analysis is based on the numerical simulation data obtained, using the hybrid finite element–discrete element (FE–DE) method. The forecasting model was compared with the numerical results and the accuracy of the model was evaluated by the root mean square (RMS) and the RMS error, the mean absolute error and the mean absolute percentage error.

The multivariate nonlinear regression model can accurately predict the nonlinear relationships between injection rate, leakoff coefficient, elastic modulus, permeability, Poisson’s ratio, pore pressure and final fracture area. The regression equations obtained from the Newton iteration of the least squares method are strong in terms of the fit to the six sensitive parameters, and the model follow essentially the same trend with the numerical simulation data, with no systematic divergence detected. Least absolutely deviation has a significantly weaker performance than the least squares method. The percentage contribution of sensitive parameters to the final fracture area is available from the simulation results and forecast model. Injection rate, leakoff coefficient, permeability, elastic modulus, pore pressure and Poisson’s ratio contribute 43.4%, −19.4%, 24.8%, −19.2%, −21.3% and 10.1% to the final fracture area, respectively, as they increased gradually. In summary, (1) the fluid injection rate has the greatest influence on the final fracture area. (2)The multivariate nonlinear regression equation was optimally obtained after 59 iterations of the least squares-based Newton method and 27 derivative evaluations, with a decidability coefficient R2 = 0.711 representing the model reliability and the regression equations fit the four parameters of leakoff coefficient, permeability, elastic modulus and pore pressure very satisfactorily. The models follow essentially the identical trend with the numerical simulation data and there is no systematic divergence. The least absolute deviation has a significantly weaker fit than the least squares method. (3)The nonlinear forecasting model of physical parameters of hydraulic fracturing established in this paper can be applied as a standard for optimizing the fracturing strategy and predicting the fracturing efficiency in situ field and numerical simulation. Its effectiveness can be trained and optimized by experimental and simulation data, and taking into account more basic data and establishing regression equations, containing more fracturing parameters will be the further research interests.

The nonlinear forecasting model of physical parameters of hydraulic fracturing established in this paper can be applied as a standard for optimizing the fracturing strategy and predicting the fracturing efficiency in situ field and numerical simulation. Its effectiveness can be trained and optimized by experimental and simulation data, and taking into account more basic data and establishing regression equations, containing more fracturing parameters will be the further research interests.

This study provides an analysis of patch repair for cracked aircraft structures. Delamination is a type of damage that affects the patch's behavior. The purpose of this study is to assess the influence of delamination on repair performance.

An analytical and numerical study using the finite element method was conducted for a cracked plate repaired with a patch containing a pre-existing delamination defect. The method for defining the contact pair surfaces and modeling the delamination interaction within the patch interface is specified using the virtual crack closure technique (VCCT) approach.

The efficiency of the repair is measured in terms of the J-integral. The effects of delamination initiation, mechanical loading, crack length and patch stacking sequences are presented. It is noted that in mode I, delamination propagation is only significant at node A. The numerical results are in good agreement with those of the analytical solution found in the literature. It is observed that the patch's behavior is strongly dependent on loading, crack size and stacking sequences in terms of reducing the structure's lifespan, especially in the presence of delamination.

The numerical modeling presented by the VCCT approach is highly valuable for studying delamination evolution. The influence of loading, crack size and stacking sequences on repair performance is discussed in this work.

This paper aims to investigate the performance of two novel numerical methods, the face-based smoothed finite element method (FS-FEM) and the edge-based smoothed finite element method (ES-FEM), which employ linear tetrahedral elements, for the purpose of strength assessment of a high-speed train hollow axle.

The calculation of stress for the wheelset, comprising an axle and two wheels, is facilitated through the application of the European axle strength design standard. This standard assists in the implementation of loading and boundary conditions and is exemplified by the typical CRH2 high-speed train wheelset. To evaluate the performance of these two methods, a hollow cylinder cantilever beam is first used as a benchmark to compare the present methods with other existing methods. Then, the strength analysis of a real wheelset model with a hollow axle is performed using different numerical methods.

The results of deflection and stress show that FS-FEM and ES-FEM offer higher accuracy and better convergence than FEM using linear tetrahedral elements. ES-FEM exhibits a superior performance to that of FS-FEM using linear tetrahedral elements, showing accuracy and convergence close to FEM using hexahedral elements.

This study channels the novel methods (FS-FEM and ES-FEM) in the static stress analysis of a railway wheelset. Based on the careful testing of FS-FEM and ES-FEM, both methods hold promise as more efficient tools for the strength analysis of complex railway structures.

Multi-well hydrofracturing is an important technology to create new fractures and expand existing fractures to increase reservoir permeability. The propagation morphology of the fracture network is affected by the disturbance between the fractures initiation sequences and spacings between adjacent wells. However, it remains unclear how well spacing and initiation sequences lead to fracture propagation, deflection and connection.

In this study, the thermal-hydro-mechanical coupling effect in the hydrofracturing process was considered, to establish a finite element-discrete element model of multistage hydrofracturing in a horizontal well. Using typical cases, the unstable propagation of hydraulic fractures in multiple horizontal wells was investigated under varying well spacing and initiation sequences. Combined with the shear stress shadow caused by in situ stress disturbed by fracture tip propagation, the quantitative indexes of fracture propagation such as length, volume, displacement vector, deflection and unstable propagation behavior of the hydrofracturing fracture network were analyzed.

The results show that the shear stress disturbance caused by multiple hydraulic fractures is a significant factor in multi-well hydrofracturing. Reducing the spacing between multiple wells increases the stress shadow area and aggravates the mutual disturbance and deflection between the fractures. The quantitative analysis results show that a decrease of well spacing reduces the total length of hydraulic fractures but increases the total volume of the fracture; compared with sequential and simultaneous fracturing, alternate fracturing can effectively reduce stress shadow area, alleviate fracture disturbance and generate larger fracture propagation length and volume.

The numerical models and results of the unstable propagation and stress evolution of the hydraulic fracture network under thermal-hydro-mechanical coupling obtained in this study can provide useful guidance for the evaluation and design of rock mass fracture networks in deep unconventional oil and gas reservoirs.

The purpose of this paper is to apply the concept of equivalent initial flaw size (EIFS) to the anisotropic nickel-based single crystal (SX) material, and to predict the fatigue life on this basis. The crack propagation law of SX material at different temperatures and the weak correlation of EIFS values verification under different loading conditions are also investigated.

A three-parameter time to crack initial (TTCI) method with multiple reference crack lengths under different loading conditions is established, which include the TTCI backstepping method and EIFS fitting method. Subsequently, the optimized EIFS distribution is obtained based on the random crack propagation rate and maximum likelihood estimation of median fatigue life. Then, an effective driving force based on anisotropic and mixed crack propagation mode is proposed to describe the crack propagation rate in the small crack stage. Finally, the fatigue life of three different temperature ESE(T) standard specimens is predicted based on the EIFS values under different survival rates.

The optimized EIFS distribution based on EIFS fitting - maximum likelihood estimation (MLE) method has the highest accuracy in predicting the total fatigue life, with the range of EIFS values being about [0.0028, 0.0875] (mm), and the mean value of EIFS being 0.0506 mm. The error between the predicted fatigue life based on the crack propagation rate and EIFS distribution for survival rates ranges from 5% to 95% and the experimental life is within two times dispersion band.

This paper systematically proposes a new anisotropic material EIFS prediction method, establishing a framework for predicting the fatigue life of SX material at different temperatures using fracture mechanics to avoid inaccurate anisotropic constitutive models and fatigue damage accumulation theory.