板钉连接CFST柱-RC梁节点梁端塑性铰区受力性能数值模拟

建立竖板-栓钉连接钢管混凝土（CFST）柱-钢筋混凝土（RC）梁节点试件（SSJD）拟静力加载试验有限元模型,并在节点损伤情况、梁端荷载-位移曲线等数值模拟结果与试验结果吻合较好的基础上,进一步开展了RC梁混凝土强度、配筋率ρs和连接竖板长度Lb及界面连接情况等对CFST柱-RC梁节点梁端塑性铰区域力学性能的影响。研究结果表明,RC梁混凝土强度对试件SSJD塑性铰区域受力性能的影响较小;适筋范围内RC梁配筋率增加可适当提高试件SSJD承载力和延性;随着连接竖板长度的增加,梁端塑性铰区域外移,梁破坏荷载增大;本研究给出的RC梁与CFST柱之间的界面抗剪承载力模拟值与计算值吻合较好,可用于界面抗剪设计。     

Seismic performance of two-story subway station structures with different isolating systems

When the sliding bearing is fixed only at the top of the middle column of the underground structure, the cracks at the side end of the middle plate should be aggravated while the seismic damage of the mid-column should be alleviated. To enhance the seismic performance of the mid-plate, a new isolation design method has been mentioned while the elastic sliding bearings are set at the top of the mid-columns and between the side end of the mid-plate and the side wall at the same time. By establishing a nonlinear finite element analysis model for the static-dynamic coupling interaction system, the seismic response characteristics of the cast-in-place station structure without a sliding bearing have been analyzed and compared with those of the station structure with the sliding bearing fixed only at the top of the middle columns, and those of the station structure with sliding bearing be fixed between the mid-plate and the sidewall at the same time. The results show that the new isolation station structures suffer fewer earthquake damages at the mid-plate and mid-columns at the same time, which can improve the overall seismic performance of the subway station structure.     

Numerical determination of pressure-dependent effective thermal conductivity in Berea sandstone

In this paper, the methods of digital rock physics are applied to determine pressure-dependent effective thermal conductivity in rock samples. Simulations are performed with an in-house three-dimensional finite volume code. In the first step, four numerical models are derived from a given tomographic scan of Berea sandstone. Consequently, simulations of the thermal conductivity at ambient conditions are performed and validated with experimental data. In a second step, a new workflow for the determination of the pressure-dependent thermal conductivity in rock samples is elaborated, tested and calibrated. Results originating from the derived workflow show very good agreement with experimental data.     

非均匀弹性介质地震波动模拟的Fourier变换方法

一、前言 地震模型技术是应用地震法研究地下储油地层分布的基础工作,它可以用来研究地,震波的传播以及各种波的特性,可以通过修改地质模型,反复进行正演,实现对实际地震剖面的最优逼近,以实现地震反问题。本文在文献[1]的基础上,应用Fourier变换方法研究了用非均匀弹性介质波动方程正演问题的解来模拟地震波的传播,探索了该方法对非均匀弹性介质地震模型的有效性及实用性。     

马17井映震灵敏的应力条件数学模拟研究

本文利用线弹性理论的二维有限元方法,计算了马17井的区域应力及井区平面与剖面的应力状态,研究该井映震灵敏的应力条件。计算结果表明井区在华北地区中强以上地震活动时具有一定的应力集中     

Equilibrium thallium isotope fractionation and its constraint on Earth’s late veneer

