优化端点条件的平面二次均匀B样条插值曲线

在利用反求法构造B样条插值曲线时，往往需要选取端点条件。 因此，可对端点条件进行优化选取，使得构造的B样条插值曲线满足特定要求。提出了一种利用曲线内能极小选取平面二次均匀B样条插值曲线端点条件的算法。首先给出了二次均匀B样条插值曲线分控制顶点与首个控制顶点（即端点条件）的递推关系式；然后给出了利用曲线内能极小优化选取首个控制顶点的算法，证明了利用该算法构造的C¹连续二次均匀B样条插值曲线为保形插值，并通过数值算例证明了算法的有效性；最后，为便于实际应用，基于MATLAB平台设计了算法所对应的图形用户界面，用户通过简单的操作即可获得光顺的C¹连续二次均匀B样条保形插值曲线。     

基于移动可变形孔洞方法的超弹性结构拓扑优化

现有隐式拓扑优化方法在进行超弹性结构拓扑优化设计时，具有设计变量多、中间设计有限元分析存在严重的收敛性和设计结果无法直接导入CAD/CAE系统等问题。为解决这些问题，提出了一种基于移动可变形孔洞的显式拓扑优化方法来进行承受大变形的超弹性结构设计，材料本构采用常用的Mooney-Rivlin模型。首先，介绍了移动可变形孔洞方法的基本思想和可变形孔洞的显式描述方法；其次，构造了基于移动可变形孔洞方法的超弹性结构拓扑优化的数学列式，给出了相应的灵敏度结果；最后，通过数值算例验证了本方法的有效性。数值结果表明，该方法可以通过较少的设计变量和非常稳健的优化过程，给出边界由B样条曲线描述且可与CAD/CAE软件无缝连接的超弹性结构设计。     

Statistical and Geometrical Way of Model Selection for a Family of Subdivision Schemes

The objective of this article is to introduce a generalized algorithm to produce the m-point n-ary approximating subdivision schemes(for any integer m, n ≥ 2). The proposed algorithm has been derived from uniform B-spline blending functions. In particular, we study statistical and geometrical/traditional methods for the model selection and assessment for selecting a subdivision curve from the proposed family of schemes to model noisy and noisy free data. Moreover, we also discuss the deviation of subdivision curves generated by proposed family of schemes from convex polygonal curve. Furthermore, visual performances of the schemes have been presented to compare numerically the Gibbs oscillations with the existing family of schemes.     

Bounds on partial derivatives of NURBS surfaces

In this paper, we estimate the partial derivative bounds for Non-Uniform Rational B-spline(NURBS) surfaces. Firstly, based on the formula of translating the product into sum of B-spline functions, discrete B-spline theory and Dir function, some derivative bounds on NURBS curves are provided. Then, the derivative bounds on the magnitudes of NURBS surfaces are proposed by regarding a rational surface as the locus of a rational curve. Finally, some numerical examples are provided to elucidate how tight the bounds are.     

最小二乘B样条函数在火箭发动机数据处理中的应用

在固体火箭发动机地面试验数据处理中,通过应用最小二乘三次B样条函数拟合的方法,进行野值辨识和修复。文中阐述了实际应用中的步骤及应注意的几点问题。为避免常规的控制序列中含有野值而影响数据辨识的问题,提出差别选取控制点循环过滤的方法,实现精度更高的野值辨识。进而采用门限判别方法,实现野值的修复。通过实测数据处理的应用表明,该算法对单个和成片野值的辨识度都很高,适用于火箭发动机地面试验数据处理。     

Computational aeroacoustics using the B-spline collocation method

One of the major problems in computational aero-acoustics is the disparity in length scales between the flow field and the acoustic field. As a result, a mapping function is normally used to achieve a non-uniform grid distribution. In this paper, a B-spline collocation method with an arbitrary grid placement capability is proposed. This capability not only allows an optimum grid distribution but also avoids the numerical complexities associated with the mapping function. The B-spline collocation method is applied to the case of spinning co-rotating vortices. The result agrees well with the matched asymptotic solution. To cite this article: R. Widjaja et al., C. R. Mecanique 333 (2005).     

THE CONSTRUCTION OF WAVELET-BASED TRUNCATED CONICAL SHELL ELEMENT USING B-SPLINE WAVELET ON THE INTERVAL

Based on B-spline wavelet on the interval (BSWI), two classes of truncated conicalshell elements were constructed to solve axisymmetric problems, i.e. BSWI thin truncated conicalshell element and BSWI moderately thick truncated conical shell element with independent slope-deformation interpolation. In the construction of wavelet-based element, instead of traditionalpolynomial interpolation, the scaling functions of BSWI were employed to form the shape functionsthrough the constructed elemental transformation matrix,and then construct BSWI element viathe variational principle. Unlike the process of direct wavelets adding in the wavelet Galerkinmethod, the elemental displacement field represented by the coefficients of wavelets expansionwas transformed into edges and internal modes via the constructed transformation matrix. BSWIelement combines the accuracy of B-spline function approximation and various wavelet-basedelements for structural analysis. Some static and dynamic numerical examples of conical shellswere studied to demonstrate the present element with higher efficiency and precision than thetraditional element.     

Full-field speckle pattern image correlation with B-Spline deformation function

A full-field speckle pattern image correlation method is presented that will determine directly the complete, two-dimensional deformation field during the image correlation process on digital images obtained using computer vision systems. In this work, a B-Spline function is used to represent the object deformation field throughout the entire image area. This is an improvement over subset-based image correlation methods by implicitly maintaining position and derivative continuity constraints among subsets up to a specified order. The control point variables within the B-Spline deformation function are optimized iteratively with the Levenberg-Marquardt method to achieve minimum disparity between the predicted and actual deformed images. Results have shown that the proposed method is computationally efficient, accurate and robust. The general framework of this method can be applied ton-dimensional image correlation systems that solve for multi-dimension vector fields.     

Nonlinear stability of large k-value shallow spherical shells with a symmetrically distributed line load

In this paper, we treat the nonlinear stability problem of shallow spherical shells with large values ofk(k=12(1–v) · 2f/h,f = shell rise,h = shell thickness) under the action of uniformly distributed line load along a circle concentric with the shell boundary. Load-deflection curves are computed at successive increments of uniformly distributed line loads by using both cubic B-spline approximations and iterative techniques. Our algorithm yields fairly good convergent results for values ofk as large as 400. The limiting case in which shells are loaded along a circle of small radius has been specially investigated and the computed critical loads are compared with those obtained with central point loads by other authors.