钢筋钢纤维混凝土梁柱节点模型化方法研究

为研究钢筋钢纤维混凝土梁柱节点抗震性能数值模拟方法,基于OpenSEES有限元平台中的梁柱节点单元(beam-column joint element, BCJE)模型,通过修正模型中剪切块和钢筋滑移弹簧的参数修正方法,提出适用于钢筋钢纤维混凝土梁柱节点的数值模型,并基于6个梁柱节点的拟静力试验结果进行模型验证分析,数值模拟结果与试验结果吻合较好,提出的数值模型能够较精确地反映节点滞回行为。在此基础上,分析了轴压比、钢纤维体积率和配箍量对梁柱节点抗震性能的影响规律,建立了节点受剪承载力的计算公式。结果表明:掺入钢纤维和增加配箍可明显改善梁柱节点的抗震性能,钢纤维体积率从0.5%增加到2.0%,极限荷载提高了18%;箍筋从1Φ8增加到3Φ8时,极限荷载提高了19.7%。     

A high‐order immersed boundary discontinuous‐Galerkin method for Poisson's equation with discontinuous coefficients and singular sources

We adopt a numerical method to solve Poisson's equation on a fixed grid with embedded boundary conditions, where we put a special focus on the accurate representation of the normal gradient on the boundary. The lack of accuracy in the gradient evaluation on the boundary is a common issue with low‐order embedded boundary methods. Whereas a direct evaluation of the gradient is preferable, one typically uses post‐processing techniques to improve the quality of the gradient. Here, we adopt a new method based on the discontinuous‐Galerkin (DG) finite element method, inspired by the recent work of [A.J. Lew and G.C. Buscaglia. A discontinuous‐Galerkin‐based immersed boundary method. International Journal for Numerical Methods in Engineering, 76:427‐454, 2008]. The method has been enhanced in two aspects: firstly, we approximate the boundary shape locally by higher‐order geometric primitives. Secondly, we employ higher‐order shape functions within intersected elements. These are derived for the various geometric features of the boundary based on analytical solutions of the underlying partial differential equation. The development includes three basic geometric features in two dimensions for the solution of Poisson's equation: a straight boundary, a circular boundary, and a boundary with a discontinuity. We demonstrate the performance of the method via analytical benchmark examples with a smooth circular boundary as well as in the presence of a singularity due to a re‐entrant corner. Results are compared to a low‐order extended finite element method as well as the DG method of [1]. We report improved accuracy of the gradient on the boundary by one order of magnitude, as well as improved convergence rates in the presence of a singular source. In principle, the method can be extended to three dimensions, more complicated boundary shapes, and other partial differential equations. Copyright © 2014 John Wiley & Sons, Ltd.     

Robin‐Neumann transmission conditions for fluid‐structure coupling: Embedded boundary implementation and parameter analysis

Partitioned procedures are appealing for solving complex fluid‐structure interaction (FSI) problems, as they allow existing computational fluid dynamics (CFD) and computational structural dynamics algorithms and solvers to be combined and reused. However, for problems involving incompressible flow and strong added‐mass effect (eg, heavy fluid and slender structure), partitioned procedures suffer from numerical instability, which typically requires additional subiterations between the fluid and structural solvers, hence significantly increasing the computational cost. This paper investigates the use of Robin‐Neumann transmission conditions to mitigate the above instability issue. Firstly, an embedded Robin boundary method is presented in the context of projection‐based incompressible CFD and finite element–based computational structural dynamics. The method utilizes operator splitting and a modified ghost fluid method to enforce the Robin transmission condition on fluid‐structure interfaces embedded in structured non–body‐conforming CFD grids. The method is demonstrated and verified using the Turek and Hron benchmark problem, which involves a slender beam undergoing large transient deformation in an unsteady vortex‐dominated channel flow. Secondly, this paper investigates the effect of the combination parameter in the Robin transmission condition, ie, α_f, on numerical stability and solution accuracy. This paper presents a numerical study using the Turek and Hron benchmark problem and an analytical study using a simplified FSI model featuring an Euler‐Bernoulli beam interacting with a two‐dimensional incompressible inviscid flow. Both studies reveal a trade‐off between stability and accuracy: smaller values of α_f tend to improve numerical stability, yet deteriorate the accuracy of the partitioned solution. Using the simplified FSI model, the critical value of α_f that optimizes this trade‐off is derived and discussed.