Concrete structures may suffer considerable restraint stresses during their hardening period. This is caused by several deformation impacts, especially temperature field changings due to hydration heat and volume chan...Concrete structures may suffer considerable restraint stresses during their hardening period. This is caused by several deformation impacts, especially temperature field changings due to hydration heat and volume changes due to autogenous shrinkage. Mainly affected are massive concrete members, but also the application of new concrete types or the erection of outstanding constructions requires further investigations in this context. 3D-FEM analyses of hydration heat induced temperature development in combination with the well known shrinkage give sufficient results for the deformation impact. The according elastic restraint stresses can be determined with consideration of the concrete’s rising elastic modulus and the restraint degree of the system. But due to duration of the heat flow process, the height of restraint stresses is strongly dependent from the viscoelasticity of the concrete. The viscoelastic effects consist of many components constituted by changing material properties influencing themselves. In practice, different simplified approaches are available for considering this in calculations. Their implementation in time step analyses is not generally admitted and requires expertise. In contrast, present research develops material models needing specific input parameters for every use case. This contribution focuses on a practicable approach considering the superposition of the viscoelastic behaviour of every stress increment in time step FEM analysis. The differentiation between the pure viscoelastic material behaviour (as it is given in the codes for idealistic conditions like creep or relaxation) and the according viscoelastic system response (addicted to the systems variable restraint degree) allows the transfer of this model into practice. As one application example of this approach, the compatibility check and the FEM-based recalculation of the monitoring program of a massive power plant slab will be presented.展开更多
We’ll study the FEM for a model for compressible miscible displacement in porous media which includes molecular diffusion and mechanical dispersion in one-dimensional space.A class of vertices-edges-elements interpol...We’ll study the FEM for a model for compressible miscible displacement in porous media which includes molecular diffusion and mechanical dispersion in one-dimensional space.A class of vertices-edges-elements interpolation operator ink is introduced.With the help of ink(not elliptic projection),the optimal error estimate in L∞(J;L2(Ω)) norm of FEM is proved.展开更多
The accurate mechanical analysis of thick-walled pressure vessel structures composed of advanced materials,such as hyperelastic and functionally graded materials(FGMs),is critical for ensuring their safety and optimiz...The accurate mechanical analysis of thick-walled pressure vessel structures composed of advanced materials,such as hyperelastic and functionally graded materials(FGMs),is critical for ensuring their safety and optimizing their design.However,conventional numerical methods can face challenges with the non-linearities inherent in hyperelasticity and the complex spatial variations in FGMs.This paper presents a novel hybrid numerical approach combining Physics-Informed Neural Networks(PINNs)with Finite Element Method(FEM)derived data for the robust analysis of thick-walled,axisymmetric,heterogeneous,hyperelastic pressure vessels with elliptical geometries.A PINN framework incorporating neo-Hookean constitutive relations is developed in MATLAB.To enhance training efficiency and accuracy,the PINN’s loss function is augmented with displacement data obtained from high-fidelity FEM simulations performed in ANSYS.The methodology is rigorously validated by comparing PINN-predicted displacement and von Mises stress fields against ANSYS benchmarks for various scenarios of FGMconfigurations(with material properties varying according to a power law)subjected to internal and external pressurization.The results demonstrate excellent agreement between the proposed hybrid PINN-FEMapproach and conventional FEMsolutions across all test cases,accurately capturing complex deformation patterns and stress concentrations.This study highlights the potential of data-augmented PINNs as an effective and accurate computational tool for tackling complex solid mechanics problems involving non-linearmaterials and significant heterogeneity,offering a promising avenue for future research in engineering design and analysis.展开更多
The dynamic evolution of fracture permeability presents a critical scientific challenge in rock masses.Understanding the mechanisms of rock mass permeability evolution is vital for engineering project design and opera...The dynamic evolution of fracture permeability presents a critical scientific challenge in rock masses.Understanding the mechanisms of rock mass permeability evolution is vital for engineering project design and operations.By integrating the discrete element method(DEM)with the finite element method(FEM),a numerical simulation framework for shear seepage in rough fractured shale has been developed to investigate the dynamic mechanisms of permeability evolution under varying confining pressures and during the shearing process.Numerical simulations were conducted on rough fractured samples under effective confining pressures ranging from 5 MPa to 20 MPa to monitor the aperture and permeability evolution of the fracture.The results of the numerical simulation are consistent with the experimental observations,indicating that both the shearing process and confining pressure significantly influence permeability.Moreover,the magnitude of the confining pressure is a crucial factor influencing the trend in permeability changes.Under a confining pressure of 5 MPa,fracture permeability initially increases significantly but decreases post-shearing.In contrast,a continuous decrease in fracture permeability is observed when the confining pressure exceeds 10 MPa.The results of the shear numerical simulation indicate that the confining pressure restricts fracture dilation during shearing,promotes the generation of rock debris,and decreases both the permeability and transmissivity of the fracture.The wear results obtained from numerical simulations are consistent with the experimental patterns and correlate with the joint roughness coefficient(JRC).This study proposed an effective numerical simulation method to reveal the evolution mechanism of fracture flow capacity,taking into account the wear of the fracture surface in shear simulations and the initial stress state of the rock in seepage simulations.This research explains the permeability evolution mechanism of fractured shale from a microscopic perspective,and the proposed numerical simulation method for shear seepage provides a powerful means to uncover the dynamic evolution mechanisms governing fracture permeability.展开更多
To explore the complex thermal-mechanical-chemical behavior in the solid-liquid cast-roll bonding(SLCRB) of Cu/Al cladding strip, numerical simulations were conducted from both macro and micro scales. In macro-scale, ...To explore the complex thermal-mechanical-chemical behavior in the solid-liquid cast-roll bonding(SLCRB) of Cu/Al cladding strip, numerical simulations were conducted from both macro and micro scales. In macro-scale, with birth and death element method, a thermo-mechanical coupled finite element model(FEM) was set up to explore the temperature and contact pressure distribution at the Cu/Al bonding interface in the SLCRB process. Taking these macro-scale simulation results as boundary conditions, we simulated the atom diffusion law of the bonding interface by molecular dynamics(MD) in micro-scale. The results indicate that the temperature in Cu/Al bonding interface deceases from 700 to 320 ℃ from the entrance to the exit of caster, and the peak of contact pressure reaches up to 140 MPa. The interfacial diffusion thickness depends on temperature and rolling reduction, higher temperature results in larger thickness, and the rolling reduction below kiss point leads to significant elongation deformation of cladding strip which yields more newborn interface with fresh metal and make the diffusion layer thinner. The surface roughness of Cu strip was found to be benefit to atoms diffusion in the Cu/Al bonding interface. Meanwhile, combined with the SEM-EDS observation on the microstructure and composition in the bonding interface of the experimental samples acquired from the castrolling bite, it is revealed that the rolling reduction and severe elongation deformation in the solid-solid contact zone below kiss point guarantee the satisfactory metallurgical bonding with thin and smooth diffusion layer. The bonding mechanisms of reactive diffusion, mechanical interlocking and crack bonding are proved to coexist in the SLCRB process.展开更多
Because insulin released by the β-cells of pancreatic islets is the main regulator of glucose levels, the quantitative modeling of their glucose-stimulated insulin secretion is of obvious interest not only to improve...Because insulin released by the β-cells of pancreatic islets is the main regulator of glucose levels, the quantitative modeling of their glucose-stimulated insulin secretion is of obvious interest not only to improve our understanding of the processes involved, but also to allow better assessment of β -cell function in diabetic patients or islet transplant recipients as well as the development of improved artificial or bioartificial pancreas devices. We have recently developed a general, local concentrations-based multiphysics computational model of insulin secretion in avascular pancreatic islets that can be used to calculate insulin secretion for arbitrary geometries of cultured, perifused, transplanted, or encapsulated islets in response to various glucose profiles. Here, experimental results obtained from two different dynamic glucose-stimulated insulin release (GSIR) perifusion studies performed by us following standard procedures are compared to those calculated by the model. Such perifusion studies allow the quantitative assessment of insulin release kinetics under fully controllable experimental conditions of varying external concentrations of glucose, oxygen, or other compounds of interest, and can provide an informative assessment of islet quality and function. The time-profile of the insulin secretion calculated by the model was in good agree- ment with the experimental results obtained with isolated human islets. Detailed spatial distributions of glucose, oxygen, and insulin were calculated and are presented to provide a quantitative visualization of various important aspects of the insulin secretion dynamics in perifused islets.展开更多
In the realization of mechanical structures, achieving stability and balance is a problem commonly encountered by engineers in the field of civil engineering, mechanics, aeronautics, biomechanics and many others. The ...In the realization of mechanical structures, achieving stability and balance is a problem commonly encountered by engineers in the field of civil engineering, mechanics, aeronautics, biomechanics and many others. The study of plate behavior is a very sensitive subject because it is part of the structural elements. The study of the dynamic behavior of free vibration structures is done by modal analysis in order to calculate natural frequencies and modal deformations. In this paper, we present the modal analysis of a thin rectangular plate simply supported. The analytical solution of the differential equation is obtained by applying the method of separating the variables. We are talking about the exact solution of the problem to the limit values. However, numerical methods such as the finite element method allow us to approximate these functions with greater accuracy. It is one of the most powerful computational methods for predicting dynamic response in a complex structure subject to arbitrary boundary conditions. The results obtained by MEF through Ansys 15.0 are then compared with those obtained by the analytical method.展开更多
The aim of this paper is to conduct experimental modal analysis and numerical simulation to verify the structural characteristics of a deployable-retractable wing for aircraft and spacecraft. A modal impact test was c...The aim of this paper is to conduct experimental modal analysis and numerical simulation to verify the structural characteristics of a deployable-retractable wing for aircraft and spacecraft. A modal impact test was conducted in order to determine the free vibration characteristics. Natural frequencies and vibration mode shapes were obtained via measurement in LMS Test. Lab. The frequency response functions were identified and computed by force and acceleration signals, and then mode shapes of this morphing wing structure were subsequently identified by PolyMAX modal parameter estimation method. FEM modal analysis was also implemented and its numerical results convincingly presented the mode shape and natural frequency characteristics were in good agreement with those obtained from experimental modal analysis. Experimental study in this paper focuses on the transverse response of morphing wing as its moveable part is deploying or retreating. Vibration response to different rotation speeds have been collected, managed and analyzed through the use of comparison methodology with each other. Evident phenomena have been discovered including the resonance on which most analysis is focused because of its potential use to generate large amplitude vibration of specific frequency or to avoid such resonant frequencies from a wide spectrum of response. Manufactured deployable-retractable wings are studied in stage of experimental modal analysis, in which some nonlinear vibration resulted should be particularly noted because such wing structure displays a low resonant frequency which is always optimal to be avoided for structural safety and stability.展开更多
文摘Concrete structures may suffer considerable restraint stresses during their hardening period. This is caused by several deformation impacts, especially temperature field changings due to hydration heat and volume changes due to autogenous shrinkage. Mainly affected are massive concrete members, but also the application of new concrete types or the erection of outstanding constructions requires further investigations in this context. 3D-FEM analyses of hydration heat induced temperature development in combination with the well known shrinkage give sufficient results for the deformation impact. The according elastic restraint stresses can be determined with consideration of the concrete’s rising elastic modulus and the restraint degree of the system. But due to duration of the heat flow process, the height of restraint stresses is strongly dependent from the viscoelasticity of the concrete. The viscoelastic effects consist of many components constituted by changing material properties influencing themselves. In practice, different simplified approaches are available for considering this in calculations. Their implementation in time step analyses is not generally admitted and requires expertise. In contrast, present research develops material models needing specific input parameters for every use case. This contribution focuses on a practicable approach considering the superposition of the viscoelastic behaviour of every stress increment in time step FEM analysis. The differentiation between the pure viscoelastic material behaviour (as it is given in the codes for idealistic conditions like creep or relaxation) and the according viscoelastic system response (addicted to the systems variable restraint degree) allows the transfer of this model into practice. As one application example of this approach, the compatibility check and the FEM-based recalculation of the monitoring program of a massive power plant slab will be presented.
基金This research is supported by the Foundation for Talents for Next Century of Shandong University
文摘We’ll study the FEM for a model for compressible miscible displacement in porous media which includes molecular diffusion and mechanical dispersion in one-dimensional space.A class of vertices-edges-elements interpolation operator ink is introduced.With the help of ink(not elliptic projection),the optimal error estimate in L∞(J;L2(Ω)) norm of FEM is proved.
文摘The accurate mechanical analysis of thick-walled pressure vessel structures composed of advanced materials,such as hyperelastic and functionally graded materials(FGMs),is critical for ensuring their safety and optimizing their design.However,conventional numerical methods can face challenges with the non-linearities inherent in hyperelasticity and the complex spatial variations in FGMs.This paper presents a novel hybrid numerical approach combining Physics-Informed Neural Networks(PINNs)with Finite Element Method(FEM)derived data for the robust analysis of thick-walled,axisymmetric,heterogeneous,hyperelastic pressure vessels with elliptical geometries.A PINN framework incorporating neo-Hookean constitutive relations is developed in MATLAB.To enhance training efficiency and accuracy,the PINN’s loss function is augmented with displacement data obtained from high-fidelity FEM simulations performed in ANSYS.The methodology is rigorously validated by comparing PINN-predicted displacement and von Mises stress fields against ANSYS benchmarks for various scenarios of FGMconfigurations(with material properties varying according to a power law)subjected to internal and external pressurization.The results demonstrate excellent agreement between the proposed hybrid PINN-FEMapproach and conventional FEMsolutions across all test cases,accurately capturing complex deformation patterns and stress concentrations.This study highlights the potential of data-augmented PINNs as an effective and accurate computational tool for tackling complex solid mechanics problems involving non-linearmaterials and significant heterogeneity,offering a promising avenue for future research in engineering design and analysis.
基金funded by the Joint Funds of the National Natural Science Foundation of China(Grant No.U23A20671)the Major Project of Inner Mongolia Science and Technology(Grant No.2021ZD0034)the Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engi-neering(Grant No.Z021003).
文摘The dynamic evolution of fracture permeability presents a critical scientific challenge in rock masses.Understanding the mechanisms of rock mass permeability evolution is vital for engineering project design and operations.By integrating the discrete element method(DEM)with the finite element method(FEM),a numerical simulation framework for shear seepage in rough fractured shale has been developed to investigate the dynamic mechanisms of permeability evolution under varying confining pressures and during the shearing process.Numerical simulations were conducted on rough fractured samples under effective confining pressures ranging from 5 MPa to 20 MPa to monitor the aperture and permeability evolution of the fracture.The results of the numerical simulation are consistent with the experimental observations,indicating that both the shearing process and confining pressure significantly influence permeability.Moreover,the magnitude of the confining pressure is a crucial factor influencing the trend in permeability changes.Under a confining pressure of 5 MPa,fracture permeability initially increases significantly but decreases post-shearing.In contrast,a continuous decrease in fracture permeability is observed when the confining pressure exceeds 10 MPa.The results of the shear numerical simulation indicate that the confining pressure restricts fracture dilation during shearing,promotes the generation of rock debris,and decreases both the permeability and transmissivity of the fracture.The wear results obtained from numerical simulations are consistent with the experimental patterns and correlate with the joint roughness coefficient(JRC).This study proposed an effective numerical simulation method to reveal the evolution mechanism of fracture flow capacity,taking into account the wear of the fracture surface in shear simulations and the initial stress state of the rock in seepage simulations.This research explains the permeability evolution mechanism of fractured shale from a microscopic perspective,and the proposed numerical simulation method for shear seepage provides a powerful means to uncover the dynamic evolution mechanisms governing fracture permeability.
基金Funded by the General Program of National Natural Science Foundation of China(Nos.51474189 and 51674222)the Excellent Youth Foundation of Hebei Scientific Committee,China(No.E2018203446)the Scientific Research Foundation of the Higher Education Institutions of Hebei Province,China(No.QN2015214)
文摘To explore the complex thermal-mechanical-chemical behavior in the solid-liquid cast-roll bonding(SLCRB) of Cu/Al cladding strip, numerical simulations were conducted from both macro and micro scales. In macro-scale, with birth and death element method, a thermo-mechanical coupled finite element model(FEM) was set up to explore the temperature and contact pressure distribution at the Cu/Al bonding interface in the SLCRB process. Taking these macro-scale simulation results as boundary conditions, we simulated the atom diffusion law of the bonding interface by molecular dynamics(MD) in micro-scale. The results indicate that the temperature in Cu/Al bonding interface deceases from 700 to 320 ℃ from the entrance to the exit of caster, and the peak of contact pressure reaches up to 140 MPa. The interfacial diffusion thickness depends on temperature and rolling reduction, higher temperature results in larger thickness, and the rolling reduction below kiss point leads to significant elongation deformation of cladding strip which yields more newborn interface with fresh metal and make the diffusion layer thinner. The surface roughness of Cu strip was found to be benefit to atoms diffusion in the Cu/Al bonding interface. Meanwhile, combined with the SEM-EDS observation on the microstructure and composition in the bonding interface of the experimental samples acquired from the castrolling bite, it is revealed that the rolling reduction and severe elongation deformation in the solid-solid contact zone below kiss point guarantee the satisfactory metallurgical bonding with thin and smooth diffusion layer. The bonding mechanisms of reactive diffusion, mechanical interlocking and crack bonding are proved to coexist in the SLCRB process.
文摘Because insulin released by the β-cells of pancreatic islets is the main regulator of glucose levels, the quantitative modeling of their glucose-stimulated insulin secretion is of obvious interest not only to improve our understanding of the processes involved, but also to allow better assessment of β -cell function in diabetic patients or islet transplant recipients as well as the development of improved artificial or bioartificial pancreas devices. We have recently developed a general, local concentrations-based multiphysics computational model of insulin secretion in avascular pancreatic islets that can be used to calculate insulin secretion for arbitrary geometries of cultured, perifused, transplanted, or encapsulated islets in response to various glucose profiles. Here, experimental results obtained from two different dynamic glucose-stimulated insulin release (GSIR) perifusion studies performed by us following standard procedures are compared to those calculated by the model. Such perifusion studies allow the quantitative assessment of insulin release kinetics under fully controllable experimental conditions of varying external concentrations of glucose, oxygen, or other compounds of interest, and can provide an informative assessment of islet quality and function. The time-profile of the insulin secretion calculated by the model was in good agree- ment with the experimental results obtained with isolated human islets. Detailed spatial distributions of glucose, oxygen, and insulin were calculated and are presented to provide a quantitative visualization of various important aspects of the insulin secretion dynamics in perifused islets.
文摘In the realization of mechanical structures, achieving stability and balance is a problem commonly encountered by engineers in the field of civil engineering, mechanics, aeronautics, biomechanics and many others. The study of plate behavior is a very sensitive subject because it is part of the structural elements. The study of the dynamic behavior of free vibration structures is done by modal analysis in order to calculate natural frequencies and modal deformations. In this paper, we present the modal analysis of a thin rectangular plate simply supported. The analytical solution of the differential equation is obtained by applying the method of separating the variables. We are talking about the exact solution of the problem to the limit values. However, numerical methods such as the finite element method allow us to approximate these functions with greater accuracy. It is one of the most powerful computational methods for predicting dynamic response in a complex structure subject to arbitrary boundary conditions. The results obtained by MEF through Ansys 15.0 are then compared with those obtained by the analytical method.
文摘The aim of this paper is to conduct experimental modal analysis and numerical simulation to verify the structural characteristics of a deployable-retractable wing for aircraft and spacecraft. A modal impact test was conducted in order to determine the free vibration characteristics. Natural frequencies and vibration mode shapes were obtained via measurement in LMS Test. Lab. The frequency response functions were identified and computed by force and acceleration signals, and then mode shapes of this morphing wing structure were subsequently identified by PolyMAX modal parameter estimation method. FEM modal analysis was also implemented and its numerical results convincingly presented the mode shape and natural frequency characteristics were in good agreement with those obtained from experimental modal analysis. Experimental study in this paper focuses on the transverse response of morphing wing as its moveable part is deploying or retreating. Vibration response to different rotation speeds have been collected, managed and analyzed through the use of comparison methodology with each other. Evident phenomena have been discovered including the resonance on which most analysis is focused because of its potential use to generate large amplitude vibration of specific frequency or to avoid such resonant frequencies from a wide spectrum of response. Manufactured deployable-retractable wings are studied in stage of experimental modal analysis, in which some nonlinear vibration resulted should be particularly noted because such wing structure displays a low resonant frequency which is always optimal to be avoided for structural safety and stability.
文摘泥石流作为高破坏性混合流体,其携带的块石对框架结构的冲击机制尚未得到充分研究。为揭示含块石泥石流冲击框架结构的动力响应与损伤机制,基于光滑粒子流体动力学-离散元法-有限元法(smoothed particle hydrodynamics-discrete element method-finite element method,SPH-DEM-FEM)耦合数值方法,构建流体-块石-结构的多耦合数值模型,模拟不同冲击速度和角度下框架结构的损伤过程,并结合两相溃坝试验来验证模型有效性。研究表明:泥石流冲击导致结构损伤经历“接触-扩散-反弹-堆积/冲击”四个阶段,当冲击速度超过6.00 m/s时,结构产生不可恢复损伤,且块石阻隔作用使流体上部冲击力大于下部,引发结构中部最先发生集中损伤;此外,冲击力峰值随速度和角度增大呈非线性增长,10.00 m/s与90°工况下框架柱底冲击力达497.17 kN,超过结构抗冲击承载力;最后,数值模拟与经验公式所得冲击力结果误差在13.95%~29.00%,数量级一致,验证了模型可靠性。研究结果可为泥石流高发区框架结构的抗冲击设计提供参考。
