Pub Date : 2026-01-12DOI: 10.1016/j.enganabound.2026.106640
Yanpeng Gong , Sishuai Li , Yue Mei , Bingbing Xu , Fei Qin , Xiaoying Zhuang , Timon Rabczuk
This study presents a finite element and virtual element (FE-VE) coupled method for thermomechanical analysis in electronic packaging structures. The approach partitions computational domains strategically, employing FEM for regular geometries to maximize computational efficiency and VEM for complex shapes to enhance geometric flexibility. Interface compatibility is maintained through coincident nodal correspondence, ensuring solution continuity across domain boundaries while reducing meshing complexity and computational overhead. Validation through electronic packaging applications demonstrates reasonable agreement with reference solutions and acceptable convergence characteristics across varying mesh densities. The method effectively captures thermal distributions and stress concentrations in multi-material systems, establishing a practical computational framework for electronic packaging analysis involving complex geometries. Source codes are available at https://github.com/yanpeng-gong/FeVeCoupled-ElectronicPackaging.
{"title":"A coupled finite element-virtual element method for thermomechanical analysis of electronic packaging structures","authors":"Yanpeng Gong , Sishuai Li , Yue Mei , Bingbing Xu , Fei Qin , Xiaoying Zhuang , Timon Rabczuk","doi":"10.1016/j.enganabound.2026.106640","DOIUrl":"10.1016/j.enganabound.2026.106640","url":null,"abstract":"<div><div>This study presents a finite element and virtual element (FE-VE) coupled method for thermomechanical analysis in electronic packaging structures. The approach partitions computational domains strategically, employing FEM for regular geometries to maximize computational efficiency and VEM for complex shapes to enhance geometric flexibility. Interface compatibility is maintained through coincident nodal correspondence, ensuring solution continuity across domain boundaries while reducing meshing complexity and computational overhead. Validation through electronic packaging applications demonstrates reasonable agreement with reference solutions and acceptable convergence characteristics across varying mesh densities. The method effectively captures thermal distributions and stress concentrations in multi-material systems, establishing a practical computational framework for electronic packaging analysis involving complex geometries. Source codes are available at <span><span>https://github.com/yanpeng-gong/FeVeCoupled-ElectronicPackaging</span><svg><path></path></svg></span>.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106640"},"PeriodicalIF":4.1,"publicationDate":"2026-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145956651","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-10DOI: 10.1016/j.enganabound.2026.106639
Junfeng Li , Yongtao Yang , Xiaodong Fu , Wenan Wu , Hangtian Song
Hydro-mechanical coupling in unsaturated media is of significant engineering relevance. To efficiently address quasi-static, long-term hydro-mechanical problems that are typically computationally expensive for conventional explicit schemes, a timestep scaling technique is integrated into the explicit numerical manifold method (NMM). A key advantage of this technique is that it effectively overcomes the strict stability limit of the mechanical solver. This approach aligns the critical timestep size of the mechanical field with that of the hydro field, bridging the significant timescale disparity between the two fields. This enables a unified and large timestep size within a staggered solution strategy, thereby drastically reducing the computational cost. Additionally, a simplified pressure-based explicit algorithm for unsaturated flow is embedded into the hydro-mechanical coupling framework. Verification and application examples confirm the substantial acceleration and high accuracy of the proposed explicit NMM with timestep scaling for hydro-mechanical coupling problems. Furthermore, its potential for extension to other multiphysics problems warrants continued investigation.
{"title":"An explicit numerical manifold method with timestep scaling for quasi-static unsaturated hydro-mechanical coupling problems","authors":"Junfeng Li , Yongtao Yang , Xiaodong Fu , Wenan Wu , Hangtian Song","doi":"10.1016/j.enganabound.2026.106639","DOIUrl":"10.1016/j.enganabound.2026.106639","url":null,"abstract":"<div><div>Hydro-mechanical coupling in unsaturated media is of significant engineering relevance. To efficiently address quasi-static, long-term hydro-mechanical problems that are typically computationally expensive for conventional explicit schemes, a timestep scaling technique is integrated into the explicit numerical manifold method (NMM). A key advantage of this technique is that it effectively overcomes the strict stability limit of the mechanical solver. This approach aligns the critical timestep size of the mechanical field with that of the hydro field, bridging the significant timescale disparity between the two fields. This enables a unified and large timestep size within a staggered solution strategy, thereby drastically reducing the computational cost. Additionally, a simplified pressure-based explicit algorithm for unsaturated flow is embedded into the hydro-mechanical coupling framework. Verification and application examples confirm the substantial acceleration and high accuracy of the proposed explicit NMM with timestep scaling for hydro-mechanical coupling problems. Furthermore, its potential for extension to other multiphysics problems warrants continued investigation.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106639"},"PeriodicalIF":4.1,"publicationDate":"2026-01-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145928768","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-09DOI: 10.1016/j.enganabound.2025.106625
Chenliang Li, Donglin Guo, Huihua Zhang
Magneto-electro-elastic (MEE) materials are pivotal in a multitude of fields. In this work, the numerical manifold method (NMM) is innovatively extended to establish 2-D numerical models for perforated MEE solids. The unique dual-cover system, namely, the mathematical cover and the physical cover, enables the NMM to discretize the physical domain with non-conforming mathematical covers straightforwardly. The governing equations and the boundary conditions for hole problems of MEE materials are firstly introduced. Then, by taking into account the governing equations, boundary conditions, and the NMM field approximations, the NMM global discrete equations are derived using the weighted residual method. Through three benchmark examples, the precision of the proposed method is verified, and it is then applied to study two more complex cases, where the effects of hole configurations, loading conditions, and polarization directions on the field quantities of perforated MEE materials are further examined.
{"title":"Modeling of magneto-electro-elastic solids with complex cutouts by the numerical manifold method","authors":"Chenliang Li, Donglin Guo, Huihua Zhang","doi":"10.1016/j.enganabound.2025.106625","DOIUrl":"10.1016/j.enganabound.2025.106625","url":null,"abstract":"<div><div>Magneto-electro-elastic (MEE) materials are pivotal in a multitude of fields. In this work, the numerical manifold method (NMM) is innovatively extended to establish 2-D numerical models for perforated MEE solids. The unique dual-cover system, namely, the mathematical cover and the physical cover, enables the NMM to discretize the physical domain with non-conforming mathematical covers straightforwardly. The governing equations and the boundary conditions for hole problems of MEE materials are firstly introduced. Then, by taking into account the governing equations, boundary conditions, and the NMM field approximations, the NMM global discrete equations are derived using the weighted residual method. Through three benchmark examples, the precision of the proposed method is verified, and it is then applied to study two more complex cases, where the effects of hole configurations, loading conditions, and polarization directions on the field quantities of perforated MEE materials are further examined.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106625"},"PeriodicalIF":4.1,"publicationDate":"2026-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145928752","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-09DOI: 10.1016/j.enganabound.2025.106621
Ziheng Li, Hong Zheng, Shouyang Huang
The original Discontinuous Deformation Analysis (DDA) faces limitations due to its linear displacement assumption and contact solution approach, making it challenging to handle nonlinear contact problems involving localized small deformations, particularly in cases like low-velocity impact on composite materials. This study presents an multi-point constraints enhanced DDA (MPC-DDA) specifically designed for modeling low-velocity impact scenarios in composites. The proposed modification incorporates two key features: (1) removal of the penalty spring mechanism preserves residual block penetration states in OC iterative results, (2) Iterative residuals are reconstructed as surface indentations through multi-point constraints, enabling full historical contact trace visualization. A novel contact algorithm utilizing virtual entrance points is developed, enabling MPC-DDA to dynamically capture the master and slave points on the contact boundary. This approach effectively represents surface indentations on composite materials through virtual entrance point displacements, significantly reducing the computational complexity associated with resolving localized small deformations characteristic of original DDA implementations. The MPC-DDA method has been implemented in a Matlab program, and its accuracy has been validated by comparing the computational results with those published in existing literature. Comparative analyses demonstrate that MPC-DDA outperforms original DDA in addressing nonlinear contact problems, particularly in terms of solution accuracy and computational efficiency.
{"title":"Multi-point constraints enhanced discontinuous deformation analysis for low-velocity impact on low composite materials","authors":"Ziheng Li, Hong Zheng, Shouyang Huang","doi":"10.1016/j.enganabound.2025.106621","DOIUrl":"10.1016/j.enganabound.2025.106621","url":null,"abstract":"<div><div>The original Discontinuous Deformation Analysis (DDA) faces limitations due to its linear displacement assumption and contact solution approach, making it challenging to handle nonlinear contact problems involving localized small deformations, particularly in cases like low-velocity impact on composite materials. This study presents an multi-point constraints enhanced DDA (MPC-DDA) specifically designed for modeling low-velocity impact scenarios in composites. The proposed modification incorporates two key features: (1) removal of the penalty spring mechanism preserves residual block penetration states in OC iterative results, (2) Iterative residuals are reconstructed as surface indentations through multi-point constraints, enabling full historical contact trace visualization. A novel contact algorithm utilizing virtual entrance points is developed, enabling MPC-DDA to dynamically capture the master and slave points on the contact boundary. This approach effectively represents surface indentations on composite materials through virtual entrance point displacements, significantly reducing the computational complexity associated with resolving localized small deformations characteristic of original DDA implementations. The MPC-DDA method has been implemented in a Matlab program, and its accuracy has been validated by comparing the computational results with those published in existing literature. Comparative analyses demonstrate that MPC-DDA outperforms original DDA in addressing nonlinear contact problems, particularly in terms of solution accuracy and computational efficiency.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106621"},"PeriodicalIF":4.1,"publicationDate":"2026-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145928753","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-09DOI: 10.1016/j.enganabound.2025.106627
Youjie Xu , Rui Yong , Lianjin Zhang , Fei Zhang , Zhenglin Mao , Xixiang Liu
The multiple multi-stage fractured horizontal wells (MFHWs) interference has been identified in tight gas reservoirs, but current method and model cannot achieve quantitative inter-well interference evaluation of non-uniform fractured region. Therefore, MFHWs mathematical model is established with consideration of multiple sub-region permeability difference. Coupling method of boundary element and source function is employed to solved the mathematical model. Inter-well interference factor is defined and used to evaluate inter-well interference degree. Total interference factor minimum value is used to determine the optimal rate ratio. The result shows that larger production time and well distance leads to small inter-well interference factor, but they have no influence optimal rate ratio. The optimal rate ratio will increase with the increasing of number of adjacent well fractures. If permeability of center region is larger than that of other region, inter-well interference factor and total interference factor will decrease and optimal rate ratio will increase. The model and method can evaluate inter-well interference degree and optimal rate ratio quantitatively, which provides guidance for horizontal well fracturing and parameter optimization design.
{"title":"An inter-well interference quantitative evaluation approach of multiple multi-stage fractured horizontal well with non-uniform simulated reservoirs volume in tight gas reservoirs","authors":"Youjie Xu , Rui Yong , Lianjin Zhang , Fei Zhang , Zhenglin Mao , Xixiang Liu","doi":"10.1016/j.enganabound.2025.106627","DOIUrl":"10.1016/j.enganabound.2025.106627","url":null,"abstract":"<div><div>The multiple multi-stage fractured horizontal wells (MFHWs) interference has been identified in tight gas reservoirs, but current method and model cannot achieve quantitative inter-well interference evaluation of non-uniform fractured region. Therefore, MFHWs mathematical model is established with consideration of multiple sub-region permeability difference. Coupling method of boundary element and source function is employed to solved the mathematical model. Inter-well interference factor is defined and used to evaluate inter-well interference degree. Total interference factor minimum value is used to determine the optimal rate ratio. The result shows that larger production time and well distance leads to small inter-well interference factor, but they have no influence optimal rate ratio. The optimal rate ratio will increase with the increasing of number of adjacent well fractures. If permeability of center region is larger than that of other region, inter-well interference factor and total interference factor will decrease and optimal rate ratio will increase. The model and method can evaluate inter-well interference degree and optimal rate ratio quantitatively, which provides guidance for horizontal well fracturing and parameter optimization design.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106627"},"PeriodicalIF":4.1,"publicationDate":"2026-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145928754","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-07DOI: 10.1016/j.enganabound.2025.106628
Yonghui Li , Mengyao Ma , Dinghao Zhang , Jiawei He
Soil seepage erosion caused by leaks in underground drainage pipelines is one of the key factors leading to ground subsidence. This study establishes a numerical model coupling computational fluid dynamics with the discrete element method (CFD-DEM) to conduct a microscopic investigation of seepage erosion phenomena in unsaturated silt. This numerical model employs the Richards equation to describe seepage processes in unsaturated silt and combines it with the VG model to determine parameters such as the soil's water capacity, saturation, and hydraulic conductivity. Furthermore, based on existing formulas for calculating the shear strength of unsaturated silt, a dynamic updating algorithm for particle cohesion strength was developed to account for saturation effects, thereby reproducing the moisture-induced degradation characteristics of silt. By comparing with laboratory tests, the numerical model was validated for its accuracy in simulating pore water pressure evolution, wetting front propagation rates, bonding strength updates, and soil particle loss processes. Finally, based on numerical simulation results, this study analyzed changes in formation water pressure, contact force chains, and pipeline stress conditions during the erosion-loss process. It revealed the intrinsic mechanism by which leakage-induced erosion of underground drainage pipelines triggers ground collapse at the microscopic level.
{"title":"A CFD-DEM method for simulating the erosion evolution of unsaturated silt induced by leakage through underground pipeline defect","authors":"Yonghui Li , Mengyao Ma , Dinghao Zhang , Jiawei He","doi":"10.1016/j.enganabound.2025.106628","DOIUrl":"10.1016/j.enganabound.2025.106628","url":null,"abstract":"<div><div>Soil seepage erosion caused by leaks in underground drainage pipelines is one of the key factors leading to ground subsidence. This study establishes a numerical model coupling computational fluid dynamics with the discrete element method (CFD-DEM) to conduct a microscopic investigation of seepage erosion phenomena in unsaturated silt. This numerical model employs the Richards equation to describe seepage processes in unsaturated silt and combines it with the VG model to determine parameters such as the soil's water capacity, saturation, and hydraulic conductivity. Furthermore, based on existing formulas for calculating the shear strength of unsaturated silt, a dynamic updating algorithm for particle cohesion strength was developed to account for saturation effects, thereby reproducing the moisture-induced degradation characteristics of silt. By comparing with laboratory tests, the numerical model was validated for its accuracy in simulating pore water pressure evolution, wetting front propagation rates, bonding strength updates, and soil particle loss processes. Finally, based on numerical simulation results, this study analyzed changes in formation water pressure, contact force chains, and pipeline stress conditions during the erosion-loss process. It revealed the intrinsic mechanism by which leakage-induced erosion of underground drainage pipelines triggers ground collapse at the microscopic level.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106628"},"PeriodicalIF":4.1,"publicationDate":"2026-01-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145928767","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-07DOI: 10.1016/j.enganabound.2025.106629
Ahmed Mostafa Shaaban , Elena Atroshchenko , Steffen Marburg
A discontinuous isogeometric boundary element method is proposed for the acoustic model, which is based on the Helmholtz time-harmonic wave propagation equation, as part of acoustic–structural interaction problems. The discontinuous formulation is combined with an offset-based collocation scheme, which enhances accuracy over the continuous approach while simplifying integration, the computation of free terms, and the representation of highly distorted elements in pole-based models. The dynamic structural model, representing the second component of the interaction problem, is formulated using isogeometric Reissner–Mindlin shell theory, which is particularly effective for modeling thin curved structures. The acoustic and structural models are directly coupled through conforming numerical meshes on the interaction surface. Isogeometric analysis is applied for both the acoustic and structural formulations, as well as to the coupling scheme, due to its ability to represent exact geometries and its superior accuracy compared to conventional numerical approaches.
Numerical examples are presented to assess the performance of the proposed solution, with results compared against available analytical solutions and previously reported outcomes based on alternative coupling strategies.
{"title":"Discontinuous isogeometric boundary elements for direct acoustic–structural coupling with Reissner–Mindlin shells","authors":"Ahmed Mostafa Shaaban , Elena Atroshchenko , Steffen Marburg","doi":"10.1016/j.enganabound.2025.106629","DOIUrl":"10.1016/j.enganabound.2025.106629","url":null,"abstract":"<div><div>A discontinuous isogeometric boundary element method is proposed for the acoustic model, which is based on the Helmholtz time-harmonic wave propagation equation, as part of acoustic–structural interaction problems. The discontinuous formulation is combined with an offset-based collocation scheme, which enhances accuracy over the continuous approach while simplifying integration, the computation of free terms, and the representation of highly distorted elements in pole-based models. The dynamic structural model, representing the second component of the interaction problem, is formulated using isogeometric Reissner–Mindlin shell theory, which is particularly effective for modeling thin curved structures. The acoustic and structural models are directly coupled through conforming numerical meshes on the interaction surface. Isogeometric analysis is applied for both the acoustic and structural formulations, as well as to the coupling scheme, due to its ability to represent exact geometries and its superior accuracy compared to conventional numerical approaches.</div><div>Numerical examples are presented to assess the performance of the proposed solution, with results compared against available analytical solutions and previously reported outcomes based on alternative coupling strategies.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106629"},"PeriodicalIF":4.1,"publicationDate":"2026-01-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145904247","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-06DOI: 10.1016/j.enganabound.2025.106612
Jiaxin Li , Shuihuai Yang , Xing Wei , Linlin Sun , Yue Yu
This paper presents the first attempt to apply the generalized finite difference method (GFDM) for the plane strain analysis of magneto-electro-elastic (MEE) materials. The computational domain is discretized into a set of nodes, where each node is associated with a local influence domain comprising its neighboring nodes. The variables at the nodes within each local influence domain are approximated using Taylor series expansion. A moving least squares method is then employed to establish a relationship between the partial derivatives of the variables and their values at the nodes in the local domain. Finally, the governing equations and boundary conditions are transformed into a system of linear equations. To improve the property of the coefficient matrix, each linear equation is divided by the maximum material parameter present in that equation. The feasibility and accuracy of the GFDM are validated through three test cases, with results compared against those from the meshless local Petrov–Galerkin, radial point interpolation, local radial basis function, boundary element, and finite element methods.
{"title":"Analysis of plane problems in magneto-electro-elastic media using the generalized finite difference method","authors":"Jiaxin Li , Shuihuai Yang , Xing Wei , Linlin Sun , Yue Yu","doi":"10.1016/j.enganabound.2025.106612","DOIUrl":"10.1016/j.enganabound.2025.106612","url":null,"abstract":"<div><div>This paper presents the first attempt to apply the generalized finite difference method (GFDM) for the plane strain analysis of magneto-electro-elastic (MEE) materials. The computational domain is discretized into a set of nodes, where each node is associated with a local influence domain comprising its neighboring nodes. The variables at the nodes within each local influence domain are approximated using Taylor series expansion. A moving least squares method is then employed to establish a relationship between the partial derivatives of the variables and their values at the nodes in the local domain. Finally, the governing equations and boundary conditions are transformed into a system of linear equations. To improve the property of the coefficient matrix, each linear equation is divided by the maximum material parameter present in that equation. The feasibility and accuracy of the GFDM are validated through three test cases, with results compared against those from the meshless local Petrov–Galerkin, radial point interpolation, local radial basis function, boundary element, and finite element methods.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106612"},"PeriodicalIF":4.1,"publicationDate":"2026-01-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145904246","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-03DOI: 10.1016/j.enganabound.2025.106626
Yifei Zhang
In this work, the linear-gradient smoothing meshfree Galerkin method (LGSM) is developed to analyze 2D static fracture problems of multiple materials, including isotropic materials, orthotropic materials, and functionally graded materials (FGMs). The discontinuous displacement fields and singular stress fields can be characterized by combining the diffraction method with intrinsic enrichment in the meshfree Galerkin method. However, the meshfree approximation based on intrinsic enrichment complicates numerical integration due to the non-polynomial nature of the shape function. Typically, high-order Gaussian integration is required for domain integration in the meshfree Galerkin method, but this is time-consuming. To address this, we introduce the linear-gradient smoothed integral scheme (LGSI) into the Galerkin weak form to construct the system matrix. LGSI uses only two Gaussian integration points at each boundary of the smoothing domain for the numerical integration. This process can significantly reduce the number of integral points and improves computational efficiency. Additionally, LGSI generates linear smoothed strains using the smoothed gradient technique, thereby providing accurate numerical results. The static stress intensity factors (SIFs) are determined using interaction integration. Through several numerical crack problems, the accuracy and efficiency of the LGSM are validated by comparing the results with those obtained from different methods in existing literature.
{"title":"A linear-gradient smoothing meshfree Galerkin method for 2D static fracture analysis of multiple materials","authors":"Yifei Zhang","doi":"10.1016/j.enganabound.2025.106626","DOIUrl":"10.1016/j.enganabound.2025.106626","url":null,"abstract":"<div><div>In this work, the linear-gradient smoothing meshfree Galerkin method (LGSM) is developed to analyze 2D static fracture problems of multiple materials, including isotropic materials, orthotropic materials, and functionally graded materials (FGMs). The discontinuous displacement fields and singular stress fields can be characterized by combining the diffraction method with intrinsic enrichment in the meshfree Galerkin method. However, the meshfree approximation based on intrinsic enrichment complicates numerical integration due to the non-polynomial nature of the shape function. Typically, high-order Gaussian integration is required for domain integration in the meshfree Galerkin method, but this is time-consuming. To address this, we introduce the linear-gradient smoothed integral scheme (LGSI) into the Galerkin weak form to construct the system matrix. LGSI uses only two Gaussian integration points at each boundary of the smoothing domain for the numerical integration. This process can significantly reduce the number of integral points and improves computational efficiency. Additionally, LGSI generates linear smoothed strains using the smoothed gradient technique, thereby providing accurate numerical results. The static stress intensity factors (SIFs) are determined using interaction integration. Through several numerical crack problems, the accuracy and efficiency of the LGSM are validated by comparing the results with those obtained from different methods in existing literature.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106626"},"PeriodicalIF":4.1,"publicationDate":"2026-01-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145884140","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-02DOI: 10.1016/j.enganabound.2025.106613
Wenjie Luo , Chaorong Li , Xudong Ling , Yilan Xue , L.O. Mubashiru
The CNN-based conditional encoder extracts prior knowledge for medical image segmentation to support the diffusion model. However, convolution operations are limited in capturing cross-channel and spatial dependencies, leading to the loss of crucial conditional information. Moreover, during the denoising UNet’s decoder, insufficient fusion of cross-layer features during up-sampling yields blurry feature maps and weak adaptability of subsequent convolutions. To tackle these challenges, we propose the Spatial and Channel Mixing Attention (SCMA) for the encoder and the Multi-Scale Feature Modulation Residual Module (MFMRM) for decoder feature fusion. MFMRM comprises the Multi-Scale Context Modulation Module (MCMM) and a Residual Dual Convolution Module (RDCM), adaptively integrating multi-resolution features to enhance representation and improve robustness to input variations. Furthermore, by introducing the Kolmogorov–Arnold Network (KAN) to optimize SCMA, we obtain KSCMA, which mitigates the curse of dimensionality and strengthens the representation of critical features. Experiments on ultrasound thyroid nodule, MRI brain tumor, and a self-constructed invasive breast cancer dataset demonstrate that our approach outperforms existing methods in segmentation accuracy. Our project is open source and available on GitHub at: https://github.com/lwj018/MSMedDiff-1.
{"title":"Enhancing cancer segmentation using conditional diffusion networks with KAN optimized attention and adaptive feature fusion","authors":"Wenjie Luo , Chaorong Li , Xudong Ling , Yilan Xue , L.O. Mubashiru","doi":"10.1016/j.enganabound.2025.106613","DOIUrl":"10.1016/j.enganabound.2025.106613","url":null,"abstract":"<div><div>The CNN-based conditional encoder extracts prior knowledge for medical image segmentation to support the diffusion model. However, convolution operations are limited in capturing cross-channel and spatial dependencies, leading to the loss of crucial conditional information. Moreover, during the denoising UNet’s decoder, insufficient fusion of cross-layer features during up-sampling yields blurry feature maps and weak adaptability of subsequent convolutions. To tackle these challenges, we propose the Spatial and Channel Mixing Attention (SCMA) for the encoder and the Multi-Scale Feature Modulation Residual Module (MFMRM) for decoder feature fusion. MFMRM comprises the Multi-Scale Context Modulation Module (MCMM) and a Residual Dual Convolution Module (RDCM), adaptively integrating multi-resolution features to enhance representation and improve robustness to input variations. Furthermore, by introducing the Kolmogorov–Arnold Network (KAN) to optimize SCMA, we obtain KSCMA, which mitigates the curse of dimensionality and strengthens the representation of critical features. Experiments on ultrasound thyroid nodule, MRI brain tumor, and a self-constructed invasive breast cancer dataset demonstrate that our approach outperforms existing methods in segmentation accuracy. Our project is open source and available on GitHub at: <span><span>https://github.com/lwj018/MSMedDiff-1</span><svg><path></path></svg></span>.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106613"},"PeriodicalIF":4.1,"publicationDate":"2026-01-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145884142","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
扫码关注我们
求助内容：
应助结果提醒方式：
京公网安备 11010802042870号