Fracture mechanics analysis of multiple cracks in anisotropic media

This paper presents a single‐domain boundary element method (BEM) for linear elastic fracture mechanics analysis in the two‐dimensional anisotropic material. In this formulation, the displacement integral equation is collocated on the un‐cracked boundary only, and the traction integral equation is collocated on one side of the crack surface only. A special crack‐tip element was introduced to capture exactly the crack‐tip behavior. A computer program with the FORTRAN language has been developed to effectively calculate the stress intensity factors of an anisotropic material. This BEM program has been verified having a good accuracy with the previous researches. Furthermore, by analyzing the different anisotropic degree cracks in a finite plate, we found that the stress intensity factors of crack tips had apparent influence by the geometry forms of cracks and media with different anisotropic degrees. Copyright © 2010 John Wiley & Sons, Ltd.     

合肥市郊低层大气的激光雷达探测研究

利用L300米散射激光雷达对合肥市郊大气边界层进行长期系统观测.分析讨论了大气边界层气溶胶消光系数与温度、湿度的关系, 大气边界层气溶胶消光系数垂直分布和时间变化的主要特征, 给出了激光雷达探测的大气边界层高度的统计特征及其与无线电气象探空仪探测大气边界层高度的比较结果.     

Application of an improved semi-Lagrangian procedure to fully-nonlinear simulation of sloshing in non-wall-sided tanks

The semi-Lagrangian procedure is widely used for updating the fully-nonlinear free surface in the time domain. However, this procedure is only available to cases when the body surface is vertical near the waterline. Present study introduces an improved semi-Lagrangian procedure which removes this ‘vertical-wall’ limitation. Coupling with the boundary element method, the improved semi-Lagrangian procedure is applied to the simulation of fully-nonlinear sloshing waves in non-wall-sided tanks. From the result comparison with the open source CFD software OpenFOAM, it is confirmed that this numerical scheme could guarantee a sufficient accuracy. Further series studies on 2D and 3D fully-nonlinear sloshing waves in wedged tanks are performed. Featured phenomena are observed which are distinct from those in wall-sided tanks.     

Multiple controls on the paleoenvironment of the Early Cambrian marine black shales in the Sichuan Basin,SW China: Geochemical and organic carbon isotopic evidence

In order to understand the paleoenvironment of the Early Cambrian black shale deposition in the western part of the Yangtze Block, geochemical and organic carbon isotopic studies have been performed on two wells that have drilled through the Qiongzhusi Formation in the central and southeastern parts of Sichuan Basin. It shows that the lowest part of the Qiongzhusi Formation has high TOC abundance, while the middle and upper parts display relative low TOC content. Redox-sensitive element (Mo) and trace elemental redox indices (e.g., Ni/Co, V/Cr, U/Th and V/(V + Ni)) suggest that the high-TOC layers were deposited under anoxic conditions, whereas the low-TOC layers under relatively dysoxic/oxic conditions. The relationship of the enrichment factors of Mo and U further shows a transition from suboxic low-TOC layers to euxinic high-TOC layers. On the basis of the Mo-TOC relationship, the Qiongzhusi Formation black shales were deposited in a basin under moderately restricted conditions. Organic carbon isotopes display temporal variations in the Qiongzhusi Formation, with a positive excursion of δ¹³C_org values in the lower part and a continuous positive shift in the middle and upper parts. All these geochemical and isotopic criteria indicate a paleoenvironmental change from bottom anoxic to middle and upper dysoxic/oxic conditions for the Qiongzhusi Formation black shales. The correlation of organic carbon isotopic data for the Lower Cambrian black shales in different regions of the Yangtze Block shows consistent positive excursion of δ¹³C_org values in the lower part for each section. This excursion can be ascribed to the widespread Early Cambrian transgression in the Yangtze Block, under which black shales were deposited.     

水平轴潮流能叶轮尖速比特性分析研究

针对设计尖速比对水平轴潮流能叶轮动力特性的影响问题开展研究工作。基于叶素-动量理论建立叶轮动力特性仿真模型,以叶片数量分别为2、3、4的叶轮为对象,考察设计尖速比对叶轮的功率系数的影响。研究结果表明:设计尖速比越大的叶轮,其最大效率也越高;叶片数量多的叶轮,其功率系数略优于叶片数量少的叶片。结合分析结果,给出了水平轴潮流能叶轮的设计尖速比的选择建议。     

