Dynamic FE simulation of four‐story steel frame modeled by solid elements and its validation using results of full‐scale shake‐table test

Dynamic finite element analyses of a four‐story steel building frame modeled as a fine mesh of solid elements are performed using E‐Simulator, which is a parallel finite element analysis software package for precisely simulating collapse behaviors of civil and building structures. E‐Simulator is under development at the National Research Institute for Earth Science and Disaster Prevention (NIED), Japan. A full‐scale shake‐table test for a four‐story frame was conducted using E‐Defense at NIED, which is the largest shaking table in the world. A mesh of the entire structure of a four‐story frame with approximately 19 million degrees of freedom is constructed using solid elements. The density of the mesh is determined by referring to the results of elastic–plastic buckling analyses of a column of the frame using meshes of different densities. Therefore, the analysis model of the frame is well verified. Seismic response analyses under 60, 100, and 115% excitations of the JR Takatori record of the 1995 Hyogoken‐Nanbu earthquake are performed. Note that the simulation does not reproduce the collapse under the 100% excitation of the Takatori record in the E‐Defense test. Therefore, simulations for the 115% case are also performed. The results obtained by E‐Simulator are compared with those obtained by the E‐Defense full‐scale test in order to validate the results obtained by E‐Simulator. The shear forces and interstory drift angles of the first story obtained by the simulation and the test are in good agreement. Both the response of the entire frame and the local deformation as a result of elastic–plastic buckling are simulated simultaneously using E‐Simulator. Copyright © 2014 John Wiley & Sons, Ltd.     

Finite‐element analysis of laminated rubber bearing of building frame under seismic excitation

Finite element analysis is carried out for a building frame supported by laminated rubber bearings to simultaneously investigate global displacement and local stress responses under seismic excitation. The frame members and the rubber bearings are discretized into hexahedral solid elements with more than 3 million degrees of freedom. The material property of rubber is represented by the Ogden model, and the frame is assumed to remain in elastic range. It is shown that the time histories of non‐uniform stress distribution and rocking behavior of the rubber bearings under a frame subjected to seismic excitation can be successfully evaluated, and detailed responses of base and frame can be evaluated through large‐scale finite element analysis. Copyright © 2015 John Wiley & Sons, Ltd.     

应力扰动对地震周期和地震矩影响的数值模拟

利用二维有限元数值模型,结合断层滑移弱化摩擦准则对断层滑动规律以及应力扰动对其影响进行了研究．数值计算结果表明,在均匀应力分布情况下, 平面断层滑动显示出典型的特征地震规律,断层面上的应力扰动对断层滑动规律产生影响,压应力增加明显延迟地震的发生时间,并增加地震释放的能量．应力扰动发生在地震破裂临界区时的影响比在震前滑移区时的影响显著．当发生在地震滑移区时,若应力扰动足够大,则压应力增大会造成地震发生时部分动力断层被暂时锁住,使得地震释放的能量变小,但可增加后续地震的能量; 而压应力减小则可导致地震规律产生更加复杂的变化,会即时触发地震．如果应力扰动发生在一个地震周期的早期,则触发的地震较小,但可导致随后的地震提前发生; 如果应力扰动发生在一个地震周期的后期,则会触发大地震．当应力扰动位于震前滑移区或破裂临界区时,小的扰动也可能产生类似的效果．应力扰动产生越晚,这种影响也越明显．应力扰动发生在破裂临界区的影响最明显．应力扰动的影响一般主要集中在应力发生扰动后的1—2个地震周期内．后续地震基本恢复无应力扰动时的特征地震规律．      

three-dimensional sensitivity coefficients of apparent resistivity and preliminary application

We calculate three-dimensional sensitivity coefficients distribution of apparent resistivity observation when Schlumberger array is used by using finite element method. Analysis results suggest that for the situation of one-dimensional positive or minus coefficient of surface medium, three-dimensional sensitivity coefficients distribution at surface shows similar patterns, and sensitivity coefficients distributions of different layered electric structures are also similar. There are two approximate ellipses at the two-dimensional surface plane between current electrodes and potential electrodes, where sensitivity coefficients are minus, and sensitivity coefficients at other areas are positive. Sensitivity coefficients at two approximate ellipses between current electrodes and potential electrodes are minus at the vertical section along monitoring line, while others are positive. From the three-dimensional view, minus sensitivity coefficients are at the two approximate half ellipsoids between current electrodes and potential electrodes when arrays are applied at surface. And coefficients near the electrodes are much greater than other areas. When resistivity of local areas at surface changes, we can qualitatively analyze the disturbing effects caused by the areas using three dimensional sensitivity coefficients distribution, and the analysis result can serve as reference for further experiment and numerical model quantitative analysis.     

常量元素氧化物含量及其比值揭示的中晚全新世以来玛曲高原的环境演变

