Automatic extraction method of force chain information and its application in the flow photoelastic experiment of granular matter

The force chain is the core of the multi-scale analysis of granular matter. Accurately extracting the force chain information among particles is of great significance to the study of particle mechanics and geological hazards caused by particle flow. However, in the photoelastic experiment, the precise identification of the branching points of force chains has not been effectively realized. Therefore, this study proposes an automatic extraction method of force chain key information. First, based on the Hough transform and the Euclidean distance, a particle geometric information identification model is established and geometric information such as particle circle center coordinates, radius, contact point location, and contact angle is extracted. Then, a particle contact force information identification model is established following the color gradient mean square method. The model realizes the rapid calibration and extraction of a large number of particle media contact force information. Next, combined with the force chain composition criterion and its quasilinear feature, an automatic extraction method of force chain information is established, which solves the problem of the accurate identification of the force chain branch points. Finally, in the photoelastic experiment of ore drawing from a single drawpoint, the automatic extraction method of force chain information is verified. The results show that the macroscopic distribution of force chains during ore drawing from a single drawpoint is left–right symmetrical. Strong force chains are mostly located on the two sides of the model but in small numbers and they mainly develop vertically. Additionally, the ends are mostly in a combination of Y and inverted Y shapes, while the middle is mostly quasilinear. Weak force chains are abundant and mostly distributed in the middle of the model, and develop in different directions. The proposed extraction method accurately extracts the force chain network from the photoelastic experiment images and dynamically characterizes the force chains of granular matter, which has significant advantages in particle geometry information extraction, force chain branch point discrimination, force chain retrieval, and force chain distribution and its azimuthal characterization. The results provide a scientific basis for studying the macroscopic and microscopic mechanical parameters of granular matter.     

Development and assessment of a single-image fringe projection method for dynamic applications

A single-image fringe projection profiling method suitable for dynamic applications was developed by combining an accurate camera calibration procedure and improved phase extraction procedures. The improved phase extraction process used a modified Hilbert transform with Laplacian pyramid algorithms to improve measurement accuracy. The camera calibration method used an accurate pinhole camera model and pixel-by-pixel calibration of the phase-height relationship. Numerical simulations and controlled baseline experiments were performed to quantify key error sources in the measurement process and verify the accuracy of the approach. Results from numerical simulations indicate that the resulting phase error can be reduced to less than 0.02 radians provided that parameters such as fringe spacing, random measured intensity noise, fringe contrast and frequency of spatial intensity noise are carefully controlled. Experimental results show that the effects of random temporal and spatial noise in typical CCD cameras for single fringe images limits the accuracy of the method to 0.04 radians in most applications. Quantitative results from application of the fringe projection method are in very good agreement with numerical predictions, demonstrating that it is possible to design both a fringe projection system and a measurement process to achieve a prespecified accuracy and resolution in the point-to-point measurement of the spatial (X, Y, Z) positions.     

粘性流基于特征线的四阶Runge-Kutta有限元法

对于二维不可压缩粘性流,通过沿流线方向的坐标变换,推导了无对流项的二维N-S（Navier-Stokes）方程。采用四阶Runge-Kutta法对N-S方程进行时间离散,并沿流线进行Taylor展开,得到显式的时间离散格式,然后利用Galerkin法对其进行空间离散,得到了高精度的有限元算法。利用本文算法对方腔驱动流和圆柱绕流进行了数值计算,通过对时间步长、网格尺寸和流场区域的计算分析,进一步验证了本文算法相比经典CBS法在时间步长、收敛性、耗散性和计算精度方面更具有优势。     

Two-grid partition of unity method for second order elliptic problems

A two-grid partition of unity method for second order elliptic problems is proposed and analyzed. The standard two-grid method is a local and parallel method usually leading to a discontinuous solution in the entire computational domain. Partition of unity method is employed to glue all the local solutions together to get the global continuous one, which is optimal in HI-norm. Furthermore, it is shown that the L^2 error can be improved by using the coarse grid correction. Numerical experiments are reported to support the theoretical results.     

裂纹问题的一致性高阶无网格法

一致性高阶无网格法能高效精确地求解连续体问题,尤其是能得到高精度的应力场。本文将该方法拓展到应力解析精度至关重要的裂纹问题（即非连续体问题）的数值分析。采用背景积分网格描述裂纹几何,基于无需增加节点额外自由度的虚拟节点法描述裂纹处位移场的间断,提出了虚拟节点的引入算法和断裂单元的数值积分方法。为进一步模拟裂纹扩展,采用相互作用积分方法计算应力强度因子,裂纹的扩展方向由最大周向应力准则确定。数值结果表明,本文发展方法能够精确地通过间断分片试验;相较于标准的高阶无网格法和低阶一致性无网格法,本文的一致性高阶无网格法显著改善了应力强度因子的计算精度,能够准确预测裂纹扩展路径。     

