首页 > 最新文献

Computer Methods in Applied Mechanics and Engineering最新文献

英文 中文
Reliability-based topology optimization using LRPIM surrogate model considering local stress and displacement constraints 利用考虑局部应力和位移约束的 LRPIM 代用模型进行基于可靠性的拓扑优化
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-16 DOI: 10.1016/j.cma.2024.117460
This paper presents a novel decoupled framework for reliability-based topology optimization (RBTO) that aims to find optimal material configurations while meeting local stiffness and strength constraints. To effectively address the nonlinear displacement and stress reliability constraints, the proposed framework replaces the conventional first-order reliability method (FORM) with the more accurate Local Radial Point Interpolation Method (LRPIM). This substitution overcomes the limitations of FORM in approximating high-dimensional nonlinear problems. The framework includes the qp-relaxation criterion and a global stress aggregation technique to avoid stress singularities. For multi-constrained optimization, the adjoint vector method is used for design sensitivity analysis, followed by a gradient-based algorithm to solve the structural optimization problem. Numerical examples are presented to validate the effectiveness of the proposed RBTO methodology, demonstrating its superiority in both accuracy and reliability compared to the Sequential Optimization and Reliability Assessment (SORA) method. The comparative analysis highlights the efficiency and precision of the proposed method across different reliability approaches, making it a robust tool for addressing complex engineering challenges.
本文提出了一种新颖的基于可靠性的拓扑优化(RBTO)解耦框架,旨在找到最佳材料配置,同时满足局部刚度和强度约束。为有效解决非线性位移和应力可靠性约束,本文提出的框架用更精确的局部径向点插值法(LRPIM)取代了传统的一阶可靠性方法(FORM)。这种替代方法克服了 FORM 在逼近高维非线性问题时的局限性。该框架包括 qp 松弛准则和全局应力聚集技术,以避免应力奇点。在多约束优化方面,使用邻接向量法进行设计敏感性分析,然后使用基于梯度的算法解决结构优化问题。通过数值示例验证了所提出的 RBTO 方法的有效性,证明其与顺序优化和可靠性评估(SORA)方法相比,在准确性和可靠性方面都更胜一筹。对比分析凸显了拟议方法在不同可靠性方法中的效率和精确性,使其成为应对复杂工程挑战的强大工具。
{"title":"Reliability-based topology optimization using LRPIM surrogate model considering local stress and displacement constraints","authors":"","doi":"10.1016/j.cma.2024.117460","DOIUrl":"10.1016/j.cma.2024.117460","url":null,"abstract":"<div><div>This paper presents a novel decoupled framework for reliability-based topology optimization (RBTO) that aims to find optimal material configurations while meeting local stiffness and strength constraints. To effectively address the nonlinear displacement and stress reliability constraints, the proposed framework replaces the conventional first-order reliability method (FORM) with the more accurate Local Radial Point Interpolation Method (LRPIM). This substitution overcomes the limitations of FORM in approximating high-dimensional nonlinear problems. The framework includes the qp-relaxation criterion and a global stress aggregation technique to avoid stress singularities. For multi-constrained optimization, the adjoint vector method is used for design sensitivity analysis, followed by a gradient-based algorithm to solve the structural optimization problem. Numerical examples are presented to validate the effectiveness of the proposed RBTO methodology, demonstrating its superiority in both accuracy and reliability compared to the Sequential Optimization and Reliability Assessment (SORA) method. The comparative analysis highlights the efficiency and precision of the proposed method across different reliability approaches, making it a robust tool for addressing complex engineering challenges.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142441963","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
A novel mesh regularization approach based on finite element distortion potentials: Application to material expansion processes with extreme volume change 基于有限元畸变势的新型网格正则化方法:应用于具有极端体积变化的材料膨胀过程
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-15 DOI: 10.1016/j.cma.2024.117444
The accuracy of finite element solutions is closely tied to the mesh quality. In particular, geometrically nonlinear problems involving large and strongly localized deformations often result in prohibitively large element distortions. In this work, we propose a novel mesh regularization approach allowing to restore a non-distorted high-quality mesh in an adaptive manner without the need for expensive re-meshing procedures. The core idea of this approach lies in the definition of a finite element distortion potential considering contributions from different distortion modes such as skewness and aspect ratio of the elements. The regularized mesh is found by minimization of this potential. Moreover, based on the concept of spatial localization functions, the method allows to specify tailored requirements on mesh resolution and quality for regions with strongly localized mechanical deformation and mesh distortion. In addition, while existing mesh regularization schemes often keep the boundary nodes of the discretization fixed, we propose a mesh-sliding algorithm based on variationally consistent mortar methods allowing for an unrestricted tangential motion of nodes along the problem boundary. Especially for problems involving significant surface deformation (e.g., frictional contact), this approach allows for an improved mesh relaxation as compared to schemes with fixed boundary nodes. To transfer data such as tensor-valued history variables of the material model from the old (distorted) to the new (regularized) mesh, a structure-preserving invariant interpolation scheme for second-order tensors is employed, which has been proposed in our previous work and is designed to preserve important properties of tensor-valued data such as objectivity and positive definiteness. As a practically relevant application scenario, we consider the thermo-mechanical expansion of materials such as foams involving extreme volume changes by up to two orders of magnitude along with large and strongly localized strains as well as thermo-mechanical contact interaction. For this scenario, it is demonstrated that the proposed regularization approach preserves a high mesh quality at small computational costs. In contrast, simulations without mesh adaption are shown to lead to significant mesh distortion, deteriorating result quality, and, eventually, to non-convergence of the numerical solution scheme.
有限元求解的精度与网格质量密切相关。特别是涉及大变形和强局部变形的几何非线性问题,往往会导致令人望而却步的大面积元素失真。在这项工作中,我们提出了一种新颖的网格正则化方法,它能以自适应的方式恢复无扭曲的高质量网格,而无需昂贵的重新网格化程序。这种方法的核心思想在于定义有限元畸变势能,同时考虑不同畸变模式的贡献,如元素的倾斜度和长宽比。正则化网格是通过最小化该变形势来实现的。此外,基于空间局部化函数的概念,该方法可针对具有强烈局部机械变形和网格畸变的区域,对网格分辨率和质量提出量身定制的要求。此外,现有的网格正则化方案通常会固定离散化的边界节点,而我们提出的网格滑动算法基于变化一致的迫击炮方法,允许节点沿问题边界无限制地切向运动。特别是对于涉及重大表面变形(如摩擦接触)的问题,与采用固定边界节点的方案相比,这种方法可以改善网格松弛。为了将材料模型的张量值历史变量等数据从旧的(扭曲的)网格转移到新的(正则化的)网格,我们采用了二阶张量的结构保持不变插值方案,该方案已在我们之前的工作中提出,旨在保持张量值数据的重要特性,如客观性和正确定性。作为一种实际应用场景,我们考虑了泡沫等材料的热机械膨胀,其中涉及高达两个数量级的极端体积变化、大而强的局部应变以及热机械接触相互作用。在这种情况下,所提出的正则化方法以较小的计算成本保持了较高的网格质量。相比之下,不进行网格自适应的模拟则会导致严重的网格畸变、结果质量下降,最终导致数值求解方案不收敛。
{"title":"A novel mesh regularization approach based on finite element distortion potentials: Application to material expansion processes with extreme volume change","authors":"","doi":"10.1016/j.cma.2024.117444","DOIUrl":"10.1016/j.cma.2024.117444","url":null,"abstract":"<div><div>The accuracy of finite element solutions is closely tied to the mesh quality. In particular, geometrically nonlinear problems involving large and strongly localized deformations often result in prohibitively large element distortions. In this work, we propose a novel mesh regularization approach allowing to restore a non-distorted high-quality mesh in an adaptive manner without the need for expensive re-meshing procedures. The core idea of this approach lies in the definition of a finite element distortion potential considering contributions from different distortion modes such as skewness and aspect ratio of the elements. The regularized mesh is found by minimization of this potential. Moreover, based on the concept of spatial localization functions, the method allows to specify tailored requirements on mesh resolution and quality for regions with strongly localized mechanical deformation and mesh distortion. In addition, while existing mesh regularization schemes often keep the boundary nodes of the discretization fixed, we propose a mesh-sliding algorithm based on variationally consistent mortar methods allowing for an unrestricted tangential motion of nodes along the problem boundary. Especially for problems involving significant surface deformation (e.g., frictional contact), this approach allows for an improved mesh relaxation as compared to schemes with fixed boundary nodes. To transfer data such as tensor-valued history variables of the material model from the old (distorted) to the new (regularized) mesh, a structure-preserving invariant interpolation scheme for second-order tensors is employed, which has been proposed in our previous work and is designed to preserve important properties of tensor-valued data such as objectivity and positive definiteness. As a practically relevant application scenario, we consider the thermo-mechanical expansion of materials such as foams involving extreme volume changes by up to two orders of magnitude along with large and strongly localized strains as well as thermo-mechanical contact interaction. For this scenario, it is demonstrated that the proposed regularization approach preserves a high mesh quality at small computational costs. In contrast, simulations without mesh adaption are shown to lead to significant mesh distortion, deteriorating result quality, and, eventually, to non-convergence of the numerical solution scheme.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142441961","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
PTPI-DL-ROMs: Pre-trained physics-informed deep learning-based reduced order models for nonlinear parametrized PDEs PTPI-DL-ROMs:预先训练的基于物理信息深度学习的非线性参数化 PDE 减阶模型
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-15 DOI: 10.1016/j.cma.2024.117404
Among several recently proposed data-driven Reduced Order Models (ROMs), the coupling of Proper Orthogonal Decompositions (POD) and deep learning-based ROMs (DL-ROMs) has proved to be a successful strategy to construct non-intrusive, highly accurate, surrogates for the real time solution of parametric nonlinear time-dependent PDEs. Inexpensive to evaluate, POD-DL-ROMs are also relatively fast to train, thanks to their limited complexity. However, POD-DL-ROMs account for the physical laws governing the problem at hand only through the training data, that are usually obtained through a full order model (FOM) relying on a high-fidelity discretization of the underlying equations. Moreover, the accuracy of POD-DL-ROMs strongly depends on the amount of available data. In this paper, we consider a major extension of POD-DL-ROMs by enforcing the fulfillment of the governing physical laws in the training process – that is, by making them physics-informed – to compensate for possible scarce and/or unavailable data and improve the overall reliability. To do that, we first complement POD-DL-ROMs with a trunk net architecture, endowing them with the ability to compute the problem’s solution at every point in the spatial domain, and ultimately enabling a seamless computation of the physics-based loss by means of the strong continuous formulation. Then, we introduce an efficient training strategy that limits the notorious computational burden entailed by a physics-informed training phase. In particular, we take advantage of the few available data to develop a low-cost pre-training procedure; then, we fine-tune the architecture in order to further improve the prediction reliability. Accuracy and efficiency of the resulting pre-trained physics-informed DL-ROMs (PTPI-DL-ROMs) are then assessed on a set of test cases ranging from non-affinely parametrized advection–diffusion–reaction equations, to nonlinear problems like the Navier–Stokes equations for fluid flows.
在最近提出的几种数据驱动的还原阶模型(ROM)中,适当正交分解(POD)与基于深度学习的还原阶模型(DL-ROM)的耦合已被证明是一种成功的策略,可以构建非侵入、高精度的代用模型,用于实时求解参数非线性时变 PDE。POD-DL-ROM 的评估成本低,由于其复杂性有限,因此训练速度也相对较快。不过,POD-DL-ROM 只能通过训练数据来解释当前问题的物理规律,而这些数据通常是通过全阶模型(FOM)获得的,依赖于对基础方程的高保真离散化。此外,POD-DL-ROM 的准确性在很大程度上取决于可用数据的数量。在本文中,我们将考虑对 POD-DL-ROMs 进行重大扩展,在训练过程中强制实现管理物理定律,即使其具有物理信息,以弥补可能的数据稀缺和/或不可用数据,并提高整体可靠性。为此,我们首先利用主干网架构对 POD-DL-ROM 进行补充,使其具备在空间域的每个点计算问题解决方案的能力,最终通过强连续公式实现基于物理损失的无缝计算。然后,我们引入了一种高效的训练策略,以限制物理信息训练阶段带来的众所周知的计算负担。特别是,我们利用为数不多的可用数据,开发了一种低成本的预训练程序;然后,我们对架构进行了微调,以进一步提高预测的可靠性。然后,我们在一组测试案例中评估了预训练的物理信息 DL-ROM (PTPI-DL-ROM)的准确性和效率,这些案例包括非参数化的平流-扩散-反应方程,以及流体流动的纳维-斯托克斯方程等非线性问题。
{"title":"PTPI-DL-ROMs: Pre-trained physics-informed deep learning-based reduced order models for nonlinear parametrized PDEs","authors":"","doi":"10.1016/j.cma.2024.117404","DOIUrl":"10.1016/j.cma.2024.117404","url":null,"abstract":"<div><div>Among several recently proposed data-driven Reduced Order Models (ROMs), the coupling of Proper Orthogonal Decompositions (POD) and deep learning-based ROMs (DL-ROMs) has proved to be a successful strategy to construct non-intrusive, highly accurate, surrogates for the <em>real time</em> solution of parametric nonlinear time-dependent PDEs. Inexpensive to evaluate, POD-DL-ROMs are also relatively fast to train, thanks to their limited complexity. However, POD-DL-ROMs account for the physical laws governing the problem at hand only through the training data, that are usually obtained through a full order model (FOM) relying on a high-fidelity discretization of the underlying equations. Moreover, the accuracy of POD-DL-ROMs strongly depends on the amount of available data. In this paper, we consider a major extension of POD-DL-ROMs by enforcing the fulfillment of the governing physical laws in the training process – that is, by making them physics-informed – to compensate for possible scarce and/or unavailable data and improve the overall reliability. To do that, we first complement POD-DL-ROMs with a <em>trunk net</em> architecture, endowing them with the ability to compute the problem’s solution at every point in the spatial domain, and ultimately enabling a seamless computation of the physics-based loss by means of the strong continuous formulation. Then, we introduce an efficient training strategy that limits the notorious computational burden entailed by a physics-informed training phase. In particular, we take advantage of the few available data to develop a low-cost pre-training procedure; then, we fine-tune the architecture in order to further improve the prediction reliability. Accuracy and efficiency of the resulting pre-trained physics-informed DL-ROMs (PTPI-DL-ROMs) are then assessed on a set of test cases ranging from non-affinely parametrized advection–diffusion–reaction equations, to nonlinear problems like the Navier–Stokes equations for fluid flows.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142433311","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
High-order multiscale method for elastic deformation of complex geometries 复杂几何体弹性变形的高阶多尺度方法
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-14 DOI: 10.1016/j.cma.2024.117436
Computational methods, such as finite elements, are indispensable for modeling the mechanical compliance of elastic solids. However, as the size and geometric complexity of the domain increases, the cost of simulations becomes prohibitive. One example is the microstructure of a porous material, such as a piece of rock or bone sample, captured by an X-ray μCT image. The solid geometry consists of numerous grains, cavities, and/or channels, with the domain large enough to allow inferring statistically converged macroscale properties. The pore-level multiscale method (PLMM) was recently proposed by the authors to reduce the associated computational cost through a divide-and-conquer strategy. The domain is decomposed into subdomains via watershed segmentation, and local basis and correction functions are built numerically, then assembled to obtain an approximate solution. However, PLMM is limited to domains corresponding to microscale porous media, incurs large errors when modeling loading conditions that generate significant bending/twisting moments locally, and it is equipped with only one mechanism to control approximation errors during a simulation. Here, we generalize PLMM into a high-order variant, called hPLMM, that removes these drawbacks. In hPLMM, appropriate mortar spaces are defined at subdomain interfaces that allow improving the boundary conditions used to solve local problems on the subdomains, thus the accuracy of the approximation. Moreover, errors can be reduced by a second mechanism wherein an upfront cost is paid prior to a simulation, useful if basis functions can be reused many times, e.g., across loading steps. Finally, the method applies not just to pore-scale, but also Darcy-scale and non-porous domains. We validate hPLMM against a range of complex 2D/3D geometries and discuss its convergence and algorithmic complexity. Implications for modeling failure and nonlinear problems are discussed.
有限元等计算方法是模拟弹性固体机械顺应性不可或缺的方法。然而,随着领域大小和几何复杂性的增加,模拟成本也变得高昂。一个例子是 X 射线 μCT 图像捕捉到的多孔材料(如岩石或骨骼样本)的微观结构。固体几何结构由许多晶粒、空腔和/或通道组成,其域大到足以推断出统计收敛的宏观尺度属性。作者最近提出了孔隙级多尺度方法(PLMM),通过分而治之的策略降低相关计算成本。该方法通过分水岭分割将域分解为子域,并通过数值方法建立局部基函数和校正函数,然后进行组合以获得近似解。然而,PLMM 仅限于微尺度多孔介质对应的域,在模拟局部产生较大弯曲/扭转力矩的加载条件时会产生较大误差,而且在模拟过程中只有一种控制近似误差的机制。在此,我们将 PLMM 推广为高阶变体,称为 hPLMM,以消除这些缺点。在 hPLMM 中,子域界面上定义了适当的迫击炮空间,可以改善用于解决子域局部问题的边界条件,从而提高近似的精度。此外,还可以通过第二种机制来减少误差,即在模拟前支付预付费用,这在基函数可以多次重复使用(如跨加载步骤)的情况下非常有用。最后,该方法不仅适用于孔隙尺度,也适用于达西尺度和非孔隙域。我们针对一系列复杂的二维/三维几何图形验证了 hPLMM,并讨论了其收敛性和算法复杂性。我们还讨论了失效和非线性问题建模的意义。
{"title":"High-order multiscale method for elastic deformation of complex geometries","authors":"","doi":"10.1016/j.cma.2024.117436","DOIUrl":"10.1016/j.cma.2024.117436","url":null,"abstract":"<div><div>Computational methods, such as finite elements, are indispensable for modeling the mechanical compliance of elastic solids. However, as the size and geometric complexity of the domain increases, the cost of simulations becomes prohibitive. One example is the microstructure of a porous material, such as a piece of rock or bone sample, captured by an X-ray <span><math><mi>μ</mi></math></span>CT image. The solid geometry consists of numerous grains, cavities, and/or channels, with the domain large enough to allow inferring statistically converged macroscale properties. The pore-level multiscale method (PLMM) was recently proposed by the authors to reduce the associated computational cost through a divide-and-conquer strategy. The domain is decomposed into subdomains via watershed segmentation, and local basis and correction functions are built numerically, then assembled to obtain an approximate solution. However, PLMM is limited to domains corresponding to microscale porous media, incurs large errors when modeling loading conditions that generate significant bending/twisting moments locally, and it is equipped with only one mechanism to control approximation errors <em>during</em> a simulation. Here, we generalize PLMM into a high-order variant, called hPLMM, that removes these drawbacks. In hPLMM, appropriate mortar spaces are defined at subdomain interfaces that allow improving the boundary conditions used to solve local problems on the subdomains, thus the accuracy of the approximation. Moreover, errors can be reduced by a second mechanism wherein an upfront cost is paid <em>prior to</em> a simulation, useful if basis functions can be reused many times, e.g., across loading steps. Finally, the method applies not just to pore-scale, but also Darcy-scale and non-porous domains. We validate hPLMM against a range of complex 2D/3D geometries and discuss its convergence and algorithmic complexity. Implications for modeling failure and nonlinear problems are discussed.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-10-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142433309","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Machine-learning-enabled discrete element method: The extension to three dimensions and computational issues 机器学习离散元素法:向三维扩展和计算问题
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-14 DOI: 10.1016/j.cma.2024.117445
The detection and resolution of contacts among irregular-shaped particles pose significant challenges in the discrete element method (DEM) and recent advancements have introduced a machine learning-enabled approach specifically tailored for contact detection and resolution in two dimensions. Building upon this progress, this paper extends the application of machine learning-enabled DEM to encompass the more complex and realistic three-dimensional (3D) scenario. Particles are modeled using a polyhedral representation with arbitrary shapes, and contact behavior is governed by an energy-conserving contact model based on contact volumes. The efficacy of the machine learning-enabled 3D DEM is evaluated through comparative analyses of computational time and simulation results across individual contact as well as whole DEM simulations against those obtained from the conventional DEM. The findings indicate that the machine learning-enabled approach adeptly identifies and resolves contacts among 3D irregular-shaped particles while accurately reproducing the mechanical characteristics of densely contacting particle assemblies. The computational issues and challenges associated with the machine learning-enabled DEM are also discussed. The study highlights that the machine learning-enabled approach significantly enhances computational efficiency, showcasing its potential to advance complex DEM simulations in a more efficient manner.
在离散元素法(DEM)中,不规则形状颗粒间接触的检测和解析是一项重大挑战,而最近的进展则引入了一种专门针对二维接触检测和解析的机器学习方法。在这一进展的基础上,本文扩展了机器学习离散元素法的应用范围,以涵盖更复杂、更现实的三维(3D)场景。粒子采用任意形状的多面体表示法建模,接触行为由基于接触体积的能量守恒接触模型控制。通过对单个接触和整个 DEM 仿真的计算时间和仿真结果与传统 DEM 仿真结果的比较分析,评估了机器学习 3D DEM 的功效。研究结果表明,支持机器学习的方法能够巧妙地识别和解决三维不规则形状颗粒之间的接触,同时准确地再现了密集接触颗粒组件的机械特性。研究还讨论了与机器学习 DEM 相关的计算问题和挑战。研究强调,机器学习支持方法显著提高了计算效率,展示了以更高效的方式推进复杂 DEM 仿真的潜力。
{"title":"Machine-learning-enabled discrete element method: The extension to three dimensions and computational issues","authors":"","doi":"10.1016/j.cma.2024.117445","DOIUrl":"10.1016/j.cma.2024.117445","url":null,"abstract":"<div><div>The detection and resolution of contacts among irregular-shaped particles pose significant challenges in the discrete element method (DEM) and recent advancements have introduced a machine learning-enabled approach specifically tailored for contact detection and resolution in two dimensions. Building upon this progress, this paper extends the application of machine learning-enabled DEM to encompass the more complex and realistic three-dimensional (3D) scenario. Particles are modeled using a polyhedral representation with arbitrary shapes, and contact behavior is governed by an energy-conserving contact model based on contact volumes. The efficacy of the machine learning-enabled 3D DEM is evaluated through comparative analyses of computational time and simulation results across individual contact as well as whole DEM simulations against those obtained from the conventional DEM. The findings indicate that the machine learning-enabled approach adeptly identifies and resolves contacts among 3D irregular-shaped particles while accurately reproducing the mechanical characteristics of densely contacting particle assemblies. The computational issues and challenges associated with the machine learning-enabled DEM are also discussed. The study highlights that the machine learning-enabled approach significantly enhances computational efficiency, showcasing its potential to advance complex DEM simulations in a more efficient manner.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-10-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142433310","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Maximum a posteriori estimation for linear structural dynamics models using Bayesian optimization with rational polynomial chaos expansions 利用贝叶斯优化与有理多项式混沌扩展对线性结构动力学模型进行最大后验估计
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-12 DOI: 10.1016/j.cma.2024.117418
Bayesian analysis enables combining prior knowledge with measurement data to learn model parameters. Commonly, one resorts to computing the maximum a posteriori (MAP) estimate, when only a point estimate of the parameters is of interest. We apply MAP estimation in the context of structural dynamic models, where the system response can be described by the frequency response function. To alleviate high computational demands from repeated expensive model calls, we utilize a rational polynomial chaos expansion (RPCE) surrogate model that expresses the system frequency response as a rational of two polynomials with complex coefficients. We propose an extension to an existing sparse Bayesian learning approach for RPCE based on Laplace’s approximation for the posterior distribution of the denominator coefficients. Furthermore, we introduce a Bayesian optimization approach, which allows to adaptively enrich the experimental design throughout the optimization process of MAP estimation. Thereby, we utilize the expected improvement acquisition function as a means to identify sample points in the input space that are possibly associated with large objective function values. The acquisition function is estimated through Monte Carlo sampling based on the posterior distribution of the expansion coefficients identified in the sparse Bayesian learning process. By combining the sparsity-inducing learning procedure with the sequential experimental design, we effectively reduce the number of model evaluations in the MAP estimation problem. We demonstrate the applicability of the presented methods on the parameter updating problem of an algebraic two-degree-of-freedom system and the finite element model of a cross-laminated timber plate.
贝叶斯分析法能够将先验知识与测量数据相结合,从而学习模型参数。通常情况下,当我们只对参数的点估计感兴趣时,就会求助于计算最大后验(MAP)估计值。我们将 MAP 估计应用于结构动态模型,其中系统响应可以用频率响应函数来描述。为了减轻重复调用昂贵模型所带来的高计算需求,我们利用有理多项式混沌扩展(RPCE)替代模型,将系统频率响应表示为两个具有复系数的多项式的有理数。我们根据分母系数后验分布的拉普拉斯近似值,为 RPCE 提出了对现有稀疏贝叶斯学习方法的扩展。此外,我们还引入了一种贝叶斯优化方法,可以在 MAP 估计的整个优化过程中自适应地丰富实验设计。因此,我们利用预期改进获取函数来识别输入空间中可能与较大目标函数值相关的样本点。获取函数是根据稀疏贝叶斯学习过程中确定的扩展系数的后验分布,通过蒙特卡罗抽样进行估计的。通过将稀疏贝叶斯学习过程与顺序实验设计相结合,我们有效地减少了 MAP 估计问题中的模型评估次数。我们在代数二自由度系统的参数更新问题和交叉层压木板的有限元模型上演示了所提出方法的适用性。
{"title":"Maximum a posteriori estimation for linear structural dynamics models using Bayesian optimization with rational polynomial chaos expansions","authors":"","doi":"10.1016/j.cma.2024.117418","DOIUrl":"10.1016/j.cma.2024.117418","url":null,"abstract":"<div><div>Bayesian analysis enables combining prior knowledge with measurement data to learn model parameters. Commonly, one resorts to computing the maximum a posteriori (MAP) estimate, when only a point estimate of the parameters is of interest. We apply MAP estimation in the context of structural dynamic models, where the system response can be described by the frequency response function. To alleviate high computational demands from repeated expensive model calls, we utilize a rational polynomial chaos expansion (RPCE) surrogate model that expresses the system frequency response as a rational of two polynomials with complex coefficients. We propose an extension to an existing sparse Bayesian learning approach for RPCE based on Laplace’s approximation for the posterior distribution of the denominator coefficients. Furthermore, we introduce a Bayesian optimization approach, which allows to adaptively enrich the experimental design throughout the optimization process of MAP estimation. Thereby, we utilize the expected improvement acquisition function as a means to identify sample points in the input space that are possibly associated with large objective function values. The acquisition function is estimated through Monte Carlo sampling based on the posterior distribution of the expansion coefficients identified in the sparse Bayesian learning process. By combining the sparsity-inducing learning procedure with the sequential experimental design, we effectively reduce the number of model evaluations in the MAP estimation problem. We demonstrate the applicability of the presented methods on the parameter updating problem of an algebraic two-degree-of-freedom system and the finite element model of a cross-laminated timber plate.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-10-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142419867","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Optimal control of constrained mechanical systems in redundant coordinates: Formulation and structure-preserving discretization 冗余坐标中受约束机械系统的优化控制:公式化和保结构离散化
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-11 DOI: 10.1016/j.cma.2024.117443
This work deals with optimal control problems for constrained mechanical systems whose motion is governed by differential algebraic equations (DAEs). Both index-3 DAEs and stabilized index-2 DAEs are considered. Two alternative formulations of the optimal control problem are compared to each other. It is shown that symmetries of the optimal control problem lead to the conservation of generalized momentum maps. These generalized momentum maps are related to quadratic invariants of the optimal control problem. A direct discretization approach is newly proposed which is (i) capable to conserve the quadratic invariants, and (ii) equivalent to the indirect approach to the optimal control problem. Numerical examples are presented to access the properties of the newly developed schemes.
本研究涉及受约束机械系统的最优控制问题,该系统的运动受微分代数方程(DAE)控制。既考虑了指数-3 DAE,也考虑了稳定的指数-2 DAE。比较了最优控制问题的两种不同表述。结果表明,最优控制问题的对称性会导致广义动量图的守恒。这些广义动量图与最优控制问题的二次不变量相关。新提出的直接离散化方法(i) 能够保持二次方不变量,(ii) 等同于最优控制问题的间接方法。通过数值示例可以了解新开发方案的特性。
{"title":"Optimal control of constrained mechanical systems in redundant coordinates: Formulation and structure-preserving discretization","authors":"","doi":"10.1016/j.cma.2024.117443","DOIUrl":"10.1016/j.cma.2024.117443","url":null,"abstract":"<div><div>This work deals with optimal control problems for constrained mechanical systems whose motion is governed by differential algebraic equations (DAEs). Both index-3 DAEs and stabilized index-2 DAEs are considered. Two alternative formulations of the optimal control problem are compared to each other. It is shown that symmetries of the optimal control problem lead to the conservation of generalized momentum maps. These generalized momentum maps are related to quadratic invariants of the optimal control problem. A direct discretization approach is newly proposed which is (i) capable to conserve the quadratic invariants, and (ii) equivalent to the indirect approach to the optimal control problem. Numerical examples are presented to access the properties of the newly developed schemes.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-10-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142419877","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Downwind and upwind approximations for primal and dual problems of elasto-plasticity with Prandtl–Reuss type material laws 具有普朗特-罗伊斯型材料定律的弹塑性原始问题和对偶问题的下风和上风近似法
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-10 DOI: 10.1016/j.cma.2024.117277
Modern adaptive finite element (FE) algorithms for solution of initial-boundary-value-problems (IBVP) employ goal-oriented error measures in order to assess the quality of computational results for the physical event under investigation. However, traditional time-stepping algorithms for solution of the corresponding dual problem run backwards-in-time, which due to additional storage requirements might become a serious drawback when an extensive number of time steps for the FEM simulation arises. In this paper, we take advantage of an end-boundary-value-problem (EBVP) associated to IBVP, with corresponding dual-problem running forwards in time. In order to obtain a unified framework for numerical approximation of primal and dual weak forms for both, IBVP and EBVP, respectively, we apply the concept of downwind and upwind approximations not only to the trial functions but also to the test functions. This results into eight different integration schemes. On this basis, as a main result of this contribution, a time-stepping algorithm is obtained, which runs forwards-in-time for the dual problem and therefore avoids the additional storage requirements of the traditional backwards-in-time stepping procedures. The presented algorithm is numerically tested and validated for a CT-specimen for elastic and elasto-plastic behavior, where the constitutive equations are written in Prandtl–Reuss type format.
用于求解初始边界值问题(IBVP)的现代自适应有限元(FE)算法采用了目标导向误差测量方法,以评估所研究物理事件的计算结果质量。然而,用于求解相应对偶问题的传统时间步进算法在时间上向后运行,当有限元模拟出现大量时间步进时,由于额外的存储要求,这可能成为一个严重的缺点。在本文中,我们利用了与 IBVP 相关的终界值问题(EBVP),相应的对偶问题在时间上向前运行。为了分别获得 IBVP 和 EBVP 原始和对偶弱形式数值近似的统一框架,我们不仅对试验函数,而且对测试函数应用了顺风和逆风近似的概念。这就产生了八种不同的积分方案。在此基础上,作为本论文的主要成果,我们获得了一种时间步进算法,该算法在对偶问题上向前-时间运行,因此避免了传统的向后-时间步进程序的额外存储要求。所提出的算法经过数值测试和验证,适用于弹性和弹塑性行为的 CT 试样,其中的构成方程是以普朗特-罗伊斯类型格式编写的。
{"title":"Downwind and upwind approximations for primal and dual problems of elasto-plasticity with Prandtl–Reuss type material laws","authors":"","doi":"10.1016/j.cma.2024.117277","DOIUrl":"10.1016/j.cma.2024.117277","url":null,"abstract":"<div><div>Modern adaptive finite element (FE) algorithms for solution of initial-boundary-value-problems (IBVP) employ goal-oriented error measures in order to assess the quality of computational results for the physical event under investigation. However, traditional time-stepping algorithms for solution of the corresponding dual problem run backwards-in-time, which due to additional storage requirements might become a serious drawback when an extensive number of time steps for the FEM simulation arises. In this paper, we take advantage of an end-boundary-value-problem (EBVP) associated to IBVP, with corresponding dual-problem running forwards in time. In order to obtain a unified framework for numerical approximation of primal and dual weak forms for both, IBVP and EBVP, respectively, we apply the concept of downwind and upwind approximations not only to the trial functions but also to the test functions. This results into eight different integration schemes. On this basis, as a main result of this contribution, a time-stepping algorithm is obtained, which runs forwards-in-time for the dual problem and therefore avoids the additional storage requirements of the traditional backwards-in-time stepping procedures. The presented algorithm is numerically tested and validated for a CT-specimen for elastic and elasto-plastic behavior, where the constitutive equations are written in Prandtl–Reuss type format.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-10-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142419879","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Dynamic topology optimization for structures exhibiting frequency-dependent material properties with prescribed frequency forbidden band 针对具有规定频率禁带的频率相关材料特性结构进行动态拓扑优化
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-09 DOI: 10.1016/j.cma.2024.117439
In dynamic vibration reduction design, the frequency-dependent material properties are crucial for the optimal configuration, especially in the problem of prescribed frequency forbidden band. In this paper, a new dynamic topology optimization method for structures with frequency-dependent material properties is proposed to achieve the vibration reduction design in the prescribed frequency forbidden band. First, a dynamic topology optimization model is established for the problem studied in this paper. This model integrates the solution method for frequency-dependent problem, dynamic isolated structures elimination method and the formulation of prescribed frequency forbidden band constraints, which are based on the research results previously developed by the authors. Additionally, different interpolation schemes are used for different number of material designs. The above optimization model is intended to consider nonlinear terms and design several frequency-dependent structures with prescribed frequency forbidden bands that are more in line with practical engineering problems, so that they can accurately avoid the operating frequency range, thus improving the service life of engineering equipment. Finally, to address common numerical problems, the "bound formulation" and "robust formulation" are employed, enhancing the applicability and robustness of the method for the application in topology optimization. The effectiveness of the developed method is supported by two types optimization problems, including single-material and bi-material examples. The cross-check results reveal that when considering frequency-dependent terms, the design results are better and closer to the practical engineering problem compared to linear structures.
在动态减振设计中,与频率相关的材料特性对优化配置至关重要,尤其是在规定频率禁带问题上。本文针对具有频率相关材料特性的结构,提出了一种新的动态拓扑优化方法,以实现规定频率禁带内的减振设计。首先,针对本文研究的问题建立了一个动态拓扑优化模型。该模型综合了频率相关问题的求解方法、动态孤立结构消除方法和规定频率禁带约束条件的表述,这些都是基于作者之前的研究成果。此外,还针对不同数量的材料设计采用了不同的插值方案。上述优化模型旨在考虑非线性项,设计出几种具有规定频率禁带、更符合实际工程问题的频率相关结构,使其能够准确避开工作频率范围,从而提高工程设备的使用寿命。最后,针对常见的数值问题,采用了 "约束公式 "和 "鲁棒公式",增强了该方法在拓扑优化应用中的适用性和鲁棒性。通过两类优化问题,包括单材料和双材料实例,证明了所开发方法的有效性。交叉检验结果表明,当考虑频率相关项时,设计结果比线性结构更好,更接近实际工程问题。
{"title":"Dynamic topology optimization for structures exhibiting frequency-dependent material properties with prescribed frequency forbidden band","authors":"","doi":"10.1016/j.cma.2024.117439","DOIUrl":"10.1016/j.cma.2024.117439","url":null,"abstract":"<div><div>In dynamic vibration reduction design, the frequency-dependent material properties are crucial for the optimal configuration, especially in the problem of prescribed frequency forbidden band. In this paper, a new dynamic topology optimization method for structures with frequency-dependent material properties is proposed to achieve the vibration reduction design in the prescribed frequency forbidden band. First, a dynamic topology optimization model is established for the problem studied in this paper. This model integrates the solution method for frequency-dependent problem, dynamic isolated structures elimination method and the formulation of prescribed frequency forbidden band constraints, which are based on the research results previously developed by the authors. Additionally, different interpolation schemes are used for different number of material designs. The above optimization model is intended to consider nonlinear terms and design several frequency-dependent structures with prescribed frequency forbidden bands that are more in line with practical engineering problems, so that they can accurately avoid the operating frequency range, thus improving the service life of engineering equipment. Finally, to address common numerical problems, the \"bound formulation\" and \"robust formulation\" are employed, enhancing the applicability and robustness of the method for the application in topology optimization. The effectiveness of the developed method is supported by two types optimization problems, including single-material and bi-material examples. The cross-check results reveal that when considering frequency-dependent terms, the design results are better and closer to the practical engineering problem compared to linear structures.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-10-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142419730","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
A hybrid upwind scheme for two-phase flow in fractured porous media 裂隙多孔介质中两相流的混合上风方案
IF 6.9 1区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2024-10-09 DOI: 10.1016/j.cma.2024.117437
Simulating the flow of two fluid phases in porous media is a challenging task, especially when fractures are included in the simulation. Fractures may have highly heterogeneous properties compared to the surrounding rock matrix, significantly affecting fluid flow, and at the same time hydraulic apertures that are much smaller than any other characteristic sizes in the domain. Generally, flow simulators face difficulties with counter-current flow, generated by gravity and pressure gradients, which hinders the convergence of non-linear solvers (Newton).
In this work, we model the fracture geometry with a mixed-dimensional discrete fracture network, thus lightening the computational burden associated to an equi-dimensional representation. We address the issue of counter-current flows with appropriate spatial discretization of the advective fluid fluxes, with the aim of improving the convergence speed of the non-linear solver. In particular, the extension of the hybrid upwinding to the mixed-dimensional framework, with the use of a phase potential upstreaming at the interfaces of subdomains.
We test the method across several cases with different flow regimes and fracture network geometries. Results show robustness of the chosen discretization and a consistent improvements, in terms of Newton iterations, compared to using phase potential upstreaming everywhere.
模拟多孔介质中两相流体的流动是一项具有挑战性的任务,尤其是当模拟中包含裂缝时。与周围的岩石基质相比,裂缝可能具有高度异质性,从而对流体流动产生重大影响,同时,裂缝的水力孔径也远小于域中其他特征尺寸。一般来说,流动模拟器在处理由重力和压力梯度产生的逆流时会遇到困难,这阻碍了非线性求解器(牛顿)的收敛。在这项工作中,我们用混合维离散断裂网络来模拟断裂几何形状,从而减轻了等维表示的计算负担。我们通过对平流流体通量进行适当的空间离散化来解决逆流问题,目的是提高非线性求解器的收敛速度。我们在不同流态和断裂网络几何形状的多个案例中测试了该方法。结果表明,所选离散方法具有稳健性,与在所有地方都使用相位逆流法相比,牛顿迭代的改进效果一致。
{"title":"A hybrid upwind scheme for two-phase flow in fractured porous media","authors":"","doi":"10.1016/j.cma.2024.117437","DOIUrl":"10.1016/j.cma.2024.117437","url":null,"abstract":"<div><div>Simulating the flow of two fluid phases in porous media is a challenging task, especially when fractures are included in the simulation. Fractures may have highly heterogeneous properties compared to the surrounding rock matrix, significantly affecting fluid flow, and at the same time hydraulic apertures that are much smaller than any other characteristic sizes in the domain. Generally, flow simulators face difficulties with counter-current flow, generated by gravity and pressure gradients, which hinders the convergence of non-linear solvers (Newton).</div><div>In this work, we model the fracture geometry with a mixed-dimensional discrete fracture network, thus lightening the computational burden associated to an equi-dimensional representation. We address the issue of counter-current flows with appropriate spatial discretization of the advective fluid fluxes, with the aim of improving the convergence speed of the non-linear solver. In particular, the extension of the hybrid upwinding to the mixed-dimensional framework, with the use of a phase potential upstreaming at the interfaces of subdomains.</div><div>We test the method across several cases with different flow regimes and fracture network geometries. Results show robustness of the chosen discretization and a consistent improvements, in terms of Newton iterations, compared to using phase potential upstreaming everywhere.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-10-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142419664","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
期刊
Computer Methods in Applied Mechanics and Engineering
全部 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