首页 > 最新文献

Finite Elements in Analysis and Design最新文献

英文 中文
Enhancing data representation in forging processes: Investigating discretization and R-adaptivity strategies with Proper Orthogonal Decomposition reduction 加强锻造过程中的数据表示:利用适当正交分解还原法研究离散化和 R-自适应策略
IF 3.5 3区 工程技术 Q1 MATHEMATICS, APPLIED Pub Date : 2024-11-07 DOI: 10.1016/j.finel.2024.104276
David Uribe, Camille Durand, Cyrille Baudouin, Régis Bigot
Effective data reduction techniques are crucial for enhancing computational efficiency in complex industrial processes such as forging. In this study, we investigate various discretization and mesh adaptivity strategies using Proper Orthogonal Decomposition (POD) to optimize data reduction fidelity in forging simulations. We focus particularly on r-adaptivity techniques, which ensure a consistent number of elements throughout the field representation, filling a gap in existing research that predominantly concentrates on h-adaptivity. Our investigation compares isotropic mesh approaches with anisotropic mesh adaptations, including gradient-based, isolines-based, and spring-energy-based methods. Through numerical simulations and analysis, we demonstrate that these anisotropic techniques provide superior fidelity in representing deformation fields compared to isotropic meshes. These improvements are achieved while maintaining a similar level of model reduction efficiency. This enhancement in representation leads to improved data reduction quality, forming the foundation for data-driven models. This research contributes to advancing the understanding of mesh adaptivity approaches and their potential applications in data-driven modeling across various industrial domains.
有效的数据缩减技术对于提高锻造等复杂工业流程的计算效率至关重要。在本研究中,我们利用适当正交分解(POD)研究了各种离散化和网格自适应策略,以优化锻造模拟中的数据缩减保真度。我们特别关注 r 自适应技术,它能确保整个场表示中元素数量的一致性,填补了现有研究中主要集中在 h 自适应方面的空白。我们的研究比较了各向同性网格方法和各向异性网格适应性,包括基于梯度的方法、基于孤立线的方法和基于弹簧能量的方法。通过数值模拟和分析,我们证明了与各向同性网格相比,这些各向异性技术在表示变形场方面具有更高的保真度。在实现这些改进的同时,还保持了类似水平的模型缩减效率。这种表示方法的改进提高了数据还原质量,为数据驱动模型奠定了基础。这项研究有助于加深对网格自适应方法及其在各工业领域数据驱动建模中的潜在应用的理解。
{"title":"Enhancing data representation in forging processes: Investigating discretization and R-adaptivity strategies with Proper Orthogonal Decomposition reduction","authors":"David Uribe,&nbsp;Camille Durand,&nbsp;Cyrille Baudouin,&nbsp;Régis Bigot","doi":"10.1016/j.finel.2024.104276","DOIUrl":"10.1016/j.finel.2024.104276","url":null,"abstract":"<div><div>Effective data reduction techniques are crucial for enhancing computational efficiency in complex industrial processes such as forging. In this study, we investigate various discretization and mesh adaptivity strategies using Proper Orthogonal Decomposition (POD) to optimize data reduction fidelity in forging simulations. We focus particularly on r-adaptivity techniques, which ensure a consistent number of elements throughout the field representation, filling a gap in existing research that predominantly concentrates on h-adaptivity. Our investigation compares isotropic mesh approaches with anisotropic mesh adaptations, including gradient-based, isolines-based, and spring-energy-based methods. Through numerical simulations and analysis, we demonstrate that these anisotropic techniques provide superior fidelity in representing deformation fields compared to isotropic meshes. These improvements are achieved while maintaining a similar level of model reduction efficiency. This enhancement in representation leads to improved data reduction quality, forming the foundation for data-driven models. This research contributes to advancing the understanding of mesh adaptivity approaches and their potential applications in data-driven modeling across various industrial domains.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"242 ","pages":"Article 104276"},"PeriodicalIF":3.5,"publicationDate":"2024-11-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142655623","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Impact of surface roughness on the formation of necking instabilities in additive manufactured porous metal plates subjected to dynamic plane strain stretching 表面粗糙度对受到动态平面应变拉伸的添加剂制造多孔金属板形成颈缩不稳定性的影响
IF 3.5 3区 工程技术 Q1 MATHEMATICS, APPLIED Pub Date : 2024-11-06 DOI: 10.1016/j.finel.2024.104275
M. Anil Kumar , J.C. Nieto-Fuentes , J.A. Rodríguez-Martínez
This paper investigates the influence of surface roughness on multiple necking formation in additive manufactured porous ductile plates subjected to dynamic plane strain stretching. For this purpose, we have developed a computational model in ABAQUS/Explicit which includes surface texture and discrete voids measured from 3D-printed metallic specimens using optical profilometry and X-ray tomography analysis, respectively. The mechanical behavior of the material is described using an elastic–plastic constitutive model, with yielding defined by the isotropic von Mises criterion, an associated flow rule, and a power-law function for the yield stress evolution which depends on plastic strain, plastic strain rate, and temperature. The finite element calculations have been conducted across a broad range of strain rates, from 5000s1 to 50000s1, to explore the interactions among inertia, surface roughness, and porosity in determining the necking pattern that emerges in the plates at large strains. The finite element results show that surface roughness induces perturbations in the deformation field of the specimen, which lead to early necking localization, while the location and number of necks formed are primarily controlled by the porous microstructure and the loading rate. The results for the neck spacing have shown quantitative agreement with the analytical stability analysis predictions and the unit-cell finite element calculations reported by Rodríguez-Martínez et al. [1]. Moreover, integrating discrete voids into simulations that already account for surface roughness results in a minor reduction in necking strain: surface roughness and porosity demonstrate similar quantitative impacts on necking ductility, which is primarily influenced by inertia effects at the highest strain rates studied. To the best of the authors’ knowledge, this paper presents the first calculations that explore dynamic plastic localization in additive manufactured metals, incorporating actual surface roughness and explicit void representation derived from experimental measurements. This work marks progress in the analysis of 3D-printed structures under impact loading, aiming to understand and predict the mechanics influencing their energy absorption capacity at high strain rates.
本文研究了表面粗糙度对受到动态平面应变拉伸的添加剂制造的多孔韧性板中多颈形成的影响。为此,我们在 ABAQUS/Explicit 中开发了一个计算模型,其中包括分别使用光学轮廓仪和 X 射线断层扫描分析法从三维打印金属试样中测量的表面纹理和离散空隙。材料的机械行为采用弹塑性构成模型进行描述,屈服由各向同性冯-米塞斯准则、相关流动规则和屈服应力演变的幂律函数(取决于塑性应变、塑性应变率和温度)定义。有限元计算的应变率范围很广,从 5000s-1 到 50000s-1,以探索惯性、表面粗糙度和孔隙率之间的相互作用,从而确定板材在大应变下出现的颈缩模式。有限元结果表明,表面粗糙度会引起试样变形场的扰动,从而导致早期颈缩局部化,而形成颈缩的位置和数量主要受多孔微结构和加载速率的控制。颈部间距的计算结果与 Rodríguez-Martínez 等人[1]报告的稳定性分析预测和单元有限元计算结果在数量上一致。此外,将离散空隙整合到已考虑表面粗糙度的模拟中会导致颈部应变的轻微降低:表面粗糙度和孔隙率对颈部延展性的定量影响相似,在研究的最高应变速率下,颈部延展性主要受惯性效应的影响。据作者所知,本文首次提出了在增材制造金属中探索动态塑性定位的计算方法,并结合了实际表面粗糙度和实验测量得出的明确空隙表示。这项工作标志着三维打印结构在冲击载荷下的分析取得了进展,其目的是了解和预测在高应变速率下影响其能量吸收能力的力学原理。
{"title":"Impact of surface roughness on the formation of necking instabilities in additive manufactured porous metal plates subjected to dynamic plane strain stretching","authors":"M. Anil Kumar ,&nbsp;J.C. Nieto-Fuentes ,&nbsp;J.A. Rodríguez-Martínez","doi":"10.1016/j.finel.2024.104275","DOIUrl":"10.1016/j.finel.2024.104275","url":null,"abstract":"<div><div>This paper investigates the influence of surface roughness on multiple necking formation in additive manufactured porous ductile plates subjected to dynamic plane strain stretching. For this purpose, we have developed a computational model in ABAQUS/Explicit which includes surface texture and discrete voids measured from 3D-printed metallic specimens using optical profilometry and X-ray tomography analysis, respectively. The mechanical behavior of the material is described using an elastic–plastic constitutive model, with yielding defined by the isotropic von Mises criterion, an associated flow rule, and a power-law function for the yield stress evolution which depends on plastic strain, plastic strain rate, and temperature. The finite element calculations have been conducted across a broad range of strain rates, from <span><math><mrow><mn>5000</mn><mspace></mspace><msup><mrow><mtext>s</mtext></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow></math></span> to <span><math><mrow><mn>50000</mn><mspace></mspace><msup><mrow><mtext>s</mtext></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow></math></span>, to explore the interactions among inertia, surface roughness, and porosity in determining the necking pattern that emerges in the plates at large strains. The finite element results show that surface roughness induces perturbations in the deformation field of the specimen, which lead to early necking localization, while the location and number of necks formed are primarily controlled by the porous microstructure and the loading rate. The results for the neck spacing have shown quantitative agreement with the analytical stability analysis predictions and the unit-cell finite element calculations reported by Rodríguez-Martínez et al. <span><span>[1]</span></span>. Moreover, integrating discrete voids into simulations that already account for surface roughness results in a minor reduction in necking strain: surface roughness and porosity demonstrate similar quantitative impacts on necking ductility, which is primarily influenced by inertia effects at the highest strain rates studied. To the best of the authors’ knowledge, this paper presents the first calculations that explore dynamic plastic localization in additive manufactured metals, incorporating actual surface roughness and explicit void representation derived from experimental measurements. This work marks progress in the analysis of 3D-printed structures under impact loading, aiming to understand and predict the mechanics influencing their energy absorption capacity at high strain rates.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"242 ","pages":"Article 104275"},"PeriodicalIF":3.5,"publicationDate":"2024-11-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142593691","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Investigation of nonlinear buckling of FGM shells using a high-order finite continuation approach 利用高阶有限延续方法研究 FGM 壳体的非线性屈曲
IF 3.5 3区 工程技术 Q1 MATHEMATICS, APPLIED Pub Date : 2024-11-06 DOI: 10.1016/j.finel.2024.104273
Oussama Elmhaia , Omar Askour , Yassir Sitli , Said Mesmoudi , Mohammed Rammane , Oussama Bourihane , Youssef Hilali
This study investigates the buckling behavior of cylindrical shells composed of Functionally Graded Materials (FGMs) when subjected to axial compression, challenging conventional assumptions regarding the influence of Poisson’s effect in homogeneous materials. To address this, we utilize a numerical approach employing the Asymptotic Numerical Method (ANM). Contrary to the expected linear pre-buckling behavior associated with a zero Poisson’s ratio, our findings reveal significant non-linearity in the response of FGM structures, emphasizing the influence of additional non-linear factors inherent in the behavior of advanced composites. Through an extensive numerical analysis conducted using a customized Matlab code, we examine the buckling and post-buckling characteristics of FGM shells with varying surface compositions, particularly focusing on configurations incorporating Al2O3 and Al on the upper surface. To elucidate our findings, we present numerical examples comparing two FGM scenarios (Al2O3/Al and Al/Al2O3) in terms of critical buckling and FGM distribution. Additionally, we validate our results by employing the commercial software Abaqus with Riks-based finite element method and Newton–Raphson solver.
本研究探讨了由功能分级材料(FGM)构成的圆柱形壳体在受到轴向压缩时的屈曲行为,挑战了有关均质材料中泊松效应影响的传统假设。为了解决这个问题,我们采用了渐近数值方法(ANM)。与预期的与零泊松比相关的线性预屈曲行为相反,我们的研究结果揭示了 FGM 结构响应中的显著非线性,强调了先进复合材料行为中固有的其他非线性因素的影响。通过使用定制的 Matlab 代码进行广泛的数值分析,我们研究了不同表面成分的 FGM 壳体的屈曲和屈曲后特性,尤其侧重于上表面含有 Al2O3 和 Al 的配置。为了阐明我们的发现,我们通过数值示例比较了两种 FGM 方案(Al2O3/Al 和 Al/Al2O3)的临界屈曲和 FGM 分布。此外,我们还利用基于 Riks 的有限元法和牛顿-拉斐逊求解器的商业软件 Abaqus 验证了我们的结果。
{"title":"Investigation of nonlinear buckling of FGM shells using a high-order finite continuation approach","authors":"Oussama Elmhaia ,&nbsp;Omar Askour ,&nbsp;Yassir Sitli ,&nbsp;Said Mesmoudi ,&nbsp;Mohammed Rammane ,&nbsp;Oussama Bourihane ,&nbsp;Youssef Hilali","doi":"10.1016/j.finel.2024.104273","DOIUrl":"10.1016/j.finel.2024.104273","url":null,"abstract":"<div><div>This study investigates the buckling behavior of cylindrical shells composed of Functionally Graded Materials (FGMs) when subjected to axial compression, challenging conventional assumptions regarding the influence of Poisson’s effect in homogeneous materials. To address this, we utilize a numerical approach employing the Asymptotic Numerical Method (ANM). Contrary to the expected linear pre-buckling behavior associated with a zero Poisson’s ratio, our findings reveal significant non-linearity in the response of FGM structures, emphasizing the influence of additional non-linear factors inherent in the behavior of advanced composites. Through an extensive numerical analysis conducted using a customized Matlab code, we examine the buckling and post-buckling characteristics of FGM shells with varying surface compositions, particularly focusing on configurations incorporating <span><math><mrow><msub><mrow><mi>Al</mi></mrow><mrow><mn>2</mn></mrow></msub><msub><mrow><mi>O</mi></mrow><mrow><mn>3</mn></mrow></msub></mrow></math></span> and <span><math><mi>Al</mi></math></span> on the upper surface. To elucidate our findings, we present numerical examples comparing two FGM scenarios (<span><math><mrow><msub><mrow><mi>Al</mi></mrow><mrow><mn>2</mn></mrow></msub><msub><mrow><mi>O</mi></mrow><mrow><mn>3</mn></mrow></msub><mo>/</mo><mi>Al</mi></mrow></math></span> and <span><math><mrow><mi>Al</mi><mo>/</mo><msub><mrow><mi>Al</mi></mrow><mrow><mn>2</mn></mrow></msub><msub><mrow><mi>O</mi></mrow><mrow><mn>3</mn></mrow></msub></mrow></math></span>) in terms of critical buckling and FGM distribution. Additionally, we validate our results by employing the commercial software Abaqus with Riks-based finite element method and Newton–Raphson solver.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"242 ","pages":"Article 104273"},"PeriodicalIF":3.5,"publicationDate":"2024-11-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142593690","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Dual failure analysis of 3D structures under cyclic loads using bFS-FEM based numerical approaches 使用基于 bFS-FEM 的数值方法对循环载荷下的 3D 结构进行双重失效分析
IF 3.5 3区 工程技术 Q1 MATHEMATICS, APPLIED Pub Date : 2024-11-02 DOI: 10.1016/j.finel.2024.104272
Phuc L.H. Ho , Canh V. Le , Changkye Lee , Dung T. Tran , Phuong H. Nguyen , Jurng-Jae Yee
Failure mechanism of 3D structures cannot always be produced by the low-order finite elements due to the so-called volumetric locking effect. In this paper, dual numerical approaches based on the bubble face-based smoothed finite element method (bFS-FEM) are developed, ensuring that the locking problem is prevented and accurate load factors of elastic-perfectly plastic structures under cyclic actions are achieved. The failure mechanisms, in terms of plastic dissipation, are realized as incremental or alternative plastic failure modes, enabling different treatments in engineering practices. Moreover, the pseudo-static approach is capable of providing three-dimensional stress fields at the failure state, which is crucial for structural design. Interaction diagrams associated with various load-types and-ranges are illustrated in numerical experiments, showing that the bearing capacity envelopes of structures under cyclic loads are evidently smaller than that of proportional loads.
由于所谓的体积锁定效应,三维结构的失效机制并不总是能由低阶有限元生成。本文开发了基于气泡面平滑有限元法(bFS-FEM)的双重数值方法,确保防止锁定问题,并实现循环作用下弹性全塑结构的精确载荷系数。从塑性耗散的角度来看,失效机理可以实现为增量或替代性塑性失效模式,从而在工程实践中实现不同的处理方法。此外,伪静态方法能够提供失效状态下的三维应力场,这对结构设计至关重要。数值实验展示了与各种荷载类型和范围相关的相互作用图,表明循环荷载下结构的承载能力包络明显小于比例荷载下的承载能力包络。
{"title":"Dual failure analysis of 3D structures under cyclic loads using bFS-FEM based numerical approaches","authors":"Phuc L.H. Ho ,&nbsp;Canh V. Le ,&nbsp;Changkye Lee ,&nbsp;Dung T. Tran ,&nbsp;Phuong H. Nguyen ,&nbsp;Jurng-Jae Yee","doi":"10.1016/j.finel.2024.104272","DOIUrl":"10.1016/j.finel.2024.104272","url":null,"abstract":"<div><div>Failure mechanism of 3D structures cannot always be produced by the low-order finite elements due to the so-called volumetric locking effect. In this paper, dual numerical approaches based on the bubble face-based smoothed finite element method (bFS-FEM) are developed, ensuring that the locking problem is prevented and accurate load factors of elastic-perfectly plastic structures under cyclic actions are achieved. The failure mechanisms, in terms of plastic dissipation, are realized as incremental or alternative plastic failure modes, enabling different treatments in engineering practices. Moreover, the pseudo-static approach is capable of providing three-dimensional stress fields at the failure state, which is crucial for structural design. Interaction diagrams associated with various load-types and-ranges are illustrated in numerical experiments, showing that the bearing capacity envelopes of structures under cyclic loads are evidently smaller than that of proportional loads.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"242 ","pages":"Article 104272"},"PeriodicalIF":3.5,"publicationDate":"2024-11-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142573362","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
3D analysis of reinforced concrete structural components using a multi-surface elasto-plastic-anisotropic-damage material model 利用多面弹塑性-各向异性-损伤材料模型对钢筋混凝土结构部件进行三维分析
IF 3.5 3区 工程技术 Q1 MATHEMATICS, APPLIED Pub Date : 2024-10-30 DOI: 10.1016/j.finel.2024.104271
A. Torabizadeh , A. Sarikaya , R.E. Erkmen
Elastic-Plastic-Damage material models are widely adopted for the numerical modelling of concrete because of their capability of representing pressure sensitive 3D material behaviour considering permanent inelastic deformations as well as degradation of material moduli beyond the elastic range. In this paper, we develop a non-associative multi-surface plastic-damage material model for the 3D solid element based finite element analysis of reinforced concrete structural components. For the non-associative plastic flow, a linear potential function is adopted, while Menetrey–Willam and Rankine surfaces are adopted as the yield surfaces in compression and tension regimes, respectively. The degradation in the material stiffness under cyclic loading is incorporated by the damage component of the material model, which is generally anisotropic and assumed to be directly dependent on the evolution of the plastic strains. This assumption leads to a computationally efficient algorithm in terms of circumventing iterations to equate the stresses between the coupled damage and plasticity components of the material model. The rigorous details of the developed return-mapping methodology considering both the Cutting-Plane as well as the Closest-Point-Projection algorithms are provided. The material model is employed for the structural level analysis, in which case the concrete bulk is modelled by using an Eight-Node, Six-Degrees-Of-Freedom per-node solid element, and the reinforcement bars and stirrups are modelled by using the conventional Two-Node, Six-Degrees-Of-Freedom per-node Euler–Bernoulli beam-bar element. The inelastic behaviour of the reinforcements is determined by using a simpler elasto-plastic-damage based material model under the assumption of uni-axial stress-strain relations. An in-house fortran software is developed for the computer implementation. Comparisons with results from literature are shown for validation purposes. The validation cases include static analyses of a beam and a column under monotonic loading as well as a shear-wall under cyclic loading.
弹性-塑性-损伤材料模型被广泛用于混凝土的数值建模,因为这些模型能够表示对压力敏感的三维材料行为,并考虑到永久非弹性变形以及超出弹性范围的材料模量退化。在本文中,我们为基于三维实体元素的钢筋混凝土结构组件有限元分析开发了一种非关联多表面塑性损伤材料模型。在非关联塑性流动中,采用了线性势函数,而在压缩和拉伸状态下,分别采用 Menetrey-Willam 和 Rankine 表面作为屈服面。材料模型的损伤部分包含了循环加载下材料刚度的衰减,损伤部分通常是各向异性的,并假定直接取决于塑性应变的演变。这一假设导致了一种计算高效的算法,即通过规避迭代来平衡材料模型的耦合损伤和塑性成分之间的应力。本文提供了所开发的返回映射方法的严格细节,其中既考虑了切削平面算法,也考虑了闭合点投影算法。该材料模型用于结构层面的分析,在这种情况下,混凝土体采用八节点、六自由度、每节点实体元素建模,钢筋和箍筋采用传统的两节点、六自由度、每节点欧拉-伯努利梁-杆元素建模。钢筋的非弹性行为是在单轴应力-应变关系的假设下,使用较简单的弹塑性损伤材料模型确定的。为计算机实施开发了一个内部 fortran 软件。与文献中的结果进行了比较,以进行验证。验证案例包括单调荷载下梁和柱的静态分析,以及循环荷载下剪力墙的静态分析。
{"title":"3D analysis of reinforced concrete structural components using a multi-surface elasto-plastic-anisotropic-damage material model","authors":"A. Torabizadeh ,&nbsp;A. Sarikaya ,&nbsp;R.E. Erkmen","doi":"10.1016/j.finel.2024.104271","DOIUrl":"10.1016/j.finel.2024.104271","url":null,"abstract":"<div><div>Elastic-Plastic-Damage material models are widely adopted for the numerical modelling of concrete because of their capability of representing pressure sensitive 3D material behaviour considering permanent inelastic deformations as well as degradation of material moduli beyond the elastic range. In this paper, we develop a non-associative multi-surface plastic-damage material model for the 3D solid element based finite element analysis of reinforced concrete structural components. For the non-associative plastic flow, a linear potential function is adopted, while Menetrey–Willam and Rankine surfaces are adopted as the yield surfaces in compression and tension regimes, respectively. The degradation in the material stiffness under cyclic loading is incorporated by the damage component of the material model, which is generally anisotropic and assumed to be directly dependent on the evolution of the plastic strains. This assumption leads to a computationally efficient algorithm in terms of circumventing iterations to equate the stresses between the coupled damage and plasticity components of the material model. The rigorous details of the developed return-mapping methodology considering both the Cutting-Plane as well as the Closest-Point-Projection algorithms are provided. The material model is employed for the structural level analysis, in which case the concrete bulk is modelled by using an Eight-Node, Six-Degrees-Of-Freedom per-node solid element, and the reinforcement bars and stirrups are modelled by using the conventional Two-Node, Six-Degrees-Of-Freedom per-node Euler–Bernoulli beam-bar element. The inelastic behaviour of the reinforcements is determined by using a simpler elasto-plastic-damage based material model under the assumption of uni-axial stress-strain relations. An in-house fortran software is developed for the computer implementation. Comparisons with results from literature are shown for validation purposes. The validation cases include static analyses of a beam and a column under monotonic loading as well as a shear-wall under cyclic loading.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"242 ","pages":"Article 104271"},"PeriodicalIF":3.5,"publicationDate":"2024-10-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142554598","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Efficient thermal modeling of laser directed energy deposition using the forward Euler scheme: Methodology, merits and limitations 使用前向欧拉方案对激光定向能沉积进行高效热建模:方法、优点和局限性
IF 3.5 3区 工程技术 Q1 MATHEMATICS, APPLIED Pub Date : 2024-10-21 DOI: 10.1016/j.finel.2024.104270
Simon Essongue , Vaibhav Nain , Muriel Carin
This paper explores mesoscale conduction-based modeling of Laser Directed Energy Deposition (LDED) for metallic materials. We benchmark the forward Euler (explicit) time integration strategy against the backward Euler (implicit) scheme using two experimentally validated simulations. Our results demonstrate the explicit scheme’s faster computational speed. Additionally, we identify previously overlooked flaws associated with its application in additive manufacturing. However, we also demonstrate that it encounters limitations when applied to LDED and highlight the need for a more stable explicit scheme.
本文探讨了基于中尺度传导的金属材料激光直接能量沉积(LDED)建模。我们使用两个经过实验验证的模拟,对前向(显式)欧拉时间积分策略和后向欧拉(隐式)方案进行了基准测试。结果表明,显式方案的计算速度更快。此外,我们还发现了该方案在增材制造中应用时以前被忽视的缺陷。不过,我们也证明了它在应用于 LDED 时遇到的限制,并强调了对更稳定的显式方案的需求。
{"title":"Efficient thermal modeling of laser directed energy deposition using the forward Euler scheme: Methodology, merits and limitations","authors":"Simon Essongue ,&nbsp;Vaibhav Nain ,&nbsp;Muriel Carin","doi":"10.1016/j.finel.2024.104270","DOIUrl":"10.1016/j.finel.2024.104270","url":null,"abstract":"<div><div>This paper explores mesoscale conduction-based modeling of Laser Directed Energy Deposition (LDED) for metallic materials. We benchmark the forward Euler (explicit) time integration strategy against the backward Euler (implicit) scheme using two experimentally validated simulations. Our results demonstrate the explicit scheme’s faster computational speed. Additionally, we identify previously overlooked flaws associated with its application in additive manufacturing. However, we also demonstrate that it encounters limitations when applied to LDED and highlight the need for a more stable explicit scheme.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"242 ","pages":"Article 104270"},"PeriodicalIF":3.5,"publicationDate":"2024-10-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142530498","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Optimum thickness design method for micro-shell structure embedded in 3D macrostructure 三维宏观结构中嵌入微壳结构的最佳厚度设计方法
IF 3.5 3区 工程技术 Q1 MATHEMATICS, APPLIED Pub Date : 2024-10-17 DOI: 10.1016/j.finel.2024.104266
Rina Nagai , Masatoshi Shimoda , Musaddiq Al Ali
In this study, we propose a multiscale thickness optimization method for designing micro-shell structure assuming that the macrostructure consists of multiple micro-shell structures. The micro-shell structures are connected to the macrostructure using the NIAH (Novel numerical implementation of asymptotic homogenization) method. The distributed thickness of the micro-shell structures is used as design variable. A squared error norm between actual and target displacements is minimized for controlling the displacements at arbitrary points of the macrostructure to the target values under the total volume constraint including the volume of the micro-shell structures. This design is formulated as a distributed optimization problem, and the thickness gradient function is theoretically derived. The derived sensitivity function is applied to the scalar-type H1 gradient method to efficiently obtain the optimal thickness distribution of the micro-shell structures. Numerical examples demonstrate the effectiveness of the proposed method to optimize the thickness distribution of complex micro-shell structures.
在本研究中,我们假设宏观结构由多个微壳结构组成,提出了一种用于设计微壳结构的多尺度厚度优化方法。使用 NIAH(渐进均质化的新颖数值实现)方法将微壳结构与宏观结构连接起来。微壳结构的分布厚度被用作设计变量。在包括微壳结构体积在内的总体积约束条件下,最小化实际位移与目标位移之间的平方误差规范,以控制宏观结构任意点的位移达到目标值。该设计被表述为分布式优化问题,并从理论上推导出厚度梯度函数。将推导出的灵敏度函数应用于标量型 H1 梯度法,可有效地获得微壳结构的最佳厚度分布。数值实例证明了所提方法在优化复杂微壳结构厚度分布方面的有效性。
{"title":"Optimum thickness design method for micro-shell structure embedded in 3D macrostructure","authors":"Rina Nagai ,&nbsp;Masatoshi Shimoda ,&nbsp;Musaddiq Al Ali","doi":"10.1016/j.finel.2024.104266","DOIUrl":"10.1016/j.finel.2024.104266","url":null,"abstract":"<div><div>In this study, we propose a multiscale thickness optimization method for designing micro-shell structure assuming that the macrostructure consists of multiple micro-shell structures. The micro-shell structures are connected to the macrostructure using the NIAH (Novel numerical implementation of asymptotic homogenization) method. The distributed thickness of the micro-shell structures is used as design variable. A squared error norm between actual and target displacements is minimized for controlling the displacements at arbitrary points of the macrostructure to the target values under the total volume constraint including the volume of the micro-shell structures. This design is formulated as a distributed optimization problem, and the thickness gradient function is theoretically derived. The derived sensitivity function is applied to the scalar-type H<sup>1</sup> gradient method to efficiently obtain the optimal thickness distribution of the micro-shell structures. Numerical examples demonstrate the effectiveness of the proposed method to optimize the thickness distribution of complex micro-shell structures.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"242 ","pages":"Article 104266"},"PeriodicalIF":3.5,"publicationDate":"2024-10-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142445856","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Adaptive stopping criterion of iterative solvers for efficient computational cost reduction: Application to Navier–Stokes with thermal coupling 有效降低计算成本的迭代求解器自适应停止准则:热耦合纳维-斯托克斯应用
IF 3.5 3区 工程技术 Q1 MATHEMATICS, APPLIED Pub Date : 2024-10-15 DOI: 10.1016/j.finel.2024.104263
Ghaniyya Medghoul, Gabriel Manzinali, Elie Hachem, Aurélien Larcher
In this article, a strategy for efficient computational cost reduction of numerical simulations for complex industrial applications is developed and evaluated on multiphysics problems. The approach is based on the adaptive stopping criterion for iterative linear solvers previously implemented for elliptic partial differential equations and the convection–diffusion equation. Control of the convergence of iterative linear solvers is inferred from a posteriori error estimators used for anisotropic mesh adaptation. Provided that the computed error indicator provides an equivalent control on the discretization error, it is a suitable ingredient to assess when enough accuracy has been reached so that iterations of algebraic solvers can be stopped. In practice the iterative solution is stopped when the algebraic error is lower than a percentage of the estimated discretization error. The proposed method proves to be an effective cost-free strategy to reduce the number of iterations needed without degrading the accuracy of the solution. The discretization in the current work is based on stabilized finite elements, while the Generalized Minimal Residual method (GMRES) is used as iterative linear solver. Numerical experiments are performed of increasing complexity, from manufactured solutions to industrial configurations to evaluate the efficiency and the strengths of the proposed adaptive method.
本文开发了一种有效降低复杂工业应用数值模拟计算成本的策略,并对多物理场问题进行了评估。该方法基于之前针对椭圆偏微分方程和对流扩散方程实施的迭代线性求解器的自适应停止准则。对迭代线性求解器收敛性的控制是通过用于各向异性网格适应的后验误差估计来推断的。如果计算出的误差指标能对离散化误差进行等效控制,那么它就可以用来评估何时达到足够的精度,从而停止代数求解器的迭代。实际上,当代数误差低于估计离散化误差的某个百分比时,迭代求解就会停止。事实证明,所提出的方法是一种有效的无成本策略,可以在不降低求解精度的情况下减少所需的迭代次数。当前工作中的离散化基于稳定有限元,而广义最小残差法(GMRES)被用作迭代线性求解器。为了评估所提出的自适应方法的效率和优势,我们进行了复杂度不断增加的数值实验,从制造解决方案到工业配置。
{"title":"Adaptive stopping criterion of iterative solvers for efficient computational cost reduction: Application to Navier–Stokes with thermal coupling","authors":"Ghaniyya Medghoul,&nbsp;Gabriel Manzinali,&nbsp;Elie Hachem,&nbsp;Aurélien Larcher","doi":"10.1016/j.finel.2024.104263","DOIUrl":"10.1016/j.finel.2024.104263","url":null,"abstract":"<div><div>In this article, a strategy for efficient computational cost reduction of numerical simulations for complex industrial applications is developed and evaluated on multiphysics problems. The approach is based on the adaptive stopping criterion for iterative linear solvers previously implemented for elliptic partial differential equations and the convection–diffusion equation. Control of the convergence of iterative linear solvers is inferred from <em>a posteriori</em> error estimators used for anisotropic mesh adaptation. Provided that the computed error indicator provides an equivalent control on the discretization error, it is a suitable ingredient to assess when enough accuracy has been reached so that iterations of algebraic solvers can be stopped. In practice the iterative solution is stopped when the algebraic error is lower than a percentage of the estimated discretization error. The proposed method proves to be an effective cost-free strategy to reduce the number of iterations needed without degrading the accuracy of the solution. The discretization in the current work is based on stabilized finite elements, while the Generalized Minimal Residual method (GMRES) is used as iterative linear solver. Numerical experiments are performed of increasing complexity, from manufactured solutions to industrial configurations to evaluate the efficiency and the strengths of the proposed adaptive method.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"242 ","pages":"Article 104263"},"PeriodicalIF":3.5,"publicationDate":"2024-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142434482","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Multi-objective topological design considering functionally graded materials and coated fiber reinforcement 考虑功能分级材料和涂层纤维加固的多目标拓扑设计
IF 3.5 3区 工程技术 Q1 MATHEMATICS, APPLIED Pub Date : 2024-10-09 DOI: 10.1016/j.finel.2024.104269
Hyunseung Ryu , Jeonghoon Yoo
This study presents a multi-objective topology optimization method tailored to structures fabricated from functionally graded materials (FGMs), coated FGMs, and coated fiber-reinforced composite materials (FRCMs) with fixed fiber thickness. The design objective is the simultaneous minimization of elastic and thermal compliance. The material properties of these composite materials were derived to generate datasets using the representative volume element method under periodic boundary conditions. Subsequently, machine learning modules were developed based on the datasets to combine with the design process. The multi-objective optimization problem was addressed using the weighted sum method ensuring the generation of the Pareto front. The adaptive weighting strategy is employed to avoid biased results toward a single objective function. To define the coated boundaries within the design domain, image post-processing techniques such as convolution filters, interpolation schemes, and erosion methods were employed on the material layout information of the optimized FGM structures. Through numerical examples, optimized material layouts for coated assemblies incorporating FGMs and FRCMs are presented, with the performance verified through objective function values.
本研究提出了一种多目标拓扑优化方法,适用于由功能分级材料(FGM)、涂层 FGM 和具有固定纤维厚度的涂层纤维增强复合材料(FRCM)制造的结构。设计目标是同时最小化弹性和热顺应性。在周期性边界条件下,使用代表性体积元素法推导出这些复合材料的材料特性,并生成数据集。随后,根据数据集开发了机器学习模块,以便与设计过程相结合。多目标优化问题采用加权和法进行处理,以确保生成帕累托前沿。采用了自适应加权策略,以避免结果偏向单一目标函数。为了在设计域内定义涂层边界,在优化的 FGM 结构的材料布局信息上采用了卷积滤波器、插值方案和侵蚀方法等图像后处理技术。通过数值示例,介绍了包含 FGM 和 FRCM 的涂层组件的优化材料布局,并通过目标函数值验证了其性能。
{"title":"Multi-objective topological design considering functionally graded materials and coated fiber reinforcement","authors":"Hyunseung Ryu ,&nbsp;Jeonghoon Yoo","doi":"10.1016/j.finel.2024.104269","DOIUrl":"10.1016/j.finel.2024.104269","url":null,"abstract":"<div><div>This study presents a multi-objective topology optimization method tailored to structures fabricated from functionally graded materials (FGMs), coated FGMs, and coated fiber-reinforced composite materials (FRCMs) with fixed fiber thickness. The design objective is the simultaneous minimization of elastic and thermal compliance. The material properties of these composite materials were derived to generate datasets using the representative volume element method under periodic boundary conditions. Subsequently, machine learning modules were developed based on the datasets to combine with the design process. The multi-objective optimization problem was addressed using the weighted sum method ensuring the generation of the Pareto front. The adaptive weighting strategy is employed to avoid biased results toward a single objective function. To define the coated boundaries within the design domain, image post-processing techniques such as convolution filters, interpolation schemes, and erosion methods were employed on the material layout information of the optimized FGM structures. Through numerical examples, optimized material layouts for coated assemblies incorporating FGMs and FRCMs are presented, with the performance verified through objective function values.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"242 ","pages":"Article 104269"},"PeriodicalIF":3.5,"publicationDate":"2024-10-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142421814","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
The influence of anomalies in supporting structures on the validation of finite-element blade bearing models 支撑结构异常对有限元叶片轴承模型验证的影响
IF 3.5 3区 工程技术 Q1 MATHEMATICS, APPLIED Pub Date : 2024-10-09 DOI: 10.1016/j.finel.2024.104268
Matthis Graßmann, Matthias Stammler, Oliver Menck, Florian Schleich
Finite-element analysis is the only means to determine the load distribution of large slewing bearings considering flexible bearing rings and supporting structures. For reliable results, the plausibility of the models need to be validated. Previous attempts on validating a finite-element model of a slewing bearing against measurement results have indicated a huge dependence of the deformation on tolerances in the supporting structures. This dependence has not yet been explored in research in favor of a focus on tolerances of the bearing itself. The present work explores different irregularities of the flange that connects to the outer ring of the bearing and their effects on bearing deformation. The results show that single dents or bulges on the flange and inclined flanges of the adapter ring significantly change the load distribution and contact angles of the bearing. They also aggravate the risk of truncation. For the calculated fatigue life however, the bearings seem to be robust to these uncertainties for the shown load cases. The dimensions of the investigated tolerances are verified by comparing the resulting deformations of the bearing outer ring against experimental data.
考虑到柔性轴承套圈和支撑结构,有限元分析是确定大型回转支承载荷分布的唯一方法。为了获得可靠的结果,需要对模型的合理性进行验证。之前根据测量结果验证回转支承有限元模型的尝试表明,变形与支承结构的公差有很大关系。由于研究重点是轴承本身的公差,因此尚未对这种依赖性进行探讨。本研究探讨了连接轴承外圈的凸缘的各种不规则情况及其对轴承变形的影响。结果表明,法兰上的单个凹痕或凸起以及适配器外圈的倾斜法兰会显著改变轴承的载荷分布和接触角。它们还增加了截断的风险。然而,就计算的疲劳寿命而言,轴承在所示载荷情况下对这些不确定性似乎是稳健的。通过将轴承外圈的变形结果与实验数据进行比较,验证了所研究的公差尺寸。
{"title":"The influence of anomalies in supporting structures on the validation of finite-element blade bearing models","authors":"Matthis Graßmann,&nbsp;Matthias Stammler,&nbsp;Oliver Menck,&nbsp;Florian Schleich","doi":"10.1016/j.finel.2024.104268","DOIUrl":"10.1016/j.finel.2024.104268","url":null,"abstract":"<div><div>Finite-element analysis is the only means to determine the load distribution of large slewing bearings considering flexible bearing rings and supporting structures. For reliable results, the plausibility of the models need to be validated. Previous attempts on validating a finite-element model of a slewing bearing against measurement results have indicated a huge dependence of the deformation on tolerances in the supporting structures. This dependence has not yet been explored in research in favor of a focus on tolerances of the bearing itself. The present work explores different irregularities of the flange that connects to the outer ring of the bearing and their effects on bearing deformation. The results show that single dents or bulges on the flange and inclined flanges of the adapter ring significantly change the load distribution and contact angles of the bearing. They also aggravate the risk of truncation. For the calculated fatigue life however, the bearings seem to be robust to these uncertainties for the shown load cases. The dimensions of the investigated tolerances are verified by comparing the resulting deformations of the bearing outer ring against experimental data.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"242 ","pages":"Article 104268"},"PeriodicalIF":3.5,"publicationDate":"2024-10-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142421786","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
期刊
Finite Elements in Analysis and Design
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1