The tunnel subjected to strike-slip fault dislocation exhibits severe and catastrophic damage.The existing analysis models frequently assume uniform fault displacement and fixed fault plane position.In contrast,post-e...The tunnel subjected to strike-slip fault dislocation exhibits severe and catastrophic damage.The existing analysis models frequently assume uniform fault displacement and fixed fault plane position.In contrast,post-earthquake observations indicate that the displacement near the fault zone is typically nonuniform,and the fault plane position is uncertain.In this study,we first established a series of improved governing equations to analyze the mechanical response of tunnels under strike-slip fault dislocation.The proposed methodology incorporated key factors such as nonuniform fault displacement and uncertain fault plane position into the governing equations,thereby significantly enhancing the applicability range and accuracy of the model.In contrast to previous analytical models,the maximum computational error has decreased from 57.1%to 1.1%.Subsequently,we conducted a rigorous validation of the proposed methodology by undertaking a comparative analysis with a 3D finite element numerical model,and the results from both approaches exhibited a high degree of qualitative and quantitative agreement with a maximum error of 9.9%.Finally,the proposed methodology was utilized to perform a parametric analysis to explore the effects of various parameters,such as fault displacement,fault zone width,fault zone strength,the ratio of maximum fault displacement of the hanging wall to the footwall,and fault plane position,on the response of tunnels subjected to strike-slip fault dislocation.The findings indicate a progressive increase in the peak internal forces of the tunnel with the rise in fault displacement and fault zone strength.Conversely,an augmentation in fault zone width is found to contribute to a decrease in the peak internal forces.For example,for a fault zone width of 10 m,the peak values of bending moment,shear force,and axial force are approximately 46.9%,102.4%,and 28.7% higher,respectively,compared to those observed for a fault zone width of 50 m.Furthermore,the position of the peak internal forces is influenced by variations in the ratio of maximum fault displacement of the hanging wall to footwall and the fault plane location,while the peak values of shear force and axial force always align with the fault plane.The maximum peak internal forces are observed when the footwall exclusively bears the entirety of the fault displacement,corresponding to a ratio of 0:1.The peak values of bending moment,shear force,and axial force for the ratio of 0:1 amount to approximately 123.8%,148.6%,and 111.1% of those for the ratio of 0.5:0.5,respectively.展开更多
准确计算户内变电站大型、复杂的噪声场分布,进而评价可采用降噪措施的减噪效果,是解决户内变电站噪声污染的关键问题。为此,综合声学有限元法(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%。展开更多
Longitudinal cracks on the tunnel lining significantly influence the performance of tunnels in operation.In this study,we propose a semi-analytical method that provides a simple and effective way to calculate the inte...Longitudinal cracks on the tunnel lining significantly influence the performance of tunnels in operation.In this study,we propose a semi-analytical method that provides a simple and effective way to calculate the internal forces of tunnel linings with multiple cracks.The semi-analytical solution is obtained using structural analysis considering the flexural rigidity for the cracked longitudinal section of the tunnel lining.Then the proposed solution is verified numerically.Using the proposed method,the influences of the crack depth and the number of cracks on the bending moment and modified crack tip stress are investigated.With the increase in crack depth,the bending moment of lining scetion adjacent to the crack decreases,while the bending moment of lining scetion far away from the crack increases slightly.The more the number of cracks in a tunnel lining,the easier the new cracks initiated.展开更多
The calculation of frost heaving with ice lens formation is still not standard for construction projects using artificial ground freezing(AGF).In fine-grained material,ice lenses may initiate and lead to significant h...The calculation of frost heaving with ice lens formation is still not standard for construction projects using artificial ground freezing(AGF).In fine-grained material,ice lenses may initiate and lead to significant heaving at the ground surface,which should be considered in advance.However,the complex processes during ice lens formation are still not fully understood and difficult to capture in a simple approach.In the past,the semi-analytical approach of Konrad and Morgenstern used one soil constant,the“segregation potential(SP)”.It has been mainly and most successfully applied to the heave calculation of natural-induced soil freezing in cold regions.Its application to AGF has been so far unsuccessful.To solve this,a new semi-analytical approach is presented in this paper.It includes AGF conditions such as bottom-up freezing,temperature gradients to reach great freezing velocities,and a distinction between two freezing states.One is the freezing-up state until a certain frost body thickness is reached(thermal transient state),and the other is a holding phase where the frost body thickness is kept constant(thermal quasi-steady state).To test its ability,the results are applied to another freezing direction,the top-down freezing.The new approach is validated using two different frost-susceptible soils and,in total,50 frost heave tests.In the thermal transient region,where the SP is applicable,the two semi-analytical approaches are compared,showing improved performance of the current method by about 15%.展开更多
振动污染在国际上已被列入“七大环境公害”之一,填充沟是振动污染治理的一种常用隔振屏障。传统的填充沟受工程造价和施工工艺等限制工程中应用较困难;为克服上述问题,提出一种新型隔振屏障-内插预制芯墙的劲芯水泥土墙。对于饱和地基...振动污染在国际上已被列入“七大环境公害”之一,填充沟是振动污染治理的一种常用隔振屏障。传统的填充沟受工程造价和施工工艺等限制工程中应用较困难;为克服上述问题,提出一种新型隔振屏障-内插预制芯墙的劲芯水泥土墙。对于饱和地基中劲芯水泥土墙的隔振问题,饱和地基采用饱和土半解析边界元法(boundary element method,BEM)建模,劲芯水泥土墙采用有限元法(finite element method,FEM)建模,根据饱和地基-水泥土墙交界面的平衡和连续性条件,分别建立了劲芯水泥土墙对入射Rayleigh波(瑞利波)的远场被动隔振和对动力机器基础振动近场主动隔振的半解析BEM-FEM耦合分析方法,并对劲芯水泥土墙的隔振效果进行了计算分析。研究结果表明:饱和地基中设置劲芯水泥土墙能够起到较好的隔振作用,其隔振效果与混凝土墙基本相当且远优于纯水泥土墙;预制芯墙在劲芯水泥土墙隔振系统中起关键作用,为避免弹性波从芯墙底绕射而降低隔振效果,预制芯墙深度应与水泥土墙深度保持一致;增大芯墙厚度对隔振效果提升不大。在等深芯墙条件下,增大墙深能显著提高远场隔振效果,但对近场隔振效果提升较小;增大墙厚也可提高隔振效果,但提高幅度不大。实际工程中,建议结合工程造价和施工工艺,选择合适的水泥土墙和预制芯墙墙厚。此外,屏障距振源距离对主动隔振效果影响较小,建议根据被保护建筑与振源之间的实际情况,选择合适位置构建隔振屏障。展开更多
To design a propeller for ship power plant,the interaction between ship hull and propeller must be taken into account.The main concern is to apply the wake effect of ship stern on the propeller performance.In this pap...To design a propeller for ship power plant,the interaction between ship hull and propeller must be taken into account.The main concern is to apply the wake effect of ship stern on the propeller performance.In this paper,a coupled BEM(Boundary Element Method)/RANS(Renolds-Averaged Navier−Stokes)solver is used to simulate propeller behind the hull in the self-propulsion test.The motivation of this work is to develop a practical tool to design marine propulsion system without suffering long computational time.An unsteady boundary element method which is also known as panel method is chosen to estimate the propeller forces.Propeller wakes are treated using a time marching wake alignment method.Also,a RANS code coupled with VoF equation is developed to consider the ship motions and wake field effects in the problem.A coupling algorithm is developed to interchange ship wake field to the potential flow solver and propeller thrust to the RANS code.Based on the difference between hull resistance and the propeller thrust,a PI controller is developed to compute the propeller RPM in every time step.Verification of the solver is carried out using the towing tank test report of a 50 m oceanography research vessel.Wake factor and trust deduction coefficient are estimated numerically.Also,the wake rollup pattern of the propeller in open water is compared with the propeller in real wake field.展开更多
The semi-analytical solutions to Fredlund and Hasan＇s one-dimensional （1D） consolidation for unsaturated soils with a semi-permeable drainage boundary are pre- seated. Two variables are introduced to transform the ...The semi-analytical solutions to Fredlund and Hasan＇s one-dimensional （1D） consolidation for unsaturated soils with a semi-permeable drainage boundary are pre- seated. Two variables are introduced to transform the two coupled governing equations of pore-water and pore-air pressures into an equivalent set of partial differential equations （PDFs）, which are easily solved by the Laplace transform method. Then, the pore-water pressure, pore-air pressure, and soil settlement are obtained in the Laplace domain. The Crump method is adopted to perform the inverse Laplace transform in order to obtain the semi-analytical solutions in the time domain. It is shown that the proposed solutions are more applicable to various types of boundary conditions and agree well with the existing solutions from the literature. Several numerical examples are provided to investigate the consolidation behavior of an unsaturated single-layer soil with single, double, mixed, and semi-permeable drainage boundaries. The changes in the pore-air and pore-water pres- sures and the soil settlement with the time factor at different values of the semi-permeable drainage boundary parameters are illustrated. In addition, parametric studies are con- ducted on the pore-air and pore-water pressures at different ratios （the air permeability coefficient to the water permeability coefficient） and depths.展开更多
This paper presents general semi-analytical solutions to Fredlund and Hasan＇s one-dimensional （1D） consolidation equations for unsaturated soils subject to different initial conditions, homogeneous boundaries and t...This paper presents general semi-analytical solutions to Fredlund and Hasan＇s one-dimensional （1D） consolidation equations for unsaturated soils subject to different initial conditions, homogeneous boundaries and time-dependent loadings. Two variables are introduced to transform the two-coupled governing equations of pore-water and poreair pressures into an equivalent set of partial differential equations （PDEs）, which are solved with the Laplace transform method. The pore-water and pore-air pressures and settlement are obtained in the Laplace transform domMn. The Crump＇s method is used to perform inverse Laplace transform to obtain the solutions in the time domain. The present solutions are more general in practical applications and show good agreement with the previous solutions in the literature.展开更多
Based on the Hellinger-Reissner （H-R） mixed variational principle for piezoelectric material, a unified 4-node Hamiltonian isoparametric element of anisotropy piezoelectric material is established. A new semi-analyt...Based on the Hellinger-Reissner （H-R） mixed variational principle for piezoelectric material, a unified 4-node Hamiltonian isoparametric element of anisotropy piezoelectric material is established. A new semi-analytical solution for the natural vibration of smart laminated plates and the transient response of the laminated cantilever with piezoelectric patch is presented. The major steps of mathematical model are as follows： the piezoelectric layer and host layer of laminated plate are considered as unattached three-dimensional bodies and discretized by the Hamiltonian isoparametric elements. The control equation of whole structure is derived by considering the compatibility of generalized displacements and generalized stresses on the interface between layers. There is no restriction for the side-face geometrical boundaries, the thickness and the number of layers of plate by the use of the present isoparametric element. Present method has wide application area.展开更多
This study focuses on the consolidation behavior and mathematical interpretation of partially-saturated ground improved by impervious column inclusion.The constitutive relations for soil skeleton,pore air and pore wat...This study focuses on the consolidation behavior and mathematical interpretation of partially-saturated ground improved by impervious column inclusion.The constitutive relations for soil skeleton,pore air and pore water for partially saturated soils are proposed in the context of partially-saturated ground improved by impervious column inclusion.Settlement equation and dissipation equations of excess pore air/water pressures for a partially saturated improved ground are then derived.The semi-analytical solutions for ground settlement and pore pressure dissipation are then obtained through the Laplace transform and validated by the existing solutions for two special cases in the literature and the numerical results obtained from the finite difference method.A series of parametric studies is finally conducted to investigate the influence of some key factors on consolidation of partially saturated ground improved by impervious column inclusion.Based on the parametric study,it can be found that a higher value of the area replacement ratio or modulus of the pile results in a longer dissipation time of excess pore air pressure(PAP),a shorter dissipation time of excess pore water pressure(PWP),and a lower normalized settlement.展开更多
基金Projects(52378411,52208404)supported by the National Natural Science Foundation of China。
文摘The tunnel subjected to strike-slip fault dislocation exhibits severe and catastrophic damage.The existing analysis models frequently assume uniform fault displacement and fixed fault plane position.In contrast,post-earthquake observations indicate that the displacement near the fault zone is typically nonuniform,and the fault plane position is uncertain.In this study,we first established a series of improved governing equations to analyze the mechanical response of tunnels under strike-slip fault dislocation.The proposed methodology incorporated key factors such as nonuniform fault displacement and uncertain fault plane position into the governing equations,thereby significantly enhancing the applicability range and accuracy of the model.In contrast to previous analytical models,the maximum computational error has decreased from 57.1%to 1.1%.Subsequently,we conducted a rigorous validation of the proposed methodology by undertaking a comparative analysis with a 3D finite element numerical model,and the results from both approaches exhibited a high degree of qualitative and quantitative agreement with a maximum error of 9.9%.Finally,the proposed methodology was utilized to perform a parametric analysis to explore the effects of various parameters,such as fault displacement,fault zone width,fault zone strength,the ratio of maximum fault displacement of the hanging wall to the footwall,and fault plane position,on the response of tunnels subjected to strike-slip fault dislocation.The findings indicate a progressive increase in the peak internal forces of the tunnel with the rise in fault displacement and fault zone strength.Conversely,an augmentation in fault zone width is found to contribute to a decrease in the peak internal forces.For example,for a fault zone width of 10 m,the peak values of bending moment,shear force,and axial force are approximately 46.9%,102.4%,and 28.7% higher,respectively,compared to those observed for a fault zone width of 50 m.Furthermore,the position of the peak internal forces is influenced by variations in the ratio of maximum fault displacement of the hanging wall to footwall and the fault plane location,while the peak values of shear force and axial force always align with the fault plane.The maximum peak internal forces are observed when the footwall exclusively bears the entirety of the fault displacement,corresponding to a ratio of 0:1.The peak values of bending moment,shear force,and axial force for the ratio of 0:1 amount to approximately 123.8%,148.6%,and 111.1% of those for the ratio of 0.5:0.5,respectively.
文摘准确计算户内变电站大型、复杂的噪声场分布,进而评价可采用降噪措施的减噪效果,是解决户内变电站噪声污染的关键问题。为此,综合声学有限元法(finite element method,FEM)求解复杂声场收敛性好及精度高的优点,及声学边界元法(boundary element method,BEM)降维求解大型声场的优势,提出了一种基于声学FEM-BEM的户内变电站噪声场求解算法。首先,建立变电站内部声源声固耦合模型,采用声学FEM求解混响噪声作用下的声固耦合响应;然后,基于声学FEM-BEM耦合理论,求解内、外耦合边界处结构单元受声固耦合激励产生的位移及应力载荷;最后,根据声压及应力载荷激发的外场声波扩散模型,基于常规Gauss数值积分法,建立外部空间声域2维BEM声学积分方程,求解外部声场。该算法在湖南某110 kV户内变电站噪声场的求解分析中得到了成功应用,与实测值的相对误差为3.61%~4.87%。
基金The authors gratefully acknowledge the financial support by the Key Project of High-speed Rail Joint Fund of National Natural Science Foundation of China(Grant No.U1934210)the Natural Science Foundation of Beijing,China(Grant No.8202037).
文摘Longitudinal cracks on the tunnel lining significantly influence the performance of tunnels in operation.In this study,we propose a semi-analytical method that provides a simple and effective way to calculate the internal forces of tunnel linings with multiple cracks.The semi-analytical solution is obtained using structural analysis considering the flexural rigidity for the cracked longitudinal section of the tunnel lining.Then the proposed solution is verified numerically.Using the proposed method,the influences of the crack depth and the number of cracks on the bending moment and modified crack tip stress are investigated.With the increase in crack depth,the bending moment of lining scetion adjacent to the crack decreases,while the bending moment of lining scetion far away from the crack increases slightly.The more the number of cracks in a tunnel lining,the easier the new cracks initiated.
基金supported by the German Research Foundation(DFG)under the project“Investigation and calculation of frost heave considering specific boundary conditions of ground freezing”(Grant No.409760547).
文摘The calculation of frost heaving with ice lens formation is still not standard for construction projects using artificial ground freezing(AGF).In fine-grained material,ice lenses may initiate and lead to significant heaving at the ground surface,which should be considered in advance.However,the complex processes during ice lens formation are still not fully understood and difficult to capture in a simple approach.In the past,the semi-analytical approach of Konrad and Morgenstern used one soil constant,the“segregation potential(SP)”.It has been mainly and most successfully applied to the heave calculation of natural-induced soil freezing in cold regions.Its application to AGF has been so far unsuccessful.To solve this,a new semi-analytical approach is presented in this paper.It includes AGF conditions such as bottom-up freezing,temperature gradients to reach great freezing velocities,and a distinction between two freezing states.One is the freezing-up state until a certain frost body thickness is reached(thermal transient state),and the other is a holding phase where the frost body thickness is kept constant(thermal quasi-steady state).To test its ability,the results are applied to another freezing direction,the top-down freezing.The new approach is validated using two different frost-susceptible soils and,in total,50 frost heave tests.In the thermal transient region,where the SP is applicable,the two semi-analytical approaches are compared,showing improved performance of the current method by about 15%.
文摘振动污染在国际上已被列入“七大环境公害”之一,填充沟是振动污染治理的一种常用隔振屏障。传统的填充沟受工程造价和施工工艺等限制工程中应用较困难;为克服上述问题,提出一种新型隔振屏障-内插预制芯墙的劲芯水泥土墙。对于饱和地基中劲芯水泥土墙的隔振问题,饱和地基采用饱和土半解析边界元法(boundary element method,BEM)建模,劲芯水泥土墙采用有限元法(finite element method,FEM)建模,根据饱和地基-水泥土墙交界面的平衡和连续性条件,分别建立了劲芯水泥土墙对入射Rayleigh波(瑞利波)的远场被动隔振和对动力机器基础振动近场主动隔振的半解析BEM-FEM耦合分析方法,并对劲芯水泥土墙的隔振效果进行了计算分析。研究结果表明:饱和地基中设置劲芯水泥土墙能够起到较好的隔振作用,其隔振效果与混凝土墙基本相当且远优于纯水泥土墙;预制芯墙在劲芯水泥土墙隔振系统中起关键作用,为避免弹性波从芯墙底绕射而降低隔振效果,预制芯墙深度应与水泥土墙深度保持一致;增大芯墙厚度对隔振效果提升不大。在等深芯墙条件下,增大墙深能显著提高远场隔振效果,但对近场隔振效果提升较小;增大墙厚也可提高隔振效果,但提高幅度不大。实际工程中,建议结合工程造价和施工工艺,选择合适的水泥土墙和预制芯墙墙厚。此外,屏障距振源距离对主动隔振效果影响较小,建议根据被保护建筑与振源之间的实际情况,选择合适位置构建隔振屏障。
文摘To design a propeller for ship power plant,the interaction between ship hull and propeller must be taken into account.The main concern is to apply the wake effect of ship stern on the propeller performance.In this paper,a coupled BEM(Boundary Element Method)/RANS(Renolds-Averaged Navier−Stokes)solver is used to simulate propeller behind the hull in the self-propulsion test.The motivation of this work is to develop a practical tool to design marine propulsion system without suffering long computational time.An unsteady boundary element method which is also known as panel method is chosen to estimate the propeller forces.Propeller wakes are treated using a time marching wake alignment method.Also,a RANS code coupled with VoF equation is developed to consider the ship motions and wake field effects in the problem.A coupling algorithm is developed to interchange ship wake field to the potential flow solver and propeller thrust to the RANS code.Based on the difference between hull resistance and the propeller thrust,a PI controller is developed to compute the propeller RPM in every time step.Verification of the solver is carried out using the towing tank test report of a 50 m oceanography research vessel.Wake factor and trust deduction coefficient are estimated numerically.Also,the wake rollup pattern of the propeller in open water is compared with the propeller in real wake field.
基金Project supported by the National Natural Science Foundation of China(Nos.41630633 and11672172)
文摘The semi-analytical solutions to Fredlund and Hasan＇s one-dimensional （1D） consolidation for unsaturated soils with a semi-permeable drainage boundary are pre- seated. Two variables are introduced to transform the two coupled governing equations of pore-water and pore-air pressures into an equivalent set of partial differential equations （PDFs）, which are easily solved by the Laplace transform method. Then, the pore-water pressure, pore-air pressure, and soil settlement are obtained in the Laplace domain. The Crump method is adopted to perform the inverse Laplace transform in order to obtain the semi-analytical solutions in the time domain. It is shown that the proposed solutions are more applicable to various types of boundary conditions and agree well with the existing solutions from the literature. Several numerical examples are provided to investigate the consolidation behavior of an unsaturated single-layer soil with single, double, mixed, and semi-permeable drainage boundaries. The changes in the pore-air and pore-water pres- sures and the soil settlement with the time factor at different values of the semi-permeable drainage boundary parameters are illustrated. In addition, parametric studies are con- ducted on the pore-air and pore-water pressures at different ratios （the air permeability coefficient to the water permeability coefficient） and depths.
基金Project supported by the National Natural Science Foundation of China(Nos.41372279 and41630633)
文摘This paper presents general semi-analytical solutions to Fredlund and Hasan＇s one-dimensional （1D） consolidation equations for unsaturated soils subject to different initial conditions, homogeneous boundaries and time-dependent loadings. Two variables are introduced to transform the two-coupled governing equations of pore-water and poreair pressures into an equivalent set of partial differential equations （PDEs）, which are solved with the Laplace transform method. The pore-water and pore-air pressures and settlement are obtained in the Laplace transform domMn. The Crump＇s method is used to perform inverse Laplace transform to obtain the solutions in the time domain. The present solutions are more general in practical applications and show good agreement with the previous solutions in the literature.
基金Project supported by the National Natural Science Foundation of China (No. 10072038)
文摘Based on the Hellinger-Reissner （H-R） mixed variational principle for piezoelectric material, a unified 4-node Hamiltonian isoparametric element of anisotropy piezoelectric material is established. A new semi-analytical solution for the natural vibration of smart laminated plates and the transient response of the laminated cantilever with piezoelectric patch is presented. The major steps of mathematical model are as follows： the piezoelectric layer and host layer of laminated plate are considered as unattached three-dimensional bodies and discretized by the Hamiltonian isoparametric elements. The control equation of whole structure is derived by considering the compatibility of generalized displacements and generalized stresses on the interface between layers. There is no restriction for the side-face geometrical boundaries, the thickness and the number of layers of plate by the use of the present isoparametric element. Present method has wide application area.
基金The financial support from National Natural Science Foundation of China (Grant Nos. 12172211 and 52078021)Shanghai Key Laboratory of Rail Infrastructure Durability and System Safety, China (Grant No. R201904)
文摘This study focuses on the consolidation behavior and mathematical interpretation of partially-saturated ground improved by impervious column inclusion.The constitutive relations for soil skeleton,pore air and pore water for partially saturated soils are proposed in the context of partially-saturated ground improved by impervious column inclusion.Settlement equation and dissipation equations of excess pore air/water pressures for a partially saturated improved ground are then derived.The semi-analytical solutions for ground settlement and pore pressure dissipation are then obtained through the Laplace transform and validated by the existing solutions for two special cases in the literature and the numerical results obtained from the finite difference method.A series of parametric studies is finally conducted to investigate the influence of some key factors on consolidation of partially saturated ground improved by impervious column inclusion.Based on the parametric study,it can be found that a higher value of the area replacement ratio or modulus of the pile results in a longer dissipation time of excess pore air pressure(PAP),a shorter dissipation time of excess pore water pressure(PWP),and a lower normalized settlement.