A convenient look-up-table based method for the compensation of non-linear error in digital fringe projection

Although the structured light system that uses digital fringe projection has been widely implemented in three-dimensional surface profile measurement, the measurement system is susceptible to non-linear error. In this work, we propose a convenient look-up-table-based(LUT-based) method to compensate for the non-linear error in captured fringe patterns. Without extra calibration, this LUT-based method completely utilizes the captured fringe pattern by recording the full-field differences. Then, a phase compensation map is established to revise the measured phase. Experimental results demonstrate that this method works effectively.     

Laguerre-Gauss collocation method for initial value problems of second order ODEs

This paper proposes a new collocation method for initial value problems of second order ODEs based on the Laguerre-Gauss interpolation. It provides the global numerical solutions and possesses the spectral accuracy. Numerical results demonstrate its high efficiency.     

A high‐order alternating direction implicit method for the unsteady convection‐dominated diffusion problem

A high‐order alternating direction implicit (ADI) method for solving the unsteady convection‐dominated diffusion equation is developed. The fourth‐order Padé scheme is used for the discretization of the convection terms, while the second‐order Padé scheme is used for the diffusion terms. The Crank–Nicolson scheme and ADI factorization are applied for time integration. After ADI factorization, the two‐dimensional problem becomes a sequence of one‐dimensional problems. The solution procedure consists of multiple use of a one‐dimensional tridiagonal matrix algorithm that produces a computationally cost‐effective solver. Von Neumann stability analysis is performed to show that the method is unconditionally stable. An unsteady two‐dimensional problem concerning convection‐dominated propagation of a Gaussian pulse is studied to test its numerical accuracy and compare it to other high‐order ADI methods. The results show that the overall numerical accuracy can reach third or fourth order for the convection‐dominated diffusion equation depending on the magnitude of diffusivity, while the computational cost is much lower than other high‐order numerical methods. Copyright © 2011 John Wiley & Sons, Ltd.     

充分发展的QDS方法二阶格式及其应用

介绍了QDS方法Euler求解器的基本思路,并从以下三个方面做了改进：（1）在统计网格宏观量时,考虑了相邻网格内流动量的梯度,发展成了二阶格式,使得在QDS算法中的粒子重构和网格宏观量统计都达到二阶格式,由此为充分发展的二阶格式;（2）在进行网格宏观量统计时,提出将按网格进行宏观量统计修改为按粒子位置更新网格宏观量,减少了对粒子进行遍历的次数,大大提高了计算效率;（3）在改进后的二阶格式中引进了斜率限制器,比较了不同斜率限制器的效果。     

高阶边界元奇异积分的一种通用高效计算方法

针对三维边界元法中曲面单元上的(弱、强、超)奇异积分提出了一种通用高效的计算方法。经极坐标变换,将奇异积分转化为常规积分;采用数值方法计算Cauchy主值积分和Hadamard有限项积分系数;引入保角变换和反曲变换消除因单元畸形或因积分点靠近单元边界而引起的周向积分奇异性。该方法可以统一处理(弱、强、超)奇异积分,并且只需要知道核函数的奇异阶数和少数几个点上的被积函数值,不依赖于积分和函数的具体选取;所需的积分点少,精度高,并且受单元畸形程度影响较小,稳定性好。采用该方法计算了声学和弹性力学中的典型奇异积分,并结合二阶Nystrm方法求解了弹性力学的边界积分方程,验证了方法的高精度和高效性。本文数值积分程序可向作者索取。     

显微液下云纹干涉测试系统

本文提出了一种具有超高灵敏度的显微液下云纹干涉法。首先使激光通过折射率ｎ＝１．６５７的液体，用氦氖激光即可获得２０８ｎｍ／条的位移场；其次，用高精度的晶体相移器配合微机图像采集和处理倍增技术，可获得５２ｎｍ／条的位移场，相当于１９２００线／ｍｍ的虚光栅，本系统光路简单，易于调整     

多介质流体力学计算的谱体积方法

