法向承力锚极限抗拔力特性

法向承力锚是一种新型的适用于深海工程的系泊基础,其极限抗拔力是锚在工程应用中的关键指标。尝试用两种不同的方法评估法向承力锚的极限抗拔力,其一是基于塑性上限分析理论;其二是运用非线性有限元数值方法。与已有的经验公式相比,所建立的计算模型不仅可考虑海床土性质,还能反映锚板定位(嵌入深度及角度)以及系缆力角度对锚极限抗拔力的影响。在与已有评估方法进行比较的基础上,还特别对锚板的嵌入深度、角度以及系缆力角度变化对极限抗拔力的影响规律进行了分析,对三种方法的适用性进行了评述。     

Statistics of Wind and Waves off Tsuyazaki, Fukuoka, in the Eastern Tsushima Strait

The variability of the sea surface wind and wind waves in the coastal area of the Eastern Tsushima Strait was investigated based on the hourly data from 1990 to 1997 obtained at a station 2 km off Tsuyazaki, Fukuoka. The annual mean wind speed was 4.84 m s⁻¹, with strong northwesterly monsoon in winter and weak southwesterly wind in summer. Significant wave heights and wave periods showed similar sinusoidal seasonal cycles around their annual means of 0.608 m and 4.77 s, respectively. The seasonal variability relative to the annual mean is maximum for wave heights, medium for wind speeds, and minimum for wave periods. Significant wave heights off Tsuyazaki turned out to be bounded by a criterion, which is proportional to the square of the significant wave period corresponding to a constant steepness, irrespective of the season or the wind speed. For terms shorter than a month, the significant wave height and the wave period were found to have the same spectral form as the inshore wind velocity: white for frequencies less than 0.2 day⁻¹ and proportional to the frequency to the −5/3 power for higher frequencies, where the latter corresponds to the inertial subrange of turbulence. The spectral levels of wave heights and wave periods in that inertial range were also correlated with those of the inshore wind velocity, though the scatter was large. This revised version was published online in August 2006 with corrections to the Cover Date.     

基于矩阵Toeplitz稀疏分解的相对波阻抗反演方法

针对利用地震道进行相对波阻抗反演中遇到的横向连续性难以保持、初始子波容错度差以及随机噪声干扰影响反演结果等问题,提出了一种基于矩阵Toeplitz稀疏分解的相对波阻抗反演方法.该方法将地震数据剖面的Toeplitz稀疏分解问题分解为两个子反演问题,其一以Toeplitz子波矩阵元素为待反演的参数,用Fused Lasso方法求解,可保证子波具有紧支集且是光滑的;其二以稀疏反射系数矩阵元素为待反演参数,用基于回溯的快速萎缩阈值迭代算法求解,大大降低了目标函数中参数选择的难度.通过交替迭代求解上述两个子反演问题可将地震数据剖面因式分解为一个Toeplitz子波矩阵和一个稀疏反射系数矩阵;然后由反射系数矩阵递推反演可以得到高分辨率的相对波阻抗剖面;利用测井资料加入低频分量后,也可得到高分辨率的绝对波阻抗剖面.Marmousi2模型生成的合成记录算例和实际地震资料算例均表明:本文方法可以从带限地震数据中有效地反演相对波阻抗,反演结果分辨率高并且能够很好地保持地震数据的横向连续性;即使在初始估计子波存在误差和地震数据被随机噪声污染的情况下也能取得较好的效果.     

基于遗传算法的反倾层状岩质边坡极限分析上限解

倾倒破坏是反倾层状边坡主要的失稳类型之一,而当反倾层状边坡内存在软弱夹层、外倾结构面等因素时,传统的纯剪断或纯折断破坏将不再适用,如何有效地建立一种该类型边坡的稳定性分析方法对边坡工程设计显得尤为重要。在此基础上,本文针对反倾层状岩质边坡平动-转动组合破坏问题,利用塑性极限分析上限定理,结合相关联流动法则和变形协调条件建立容许运动场,构建考虑平动加转动的岩质边坡稳定性极限分析理论计算模型,进一步采用遗传优化算法进行全局极值搜索,得到基于遗传算法的平动加转动岩质边坡极限分析稳定性系数上限解和相应的最优折断面。采用该方法与改进的 Goodman Bray方法进行对比,计算结果表明所提方法具有一定的适用性。同时利用该方法,对反倾岩质边坡的稳定性影响因素进行了敏感性分析。     

磷化氢在海水中的转化及其影响因素研究