Corrected Formula of Bed Resistance Coefficient for Plane Numerical Simulation of Tidal Current

By means of a logarithm law for the velocity profile,a corrected formula of bed resistancecoefficient,which involves many factors such as gradient of still water depth,variation of surfaceelevation,flow direction,and so on,is derived from the 3D governing equations of tidal current by aver-aging over the whole water depth.Theoretical analysis and application have shown that the 2D plane tidalcurrent numerical model would be more reasonable and could be applied to steep bottom topographywhen the corrected bed resistance coefficient is used,therefore the results of reproduction simulation andengineering calculation would be more scientific and reasonable.     

椭圆度-凹坑双缺陷海底管道局部屈曲特性研究

海底管道在制造、埋设以及使用过程中极易产生椭圆度-凹坑双缺陷,双缺陷影响管道局部屈曲,对含椭圆度-凹坑双缺陷海底管道的局部屈曲特性研究十分必要。现行规范中采用等效椭圆度对含椭圆度-凹坑海底管道进行评估,该方法无法准确评估不同缺陷形式的屈曲特性。采用形状系数对含椭圆度-凹坑双缺陷的海底管道进行评估,运用有限元软件ABAQUS进行数值模拟,并进行试验验证。在此基础上对含有不同凹坑深度、不同椭圆度的海底管道进行局部屈曲的数值模拟,计算不同形状椭圆度、含有不同凹坑深度海底管道的形状系数,对其进行敏感性分析。计算结果表明:形状系数对海底管道椭圆度、凹坑深度、径厚比敏感性较强;对凹坑宽度敏感性较弱。     

颗粒材料数值样本的坐标排序生成技术

颗粒材料离散颗粒模型的数值模拟结果与颗粒材料的数值样本密切相关,随着离散单元在颗粒材料数值模拟领域的广泛应用,颗粒材料的数值样本生成技术日益受到重视。基于RSA模型研究如何使随机生成的颗粒材料更密实,对均匀颗粒而言亦即如何在指定区域内生成更多的颗粒,讨论了4类修正方案,并建议了一种基于坐标排序的样本生成技术。研究表明,在传统的颗粒体随机生成技术基础上,通过对随机生成的x坐标序列或y坐标序列进行排序,可使生成的颗粒材料数值样本更密实。     

An analysis of horizontal dispersion due to the drift current with an Ekman layer

An analytical method for describing horizontal matter dispersion in shear currents is presented using a tensor expression from the point of view that matter dispersion due to the shear effect should be one of the principal mixing dilution processes. Although the behavior of horizontal dispersion is considerably more complicated than common longitudinal dispersion, the present study elucidates the vertical structure of dispersion and the dispersing process from the initial to the stationary stage, besides the usual depth-averaged dispersion coefficient at the stationary stage. As one of the typical applications of horizontal dispersion, dispersion due to the pure drift current with an Ekman layer is examined theoretically using the present method. This examination reveals that the displacement of the centroid and the major axis of dispersion are twisted in the vertical direction more than the direction of the current vector forming the Ekman spiral; that the variance increases in proportion to the third power of the elapsed time; and that the dispersion coefficient at the stationary stage remains constant, independent of the depth normalized by an Ekman layer thickness. Such dependence of the dispersion coefficient in the steady current is shown to be different from that in the oscillatory current, which is inversely proportional to the depth normalized by a Stokes layer thickness. This is considered to be induced by the difference of the vertical profiles of the first order moment in both currents, that is, the shear region of the first order moment is restricted around the floor by the alternation of the current shear in the oscillatory current while it is diffused in the whole depth in the steady current.     

Numerical Modeling on Suspended Sediment Transportation in the Hangzhou Bay

In this paper the characteristics of tidal flow and seasonal variation of seidment content in theHangzhou Bay and their affecting factors are studied.Field investigations and data analysis indicate thatthe sediment movement is mainly influenced by the Yangtze estuary and the sediment of the Yangtze estua-ry is induced by wind wave and tidal flow.Owing to the variation of dynamic conditions,the instanta-neous sediment content is controlled by tidal flow,wind wave,depth of water and tidal range synthetically.A sediment content relationship formula is established with related factors.A non-equilibrium2-dimensional numerical model of suspended sediment transportation is set up,and the finite element meth-od is applied.The computation results of the model is in accordance with field data.