Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method     

Omid Ghaffaripour  Golnaz A. Esgandani  Arman Khoshghalb  Babak Shahbodaghkhan 《国际地质力学数值与分析法杂志》2019,43(11):1919-1955

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.  相似文献   

用双三次样条函数和GPS资料反演现今中国大陆构造形变场   总被引：38，自引：20，他引：18  

杨少敏  游新兆  杜瑞林  乔学军  王琪 《大地测量与地球动力学》2002,22(1):68-75

将中国大陆现今构造变动视为一种连续的地壳变形,利用双三次样条函数模拟了近期GPS测定的大陆内部及周边地区412个测站速率,反演大陆地区自洽的构造变动速度场和应变率场.模拟结果显示:印度板块与欧亚板块的碰撞、挤压是构成中国大陆内部岩石层水平形变的主要驱动力.印度板块在东喜马拉雅构造结深深插入青藏高原,造成地壳大规模的缩短和抬升.青藏高原东南部的喜马拉雅带、拉萨和羌塘地块以及青藏高原东南边的川滇地区,内部构造活动强烈,其内部的构造变形包含地壳碎片的冲断、褶皱和侧向逃逸.大陆地壳(或岩石圈)的增厚,尤其是喜马拉雅山脉南北向的快速缩短和青藏高原东西向的缓慢拉张,大约吸收了印欧板块会聚量的85%,西藏中东地区东西向的拉张速率达到了(16±2.0)mm/a,且顺时针方向扭转明显.印度板块相对欧亚板块运动的欧拉极为(29.7°N, 19.3°E, 0.392°/Ma);华南地块相对于欧亚大陆向东(102°±7.4°)南的运动速率是(11±1.54)mm/a,华南块体相对欧亚板块运动的欧拉极为(62.25°N, 126.56°E, 0.141°/Ma);塔里木地块相对较稳定,其西部运动速度高于东部运动速度,作顺时针方向旋转.总体上讲,中国大陆运动方向为北偏东呈辐射状,从西部近南北方向的运动转向东部地区东南方向的运动,绕东喜马拉雅构造结有一顺时针方向的旋转.横穿喜马拉雅构造带及青藏内部的南北向压缩速率为(19±2.0)mm/a,横穿西天山构造带的南北向压缩平均速率为(13±1.5)mm/a,横穿东天山构造带的南北向压缩平均速率为(6.0±1.4)mm/a.阿尔金断裂带的左旋走滑速率为(6±1.2)mm/a.  相似文献   

SAMPLING REMOTELY SENSED IMAGERY FOR STORAGE,RETRIEVAL, AND RECONSTRUCTION*     

P. M. Atkinson  P. J. Curran  R. Webster 《The Professional geographer》1990,42(3):345-353

A major difficulty in remote sensing is handling the many data from sensors aboard aircraft and satellites. In this paper we identify an optimal procedure for sampling remotely sensed data before their storage or on their retrieval. The procedure depends on spatial correlation in the scene and uses kriging to estimate values that have been lost. An example in which data from an airborne multispectral scanner could be diminished to only about one tenth without serious loss of precision illustrates the method.  相似文献   

重力插值方法研究     

向文 《大地测量与地球动力学》1997,(2)

为了根据离散观测数据构制连续空间重力变化图像，分析和讨论了３种数值插值方法，计算结果表明多面函数方法插值精度最高。由于逐步回归分析筛选核函数中心点的计算繁琐，文中提出根据分形理论和Ｓｈａｎｎｏｎ取样定理来确定核函数中心点。对滇西试验场进行模拟试算，插值精度可达到４～５（１０－８ｍｓ－２）。  相似文献   

Heat capacity and phase equilibria of wadeite-type K2Si4O9     

Wenjun Yong  E. Dachs  A. C. Withers  E. J. Essene 《Contributions to Mineralogy and Petrology》2008,155(2):137-146

The low-temperature heat capacity (C _p) of Si-wadeite (K₂Si₄O₉) synthesized with a piston cylinder device was measured over the range of 5–303 K using the heat capacity option of a physical properties measurement system. The entropy of Si-wadeite at standard temperature and pressure calculated from the measured heat capacity data is 253.8 ± 0.6 J mol⁻¹ K⁻¹, which is considerably larger than some of the previous estimated values. The calculated phase transition boundaries in the system K₂O–Al₂O₃–SiO₂ are generally consistent with previous experimental results. Together with our calculated phase boundaries, seven multi-anvil experiments at 1,400 K and 6.0–7.7 GPa suggest that no equilibrium stability field of kalsilite + coesite intervenes between the stability field of sanidine and that of coesite + kyanite + Si-wadeite, in contrast to previous predictions. First-order approximations were undertaken to calculate the phase diagram in the system K₂Si₄O₉ at lower pressure and temperature. Large discrepancies were shown between the calculated diagram compared with previously published versions, suggesting that further experimental or/and calorimetric work is needed to better constrain the low-pressure phase relations of the K₂Si₄O₉ polymorphs. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.  相似文献   