对玛曲高原具有代表性的黄土-古土壤-风成砂剖面(OL剖面)常量元素氧化物含量及其比值的地球化学特征分析表明:(1)常量元素氧化物总含量为91.93%且SiO₂(67.93%) >Al₂O₃(11.2%) >Fe₂O₃(3.97%) >CaO(3.76%)>K₂O(2.16%) >Na₂O(1.81%) >MgO(1.1%);(2)玛曲高原处于较弱的脱Ca、Na初级风化阶段,并伴有碳酸盐的淋失;(3)玛曲高原的气候经历了冷干-凉润-温湿-凉干的冷暖旋回。依据常量元素氧化物含量及其比值揭示的环境特征并结合¹⁴C测年将玛曲高原中晚全新世以来的环境演变划分为4.6 ka BP以前的风沙活动强烈期、4.6～4.3 ka BP的风沙活动较弱期、4.3～0.3 ka BP的固结成壤期和0.3 ka BP至今的风沙活动增强期4个阶段。     

基于要素流的长三角城市群空间布局

以长三角城市群30个地级市市区为研究单元,从要素流动的动态视角,建立要素流的数理模型,借助于ArcGIS软件,对长三角城市群内城市的经济流、交通流和信息流的强度和方向进行定量分析和评价,运用赋值法对3种要素流总量进行评分,并进行空间层级划分。得出核心区（包括上海、江苏沿江八市、环杭州湾地区等地市）和辐射影响区（位于长三角两翼的地市）。基于城市间的空间联系特征和城市群的空间层级结构,对长三角城市群的空间发展布局提出进一步引导和规划,即上海为核心,苏州、南京、杭州、合肥和徐州为区域性中心城市,构成以沪宁合、沪杭、沪徐和沪温为4个发展轴的反“K”字型的空间形态,依据城市区位和要素联系特征,将城市群进一步划分为6个城市组团,形成以中心城市为增长节点、以发展轴为联结线、以城市组团为发展域的“点―线―面”有序发展的空间互动格局。     

裂隙岩体三维柔度张量的数值试验方法

在裂隙岩体二维柔度张量数值试验的基础上,建立了裂隙岩体三维柔度张量及其表征单元体积（REV）尺度研究的简化数值试验方法。首先根据岩体裂隙的统计参数及分布规律,运用蒙特卡罗方法在研究域内获得与实际岩体裂隙同分布的三维随机裂隙网络,提取位于不同方位的岩体试件,运用二维柔度张量的数值试验方法求出各个平面方向上的二维柔度张量,然后根据二维与三维柔度张量的拓扑关系,用数学方法求解裂隙岩体的三维柔度张量。对于含3组正交全贯通裂隙的岩体,通过数值试验获得的柔度张量与理论解相比,其主对角线上各参数的误差在5%以内,表明该数值试验方法较为可靠。最后将此方法运用到小湾水电站工程中,确定左岸坝区裂隙岩体的应力REV为11 m×11 m×11 m,并获得该区域裂隙岩体的三维柔度张量。     

粒间胶结接触力学特性的三维试验研究

为实现结构性砂土离散元接触模型合理性的三维试验验证,设计了一套可用于三维半球形理想胶结颗粒成型及实现不同加载条件下的接触力学特性测试装置,制备了一定胶结尺寸的环氧树脂半球形颗粒胶结试样,在一系列辅助加载装置中初步开展了不同加载条件（拉伸、压缩、剪切、弯转、扭转）下的力学性能测试。结果表明,该装置可用于实现三维情况下胶结颗粒接触力学特性测试;不同加载条件下的实测试验结果与二维试验成果基本一致;峰值剪切、弯矩、扭矩随着法向荷载的增大呈现先增大后减小的趋势,存在一个相同的临界法向荷载。     

土体直剪力学特性颗粒尺度效应理论与试验研究

土体是一种颗粒物质,其强度与变形特性具有显著的颗粒尺度效应。根据土体颗粒间的连结性状和微重比,将土颗粒划分为基体颗粒与加强颗粒。构建了反映土体内部材料信息和颗粒特征信息的土体胞元,基于应变梯度理论建立可以描述土体颗粒尺度效应的土体胞元模型。设计一系列饱和重塑土的直接快剪试验以研究土体直剪力学特性的颗粒尺度效应,并定量计算了土体胞元模型的应变梯度和內禀尺度等微细观计算参数。试验结果表明,土体的剪切屈服应力随加强颗粒体积比和平均应变梯度的增加而增加,且与加强颗粒体积比呈近似线性关系,与平均应变梯度呈抛物线关系;加强颗粒粒径对土体的剪切屈服应力影响不明显。土体剪切屈服应力的试验结果与土体胞元模型的预测结果一致。     

用空气单元法求解渗流场的逸出边界问题

依据空气单元模拟排水孔的思路,在渗流场任意可能逸出边界外布置一层渗透性较强的单元,与其他实体单元一起进行渗流计算,可以实现该边界的自动求解。与现行方法相比,不仅省去了逸出边界迭代计算的麻烦,而且避免了由于逸出点定位不准确而可能引起的渗流计算不合理或收敛困难等缺陷。算例分析表明,空气单元法的计算精度与空气单元的相对渗透系数R有关,与厚度L无关。当选取了合适的R值后,可以较好地逼近现行方法。对于因引入空气单元而造成的强弱透水介质界面处数值拟真性差的现象,还提出了加入接触面薄层单元的改进方法。方法改进后,计算精度可进一步提高。最后,以小湾水电站22#坝段渗流场计算为例,说明该法的实用性和可行性。