针对高维及多物理耦合计算耗费大等困难,设计适合多介质流动模拟的模板紧致、易于并行、高阶精度、计算耗费小的谱体积方法。该方法是求解双曲型守恒率谱体积方法的直接推广,针对多介质流动物质界面捕捉的困难,利用拟守恒格式的思想避免物质界面处的非物理振荡。数值模拟结果表明,本方法具有高阶精度、高分辨率,且节约计算量,并且可以有效避免物质界面处非物理振荡。     

求解温盐双扩散系统的一种高精度方法

引入边界拟合坐标系来研究温盐双扩散系统。为了提高求解的精确性 ,对流项采用四阶精度的迎风格式 ,扩散项和涡量方程的浮力项则采用四阶精度的中心差 ,因此本文的方法是高精度的方法。首先针对温度占优( Rρ=0 .32 )和盐度占优 ( Rρ=1 .6 8)的情形进行了验证性计算 ,得到了与前人一致的结果。进一步 ,本文系统给出了不同的盐度通量强度下的流动形态 ,包括对称结构 ,不对称结构 ,反转结构等 ,结果与前人的吻合。     

High order multiplication perturbation method for singular perturbation problems

This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singularly perturbed two-point boundary value problems are first transformed into the singularly perturbed initial value problems. With the variable coefficient dimensional expanding, the non-homogeneous ordinary dif- ferential equations （ODEs） are transformed into the homogeneous ODEs, which are then solved by the high order multiplication perturbation method. Some linear and nonlinear numerical examples show that the proposed method has high precision.     

Coupling of high order multiplication perturbation method and reduction method for variable coefcient singular perturbation problems

Based on the precise integration method(PIM), a coupling technique of the high order multiplication perturbation method(HOMPM) and the reduction method is proposed to solve variable coefcient singularly perturbed two-point boundary value problems(TPBVPs) with one boundary layer. First, the inhomogeneous ordinary diferential equations(ODEs) are transformed into the homogeneous ODEs by variable coefcient dimensional expansion. Then, the whole interval is divided evenly, and the transfer matrix in each sub-interval is worked out through the HOMPM. Finally, a group of algebraic equations are given based on the relationship between the neighboring sub-intervals, which are solved by the reduction method. Numerical results show that the present method is highly efcient.     

Coupling of high order multiplication perturbation method and reduction method for variable coefficient singular perturbation problems

Based on the precise integration method （PIM）, a coupling technique of the high order multiplication perturbation method （HOMPM） and the reduction method is proposed to solve variable coefficient singularly perturbed two-point boundary value prob lems （TPBVPs） with one boundary layer. First, the inhomogeneous ordinary differential equations （ODEs） are transformed into the homogeneous ODEs by variable coefficient dimensional expansion. Then, the whole interval is divided evenly, and the transfer ma trix in each sub-interval is worked out through the HOMPM. Finally, a group of algebraic equations are given based on the relationship between the neighboring sub-intervals, which are solved by the reduction method. Numerical results show that the present method is highly efficient.     

The micro-connection of the second-order moiré fringes and its order determination,and the determination of the boundary strains

The micro-connection method for determining the centre lines of second-order moiré fringes presented in this paper can enhance the accuracy in measuring strain fields, and generally can determine the absolute order of the second-order moiré fringes. The strain data obtained from curved beam experiment are in good accordance with theory. The relationship between the second-order moiré fringe and the strains at the specific points of the specimen is derived. Hence a reciprocate shift method is presented for determining strain distributions in the non-overlapping region of the shifted moiré patterns, (usually in the region of specimen where no second-order moiré fringes occur is the boundary region).     

空间目标初始轨道确定吉布斯问题误差分析

现代战争中的空间武器平台捕获敌方的空间目标,确定其运行轨道后实施军事打击,所以对轨道确定参数的误差进行分析对于精确打击具有非常重要的意义。为了准确分析吉布斯三位置矢量定轨法的轨道根数的确定精度,根据吉布斯方法的轨道确定模型,采用向量求导的方法,详细推导了轨道根数对于观测位置误差的灵敏度矩阵,并给出了轨道根数误差与观测位置误差之间的关系,明确初始轨道的确定精度及误差传播规律。最后对吉布斯方法及其误差分析进行仿真分析验证,验证了方法的正确性。     

Photoelastic experimental hybrid method for fracture mechanics of anisotropic materials

In this paper, the isotropic and anisotropic photoelastic experimental hybrid methods for fracture mechanics are developed. Using the photoelastic experimental hybrid method, it is demonstrated that one can precisely obtain stress intensity factors and separate the stress components of isotropic and anisotropic plate problems from the only isochromatics.