通过室内模拟实验,研究了磷化氢在海水中的转化过程,分析了光照条件、氧化条件、通入方式等环境因子对磷化氢转化的影响。结果表明,在日光照射条件下,磷化氢在海水中可部分转化为磷酸盐,但转化率相对较低,转化过程复杂;环境因子在磷化氢转化过程中起重要作用,日光照射、紫外光照射、氧化条件等均能促进海水中磷化氢向磷酸盐和总磷的转化,在紫外光照射条件下检测到磷化氢向磷酸盐的最大转化率为66.6%。本研究对探索磷化氢在海洋环境生物地球化学循环中的作用具有重要意义。     

CONVEX MODEL FOR SEISMIC DESIGN OF STRUCTURES—I: ANALYSIS

A convex model is used to estimate the maximum response of structural systems subjected to uncertain seismic excitations. The convex model is based on the assumption that the energy of the excitation is bounded . A reduction factor, defined in the modal domain by dividing the results obtained from the convex model by those from the time history of the actual record, is used to calibrate the convex model. An average reduction factor is also defined by averaging a set of excitation-specific reduction factors. The average reduction factor can be used for unknown excitations with an assumed energy bound and certain common earthquake characteristics. These common characteristics can be defined either by a set of previous earthquakes in the region or by regional earthquake spectra. The convex model using the average reduction factor yields acceptable predictions of the maximum response.     

Assessment of radionuclide transport uncertainty in the unsaturated zone of Yucca Mountain

The present study assesses the uncertainty of flow and radionuclide transport in the unsaturated zone at Yucca Mountain using a Monte Carlo method. Matrix permeability, porosity, and sorption coefficient are considered random. Different from previous studies that assume distributions of the parameters, the distributions are determined in this study by applying comprehensive transformations and rigorous statistics to on-site measurements of the parameters. The distribution of permeability is further adjusted based on model calibration results. Correlation between matrix permeability and porosity is incorporated using the Latin Hypercube Sampling method. After conducting 200 Monte Carlo simulations of three-dimensional unsaturated flow and radionuclide transport for conservative and reactive tracers, the mean, variances, and 5th, 50th, and 95th percentiles for quantities of interest (e.g., matrix liquid saturation and water potential) are evaluated. The mean and 50th percentile are used as the mean predictions, and their associated predictive uncertainties are measured by the variances and the 5th and 95th percentiles (also known as uncertainty bounds). The mean predictions of matrix liquid saturation and water potential are in reasonable agreement with corresponding measurements. The uncertainty bounds include a large portion of the measurements, suggesting that the data variability can be partially explained by parameter uncertainty. The study illustrates propagation of predictive uncertainty of percolation flux, increasing downward from repository horizon to water table. Statistics from the breakthrough curves indicate that transport of the reactive tracer is delayed significantly by the sorption process, and prediction on the reactive tracer is of greater uncertainty than on the conservative tracer because randomness in the sorption coefficient increases the prediction uncertainty. Uncertainty in radionuclide transport is related to uncertainty in the percolation flux, suggesting that reducing the former entails reduction in the latter.     

矩阵相乘的压缩存储算法

介绍了对稀疏矩阵进行压缩存储时,稀疏矩阵相乘运算的基本思想和算法.在此基础上,探讨了避免矩阵中零元素相乘的无效操作,奠定了该算法的矩阵变换基础.     

非线性破坏准则与岩土材料地基承载力研究

根据线性Mohr-Coulomb破坏准则,前人运用极限平衡法、滑移线理论或极限分析等方法求解地基承载力问题,但实际上岩土在剪切破坏过程中破坏准则具有非线性。因此,为了研究非线性破坏准则对地基承载力的影响,基于上限定理,通过“切线法”引进变量,根据能量耗散情况,将承载力问题转变为非线性规划问题,运用“序列二次规划算法”求出地基承载力的最优解。数值计算结果表明,当非线性破坏准则转变为线性破坏准则时,非线性参数对地基承载力有重要影响。     

Stage–discharge estimation in straight compound channels using isovel contours

For various hydrological applications such as flood control projects, a knowledge of stage–discharge relationship is of particular interest to river engineers. Stage–discharge curves in compound channels cannot be easily predicted in comparison with single channels due to their 3D characteristics of flow. In this paper, the concept of cross‐sectional isovel contours is used for estimation of stage–discharge curves in compound channels. The multivariate Newton's method is applied to the difference between the observed and estimated data to optimize the exponent values of the governing parameters. The accuracy of the proposed model is tested successfully against available experimental results, which are taken from the Flood Channel Facility (FCF) laboratory. Then the results are compared with the Single and Divided Channel Methods (SCM and DCM, respectively), the Weighted Divided Channel Method (WDCM), the exchange discharge method (EDM), and the Coherence Method (COHM). The average values of mean absolute percentage error (MAPE) and normalized root mean square error (NRMSE) in discharge estimation based on each referenced section at any level for 6 sections of the experimental cases are within 3.1% and 0.023, respectively. The biggest advantage of the proposed method is its inherent simplicity, which does not need any calibration.