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.展开更多
坝体抗震设计和评估需要准确计算无限水库动力响应.基于比例边界有限元法(scaled boundary finite element method,SBFEM)力学推导技术,推导了顺河向地震激励下等横截面无限水域频域响应计算公式,利用Fourier逆变换建立了时域响应控制方...坝体抗震设计和评估需要准确计算无限水库动力响应.基于比例边界有限元法(scaled boundary finite element method,SBFEM)力学推导技术,推导了顺河向地震激励下等横截面无限水域频域响应计算公式,利用Fourier逆变换建立了时域响应控制方程,通过线性叠加推导了顺河、横河、竖直三向组合地震激励下的无限水域频域和时域响应的SBFEM计算公式.结合有限元法,建立了无限水库频域和时域响应的FEM-SBFEM耦合方程.分析了地震激励下的二维、三维等横截面无限水库频域、时域响应,数值验证了所建立计算公式的正确性.所发展的FEM-SBFEM公式体系可推广应用于库底库岸具有吸收性的、横截面有任意几何形状的无限水库谐响应及瞬态响应分析.展开更多
高位滑坡对建筑集群的冲击破坏时常导致严重的人员伤亡,基于光滑粒子流体动力学-离散元法-有限元法(smoothed particle hydrodynamics-discrete element method-finite element method,SPH-DEM-FEM)耦合的数值模型,开展了高位滑坡对框...高位滑坡对建筑集群的冲击破坏时常导致严重的人员伤亡,基于光滑粒子流体动力学-离散元法-有限元法(smoothed particle hydrodynamics-discrete element method-finite element method,SPH-DEM-FEM)耦合的数值模型,开展了高位滑坡对框架结构建筑群的冲击过程、建筑结构破坏机理、冲击力时程与框架柱关键点应力和弯矩等动力机制研究。研究结果表明:SPH-DEM-FEM耦合数值方法能够有效地模拟碎石土滑坡中土(SPH)石(DEM)混合物的抛射弹跳、爬高绕流冲击运动过程。考虑了常规建筑垂直、平行于滑坡流向的三排建筑组合布局,位于滑坡近端的纵向排列建筑表现为连续性倾倒破坏,横向排列的建筑则呈现整体倾倒破坏;因前排建筑群对滑坡冲击能量的耗散及滑坡自身摩擦耗能,位于滑坡后端建筑表现为引流面墙体和前排柱发生局部破坏,结构保持稳定,损毁程度依次为上游无建筑缓冲耗能的建筑>有横向排列的建筑>有纵向排列的建筑;纵向、横向排列的建筑冲击力衰减幅度分别31%、21%。横向框架建筑整体倾倒的损毁机制表现为框架柱的直接剪断或节点塑形铰链失效;纵向框架建筑连续性倾倒的损毁机制表现为前排框架柱的失效引起后排框架柱轴向压力和极限弯矩增加,持续冲击荷载超过其极限弯矩致使后排框架柱发生弯曲破坏,最终结构倾倒。系统能量在动能、内能和摩擦耗能间转化,其中摩擦耗能占65.5%,结构耗能占23.6%,动能快速下降与内能急剧增加是建筑破坏的关键特征。展开更多
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.展开更多
文摘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.
文摘坝体抗震设计和评估需要准确计算无限水库动力响应.基于比例边界有限元法(scaled boundary finite element method,SBFEM)力学推导技术,推导了顺河向地震激励下等横截面无限水域频域响应计算公式,利用Fourier逆变换建立了时域响应控制方程,通过线性叠加推导了顺河、横河、竖直三向组合地震激励下的无限水域频域和时域响应的SBFEM计算公式.结合有限元法,建立了无限水库频域和时域响应的FEM-SBFEM耦合方程.分析了地震激励下的二维、三维等横截面无限水库频域、时域响应,数值验证了所建立计算公式的正确性.所发展的FEM-SBFEM公式体系可推广应用于库底库岸具有吸收性的、横截面有任意几何形状的无限水库谐响应及瞬态响应分析.
文摘高位滑坡对建筑集群的冲击破坏时常导致严重的人员伤亡,基于光滑粒子流体动力学-离散元法-有限元法(smoothed particle hydrodynamics-discrete element method-finite element method,SPH-DEM-FEM)耦合的数值模型,开展了高位滑坡对框架结构建筑群的冲击过程、建筑结构破坏机理、冲击力时程与框架柱关键点应力和弯矩等动力机制研究。研究结果表明:SPH-DEM-FEM耦合数值方法能够有效地模拟碎石土滑坡中土(SPH)石(DEM)混合物的抛射弹跳、爬高绕流冲击运动过程。考虑了常规建筑垂直、平行于滑坡流向的三排建筑组合布局,位于滑坡近端的纵向排列建筑表现为连续性倾倒破坏,横向排列的建筑则呈现整体倾倒破坏;因前排建筑群对滑坡冲击能量的耗散及滑坡自身摩擦耗能,位于滑坡后端建筑表现为引流面墙体和前排柱发生局部破坏,结构保持稳定,损毁程度依次为上游无建筑缓冲耗能的建筑>有横向排列的建筑>有纵向排列的建筑;纵向、横向排列的建筑冲击力衰减幅度分别31%、21%。横向框架建筑整体倾倒的损毁机制表现为框架柱的直接剪断或节点塑形铰链失效;纵向框架建筑连续性倾倒的损毁机制表现为前排框架柱的失效引起后排框架柱轴向压力和极限弯矩增加,持续冲击荷载超过其极限弯矩致使后排框架柱发生弯曲破坏,最终结构倾倒。系统能量在动能、内能和摩擦耗能间转化,其中摩擦耗能占65.5%,结构耗能占23.6%,动能快速下降与内能急剧增加是建筑破坏的关键特征。
基金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.