Equilibrium isotope fractionation of thallium(Tl) includes the traditional mass-dependent isotope fractionation effect and the nuclear volume effect(NVE). The NVE dominates the overall isotope fractionation, especially at high temperatures. Heavy Tl isotopes tend to be enriched in oxidized Tl^3+-bearing species. Our NVE fractionation results of oxidizing Tl^+ to Tl^3+ can explain the positive enrichments observed in ferromanganese sediments. Experimental results indicate that there could be0.2–0.3 e-unit fractionation between sulfides and silicates at 1650 ℃. It is consistent with our calculation results,which are in the range of 0.17–0.38 e-unit. Importantly,Tl’s concentration in the bulk silicate Earth(BSE) can be used to constrain the amount of materials delivered to Earth during the late veneer accretion stage. Because the Tl concentration in BSE is very low and its Tl isotope composition is similar with that of chondrites, suggesting either no Tl isotope fractionation occurred during numerous evaporation events, or the Tl in current BSE was totally delivered by late veneer. If it is the latter, the Tl-contentbased estimation could challenge the magnitude of late veneer which had been constrained by the amount of highly siderophile elements in BSE. Our results show that the lateaccreted mass is at least five-times larger than the previously suggested magnitude, i.e., 0.5 wt% of current Earth’s mass. The slightly lighter 205 Tl composition of BSE relative to chondrites is probable a sign of occurrence of Tlbearing sulfides, which probably were removed from the mantle in the last accretion stage of the Earth.     

准各向异性粘弹介质地震波的数字仿真

本文论述准各向异性粘弹波动方程有限元数值解地震波响应的原理与计算方法,并给出地震衰减因子与方程中粘滞系数的计算关系式。列出了所研制的地面记录地震剖面与井中记录垂直地震剖面,以及波场时间切片图、散度与旋度图、应变能与动能的时间切片图。本方法适于任意分层形态且层间物性参数为强间断的实际地震剖面的仿真,旨在研究正问题以实现反问题。     

A survey of numerical methods for stochastic differential equations

The development of numerical methods for stochastic differential equations has intensified over the past decade. The earliest methods were usually heuristic adaptations of deterministic methods, but were found to have limited accuracy regardless of the order of the original scheme. A stochastic counterpart of the Taylor formula now provides a framework for the systematic investigation of numerical methods for stochastic differential equations. It suggests numerical schemes, which involve multiple stochastic integrals, of higher order of convergence. We shall survey the literature on these and on the earlier schemes in this paper. Our discussion will focus on diffusion processes, but we shall also indicate the extensions needed to handle processes with jump components. In particular, we shall classify the schemes according to strong or weak convergence criteria, depending on whether the approximation of the sample paths or of the probability distribution is of main interest.     

Lava flow hazard at Nyiragongo volcano,D.R.C.

The 2002 eruption of Nyiragongo volcano constitutes the most outstanding case ever of lava flow in a big town. It also represents one of the very rare cases of direct casualties from lava flows, which had high velocities of up to tens of kilometer per hour. As in the 1977 eruption, which is the only other eccentric eruption of the volcano in more than 100 years, lava flows were emitted from several vents along a N–S system of fractures extending for more than 10 km, from which they propagated mostly towards Lake Kivu and Goma, a town of about 500,000 inhabitants. We assessed the lava flow hazard on the entire volcano and in the towns of Goma (D.R.C.) and Gisenyi (Rwanda) through numerical simulations of probable lava flow paths. Lava flow paths are computed based on the steepest descent principle, modified by stochastically perturbing the topography to take into account the capability of lava flows to override topographic obstacles, fill topographic depressions, and spread over the topography. Code calibration and the definition of the expected lava flow length and vent opening probability distributions were done based on the 1977 and 2002 eruptions. The final lava flow hazard map shows that the eastern sector of Goma devastated in 2002 represents the area of highest hazard on the flanks of the volcano. The second highest hazard sector in Goma is the area of propagation of the western lava flow in 2002. The town of Gisenyi is subject to moderate to high hazard due to its proximity to the alignment of fractures active in 1977 and 2002. In a companion paper (Chirico et al., Bull Volcanol, in this issue, 2008) we use numerical simulations to investigate the possibility of reducing lava flow hazard through the construction of protective barriers, and formulate a proposal for the future development of the town of Goma.     

结构动力方程的显式级数积分格式

将Newmark-β法中常平均加速度法的基本假定与精细指数算法结合,根据指数矩阵的Taylor级数展开式,提出了动力方程的显式级数解,并设计了相应的时程积分算法.该算法的精度可根据Taylor级数展开式的项数进行灵活控制.算例的结果表明:在满足稳定性条件的前提下,随着时间步长的增加,其精度优于传统的时程积分法.通过稳定性的分析,指出其稳定性条件是显然满足的.