一种计算台风风场的方法     

陈孔沫 《热带海洋学报》1994,(2)

揭示Ｒａｎｋｉｎｅ涡风场模式和Ｊｅｌｅｓｎｉａｎｓｋｉ风场模式之间的联系，并设计了一种台风风场分布模式，它的风速分布曲线落在Ｊｅｌｅｓｎｉａｎｓｋｉ和Ｒａｎｋｉｎｅ涡两个风场模式的风速分布曲线之间，具有一个既优于Ｊｅｌｅｓｎｉａｎｓｋｉ又优于Ｒａｎｋｉｎｅ涡的风速衰减速率，因此它同时克服了Ｒａｎｋｉｎｅ涡模式计算风速偏小和Ｊｅｌｅｓｎｉａｎｓｋｉ模式计算风速偏大的缺点，以一种比较合理的变化趋势向远方衰减，成为一个比较切合实际的台风风场分布模式。同时，文中提出的移行台风风场计算方法对宫崎正卫、上野武夫和Ｊｅｌｅｓｎｉａｎｓｋｉ模式都有一定的改进。  相似文献   

Numerical and simplified methods for the calculation of the total horizontal wave force on a perforated caisson with a top cover     

Xuefeng Chen  Yucheng Li  Bin Teng 《Coastal Engineering》2007

The VOF method and the k–ε model, combined with the equation of state of air at constant temperature, have been used to calculate the total horizontal wave force caused by monochromatic waves acting on a perforated caisson with a top cover. From comparison of various parameters, such as the total horizontal force, the pressure difference on the front wall, the pressure on the back wall and the pressure on the top cover, between the numerical results and test data, it can be seen that the numerical results agree well with the test data. It is concluded that the method described in this paper can be utilized to calculate wave forces acting on perforated caissons with a top cover in the case of nonovertopping, nonbreaking waves. A simplified method to calculate the total horizontal force has been developed, based on test data, using a least-squares method. A comparison between the numerical results and the values calculated from the simplified equations shows good agreement. Therefore the simplified equations can be used in engineering applications to evaluate the total horizontal force on a perforated caisson with a top cover.  相似文献   

黑潮区台风域海-气热输送量的计算          下载免费PDF全文

林晓能 《海洋与湖沼》1993,24(3):319-324

根据1972和1974年经过黑潮区的两个台风个例,计算和分析了台风各部位在不同的海气条件下热量输送的垂直分布特点。表明,在台风前部信风带中,热量垂直输送限于700hPa以下,在1000hPa处有约167.47—209.34J/(cm~2·h)的最大值,700hPa以上热量向下输送;在台风中心区,热量输送量最大,在1000hPa处接近420J/(cm~2·h),在850hPa附近可能有热输送量的最大值,其上随高度减小;在台风后部偏南气流中,向上热输送量较小;各层最大值不超过84J/(cm~2·h)。  相似文献   

斜坡提砌石护面厚度的计算     

张怀慧  李玉成 《海洋工程》2002,20(2):14-18

介绍了斜坡式防波堤砌石护面中的插砌条石和干砌块石在波浪作用下的力学图式及作用机理 ,重点给出考虑摩擦力作用时的砌石块体稳定厚度的计算公式 ,并结合工程实例验证了公式的适用性  相似文献   

港口靠船桩工作性状计算的双参数法   总被引：8，自引：0，他引：8  

谢耀峰 《海洋工程》2002,20(2):38-42,48

港口靠船桩是承受重复性荷载或多循环荷载的一种横向承载桩。双参数法能较好地描述靠船桩 (单桩、群桩 )的桩土工作性状。用双参数法进行计算时 ,所给的地面条件必须由重复荷载或多循环荷载作用下实测出来。根据现场试桩资料 ,标定桩土参数 ,即k =mx1/n的指数l/n、桩土相对柔度系数α、桩土综合刚度EI等值 ,可以用来设计试桩附近局部地区同类条件下的其它长桩。对于没有进行试桩的某些工程 ,可以联合运用双参数法和p -y曲线法进行桩的分析  相似文献