船舶与海洋平台结构冰载荷的高性能扩展多面体离散元方法1)
HIGH-PERFORMANCE DILATED POLYHEDRAL BASED DEM FOR ICE LOADS ON SHIP AND OFFSHORE PLATFORM STRUCTURES1)
通讯作者: 2) 季顺迎,教授,主要研究方向:计算颗粒力学及其工程应用. E-mail:jisy@dlut.edu.cn
收稿日期: 2019-08-31 接受日期: 2019-10-8 网络出版日期: 2019-10-08
基金资助: |
|
Received: 2019-08-31 Accepted: 2019-10-8 Online: 2019-10-08
作者简介 About authors
船舶与海洋平台结构的冰载荷是寒区海洋工程结构物设计中的关键参数,而离散元方法是有效计算结构冰载荷的重要手段. 本文采用基于闵可夫斯基和原理的扩展多面体离散元方法模拟船舶与海洋平台结构的相互作用过程. 其中,构造扩展多面体的近似包络函数并建立了基于优化模型的快速接触搜索算法;考虑单元间粘结作用的刚度软化过程建标识码元间的粘结-破碎模型. 同时,发展了 CPU-GPU 协同异构环境下的高性能并行算法. 为分析海冰与海洋结构作用中的冰载荷,采用ISO标准验证了扩展多面体离散元分析结构冰载荷的准确性. 采用离散元方法计算了船舶结构的冰载荷,研究了船舶结构表明的线载荷分布特点,并采用船舶结构冰阻力经验公式验证了计算结果的合理性. 采用离散元方法计算了平整冰区与多桩腿平台结构的相互作用,分析各桩腿上的冰载荷特点. 针对碎冰区的海冰管理过程,采用离散元方法分析了船舶结构绕行过程中的船舶和海洋平台结构冰载荷. 本文方法可有效应用于海洋结构冰载荷分析,能为极地船舶与海洋平台结构的设计和安全运行提供科学的分析手段.
关键词:
The ice loads on ship and offshore platform structures is the key factor in structure design for cold regions. The discrete element method (DEM) is an important approach to determine the ice load on structures. According to the Minkowski sum theory, the dilated polyhedra based DEM is employed to simulate the interaction between sea ice and ship and offshore platform structures in this paper. In the dilated polyhedra based DEM, the enveloped function of the dilated polyhedron is generated to establish the fast contact detection algorithm based on the optimization model. Meanwhile, the bond-break model between elements is established by considering the stiffness softening process between bonded elements. Accordingly, the high-performance algorithm based on CPU-GPU cooperative-heterogeneous environment is developed. The ISO standard is employed to validate the ice load determined by the dilated polyhedra based DEM for better engineering applications of the interaction between sea ice and marine structures. The ice load on ship hull is calculated by the proposed method while the line load distribution on ship hull is studied. The ice resistance of ship hull is compared with the result by Lindqvist empirical formula to validate the accuracy of DEM simulations. The interaction between level ice and multi-leg platform is simulated while the ice load on each leg is analyzed. For the ice management in broken ice regions, the ice load on ship and offshore structures is simulated when the ship navigates around the offshore platform in circle. The proposed method can be effectively applied in the analysis of ice load on marine structures, and can provide a scientific approach for the design and safety operation of ship and offshore structures.
Keywords:
本文引用格式
刘璐, 尹振宇, 季顺迎.
Liu Lu, Yin Zhenyu, Ji Shunying.
引 言
海冰与船舶结构作用过程中,由于船体不同部位几何结构的差异,海冰会呈现弯曲、挤压、屈曲、剪切等不同的破坏模式. 在海冰和船艏的作用过程中, 弯曲破坏产生的冰载荷最为显著[4]. 该过程受结构形状和海冰参数的影响,海冰会发生径向或环向断裂破坏. 当海冰作用于船体舯部和船艏两侧时,由于该部分结构的垂向倾角很小,海冰呈挤压破坏状态. 此时海冰在船体表面的局部冰压力分布成为船体结构安全的关键因素[5-6]. 在海冰与海洋平台结构的相互作用研究中,主要关注结构的最大静冰力和交变动冰力[7-8]. 根据结构的几何形状,海冰与直立结构和坡面结构作用时的破坏模式和冰载荷形式区别较大. 海冰与直立结构作用主要发生挤压破碎,冰载荷可能表现为较高频率的持续性冰力[9];海冰与坡面结构作用时主标识码弯曲破坏,冰载荷通常表现为规律性的交变冰力[10].
离散元方法自 20 世纪 80 年代开始被用于海冰碰撞和漂移等研究中[11]. 近年来,基于离散元方法的海冰与船体和标识码台结构间的相互作用研究取得了较大进展. 海冰离散元方法的计算单元可采用球体[12]、圆盘[13]、多面体[14-15]和扩展多面体[16]等不同形态. 其中,球体单元通过平行粘结模型进行粘结并考虑单元之间的粘结失效,使离散的球体单元粘结为平整冰排,在海流的作用下与海标识码发生碰撞进而破碎,具有模型简单和计算效率高的特点[17]. 针对球体单元的尺寸效应带来参数不确定的问题,可通过离散元模拟海冰的单轴压缩和三点弯曲试验并结合渤海海冰的强度试验结标识码用无量纲化的单元尺寸参数分析并总结试样强度的变化规律,从而能够建立标准化的海冰宏微观强度关系并用于海冰的标识码拟中[18-19]. 在宏观尺度上,球体单元需要大量的单元粘结形成较大的冰区,而多面体单元可用较少的单元数模拟较大区域的海冰,具标识码的优势[20-21]. 多面体离散元模型主要模拟平整冰的破碎、重叠和堆积过程,可合理地描述冰块的几何形态. 针对碎冰区冰块形态呈多边形且随机分布的特点,还可采用 Voronoi 切割算法生成具有不同冰块尺寸、几何形态和密集度标识码分布[22]. 目前,采用块体离散元方法对堆积冰的受压形变过程、海冰与船体结构的相互作用过程进行了数值模拟[23-24]. 但多面体的接触检测等面临计算效率低的问题,合理构建其粘结-破碎模型也是基于多面体的海冰离散元方标识码关键问题[25-28].
最近,基于扩展多面体的快速搜索算法标识码-破碎模型发展较快,可应用于海冰与结构相互作用过程的分析中[29]. 基于闵标识码基和原理可将基本的几何单元与球体单元进行空间叠加从而构成扩展单元. 采用多面体和球体的闵可夫斯基和可构造扩展多面体单元[30]. 据此可发展基于扩展多面体标识码-破碎模型[31-32]和 DEM-LBM (Lattice Boltzmann method) 的流固耦合模型[33],并用于脆性材料水力压裂过程的模拟分析中. 由于扩展多面体单元将尖锐的角点和棱边转换为光滑的球面和圆柱面,在接触判断中避免了多面体之间难以通过几何元素直接判标识码的缺点,提高了接触搜索的效率和稳定性;还可采用精确化的Hertz接触力模型计算接触力,提高接触力计算的准确性[16,34]. 另外,根据类似的“扩展单元”方法,采用多面体扩展函数构造与扩展多面体类似的几何形状,可通过扩展函数表达式之间的优标识码求解单元之间的接触搜索问题[35-36]. 通过优化模型解决接触问题可有效避免几何元素之间复杂的几何判断,也可避免几何方法中产生多个接触点引起的物理失真,标识码计算中具有较高的执行效率[37-38]. 在扩展多面体单元的粘结-破碎模型中,考虑粘结作用的刚度软化过程可避免粘结失效时突然的能量释放,标识码高海冰断裂过程模 拟的平稳 性[29].
本文基于闵可夫斯基和原理构造扩展多面体单元,发展扩展多面体间的快速搜索算法和粘结-破碎模型,标识码CPU-GPU协同方式进 行并行计算. 采用该方法计算海冰与船舶和海洋平台结构的相互作用过程,通过ISO标准公式和现场监测数据验证分析结构冰载荷计算的合理标识码而采用将该方法应用于具体的工程分析中.
1 扩展多面体离散元方法及高性能计算
采用二维的 Voronoi 切割算法将二维平面划分为若干多边形,将多边形在垂直纸面的方向上进行拉伸即形成具有一标识码的多面体,根据闵可夫斯基和原理将该多面体和球体结合即可构造扩展多面体形态的海冰单元.
1.1 海冰的扩展多面体构造
假设空间中两个几何体代表的点集 $A$ 和 $B$,且集合 $A$ 和 $B$ 为几何封闭形态,那么集合 $A$ 和 $B$ 的闵可夫标识码可定义为
图1
图1
基于闵可夫斯基和的扩展多面体单元
Fig.1
Dilated polyhedral element based on the Minkowski sum
Voronoi 切割算法根据随机的空间点将空间划分为随机的任意多面体或多边形,可合理表述自然条件下海冰的标识码形状规律[39]. 这里采用二维 Voronoi 切割算法,将一定范围的二维平面划分为若干任意多边形,如图2(a) 所示. 其中,多边形的形状可根据随机点的分布进行控制[40]. 然后将每个多边形在垂直纸面方向进行拉伸形成多面体,根据闵可夫斯基和原理则可生成扩展多面体形态的海冰单元,如图2(b) 所示. 其中,$h$ 为海冰厚度. 若$r$为扩展半径,则 $t= h - 2r$. 对于碎冰初始场的生成,则可通过实际海冰面积和冰区范围控制碎冰密集度[22];而对于平整冰的生成,则可直接划分网格.
图2
图2
Voronoi 切割算法生成扩展多面体形态的海冰单元
Fig.2
Dilated polyhedron shaped sea ice element based on Voronoi tessellation algorithm
1.2 扩展多面体离散元的接触模型
在扩展多面体单元的接触模型中,可构造扩展多面体的近似包络函数,从而采用优化模型建立单元间的接触搜索算法. 通 过二阶多面体扩展函数与球面函数的加权求和形式可构造具有多面体特征的光滑函数[36],其归一化形式可写作
式中,$k$ 为光滑度系数;$R$ 为球面函数的半径.
图3
图3
扩展多面体的包络函数随参数$k$的变化
Fig.3
The envelope function of the dilated polyhedron under different parameter $k$
式中,$f_{A}$ 和 $f_{B}$ 对应两个单元的包络函数. 如图4 所示,两个单元接触时满足上式的解 $C_{AB} $ 同时位于两个包 络函数内部,其可被视为接触重叠区域的接触中心点.
图4
根据拉格朗日乘子法构造拉格朗日函数如下
由 $\nabla L\left( {{\pmb X},\lambda } \right) = 0$ 可将该优化问题变为四元非线性方程组,可采用牛顿迭代法进行求解该方程组. 在牛顿迭代法的初始点选取标识码先在迭代中通过调整光滑度系数 $k$ 由 1.0 渐变至较小值,从而使初始点逐渐逼近最终解附近,从而提高迭代的收敛效率[41]. 另外,还可在迭代中设置稳定性系数,提高迭代的稳定性[42]. 求得式 (3) 表示的优化模型结果 ${\pmb C}_{AB}$ 之后,由该点出发迭代搜索基本多面体上距离最近的点来确定两个颗粒上的接触点,进而可确定接触重叠量和接触法向. 由于该计算并没有复杂的判断过程,且接触中心点 ${\pmb C}_{AB}$ 一般距这两个最近点十分接近,所以该迭代过程并不会耗费较多计算资源[41].
式中,$k_{\rm n} $ 为法向接触刚度,$k_{\rm n} = 4E^\ast \sqrt {R^\ast } / 3$,$E^\ast $ 和 $R^\ast $ 是等效的弹性模量等效颗粒半径 [45];$\delta _{\rm n} $ 是两颗粒的法向接触重叠量;${\delta }”$ 是法向相对速度;$c_{\rm n} $ 为法向阻尼系数
式中,$m_{AB}$ 为等效质量[45], $\zeta _{\rm n} $ 是无量纲阻尼系数
式中,$e$ 是回弹系数,可参考不同碰撞速度下颗粒回弹系数的计算[46].
切向接触力与切向接触重叠量和摩擦力有关,一般切向力最大不超过摩擦力. 这里可采用简化的方式避免与摩擦力的频繁比较和判断[44]. 其中,切向弹性力$F_{\rm t}^{\rm e} $可写作
式中,$\mu $ 是摩擦系数;$\delta _{\rm t} $ 是切向重叠量,一般根据相对切向速度并逐步叠加的方法累计计算切向重叠量;$\delta _{\rm t}^{\max} $ 是最大切向重叠量,其主要由摩擦系数和法向重叠量决定,且 $\delta _{\rm t}^{\max} = \mu \delta _{\rm n} \left( {2 - \nu } \right) / \left( {2 - 2\nu } \right)$,$\nu$ 为泊松比. 扩展多面体单元间的切向黏滞力 $F_{\rm t}^{\rm v} $ 可写作
式中,${\delta }”$ 是切向相对速度;$c_{\rm t} $是切向阻尼系数,这里取 $c_{\rm t} = c_{\rm n} / \left( {2\left( {1 + \nu } \right)} \right)$[47]. 该切向力计算方法避免了每个时间步中与摩擦力的大小对比,有效地提高了计算效率.
由于本文的接触搜索只有一个接触点,即互相接触的两个单元只有一个接触力,故以上接触力模型无法消除单元绕标识码向的非物理转动. 为此,这里在单元的转动方程中引入绕接触法向的滚动摩阻力矩 ${\pmb M}_{\rm r}$[43]
式中,$\mu _{\rm r} $ 是滚动摩擦系数;${\pmb \omega }_{ AB} $ 是相对角速度.
1.3 扩展多面体离散元的粘结-破碎模型
1.3.1 扩展多面体单元的粘结模型
图5
式中,${\pmb d}$ 是两个粘结节点之间的距离向量;${\pmb n}$ 为交界面的法向,这里取为两个单元粘结面的外法向之差,即 ${\pmb n} ={\rm unit} ({\pmb n}_{i}- {\pmb n}_{j})$;$C_{ij}$ 为两个粘结单元之间的特征长度,与刚体有限元中的定义相同,标识码个粘结面到各自单元质心的距离之和,可写作
由于粘结力会传递到单元质心上进行合力计算,因此这里切向应变只考虑与法向垂直的一个方向应变标识码向为 ${\pmb t} ={\rm unit} ({\pmb d} - ({\pmb d} \cdot {\pmb n} ) {\pmb n} )$,切向应变可写作
根据图5(b) 所示的粘结力模型,采用三维条件下的弹性矩阵,两个粘结单元在交界面上的弹性应力 ${\pmb \sigma }^{\rm e}$ 可写作
式中,$E$ 为材料弹性模量;$\nu $ 为材料泊松比;${\pmb \sigma }^{\rm e}$ 考虑法向和切向的应力,${\pmb \sigma }^{\rm e} = \{\sigma ,\tau \}^{\rm T}$;${\pmb \varepsilon }$ 考虑粘结节点之间的法向和切向应变, 即 ${\pmb \varepsilon } = \{\varepsilon _{\rm n}, \varepsilon _{\rm t} \}^{\rm T}$.
在粘结力模型中,考虑法向和切向阻尼系数 $C_{\rm n} $ 和 $C_{\rm s} $ 表示的黏滞作用. 该阻尼效应可代表真实的物理阻尼并消耗动能,从而提高模拟的稳定性. 一般情况下,通过与刚度相关的常系数 $\beta $ 计算阻尼系数. 因此,由黏滞作用产生的应力 ${\pmb \sigma }^{\rm v}$ 可写作
式中,${\pmb \varepsilon }”$是应变率,且${\pmb \varepsilon }''=\{ \dot \varepsilon_n,\dot \varepsilon_t\}^{\rm T}$. 两个粘结节点之间的粘结力可写作
式中,$A$ 是交界面面积;$n$ 是一个粘结面上的粘结节点个数. 通过上式计算单元受到的粘结力,将其转换到单元质心上,并结合接触力等其他受力,即可得到每个单元质心受到的合力情况.
可以看出,直接在整体坐标系下将节点间的粘结力传递到单元质心上,每个单元的运动可直接在单元的质心上进行求解. 因标识码方式可高效地融入离散元方法中并形成算法和数据结构高度统一的数值方法.
1.3.2 扩展多面体单元间的粘结失效模型
在球体单元常用的平行粘结模型中,大多采用拉剪分区的方式判断粘结的失效,即法向拉伸和切向剪切分别判断标识码否达到临界值,若其中一个失效则整个粘结作用失效[50]. 显然,该方法对材料的破坏过程描述较为片面,不能合理地评估裂纹附近材料的断裂特性. 因此,本文同时考虑法向拉伸和剪切应力计算临界应力强度,在超过临界应力强度之后引入损伤的概念从而标识码度的软化过程,结合拉伸破坏和剪切破坏并采用混合的断裂能模型计算断裂能,从而根据断裂能确定粘结节点之间的临界变形. 若节点间的变形大于该临界变形,则该粘结失效.
式中,$\bar{\sigma }”$ 和 ${\tau }”$ 分别是法向和切向粘结强度; $\left\langle \right\rangle $ 表示 Macaulay 括号,且有: $\left\langle x \right\rangle = x$,若 $x \geqslant 0$;$\left\langle x \right\rangle = 0$,若 $x < 0$. 这里假设拉应力为正且压应力为负,那么 $\left\langle \sigma \right\rangle $ 即表示只考虑拉伸应力. 图6 为该强度判定准则在拉伸和剪切应力坐标系中的表示,图中曲线内形状为椭圆,椭圆之外即表示材料损伤. 另外,拉剪分区判断准则即为图中虚线与横纵轴围成的长方形边界,因此本文采用的强度准则比拉剪分区方法“更易”发生损伤或破坏.
图6
图6
拉伸和剪切应力结合的强度准则
Fig.6
Hybrid failure criterion by tensile and shear strength
式中,$C$ 是黏聚力;$\mu _{\rm b} = \tan\theta $ 是内摩擦系数,$\theta $ 是内摩擦角.
图7
图7
Mohr-Coulomb准则确定切向粘结强度
Fig.7
Shear bond strength by Mohr- Coulomb criterion
为避免采用拉剪分区方式判断粘结失效,本文的断裂准则通过等效的应力 $\sigma _{\rm m} $ 和变形 $\delta _{\rm m} $ 构造本构关系,可分别写作[54]
式中,$\delta _{\rm n} $ 和 $\delta _{\rm s}\sigma _{\rm m}^0 $ 和对应的极限变形 $\delta _{\rm m}^0 $. 分别是法向和切向变形量. 当粘结节点的应力状态满足式 (17) 时,可视为该节点的应力状态达到了极限强度,那么根据式 (19) 可计算等效的极限强度
如图8 所示,采用双线性模型定义等效应力和变形之间的关系. 这里引入复合损伤$D$描述极限强度后的刚度软化过程
图8
图8
等效应力和变形确定的双线性本构关系
Fig.8
Bi-linear constitutive relationship between equivalent stress and deformation
式中,$\delta _{\rm m}^{\max}$ 是加载历程中粘结节点上发生过的最大变形;$\delta _{\rm m}^{\rm f} $ 是粘结失效时的临界变形,即当 $\delta _{\rm m} > \delta _{\rm m}^{\rm f} $时粘结失效. 那么刚度软化阶段材料的刚度 ${k}”$ 变为
式中,$k$ 为极限强度范围内的材料刚度.
式中,$G_{\rm I}^{\rm c} $ 和 $G_{\rm II }^{\rm c} $ 分别对应于拉伸和剪切型裂纹的临界断裂能;$\eta$ 为常系数,对于脆性材料可取 $\eta = 1.75$[56]. 对应当前状态的拉伸和剪切型断裂能 $G_{\rm I} $ 和 $G_{\rm II } $ 可通过拉伸和剪切的变形计算二者的比例,写作
1.4 离散元高性能计算方法
1.4.1 球形包围盒的空间网格法
球形包围盒即是可完全包围单元的最小球体. 采用所有单元中最大的球形包围盒直径的 2 倍定义网格的尺寸,则标识码算域划分为由正方形网格构成的空间网格结构. 通过每个单元形心的空间坐标 $\left( {x,y,z} \right)$ 可计算其网格坐标 $\left( {i_x ,i_y ,i_z } \right)$
式中,$x_{\min} $, $y_{\min}$和$z_{\min} $ 分别代表计算域在空间中的最小值;$d_{\rm grid} $ 是网格尺寸;${\rm int}( )$ 是取整运算符. 若计算域的网格个数为 $N_{\rm grid} = N_x \times N_y \times N_z $,可将三维网格坐标对应到一维的网格序列上
图9
根据该包含关系,在某个单元的搜索中,只搜索其所在网格的相邻网格中的单元. 在二维空间中是以该单元所在网格为中心的3$\times $3网格范围,对应三维空间中是以该单元所在网格为中心的 3$\times $3$\times $3 的网格范围. 进一步比较球心距离和两球形包围盒半径之和判断单元是否可能接触,从而可建立每个单元的邻居列表.
GPU 对大量类型高度统一且相互无依赖的大规模数据计算具有更高的计算效率[59],而空间网格法标识码任务可分解为大量相互无依赖的数值计算任务,非常适合 GPU 的计算. 因此,这里采用 GPU 并行算法执行球形包围盒的空间网格法计算.
1.4.2 多面体包围盒的改进分离轴方法
对于邻居列表中的邻居单元,可采用多面体的方向包围盒 (oriented bounding box, OBB) 对每个其做更加细致的判断. 对于扩展多面体单元,该包围盒的形状是可包围该扩展多面体的最小多面体. 在 OBB 的接触判断中,最常采用的是分离轴方法 (separate axis algorithm, SAA) [60]. 一般来说,SAA适用于任意形状空间体的接触判断,包括凹和凸的. 为简化描述,二维条件下凹多边形的SAA方法如图10(a) 所示. 多边形$i$的边 ${ A B}$ 外法向为 ${\pmb n}_{ab}$,那么以该法向 ${\pmb n}_{ab}$ 为轴向建立分离轴 $l^{ab}$. 将多边形 $i$ 和待判断的多边形 $j$ 上的所有点投影到分离轴上. 分别计算两个多边形在分离轴上的投影距离 $l_i^{ab} $ 和 $l_j^{ab} $,判断两个投影是否有交叉. 如果存在交叉,则可能接触;若不存在交叉,则不存在接触.
图10
可以看出,分离轴方法需要对两个单元的每个角点、每个面进行计算并综合比较结果才能判断是否接触. 凸多面体具有所有几何 元素都位于其某个面单侧的特点. 因此,这里提出改进的 SAA 方法对凸多面体进行判断. 同样以二维情况为例,如图10(b) 所示. 从 多边形$i$的边 ${AB}$ 上的任意一点到待判断多边形 $j$ 上角点的向量为 ${\pmb e}_k^{ab} $,对于所有多边形 $j$ 的角点,若都能满足以下条件
那么可认为多边形 $i$ 与多边形 $j$ 不存在接触. 在实际判断两个多面体是否接触的计算中,如果在其中一个多面体中存在符合式 (26) 的面,则可判定两个多面标识码触,且不需对后面未做计算的角点和面进行计算,从而节省了计算量. 在最坏情况下,该方法与SAA方法计算量相同;而在最好情况下,该方法的计算效率会远高于 SAA 方法.
改进的 SAA 算法需要反复进行逻辑判断,适合在具有更强逻辑运算能力的 CPU 上进行. 因此,这里采用 CPU 并行算法执行改进的 SAA 算法.
2 扩展多面体离散元方法的验证
2.1 坡面结构冰载荷
图11
图12
式中,$\xi $ 是坡面几何参数
$ \xi = \left( {\sin \alpha + \mu \cos \alpha } \right) / \left( {\cos \alpha - \mu \sin \alpha } \right) $
$l_{\rm c} $是特征长度
图13 是离散元模拟和 ISO 标准计算的坡面结构冰载荷对比,图中标出了两者的对比误差. 对比表明离散元模拟结果与 ISO 标准在趋势上基本一致,在数值上也十分接近. 由于 30$^\circ$ 和 45$^\circ$ 时海冰会发生较完整的弯曲破坏,因此模拟冰载荷大小与 ISO 结果误差较小. 当斜坡角度增大到 60$^\circ$ 和 75$^\circ$ 时,海冰的破坏模式不完全为弯曲破坏,故冰载荷的误差相对较大. 一般情况下,ISO 标准定义的是最大值,因此总体上本文的离散元模拟结果符合 ISO 标准,能够说明基于扩展标识码的离散元方法具有较好的合理性.
图13
2.2 锥体平台结构冰载荷
图14
平整冰与锥体结构相互作用的模拟结果如图15 所示. 可以看出,平整冰在与锥体结构接触后发生弯曲破坏,破碎海标识码着锥体向上攀爬,从而在冰载荷时程曲线中形成规律性的峰值载荷. 渤海的现场监测数据中包含大量不同冰厚条件下的冰载荷,因此这里采用不同的冰厚模拟平整冰与锥体结构的相互作标识码过 ISO 标准和现场实测数据分析计算结果的合理性.
图15
图16 是离散元模拟、ISO 标准和渤海的现场监测数据对比. 可以看出,离散元计算结果在趋势上与 ISO 标准和现场监测结果保持一致,数值上基本处于 ISO 标准和现场监测数据之间. 当 $h < 0.15$m 时,离散元结果与 ISO 标准值更为接近;当 $h > 0.15$m 时,离散元结果明显小于 ISO 标准值. 因此,本文标识码元模拟结果与 ISO 标准值和渤海现场监测标识码具有较好的一致性,充分说明了扩展多面体离标识码法在海冰与海洋结构相互作用模拟中的合 理性.
图16
图16
离散元、ISO 标准和现场实测的锥体结构冰载荷
Fig.16
Ice load on cone by DEM, ISO standard and field data
3 船舶结构冰载荷的离散元分析
3.1 船舶结构冰载荷分析
图17
图18 所示是船体结构在平整冰区的破冰过程模拟. 从图中可以看出,平整冰在船艏的作用下发生弯曲破坏,且沿船艏两侧扩展开. 海冰在船体运行方向上发生持续的破碎,冰面上会出现与船宽相当的开阔水道. 同时,破碎后的冰块会浮在开阔水道上,部分碎冰会附着在船底并与船底摩擦,随后沿着船底运动到船艉. 平整冰破碎后会在船艏部区域发生堆叠淤积. 随着船体持续运动,该部分碎冰会与船舯、船艉发生摩擦并被逐渐排开.
图18
图18
船体结构在平整冰区的破冰过程模拟
Fig.18
Simulation of ship icebreaking in level ice region
图19 是船体破冰过程中结构上 3 个方向的冰载荷时程曲线. 从图中可看出,$x$ 和 $z$ 方向的冰载荷体现了海冰破坏模式主要为弯曲破坏. $x$ 方向冰载荷对船体在破冰过程影响最大,也是船体运动的主要推进方向. 因此,尽管在 $y$ 和 $z$ 方向上也受到了较大的冰载荷,但通常主要研究 $x$ 方向的冰载荷. 在后续的模拟中,主要分析 $x$ 方向的冰载荷.
图19
采用 Lindqvist 经验公式对比不同冰厚条件下船舶结构的冰阻力,即 $x$ 方向的冰载荷. Lindqvist 公式表示的船舶结 构冰阻力 $R_{i}$ 可写作[66]
式中,$v$ 是船冰相对速度;$L$ 是水线处的船长;$R_{\rm c} $, $R_{\rm b} $ 和 $R_{\rm s} $ 分别表示船艏挤压破坏阻力、弯曲破坏阻力和浸没阻力.
挤压破坏阻力 $R_{\rm c} $ 可写作
式中,$\sigma _{\rm f} $ 是海冰弯曲强度;$\varphi $ 为首柱倾角;$\alpha $ 为水线进角;$\psi = {\rm arctan}\left( {\tan\varphi / \sin\alpha } \right)$.
弯曲破坏阻力 $R_{\rm b} $ 写作
式中,$B$ 为船舶结构型宽.
浸没阻力 $R_{\rm s} $ 写作
式中
式中,$T$ 为船体吃水深度.
图20
图20
离散元和Lindqvist公式的船舶冰阻力对比
Fig.20
Ice force on ship hull by DEM and Lindqvist formula
图21
图22
图23
3.2 双船冰载荷分析
多艘破冰船或极区科考船在冰区协同作业时,会出现多船平行进入冰区并以平行方向在冰区破冰航行的情况. 如图24 所标识码拟两艘船前后在平行方向上以不同速度在冰区航行,初始时船一在前船二在后,且在水平方向上间隔一定距离 $d$. 另外,船一船速为 1.0m/s,船二船速为 1.5m/s.
图24
图25
图25
平整冰区两船平行航行的模拟
Fig.25
Simulation of double ships navigating parallel in level ice region
图26
图27
4 海冰管理中船舶与海洋平台结构冰载荷
采用扩展多面体的离散元方法计算多桩腿海洋平台结构与平整冰的作用过程,分析不同桩腿上的冰载荷特点. 同时考虑标识码理过程中船舶和海洋平台结构同时与海冰的作用过程,分析船舶结构在绕海洋平台航行过程中的冰载荷.
4.1 平整冰区多桩腿海洋平台结构的冰载荷分析
图28
图29
图30
4.2 海冰管理中船舶及海洋平台结构的冰载荷
海洋平台结构在冰区运行过程中可能会受到严重的冰载荷影响,在恶劣冰况下需要破冰船协助海洋平台进行破冰作标识码低海洋平台结构上的冰载荷. 采用扩展多面体离散元方法计算船舶结构在碎冰区的海洋平台结构的绕行过程,分析船舶和海洋平台结构的冰载荷. 如图31 所示,船舶结构以半径 100m 绕海洋平台结构运动,行进速度为 1.0m/s. 为更合理地分析船舶结构的冰阻力,定义船舶结构的标架坐标系. 船舶结构的冰载荷均在标架坐标系内进行统计. 另外,为生成更加随机的碎冰初始场,初始时刻给每个单元绕 $z$ 轴的随机转速. 图32 是海冰密集度为 60${\%}$ 时海冰管理过程的离散元模拟,海冰在与船舶和海洋平台作用后均形成了开阔水道.
图31
图32
图33
图33
海冰管理过程中船舶和海洋平台结构冰载荷
Fig.33
Ice load on ship and platform during ice management
5 结 论
本文采用闵可夫斯基和构造了扩展多面体单元,根据扩展多面体的包络函数发展了单元接触搜索的优化算法,同时建立标识码刚度软化过程的单元间粘结-破碎模型,由此建立了基于扩展多面体的离散元方法,并将该方法应用于船舶标识码平台结构冰载荷的计算分析中. 为验证本文方法的合理性,采用 ISO 标准和渤海现场监测冰载荷数据验证分析了本文计算结果的准确度. 采用本文方法进一步模拟了工程实际中船舶和海洋平台结构的冰载荷. 相关计算分析结果表明,扩展多面体的离散元方法在海冰与海洋结构的相互作用分析中体现出了很好的适应性,对船标识码洋结构的冰载荷分析准确性较高. 该方法可对复杂工况下的海冰与海洋结构相互作用过程进行详细地分析,为极地船舶与海洋平台结构的设计和安全运行提供了科学的分析手段.
此外,发展单元内的破碎准则从而建立单元间和单元内的多级破碎模型,可有效提高海冰与结构相互作用的分析精度标识码计算效率. 为更加全面详细地分析海冰与结构相互作用过程,还需考虑海水的影响并发展相应的DEM-CFD流固耦合模型.
参考文献
Local ice pressures measured on a strain gauge panel during the CCGS Terry Fox bergy bit impact study
,APETALA2/ETHYLENE RESPONSE FACTOR transcription factors (AP2/ERFs) play crucial roles in adaptation to stresses such as those caused by pathogens, wounding and cold. Although their name suggests a specific role in ethylene signalling, some ERF members also co-ordinate signals regulated by other key plant stress hormones such as jasmonate, abscisic acid and salicylate. We analysed a set of ERF proteins from three divergent plant species for intrinsically disorder regions containing conserved segments involved in protein-protein interaction known as Molecular Recognition Features (MoRFs). Then we correlated the MoRFs identified with a number of known functional features where these could be identified. Our analyses suggest that MoRFs, with plasticity in their disordered surroundings, are highly functional and may have been shuffled between related protein families driven by selection. A particularly important role may be played by the alpha helical component of the structured DNA binding domain to permit specificity. We also present examples of computationally identified MoRFs that have no known function and provide a valuable conceptual framework to link both disordered and ordered structural features within this family to diverse function.
抗冰平台的冰振疲劳分析
,
Fatigue-life analysis of ice-resistant platforms
,
Research trends in ice mechanics
,Seasonal acceleration of the Greenland Ice Sheet is influenced by the dynamic response of the subglacial hydrologic system to variability in meltwater delivery to the bed via crevasses and moulins (vertical conduits connecting supraglacial water to the bed of the ice sheet). As the melt season progresses, the subglacial hydrologic system drains supraglacial meltwater more efficiently, decreasing basal water pressure and moderating the ice velocity response to surface melting. However, limited direct observations of subglacial water pressure mean that the spatiotemporal evolution of the subglacial hydrologic system remains poorly understood. Here we show that ice velocity is well correlated with moulin hydraulic head but is out of phase with that of nearby (0.3-2 kilometres away) boreholes, indicating that moulins connect to an efficient, channelized component of the subglacial hydrologic system, which exerts the primary control on diurnal and multi-day changes in ice velocity. Our simultaneous measurements of moulin and borehole hydraulic head and ice velocity in the Paakitsoq region of western Greenland show that decreasing trends in ice velocity during the latter part of the melt season cannot be explained by changes in the ability of moulin-connected channels to convey supraglacial melt. Instead, these observations suggest that decreasing late-season ice velocity may be caused by changes in connectivity in unchannelized regions of the subglacial hydrologic system. Understanding this spatiotemporal variability in subglacial pressures is increasingly important because melt-season dynamics affect ice velocity beyond the conclusion of the melt season.
Effect of dynamic bending of level ice on ship’s continuous-mode icebreaking
,
Statistical analysis of ice crushing pressures on a ship’s hull during hull-ice interaction
,
Ice bound shipping is considerably more dangerous than shipping under more favourable water conditions. During winter the Baltic Sea can be covered with ice over an average of about 45% of its surface. The ice mostly reaches its maximum extent in late February or early March. The Bay of Bothnia located in the northern basin of the Gulf of Bothnia is prone to extreme frozen ice conditions. During winter the typical ice thickness in the northernmost areas of the Bay of Bothnia is about 70 cm for land-fast sea ice.
In this paper the ice crushing pressure on a ship's hull is analysed statistically. A series of tests took place over four days in the Aalto ice tank and involved two different ship models. The peak pressures caused by ice were calculated from four sensor sheets located at different positions on the hull. In the statistical part, the sum of the forces acting on a section of a ship's hull was modelled as a Poisson random process. Analysis of the cumulative distribution function (CDF) was modelled for the total force. The full scale comparisons of pressure-area values and line loads as a function of the load width were analysed by comparing the model scale data with full scale measurements onboard MS Arcturus and IB Sisu. Finally a short comparison is also conducted with the load level obtained from damage statistics gathered in the Baltic Sea. (C) 2011 Elsevier B.V.
Estimation of local ice pressure using up-crossing rate
,The variation of physical and mechanical properties of the lightweight bulk filling material with cement and expanded polystyrene (EPS) beads contents under different confining pressures is important to construction and geotechnical applications. In this study, a lightweight bulk filling material was firstly fabricated with Singapore marine clay, ordinary Portland cement and EPS. Then, the influences of EPS beads content, cement content, curing time and confining pressure on the mass density, stress-strain behavior and compressive strength of this lightweight bulk filling material were investigated by unconsolidated and undrained (UU) triaxial tests. In these tests, the mass ratios of EPS beads to dry clay (E/S) were 0%, 0.5%, 1%, 2%, and 4% and the mass ratios of cement to dry clay (C/S) were 10% and 15%. Thirdly, a series of UU triaxial tests were performed at a confining pressure of 0 kPa, 50 kPa, 100 kPa, and 150 kPa after three curing days, seven curing days, and 28 curing days. The results show that the mass density of this lightweight bulk filling material was mainly controlled by the E/S ratio. Its mass density decreased by 55.6% for the C/S ratio 10% and 54.9% for the C/S ratio 15% when the E/S ratio increased from 0% to 4% after three curing days. Shear failure more easily occurred in the specimens with higher cement content and lower confining pressure. The relationships between compressive strength and mass density or failure strain could be quantified by the power function. Increasing cement content and reducing EPS beads content will increase mass density and compressive strength of this lightweight bulk filling material. The compressive strength with curing time can be expressed by a logarithmic function with fitting correlation coefficient ranging from 0.83 to 0.97 for five confining pressures. These empirical formulae will be useful for the estimation of physical and mechanical properties of lightweight concretes in engineering application.
Quantitative analysis of ice sheet failure against an inclined plane
,
Model tests of four-legged jacket platforms in ice: Part 1. Model tests and results
,Two series of model tests were performed to investigate the dynamic ice loads on two types of four-legged jacket platforms and the consequently structural responses. The first test series was carried out with cylindrical models and the second test series with conical models at the water level. To simulate the ice induce vibrations on multi-legged structures, a new series of test equipments was designed. The ice load on each model leg, structural displacements in three directions, and foundation reactions in two directions were measured during the tests. Different ice failure modes and structural responses were observed for different ice conditions and structural stiffness values. This paper provides a description of the test techniques and the experimental results. The important test phenomena are also discussed in this paper. (C) 2013 Elsevier B.V.
Model test study of the nonsimultaneous failure of ice before wide conical structures
,
On applying granular flow theory to a deforming broken ice field
,
锥角对锥体结构抗冰性能影响的离散元分析
,
Influence of cone angle on anti-icing performance of conical structure with numerical simulations of discrete element method
Simulation of pancake ice load on a circular cylinder in a wave and current field
,APETALA2/ETHYLENE RESPONSE FACTOR transcription factors (AP2/ERFs) play crucial roles in adaptation to stresses such as those caused by pathogens, wounding and cold. Although their name suggests a specific role in ethylene signalling, some ERF members also co-ordinate signals regulated by other key plant stress hormones such as jasmonate, abscisic acid and salicylate. We analysed a set of ERF proteins from three divergent plant species for intrinsically disorder regions containing conserved segments involved in protein-protein interaction known as Molecular Recognition Features (MoRFs). Then we correlated the MoRFs identified with a number of known functional features where these could be identified. Our analyses suggest that MoRFs, with plasticity in their disordered surroundings, are highly functional and may have been shuffled between related protein families driven by selection. A particularly important role may be played by the alpha helical component of the structured DNA binding domain to permit specificity. We also present examples of computationally identified MoRFs that have no known function and provide a valuable conceptual framework to link both disordered and ordered structural features within this family to diverse function.
Discrete element analysis of ice loads on ships and structures
,In this paper, the versatility of discrete element method (DEM) in modelling ice-related problems is discussed and further demonstrated using the results from the DEM works conducted by the National Research Council's Institute for Ocean Technology (NRC-IOT) using a commercial code DECICE. These works include a wide range of ice-structure and ice-ship interaction problems of current interest, i.e. ice loads on conical structures, jamming of floes at bridge piers, modelling of the mechanical behaviour of ice rubble, pack ice stability and associated forces on offshore structures, rubble loads exerted on an inclined retaining wall, ridge keel resistance during seabed scouring, dynamic response of a moored conical drill-ship in ice and ship manoeuvring performance in ice. Representative simulations for each case are presented including load, motion and/or interaction process, whichever is appropriate. The simulations from DECICE were compared with experimental data and found satisfactory in terms of accuracy and real-time simulation capability. The accuracy is important for design and engineering of marine structures and ships, whereas the real-time simulation capability allows it to be used in marine simulators for personnel training and marine operation assessment. The performance of DECICE is also addressed and improved via the implementation of a new contact detection sequence and parallel considerations.
斜坡式防波堤前海冰堆积数值模拟
,本文将海冰在斜坡式防波堤前的堆积过程简化为二维的动力学问题,假定海冰为粘弹性体,在多边块体运动模拟系统基础上,采用离散元数值方法对冰排爬升堆积过程进行了实时仿真模拟.通过与二维理论断裂模型及物理模型试验结果对比,证明了此数值模型的合理性与可行性.
Numerical simulation of sea ice pile-up on inclined breakwater
本文将海冰在斜坡式防波堤前的堆积过程简化为二维的动力学问题,假定海冰为粘弹性体,在多边块体运动模拟系统基础上,采用离散元数值方法对冰排爬升堆积过程进行了实时仿真模拟.通过与二维理论断裂模型及物理模型试验结果对比,证明了此数值模型的合理性与可行性.
基于扩展多面体的离散单元法及其作用于圆桩的冰载荷计算
,
Dilated polyhedra based discrete element method and its application of ice load on cylindrical pile
海冰与自升式海洋平台相互作用GPU离散元模拟
,
GPU-based discrete element modelling of interaction between sea ice and jack-up platform structure
Validation of microparameters in discrete element modeling of sea ice failure process
,
水位变化对正倒锥体冰载荷影响的离散元分析
,
Influence of water level on ice load on upward-downward conical structure based on DEM analysis
Modelling the reorientation of sea-ice faults as the wind changes direction
,A discrete-element model of sea ice is used to study how a 900 change in wind direction alters the pattern of faults generated through mechanical failure of the ice. The sea-ice domain is 400 km in size and consists of polygonal floes obtained through a Voronoi tessellation. Initially the floes are frozen together through viscous-elastic joints that can break under sufficient compressive, tensile and shear deformation. A constant wind-stress gradient is applied until the initially frozen ice pack is broken into roughly diamond-shaped aggregates, with crack angles determined by wing-crack formation. Then partial refreezing of the cracks delineating the aggregates is modelled through reduction of their length by a particular fraction, the ice pack deformation is neglected and the wind stress is rotated by 90 degrees. New cracks form, delineating aggregates with a different orientation. Our results show the new crack orientation depends on the refrozen fraction of the initial faults: as this fraction increases, the new cracks gradually rotate to the new wind direction, reaching 90 degrees for fully refrozen faults. Such reorientation is determined by a competition between new cracks forming at a preferential angle determined by the wing-crack theory and at old cracks oriented at a less favourable angle but having higher stresses due to shorter contacts across the joints.
Floe formation in Arctic sea ice
,Dimethyl sulfide (DMS) production in the northern Arctic Ocean has been considered to be minimal because of high sea ice concentration and extremely low productivity. However, we found DMS concentration (1-33 nM) in melt ponds on sea ice at a very high latitude (78°N) in the central Arctic Ocean to be up to ten times that in the adjacent open ocean (&lt;3 nM). We divided melt ponds into three categories: freshwater melt ponds, brackish melt ponds, and open saline melt ponds. Melt ponds from each category had different formation mechanisms and associated DMS contents. Closed brackish ponds (salinity of &gt;20) had particularly high DMS concentration. Water in brackish ponds was mixed with open ocean water in the past via a hole at the bottom of the floe that kept the pond open to the ocean; therefore, unlike freshwater melt ponds, brackish ponds became sites of DMS accumulation. Our results suggest that continuous increase in melt pond coverage on Arctic sea ice could considerably impact future Arctic climate as well as enhancing DMS concentration in the Arctic atmosphere.
Ice load on floating structure simulated with dilated polyhedral discrete element method in broken ice field
,
Limit mechanisms for ice loads on inclined structures: Buckling
,APETALA2/ETHYLENE RESPONSE FACTOR transcription factors (AP2/ERFs) play crucial roles in adaptation to stresses such as those caused by pathogens, wounding and cold. Although their name suggests a specific role in ethylene signalling, some ERF members also co-ordinate signals regulated by other key plant stress hormones such as jasmonate, abscisic acid and salicylate. We analysed a set of ERF proteins from three divergent plant species for intrinsically disorder regions containing conserved segments involved in protein-protein interaction known as Molecular Recognition Features (MoRFs). Then we correlated the MoRFs identified with a number of known functional features where these could be identified. Our analyses suggest that MoRFs, with plasticity in their disordered surroundings, are highly functional and may have been shuffled between related protein families driven by selection. A particularly important role may be played by the alpha helical component of the structured DNA binding domain to permit specificity. We also present examples of computationally identified MoRFs that have no known function and provide a valuable conceptual framework to link both disordered and ordered structural features within this family to diverse function.
Discrete element simulation of the resistance of a ship in unconsolidated ridges
,APETALA2/ETHYLENE RESPONSE FACTOR transcription factors (AP2/ERFs) play crucial roles in adaptation to stresses such as those caused by pathogens, wounding and cold. Although their name suggests a specific role in ethylene signalling, some ERF members also co-ordinate signals regulated by other key plant stress hormones such as jasmonate, abscisic acid and salicylate. We analysed a set of ERF proteins from three divergent plant species for intrinsically disorder regions containing conserved segments involved in protein-protein interaction known as Molecular Recognition Features (MoRFs). Then we correlated the MoRFs identified with a number of known functional features where these could be identified. Our analyses suggest that MoRFs, with plasticity in their disordered surroundings, are highly functional and may have been shuffled between related protein families driven by selection. A particularly important role may be played by the alpha helical component of the structured DNA binding domain to permit specificity. We also present examples of computationally identified MoRFs that have no known function and provide a valuable conceptual framework to link both disordered and ordered structural features within this family to diverse function.
Polygon-based contact description for modeling arbitrary polyhedra in the discrete element method
,
堆石颗粒压缩破碎强度的尺寸效应
,
Size effect on the crushing strengths of rock particles
A shrunken edge algorithm for contact detection between convex polyhedral blocks
,
石英玻璃圆环高速膨胀碎裂过程的离散元模拟
,
Discrete element simulations of the high velocity expensionand fragmentation of qartz glass rings
Bond and fracture model in dilated polyhedral DEM and its application to simulate breakage of brittle materials
,
Minkowski-Voronoi diagrams as a method to generate random packings of spheropolygons for the simulation of soils
,Minkowski operators (dilation and erosion of sets in vector spaces) have been extensively used in computer graphics, image processing to analyze the structure of materials, and more recently in molecular dynamics. Here, we apply those mathematical concepts to extend the discrete element method to simulate granular materials with complex-shaped particles. The Voronoi-Minkowski diagrams are introduced to generate random packings of complex-shaped particles with tunable particle roundness. Contact forces and potentials are calculated in terms of distances instead of overlaps. By using the Verlet method to detect neighborhood, we achieve CPU times that grow linearly with the body's number of sides. Simulations of dissipative granular materials under shear demonstrate that the method maintains conservation of energy in accord with the first law of thermodynamics. A series of simulations for biaxial test, shear band formation, hysteretic behavior, and ratcheting show that the model can reproduce the main features of real granular-soil behavior.
Breaking processes in three-dimensional bonded granular materials with general shapes
,The paper presents an extension to the spheropolyhedra method for the simulation of granular materials comprising particles of general shapes with bonding. A bonding, cement, or cohesion model for particles sharing common faces is introduced. The bonding force is elastic and has a strain-based breaking threshold for modelling fracture. An initial study is conducted based on the Brazilian tensile test to check how the parameters of the proposed model affect the principal variables measured in this test. Afterwards, solid cubic blocks are then subjected to a triaxial test to explore the mathematical macroscopic failure model. It is found that the peak strength envelope is the product of the superposition of frictional and fracture failure mechanisms. The fracture failure is mainly produced by an avalanche of broken cohesive bonds. The intensity of the avalanche exhibits a power law distribution, as reported in previous studies. The method allows for random divisions of solid bodies without any pre-existing internal voids. It offers a natural, effective tool to model, simulate and study fragmentation processes in 3D. (C) 2011 Elsevier B.V.
A calibration methodology to obtain material parameters for the representation of fracture mechanics based on discrete element simulations
,
A coupled discrete element Lattice Boltzmann method for the simulation of fluid-solid interaction with particles of general shapes
,The main aim of this article is to compare various heating methods of the shrink-fit connection and to analyse its shrinkage mechanism. The shrink fit connection consists of the crank and the pin or the ring and the pin, which have negative value of the geometrical clearance. The experimental part of research included the measurement of the temperature at the various points in the surface of the crank, the measurement of the residual stress 5 mm from the crank hole and the measurement of the maximum torsional moment, at which sliding of the shrink-fit connection occurs. Conformity between data from the Finite Element Analysis and the experimental data was obtained in relation to the von Mises stress. The analysis results show that Lamé equations can be successfully used for calculation of the load capacity even for parts with the irregular shape.
尼龙粉末在SLS 预热温度下的离散元模型参数确定及其流动特性分析
,
Determination of discrete element model contact parameters of nylon powder at SLS preheating temperature and its flow characteristics
Potential particles for the modelling of interlocking media in three dimensions
,
Potential particles: A method for modelling non-circular particles in DEM
,Active development of quantum informational components such as quantum computers and quantum key distribution systems requires parameter characterization of single photon detectors. A key property of the single photon detectors is detection efficiency. One of the methods of the detection efficiency measurement, as listed in the international standard ETSI, is the reference-free twin-photon-based Klyshko method. The signal-to-noise ratio (SNR) of this method depends on the combination of the pump wavelength, the nonlinear crystal's axis angle, and the type of detector's sensitive element. When the combination is difficult, one has to deal with the low SNR of the detector counts measurement. To gain the high SNR, one has to average the long record complicated with the &quot;random telegraph signal&quot; noise. This type of noise exhibits high spectral density at a zero frequency, where simple averaging works. The heterodyne based method we have proposed is to perform averaging at the higher frequency of the modulation introduced to the standard Klyshko measurement scheme. The method was numerically simulated and experimentally tested. The 14 times improvement in SNR for the proposed method relative to the simple averaging was demonstrated by the numerical simulation and confirmed experimentally.
A new algorithm for contact detection between convex polygonal and polyhedral particles in the discrete element method
,Active development of quantum informational components such as quantum computers and quantum key distribution systems requires parameter characterization of single photon detectors. A key property of the single photon detectors is detection efficiency. One of the methods of the detection efficiency measurement, as listed in the international standard ETSI, is the reference-free twin-photon-based Klyshko method. The signal-to-noise ratio (SNR) of this method depends on the combination of the pump wavelength, the nonlinear crystal's axis angle, and the type of detector's sensitive element. When the combination is difficult, one has to deal with the low SNR of the detector counts measurement. To gain the high SNR, one has to average the long record complicated with the &quot;random telegraph signal&quot; noise. This type of noise exhibits high spectral density at a zero frequency, where simple averaging works. The heterodyne based method we have proposed is to perform averaging at the higher frequency of the modulation introduced to the standard Klyshko measurement scheme. The method was numerically simulated and experimentally tested. The 14 times improvement in SNR for the proposed method relative to the simple averaging was demonstrated by the numerical simulation and confirmed experimentally.
A new contact detection algorithm for three-dimensional non-spherical particles
,
Characterizing the size and shape of sea ice floes
,Monitoring drift ice in the Arctic and Antarctic regions directly and by remote sensing is important for the study of climate, but a unified modeling framework is lacking. Hence, interpretation of the data, as well as the decision of what to measure, represent a challenge for different fields of science. To address this point, we analyzed, using statistical physics tools, satellite images of sea ice from four different locations in both the northern and southern hemispheres, and measured the size and the elongation of ice floes (floating pieces of ice). We find that (i) floe size follows a distribution that can be characterized with good approximation by a single length scale , which we discuss in the framework of stochastic fragmentation models, and (ii) the deviation of their shape from circularity is reproduced with remarkable precision by a geometric model of coalescence by freezing, based on random Voronoi tessellations, with a single free parameter expressing the shape disorder. Although the physical interpretations remain open, this advocates the parameters and as two independent indicators of the environment in the polar regions, which are easily accessible by remote sensing.
Role of cell regularity and relative density on elasto-plastic compression response of random honeycombs generated using Voronoi diagrams
,Although a resistivity saturation (minimum conductivity) is often observed in disordered metallic solids, such phenomena in the corresponding liquids are not known. Here we report a saturation of the electrical resistivity in Zr_{64}Ni_{36} and Cu_{50}Zr_{50} liquids above a dynamical crossover temperature for the viscosity (T_{A}). The measurements were made for the levitated liquids under the microgravity conditions of the International Space Station. Based on recent molecular dynamics simulations, the saturation is likely due to the ineffectiveness of electron-phonon scattering above T_{A} when the phonon lifetime becomes too short compared to the electron relaxation time. This is different from the conventional resistivity saturation mechanisms in solids.
基于扩展多面体包络函数的快速接触搜索算法
,
A fast detection algorithm based on the envelope function of dilated polyhedron
Efficient implementation of superquadric particles in discrete element method within an open-source framework
,
A novel method based on orientation discretization for discrete element modeling of non-spherical particles
,
Flow characteristics and discharge rate of ellipsoidal particles in a flat bottom hopper
,
A theoretical model of collision between soft-spheres with Hertz elastic loading and nonlinear plastic unloading
,
Coefficient of restitution of colliding viscoelastic spheres
,We numerically investigate head-on collisions of isothermal viscoelastic spheres. We find that the restitution coefficient oscillates against the impact speed if the solid viscosity inside the sphere is small enough. We confirm that the oscillation arises from the resonance between the duration of contact and the eigenfrequencies of the sphere. This oscillation disappears if there exists the strong solid viscosity in spheres. We also find that a sinusoidal behavior of the restitution coefficient against the initial phase in the eigenmodes for collisions between a thermally activated sphere and a flat wall. As a result, the restitution coefficient can exceed unity if the impact speed of the colliding sphere is nearly equal to or slower than the thermal speed. We have confirmed the existence of the fluctuation theorem for impact processes through our simulation.
Internal parameters and regime map for soft polydispersed granular materials
,
Static and dynamic stability assessment of slopes or dam foundations using a rigid body-spring element method
,Organic matter is responsible for the generation of hydrocarbons during the thermal maturation of source rock formation. This geochemical process engenders a network of organic hosted pores that governs the flow of hydrocarbons from the organic matter to fractures created during the stimulation of production wells. Therefore, it can be reasonably assumed that predictions of potentially recoverable confined hydrocarbons depend on the geometry of this pore network. Here, we analyze mesoscale structures of three organic porous networks at different thermal maturities. We use electron tomography with subnanometric resolution to characterize their morphology and topology. Our 3D reconstructions confirm the formation of nanopores and reveal increasingly tortuous and connected pore networks in the process of thermal maturation. We then turn the binarized reconstructions into lattice models including information from atomistic simulations to derive mechanical and confined fluid transport properties. Specifically, we highlight the influence of adsorbed fluids on the elastic response. The resulting elastic energy concentrations are localized at the vicinity of macropores at low maturity whereas these concentrations present more homogeneous distributions at higher thermal maturities, due to pores' topology. The lattice models finally allow us to capture the effect of sorption on diffusion mechanisms with a sole input of network geometry. Eventually, we corroborate the dominant impact of diffusion occurring within the connected nanopores, which constitute the limiting factor of confined hydrocarbon transport in source rocks.
A rigid particle model for rock fracture following the voronoi tessellation of the grain structure: formulation and validation
,
A bonded-particle model for rock
,Organic matter is responsible for the generation of hydrocarbons during the thermal maturation of source rock formation. This geochemical process engenders a network of organic hosted pores that governs the flow of hydrocarbons from the organic matter to fractures created during the stimulation of production wells. Therefore, it can be reasonably assumed that predictions of potentially recoverable confined hydrocarbons depend on the geometry of this pore network. Here, we analyze mesoscale structures of three organic porous networks at different thermal maturities. We use electron tomography with subnanometric resolution to characterize their morphology and topology. Our 3D reconstructions confirm the formation of nanopores and reveal increasingly tortuous and connected pore networks in the process of thermal maturation. We then turn the binarized reconstructions into lattice models including information from atomistic simulations to derive mechanical and confined fluid transport properties. Specifically, we highlight the influence of adsorbed fluids on the elastic response. The resulting elastic energy concentrations are localized at the vicinity of macropores at low maturity whereas these concentrations present more homogeneous distributions at higher thermal maturities, due to pores' topology. The lattice models finally allow us to capture the effect of sorption on diffusion mechanisms with a sole input of network geometry. Eventually, we corroborate the dominant impact of diffusion occurring within the connected nanopores, which constitute the limiting factor of confined hydrocarbon transport in source rocks.
Modeling the particle breakage of rockfill materials with the cohesive crack model
,Active development of quantum informational components such as quantum computers and quantum key distribution systems requires parameter characterization of single photon detectors. A key property of the single photon detectors is detection efficiency. One of the methods of the detection efficiency measurement, as listed in the international standard ETSI, is the reference-free twin-photon-based Klyshko method. The signal-to-noise ratio (SNR) of this method depends on the combination of the pump wavelength, the nonlinear crystal's axis angle, and the type of detector's sensitive element. When the combination is difficult, one has to deal with the low SNR of the detector counts measurement. To gain the high SNR, one has to average the long record complicated with the &quot;random telegraph signal&quot; noise. This type of noise exhibits high spectral density at a zero frequency, where simple averaging works. The heterodyne based method we have proposed is to perform averaging at the higher frequency of the modulation introduced to the standard Klyshko measurement scheme. The method was numerically simulated and experimentally tested. The 14 times improvement in SNR for the proposed method relative to the simple averaging was demonstrated by the numerical simulation and confirmed experimentally.
A hybrid approach for modeling of breakable granular materials using combined finite-discrete element method
,
Numerical simulation of breakages of concrete armour units using a three-dimensional fracture model in the context of the combined finite-discrete element method
,Active development of quantum informational components such as quantum computers and quantum key distribution systems requires parameter characterization of single photon detectors. A key property of the single photon detectors is detection efficiency. One of the methods of the detection efficiency measurement, as listed in the international standard ETSI, is the reference-free twin-photon-based Klyshko method. The signal-to-noise ratio (SNR) of this method depends on the combination of the pump wavelength, the nonlinear crystal's axis angle, and the type of detector's sensitive element. When the combination is difficult, one has to deal with the low SNR of the detector counts measurement. To gain the high SNR, one has to average the long record complicated with the &quot;random telegraph signal&quot; noise. This type of noise exhibits high spectral density at a zero frequency, where simple averaging works. The heterodyne based method we have proposed is to perform averaging at the higher frequency of the modulation introduced to the standard Klyshko measurement scheme. The method was numerically simulated and experimentally tested. The 14 times improvement in SNR for the proposed method relative to the simple averaging was demonstrated by the numerical simulation and confirmed experimentally.
A unified potential-based cohesive model of mixed-mode fracture
,
Measurement of mixed-mode delamination fracture toughness of unidirectional glass/epoxy composites with mixed-mode bending apparatus
,
Numerical simulation of mixed-mode progressive delamination in composite materials
,
Collision detection of convex polyhedra on the NVIDIA GPU architecture for the discrete element method
,
Parallel-vector algorithms for particle simulations on shared-memory multiprocessors
,
Incompressible SPH (ISPH) with fast Poisson solver on a GPU
,
Dense packings of the Platonic and Archimedean solids
,
ISO 19906: 2010, Petroleum and natural gas industries-Arctic offshore structures
,
Quantitative analysis of ice sheet failure against an inclined plane
,
Numerical investigations of ship-ice interaction and maneuvering performance in level ice
,
Ice forces acting on towed ship in level ice with straight drift. Part II: Numerical simulation
,A search for narrow low-mass resonances decaying to quark-antiquark pairs is presented. The search is based on proton-proton collision events collected at 13 TeV by the CMS detector at the CERN LHC. The data sample corresponds to an integrated luminosity of 35.9 fb^{-1}, recorded in 2016. The search considers the case where the resonance has high transverse momentum due to initial-state radiation of a hard photon. To study this process, the decay products of the resonance are reconstructed as a single large-radius jet with two-pronged substructure. The signal would be identified as a localized excess in the jet invariant mass spectrum. No evidence for such a resonance is observed in the mass range 10 to 125 GeV. Upper limits at the 95% confidence level are set on the coupling strength of resonances decaying to quark pairs. The results obtained with this photon trigger strategy provide the first direct constraints on quark-antiquark resonance masses below 50 GeV obtained at a hadron collider.
A straightforward method for calculation of ice resistance of ships
//,
Parallel channels’ fracturing mechanism during ice management operations. Part II: Experiment
,APETALA2/ETHYLENE RESPONSE FACTOR transcription factors (AP2/ERFs) play crucial roles in adaptation to stresses such as those caused by pathogens, wounding and cold. Although their name suggests a specific role in ethylene signalling, some ERF members also co-ordinate signals regulated by other key plant stress hormones such as jasmonate, abscisic acid and salicylate. We analysed a set of ERF proteins from three divergent plant species for intrinsically disorder regions containing conserved segments involved in protein-protein interaction known as Molecular Recognition Features (MoRFs). Then we correlated the MoRFs identified with a number of known functional features where these could be identified. Our analyses suggest that MoRFs, with plasticity in their disordered surroundings, are highly functional and may have been shuffled between related protein families driven by selection. A particularly important role may be played by the alpha helical component of the structured DNA binding domain to permit specificity. We also present examples of computationally identified MoRFs that have no known function and provide a valuable conceptual framework to link both disordered and ordered structural features within this family to diverse function.
Parallel channels’ fracturing mechanism during ice management operations. Part I: Theory
,APETALA2/ETHYLENE RESPONSE FACTOR transcription factors (AP2/ERFs) play crucial roles in adaptation to stresses such as those caused by pathogens, wounding and cold. Although their name suggests a specific role in ethylene signalling, some ERF members also co-ordinate signals regulated by other key plant stress hormones such as jasmonate, abscisic acid and salicylate. We analysed a set of ERF proteins from three divergent plant species for intrinsically disorder regions containing conserved segments involved in protein-protein interaction known as Molecular Recognition Features (MoRFs). Then we correlated the MoRFs identified with a number of known functional features where these could be identified. Our analyses suggest that MoRFs, with plasticity in their disordered surroundings, are highly functional and may have been shuffled between related protein families driven by selection. A particularly important role may be played by the alpha helical component of the structured DNA binding domain to permit specificity. We also present examples of computationally identified MoRFs that have no known function and provide a valuable conceptual framework to link both disordered and ordered structural features within this family to diverse function.
/
〈 | 〉 |