针对飞机典型部位在遭到高速破片攻击后结构整体的战伤状态及破片的剩余行为开展数值模拟。应用LS-DYNA软件,结合有限单元方法(finite element method,FEM)和光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)两者的优势,建立...针对飞机典型部位在遭到高速破片攻击后结构整体的战伤状态及破片的剩余行为开展数值模拟。应用LS-DYNA软件,结合有限单元方法(finite element method,FEM)和光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)两者的优势,建立自适应的FEM-SPH耦合模拟方法,并构建2种飞机典型部位的计算模型,采用六面体网格局部细化方法实现了核心位置的精确模拟,并进行试验来验证数值模型;开展了一系列高速冲击战伤模拟,对比了不同工况下破片高速冲击结构后形成的碎片云和破口形貌,并对破片的剩余速度和质量进行分析,确定了破片在结构蒙皮上的临界跳飞角。结果表明:自适应FEM-SPH耦合算法的计算结果与试验结果吻合良好,能够对破片高速冲击战伤进行有效准确模拟;碎片云分布形状随破片速度增加变得狭长,冲击角度会改变碎片云和结构破口形状朝向;碎片云高度和扩散速度随破片速度或角度的变化趋势基本一致并都呈线性关系;破片的速度减少量不随初始速度变化,质量减少量则与冲击速度正相关,两者与冲击角度都负相关;破片临界跳飞角与冲击速度大小基本呈线性关系。展开更多
泥石流是我国西南山区常见的地质灾害。架空输电杆塔在泥石流的冲击下往往发生基础破坏甚至会造成杆塔倒塌。首先采用光滑粒子流体动力学(smoothed particle hydrodynamics,简称SPH)方法和有限元方法(finite element method,简称FEM)相...泥石流是我国西南山区常见的地质灾害。架空输电杆塔在泥石流的冲击下往往发生基础破坏甚至会造成杆塔倒塌。首先采用光滑粒子流体动力学(smoothed particle hydrodynamics,简称SPH)方法和有限元方法(finite element method,简称FEM)相耦合的三维数值方法模拟了泥石流对杆塔基础的冲击作用;在与相关模型试验结果验证的基础上,开展了不同泥石流密度、黏度系数及初始速度条件下对输电塔基础的冲击力作用的参数分析;研究结果表明:随着泥石流初始速度的增加,冲击力峰值会随之增大;前排基础的冲击力峰值均大于后排基础;泥石流冲击过程特性受到泥石流密度和黏度系数影响。与稀性泥石流相比:黏性泥石流冲击基础后,基础下游真空区相对要小;此外,将数值模拟结果与Kwan冲击力公式及铁二院推荐的冲击压力设计公式预测值进行对比分析可以发现:Kwan冲击力公式能较好地预测出基础所受泥石流冲击力的平均趋势,最大预测误差低于30%,铁二院公式预测的稀性和黏性泥石流的冲击压力平均偏低分别约17%和28%。相关研究结果有望为泥石流频发区域输电塔基础的设计和风险评估提供一定的参考依据。展开更多
通过设计开展静爆实验,对土中浅埋爆炸成坑效应及影响因素进行研究。基于所获取的静爆试验数据,建立SPH-FEM数值模型,对土中浅埋爆炸过程及爆坑参数进行分析,验证SPH-FEM方法模拟土中浅埋爆炸的有效性。进一步通过数值计算对比,讨论了SP...通过设计开展静爆实验,对土中浅埋爆炸成坑效应及影响因素进行研究。基于所获取的静爆试验数据,建立SPH-FEM数值模型,对土中浅埋爆炸过程及爆坑参数进行分析,验证SPH-FEM方法模拟土中浅埋爆炸的有效性。进一步通过数值计算对比,讨论了SPH-FEM方法与ALE方法在土中爆炸数值模拟过程中的异同、爆坑半径与炸药埋设深度的关系以及炸药量与爆坑体积的关系,并通过数学拟合得到了9.35 kg TNT条件下,炸药埋设深度与爆坑半径的经验曲线公式。研究表明:SPH-FEM方法较ALE方法模拟土中浅埋爆炸更加精细且计算效率更高。所建立的炸药埋设深度与爆坑半径的经验曲线公式可以有效预测浅埋爆炸的爆坑半径,且随着炸药埋设深度的增加,爆坑半径将先随之增大到峰值,之后爆坑半径将会随炸药埋设深度的增加而迅速减小。展开更多
准确计算户内变电站大型、复杂的噪声场分布,进而评价可采用降噪措施的减噪效果,是解决户内变电站噪声污染的关键问题。为此,综合声学有限元法(finite element method,FEM)求解复杂声场收敛性好及精度高的优点,及声学边界元法(boundary ...准确计算户内变电站大型、复杂的噪声场分布,进而评价可采用降噪措施的减噪效果,是解决户内变电站噪声污染的关键问题。为此,综合声学有限元法(finite element method,FEM)求解复杂声场收敛性好及精度高的优点,及声学边界元法(boundary element method,BEM)降维求解大型声场的优势,提出了一种基于声学FEM-BEM的户内变电站噪声场求解算法。首先,建立变电站内部声源声固耦合模型,采用声学FEM求解混响噪声作用下的声固耦合响应;然后,基于声学FEM-BEM耦合理论,求解内、外耦合边界处结构单元受声固耦合激励产生的位移及应力载荷;最后,根据声压及应力载荷激发的外场声波扩散模型,基于常规Gauss数值积分法,建立外部空间声域2维BEM声学积分方程,求解外部声场。该算法在湖南某110 kV户内变电站噪声场的求解分析中得到了成功应用,与实测值的相对误差为3.61%~4.87%。展开更多
碎屑流是我国山区最危险的地质灾害之一,山区桥墩常受到碎屑流冲击而开裂、倾斜甚至倒塌,给山区桥梁建设、运营带来严重的安全隐患。采用离散元方法(discrete element method,DEM)和有限元方法(finite element method,FEM)耦合的三维数...碎屑流是我国山区最危险的地质灾害之一,山区桥墩常受到碎屑流冲击而开裂、倾斜甚至倒塌,给山区桥梁建设、运营带来严重的安全隐患。采用离散元方法(discrete element method,DEM)和有限元方法(finite element method,FEM)耦合的三维数值模拟方法模拟了碎屑流对双柱式桥墩的冲击效应,并结合斜槽试验,验证了耦合方法的准确性,进一步分析了碎屑流冲击坡度、距离和体积密度对桥墩冲击力的影响规律。结果表明,最大冲击力与碎屑流冲击坡度、距离和体积密度分别呈幂函数(指数大于1)、幂函数(指数小于1)和线性正相关。冲击坡度、距离和体积密度对最大冲击力的敏感度值分别为3.012、0.202、0.804,在桥梁碎屑流灾害防治时需重视冲击坡度和体积密度的影响。将冲击力的数值模拟值与流体动力学模型预测值对比分析表明,流体动力学模型理论公式能较好地预测桥墩所受的最大冲击力,最大预测误差低于23.6%。相关研究结果可为山区桥梁碎屑流灾害防治与设计提供一定的参考依据。展开更多
A numerical model using the coupled smoothed particle hydrodynamics-finite element method (SPH-FEM) approach is presented for analysis of structures under blast loads. The analyses on two numerical cases, one for free...A numerical model using the coupled smoothed particle hydrodynamics-finite element method (SPH-FEM) approach is presented for analysis of structures under blast loads. The analyses on two numerical cases, one for free field explosive and the other for structural response under blast loads, are performed to model the whole processes from the propagation of the pressure wave to the response of structures. Based on the simulation, it is concluded that this model can be used for reasonably accurate explosive analysis of structures. The resulting information would be valuable for protecting structures under blast loads.展开更多
CO_(2) drilling is a promising underbalance drilling technology with great advantages,such as lower cutting force,intense cooling and excellent lubrication.However,in the underbalance drilling,the mechanism of the cou...CO_(2) drilling is a promising underbalance drilling technology with great advantages,such as lower cutting force,intense cooling and excellent lubrication.However,in the underbalance drilling,the mechanism of the coupling CO_(2) jet and polycrystalline-diamond-compact(PDC)cutter are still unclear.Whereby,we established a coupled smoothed particle hydrodynamics/finite element method(SPH/FEM)model to simulate the composite rock-breaking of high-pressure CO_(2) jet&PDC cutter.Combined with the experimental research results,the mechanism of composite rock-breaking is studied from the perspectives of rock stress field,cutting force and jet field.The results show that the composite rock-breaking can effectively relieve the influence of vibration and shock on PDC cutter.Meanwhile,the high-pressure CO_(2) jet has a positive effect on carrying rock debris,which can effectively reduce the temperature rising and the thermal wear of the PDC cutter.In addition,the effects of CO_(2) jet parameters on composite rock-breaking were studied,such as jet impact velocity,nozzle diameter,jet injection angle and impact distance.The studies show that when the impact velocity of the CO_(2) jet is greater than 250 m/s,the CO_(2) jet could quickly break the rock.It is found that the optimal range of nozzle diameter is 1.5–2.5 mm,the best injection angle of CO_(2) jet is 60,the optimal impact distance is 10 times the nozzle diameter.The above studies could provide theoretical supports and technical guidance for composite rock-breaking,which is useful for the CO_(2) underbalance drilling and drill bit design.展开更多
The coupling iteration (CI) of the finite element method(FEM) is used to simulate the magnetic and mechanical characteristics for a GMM actuator. The convergent ability under different prestress and different load typ...The coupling iteration (CI) of the finite element method(FEM) is used to simulate the magnetic and mechanical characteristics for a GMM actuator. The convergent ability under different prestress and different load types is investigated. Then the calculated deformations are compared with the experimental values. The results convince that the CI of FEM is suitable for the simulation of energy coupling and transformation mechanism of the GMM. At last, the output deformation properties are studied under different input currents, showing that there is a good compromise between good linearity and large strain under the prestress 6 MPa.展开更多
The model of pressure solution for granular aggregate was introduced into the FEM code for analysis of thermo-hydro-mechanical (T-H-M) coupling in porous medium. Aiming at a hypothetical nuclear waste repository in an...The model of pressure solution for granular aggregate was introduced into the FEM code for analysis of thermo-hydro-mechanical (T-H-M) coupling in porous medium. Aiming at a hypothetical nuclear waste repository in an unsaturated quartz rock mass, two computation conditions were designed:1) the porosity and the permeability of rock mass are functions of pressure solution; 2) the porosity and the permeability are constants. Then the corresponding numerical simulations for a disposal period of 4 a were carried out, and the states of temperatures, porosities and permeabilities, pore pressures, flow velocities and stresses in the rock mass were investigated. The results show that at the end of the calculation in Case 1, pressure solution makes the porosities and the permeabilities decrease to 10%-45% and 0.05%-1.4% of their initial values, respectively. Under the action of the release heat of nuclear waste, the negative pore pressures both in Case 1 and Case 2 are 1.2-1.4 and 1.01-1.06 times of the initial values, respectively. So, the former represents an obvious effect of pressure solution. The magnitudes and distributions of stresses within the rock mass in the two calculation cases are the same.展开更多
In this study, the magnetohydrodynamic (MHD) flow through a circular pipe under the influence of a transverse mag- netic field when the outside medium is also electrically conducting is solved numerically by using FEM...In this study, the magnetohydrodynamic (MHD) flow through a circular pipe under the influence of a transverse mag- netic field when the outside medium is also electrically conducting is solved numerically by using FEM-BEM coupling approach. The coupled partial differential equations defined for the interior medium are transformed into homogenous modified Helmholtz equations. For the exterior medium on an infinite region, the Laplace equation is considered for the exterior magnetic field. Unknowns in the equations are also related with the corresponding Dirichlet and Neumann type coupled boundary conditions. Unknown values of the magnetic field on the boundary and for the exterior region are obtained by using BEM, and the unknown velocity and magnetic field inside the pipe are obtained by using SUPG type stabilized FEM. Computations are carried for very high values of magnetic Reynolds numbers Rm1, Reynolds number Re and magnetic pressure Rh of the fluid. The results show that using stabilized method enables us to get stable and accurate numerical approximations consistent with the physical configuration of the problem over rough mesh which also results a cheap computational cost.展开更多
文摘针对飞机典型部位在遭到高速破片攻击后结构整体的战伤状态及破片的剩余行为开展数值模拟。应用LS-DYNA软件,结合有限单元方法(finite element method,FEM)和光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)两者的优势,建立自适应的FEM-SPH耦合模拟方法,并构建2种飞机典型部位的计算模型,采用六面体网格局部细化方法实现了核心位置的精确模拟,并进行试验来验证数值模型;开展了一系列高速冲击战伤模拟,对比了不同工况下破片高速冲击结构后形成的碎片云和破口形貌,并对破片的剩余速度和质量进行分析,确定了破片在结构蒙皮上的临界跳飞角。结果表明:自适应FEM-SPH耦合算法的计算结果与试验结果吻合良好,能够对破片高速冲击战伤进行有效准确模拟;碎片云分布形状随破片速度增加变得狭长,冲击角度会改变碎片云和结构破口形状朝向;碎片云高度和扩散速度随破片速度或角度的变化趋势基本一致并都呈线性关系;破片的速度减少量不随初始速度变化,质量减少量则与冲击速度正相关,两者与冲击角度都负相关;破片临界跳飞角与冲击速度大小基本呈线性关系。
文摘泥石流是我国西南山区常见的地质灾害。架空输电杆塔在泥石流的冲击下往往发生基础破坏甚至会造成杆塔倒塌。首先采用光滑粒子流体动力学(smoothed particle hydrodynamics,简称SPH)方法和有限元方法(finite element method,简称FEM)相耦合的三维数值方法模拟了泥石流对杆塔基础的冲击作用;在与相关模型试验结果验证的基础上,开展了不同泥石流密度、黏度系数及初始速度条件下对输电塔基础的冲击力作用的参数分析;研究结果表明:随着泥石流初始速度的增加,冲击力峰值会随之增大;前排基础的冲击力峰值均大于后排基础;泥石流冲击过程特性受到泥石流密度和黏度系数影响。与稀性泥石流相比:黏性泥石流冲击基础后,基础下游真空区相对要小;此外,将数值模拟结果与Kwan冲击力公式及铁二院推荐的冲击压力设计公式预测值进行对比分析可以发现:Kwan冲击力公式能较好地预测出基础所受泥石流冲击力的平均趋势,最大预测误差低于30%,铁二院公式预测的稀性和黏性泥石流的冲击压力平均偏低分别约17%和28%。相关研究结果有望为泥石流频发区域输电塔基础的设计和风险评估提供一定的参考依据。
文摘通过设计开展静爆实验,对土中浅埋爆炸成坑效应及影响因素进行研究。基于所获取的静爆试验数据,建立SPH-FEM数值模型,对土中浅埋爆炸过程及爆坑参数进行分析,验证SPH-FEM方法模拟土中浅埋爆炸的有效性。进一步通过数值计算对比,讨论了SPH-FEM方法与ALE方法在土中爆炸数值模拟过程中的异同、爆坑半径与炸药埋设深度的关系以及炸药量与爆坑体积的关系,并通过数学拟合得到了9.35 kg TNT条件下,炸药埋设深度与爆坑半径的经验曲线公式。研究表明:SPH-FEM方法较ALE方法模拟土中浅埋爆炸更加精细且计算效率更高。所建立的炸药埋设深度与爆坑半径的经验曲线公式可以有效预测浅埋爆炸的爆坑半径,且随着炸药埋设深度的增加,爆坑半径将先随之增大到峰值,之后爆坑半径将会随炸药埋设深度的增加而迅速减小。
文摘准确计算户内变电站大型、复杂的噪声场分布,进而评价可采用降噪措施的减噪效果,是解决户内变电站噪声污染的关键问题。为此,综合声学有限元法(finite element method,FEM)求解复杂声场收敛性好及精度高的优点,及声学边界元法(boundary element method,BEM)降维求解大型声场的优势,提出了一种基于声学FEM-BEM的户内变电站噪声场求解算法。首先,建立变电站内部声源声固耦合模型,采用声学FEM求解混响噪声作用下的声固耦合响应;然后,基于声学FEM-BEM耦合理论,求解内、外耦合边界处结构单元受声固耦合激励产生的位移及应力载荷;最后,根据声压及应力载荷激发的外场声波扩散模型,基于常规Gauss数值积分法,建立外部空间声域2维BEM声学积分方程,求解外部声场。该算法在湖南某110 kV户内变电站噪声场的求解分析中得到了成功应用,与实测值的相对误差为3.61%~4.87%。
文摘碎屑流是我国山区最危险的地质灾害之一,山区桥墩常受到碎屑流冲击而开裂、倾斜甚至倒塌,给山区桥梁建设、运营带来严重的安全隐患。采用离散元方法(discrete element method,DEM)和有限元方法(finite element method,FEM)耦合的三维数值模拟方法模拟了碎屑流对双柱式桥墩的冲击效应,并结合斜槽试验,验证了耦合方法的准确性,进一步分析了碎屑流冲击坡度、距离和体积密度对桥墩冲击力的影响规律。结果表明,最大冲击力与碎屑流冲击坡度、距离和体积密度分别呈幂函数(指数大于1)、幂函数(指数小于1)和线性正相关。冲击坡度、距离和体积密度对最大冲击力的敏感度值分别为3.012、0.202、0.804,在桥梁碎屑流灾害防治时需重视冲击坡度和体积密度的影响。将冲击力的数值模拟值与流体动力学模型预测值对比分析表明,流体动力学模型理论公式能较好地预测桥墩所受的最大冲击力,最大预测误差低于23.6%。相关研究结果可为山区桥梁碎屑流灾害防治与设计提供一定的参考依据。
基金National Basic Research Program (973) of China (No. 2002CB412709)the National Natural Science Foun-dation of China (No. 50378054)
文摘A numerical model using the coupled smoothed particle hydrodynamics-finite element method (SPH-FEM) approach is presented for analysis of structures under blast loads. The analyses on two numerical cases, one for free field explosive and the other for structural response under blast loads, are performed to model the whole processes from the propagation of the pressure wave to the response of structures. Based on the simulation, it is concluded that this model can be used for reasonably accurate explosive analysis of structures. The resulting information would be valuable for protecting structures under blast loads.
基金This work was supported by the National Natural Science Foundation of China(No.52004236)Sichuan Science and Technology Program(No.2021JDRC0114)+4 种基金the Starting Project of Southwest Petroleum University(No.2019QHZ009)the China Postdoctoral Science Foundation(No.2020M673285)the Open Project Program of Key Laboratory of Groundwater Resources and Environment(Jilin University)Ministry of Education(No.202005009KF)the Chinese Scholarship Council(CSC)funding(No.202008515107).
文摘CO_(2) drilling is a promising underbalance drilling technology with great advantages,such as lower cutting force,intense cooling and excellent lubrication.However,in the underbalance drilling,the mechanism of the coupling CO_(2) jet and polycrystalline-diamond-compact(PDC)cutter are still unclear.Whereby,we established a coupled smoothed particle hydrodynamics/finite element method(SPH/FEM)model to simulate the composite rock-breaking of high-pressure CO_(2) jet&PDC cutter.Combined with the experimental research results,the mechanism of composite rock-breaking is studied from the perspectives of rock stress field,cutting force and jet field.The results show that the composite rock-breaking can effectively relieve the influence of vibration and shock on PDC cutter.Meanwhile,the high-pressure CO_(2) jet has a positive effect on carrying rock debris,which can effectively reduce the temperature rising and the thermal wear of the PDC cutter.In addition,the effects of CO_(2) jet parameters on composite rock-breaking were studied,such as jet impact velocity,nozzle diameter,jet injection angle and impact distance.The studies show that when the impact velocity of the CO_(2) jet is greater than 250 m/s,the CO_(2) jet could quickly break the rock.It is found that the optimal range of nozzle diameter is 1.5–2.5 mm,the best injection angle of CO_(2) jet is 60,the optimal impact distance is 10 times the nozzle diameter.The above studies could provide theoretical supports and technical guidance for composite rock-breaking,which is useful for the CO_(2) underbalance drilling and drill bit design.
基金This project is supported by National Natural Science Foundation of China (No.50077019).
文摘The coupling iteration (CI) of the finite element method(FEM) is used to simulate the magnetic and mechanical characteristics for a GMM actuator. The convergent ability under different prestress and different load types is investigated. Then the calculated deformations are compared with the experimental values. The results convince that the CI of FEM is suitable for the simulation of energy coupling and transformation mechanism of the GMM. At last, the output deformation properties are studied under different input currents, showing that there is a good compromise between good linearity and large strain under the prestress 6 MPa.
基金Project(2010CB732101) supported by the National Basic Research Program of ChinaProject(51079145) supported by the National Natural Science Foundation of ChinaProject(2009BAK53B03) supported by the National Key Technology R&D Program of China
文摘The model of pressure solution for granular aggregate was introduced into the FEM code for analysis of thermo-hydro-mechanical (T-H-M) coupling in porous medium. Aiming at a hypothetical nuclear waste repository in an unsaturated quartz rock mass, two computation conditions were designed:1) the porosity and the permeability of rock mass are functions of pressure solution; 2) the porosity and the permeability are constants. Then the corresponding numerical simulations for a disposal period of 4 a were carried out, and the states of temperatures, porosities and permeabilities, pore pressures, flow velocities and stresses in the rock mass were investigated. The results show that at the end of the calculation in Case 1, pressure solution makes the porosities and the permeabilities decrease to 10%-45% and 0.05%-1.4% of their initial values, respectively. Under the action of the release heat of nuclear waste, the negative pore pressures both in Case 1 and Case 2 are 1.2-1.4 and 1.01-1.06 times of the initial values, respectively. So, the former represents an obvious effect of pressure solution. The magnitudes and distributions of stresses within the rock mass in the two calculation cases are the same.
文摘In this study, the magnetohydrodynamic (MHD) flow through a circular pipe under the influence of a transverse mag- netic field when the outside medium is also electrically conducting is solved numerically by using FEM-BEM coupling approach. The coupled partial differential equations defined for the interior medium are transformed into homogenous modified Helmholtz equations. For the exterior medium on an infinite region, the Laplace equation is considered for the exterior magnetic field. Unknowns in the equations are also related with the corresponding Dirichlet and Neumann type coupled boundary conditions. Unknown values of the magnetic field on the boundary and for the exterior region are obtained by using BEM, and the unknown velocity and magnetic field inside the pipe are obtained by using SUPG type stabilized FEM. Computations are carried for very high values of magnetic Reynolds numbers Rm1, Reynolds number Re and magnetic pressure Rh of the fluid. The results show that using stabilized method enables us to get stable and accurate numerical approximations consistent with the physical configuration of the problem over rough mesh which also results a cheap computational cost.