自由单元法及其在结构分析中的应用

通过吸收有限元与无网格法的优点,提出了一种新的数值方法------自由单元法.此方法在离散方面,采用有限元法中的等参单元,表征几何形状和进行物理量的插值;在算法方面,采用单元配点技术,逐点产生系统方程.主要特点是,在每个配置点只需要一个和周围自由选择的节点而形成的一个独立的等参单元,因而不需要考虑物理量在单元之间的相互连接关系与导数连续性问题. 本文介绍强形式与弱形式两种自由单元法,前者直接由控制方程和边界条件直接产生系统方程,后者通过在自由单元上建立控制方程的加权余量式产生弱形式积分式,并通过像传统有限元法中的积分过程建立系统方程组.本文提出的方法是一种单元配点法,对于域内点为了获得较高的导数精度,需要采用至少具有一个内部点的等参单元,为此除了可使用各阶次的拉格朗日四边形单元外, 还 给出了七节点三角形等参单元,用于模拟较为复杂的几何形状问题.      

New regularization method and iteratively reweighted algorithm for sparse vector recovery

Motivated by the study of regularization for sparse problems,we propose a new regularization method for sparse vector recovery.We derive sufficient conditions on the well-posedness of the new regularization,and design an iterative algorithm,namely the iteratively reweighted algorithm(IR-algorithm),for efficiently computing the sparse solutions to the proposed regularization model.The convergence of the IR-algorithm and the setting of the regularization parameters are analyzed at length.Finally,we present numerical examples to illustrate the features of the new regularization and algorithm.     

Mesh regularization for an ALE code based on the limitation of the Lagrangian mesh velocity

The Lagrangian approach is usually used for the simulation of flow with strong shock waves. Moreover, this approach is particularly well suited to treatment of material interfaces in the case of multimaterial flows.Unfortunately, this formulation leads to very large deformations in the mesh. The arbitrary Lagrangian‐Eulerian method overcomes this drawback by using a mesh regularization that is based on an analysis of cell geometry. The regularization step may be considered as a method used to correct the nonconvex and potentially tangled cells that constitute the mesh. In this paper, we present a new approach to mesh regularization. Instead of using a purely geometric criterion, we propose that the mesh evolution is computed on the basis of the flow vorticity. This approach is called the large Eddy limitation method, and it is aimed here to be used in finite volume direct arbitrary Lagrangian‐Eulerian methods. The large Eddy limitation method is general, which means that it is not restricted to applications in the finite volume framework dedicated to fluid flow simulation; for instance, it could also be naturally applied to the finite element framework.     

Stress state of flexible plates with a hole that are subject to bending

This paper proposes and mathematically proves a method to regularize singular iterations used to solve geometrically nonlinear problems of bending of thin plates with a hole. It is established that the new approximations for deflections and stress functions satisfy the boundary conditions exactly and Karman’s equations asymptotically exactly. The regularization method and the approximations for solutions are based on new methods of summing functional series and sequences __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 1, pp. 100–116, January 2007.     

热结构瞬态响应的耦合灵敏度分析方法与优化设计

研究结构瞬态热变形和热应力的灵敏度分析方法及其优化设计,灵敏度计算给出了直接法和伴随法两种算法.考虑了温度场的耦合作用,在直接法中需要计算温度场对设计变量的导数,在伴随法中需要计算热载荷对温度场的导数.数值算例验证了该方法的精度.伴随法在应用程序中的实现,为大型结构优化提供了高效率的灵敏度计算方法.     

扩散抛物化Navier-Stokes方程数值解法评述

20世纪60年代末期在边界层理论基础上发展起来的各种简化Navier-Stokes (N-S)方程(统称为扩散抛物化N-S方程)及其算法, 较为彻底地解决了无黏流及黏流的相互干扰问题, 并为高雷诺数大型复杂黏性流场的数值模拟开辟了新的途径. 本文将系统地评述这一领域的主要成果, 包括各种简化N-S模型的优缺点; 数学奇性及正则化方法; 代表性的数值解法以及最近几年的新进展.      

Investigation of relaxation properties of polymer melts by comparison of relaxation time spectra calculated with different algorithms

In this paper a recently introduced algorithm (Brabec and Schausberger, 1995) for the calculation of relaxation time spectra is compared with two standard methods, i.e., Weese's regularization, and Baumgaertel's and Winter's regression algorithm. A reasonable agreement between those three algorithms is found for the relaxation properties of mono-, polydisperse, bi-, and multimodal polystyrene samples. All three numerical methods reproduce the relaxation properties for long and medium times correctly, but they show some disagreement at short times because of sparse experimental data. The high numerical accuracy opens the possibility to test and improve the physical models which underlie the calculations. The good agreement of the different algorithms suggests that small inconsistencies to physical models are not due to a failure of the numerical methods, but due to an insufficiency of the generalized Maxwell model.     

An Uncertainty Estimator for Slitting Method Residual Stress Measurements Including the Influence of Regularization

This paper describes the development of a new uncertainty estimator for slitting method residual stress measurements. The new uncertainty estimator accounts for uncertainty in the regularization-based smoothing included in the residual stress calculation procedure, which is called regularization uncertainty. The work describes a means to quantify regularization uncertainty and then, in the context of a numerical experiment, compares estimated uncertainty to known errors. The paper further compares a first-order uncertainty estimate, established by a repeatability experiment, to the new uncertainty estimator and finds good correlation between the two estimates of precision. Furthermore, the work establishes a procedure for automated determination of the regularization parameter value that minimizes total uncertainty. In summary, the work shows that uncertainty in the regularization parameter is a significant contributor to the total uncertainty in slitting method measurements and that the new uncertainty estimator provides a reasonable estimate of single measurement uncertainty.

     

热结构稳态响应的耦合灵敏度分析方法

研究结构稳态热变形和热应力的灵敏度分析方法，给出了直接法和伴随法两种算法。考虑了温度场的耦合作用，在直接法中需要计算温度场对设计变量的导数，在伴随法中需要计算热载荷对温度场的导数。对尺寸和形状两类设计变量的灵敏度分析算例，验证了本文方法的精度。伴随法在应用程序中的实现，为大型结构优化提供了高效率的灵敏度计算方法。     

An unfitted interior penalty discontinuous Galerkin method for incompressible Navier–Stokes two‐phase flow

A discontinuous Galerkin method for the solution of the immiscible and incompressible two‐phase flow problem based on the nonsymmetric interior penalty method is presented. Therefore, the incompressible Navier–Stokes equation is solved for a domain decomposed into two subdomains with different values of viscosity and density as well as a singular surface tension force. On the basis of a piecewise linear approximation of the interface, meshes for both phases are cut out of a structured mesh. The discontinuous finite elements are defined on the resulting Cartesian cut‐cell mesh and may therefore approximate the discontinuities of the pressure and the velocity derivatives across the interface with high accuracy. As the mesh resolves the interface, regularization of the density and viscosity jumps across the interface is not required. This preserves the local conservation property of the velocity field even in the vicinity of the interface and constitutes a significant advantage compared with standard methods that require regularization of these discontinuities and cannot represent the jumps and kinks in pressure and velocity. A powerful subtessellation algorithm is incorporated to allow the usage of standard time integrators (such as Crank–Nicholson) on the time‐dependent mesh. The presented discretization is applicable to both the two‐dimensional and three‐dimensional cases. The performance of our approach is demonstrated by application to a two‐dimensional benchmark problem, allowing for a thorough comparison with other numerical methods. Copyright © 2012 John Wiley & Sons, Ltd.     

MESHLESS METHOD FOR NUMERICAL SIMULATION OF SHOCK-INDUCED COMBUSTION

发展了基于无网格方法的激波诱导燃烧流场数值模拟算法. 该算法采用二维多组分Euler方程,在点云离散的基础上采用曲面逼近计算空间导数,引入多组分HLLC （Harten-Lax-van Leer-contact） 格式计算无黏通量,运用四阶Runge-Kutta 法进行时间显式推进,化学动力学采用有限速率反应模型. 对不同预混气体中的激波诱导燃烧流场进行了数值模拟,结果同相关文献吻合较好,验证了算法的正确性.     

Slender-body theory for viscous flow via dimensional reduction and hyperviscous regularization

A new slender-body theory for viscous flow, based on the concepts of dimensional reduction and hyperviscous regularization, is presented. The geometry of flat, elongated, or point-like rigid bodies immersed in a viscous fluid is approximated by lower-dimensional objects, and a hyperviscous term is added to the flow equation. The hyperviscosity is given by the product of the ordinary viscosity with the square of a length that is shown to play the role of effective thickness of any lower-dimensional object. Explicit solutions of simple problems illustrate how the proposed method is able to represent with good approximation both the velocity field and the drag forces generated by rigid motions of the immersed bodies, in analogy with classical slender-body theories. This approach has the potential to open up the way to more effective computational techniques, since geometrical complexities can be significantly reduced. This, however, is achieved at the expense of involving higher-order derivatives of the velocity field. Importantly, both the dimensional reduction and the hyperviscous regularization, combined with suitable numerical schemes, can be used also in situations where inertia is not negligible.     

Force reconstruction based on the union of singular spectral analysis and regularization method

Abstract

To overcome the problem of some ill-posed inverse problem of force reconstruction, which is caused by the noise in the measured responses and small singular values of the structure, a technology of force reconstruction based on a hybrid method of singular spectral analysis (SSA) and the Landweber regularization method is proposed in this study for the first time. The SSA is used to filter the structural response before using Landweber regularization. A new choice method of phase space reconstruction dimension is theoretically proposed, and the minimum embedding dimension is determined by the concept of optimizing difference spectrum theory. The feasibility of this method was demonstrated through three kinds of force reconstructions. The numerical simulation results and an acoustic vibration experiment demonstrated that the proposed method is more effective than the traditional method.

Communicated by Wei-Chau Xie.     

求解非线性反问题的大范围收敛梯度正则化算法

基于同伦映射的思想,改进了求解非线性反问题的梯度正则化算法。通过路径跟踪有效地拓宽了梯度正则化算法求解的收敛范围。对于正则化参数的修正,通过引入拟Sigmoid函数,提出了一种下降速率可调的连续化参数修正方法,在保证迭代稳定的条件下,得到较好的计算效率,同时保证该算法具有很好的抵抗观测噪声能力。实际算例表明,该方法收敛范围宽,计算效率高,在存在较强观测噪声的条件下也能得到很好的反演结果。     

A new implicit surface tension implementation for interfacial flows

A new implementation of surface tension effects in interfacial flow codes is proposed which is both fully implicit in space, that is the interface never has to be reconstructed, and also semi‐implicit in time, with semi‐implicit referring to the time integration of the surface tension forces. The main idea is to combine two previously separate techniques to yield a new expression for the capillary forces. The first is the continuum surface force (CSF) method, which is used to regularize the discontinuous surface tension force term. The regularization can be elegantly implemented with the use of distance functions, which makes the level set method a suitable choice for the interface‐tracking algorithm. The second is to use a finite element discretization together with the Laplace–Beltrami operator, which enables simple reformulation of the surface tension term into its semi‐implicit equivalent. The performance of the new method is benchmarked against standard explicit methods, where it is shown that the new method is significantly more robust for the chosen test problems when the time steps exceed the numerical capillary time step restriction. Some improvements are also found in the average number of nonlinear iterations and linear multigrid steps taken while solving the momentum equations. Copyright © 2005 John Wiley & Sons, Ltd.     

面载荷识别的TSVD正则化方法

线性系统的响应可以表示为单位脉冲响应函数与动态载荷的卷积分,经离散后得到载荷识别反问题。引入TSVD正则化方法处理反问题的不适定性,对不同噪声的实测响应进行了面载荷的反演求解,并对识别精度做出了分析。算例表明,TSVD正则化方法对响应误差的适应性强,识别精度高。     

基于模态参数灵敏度的损伤方程组求解正则化方法研究

基于模态参数的结构损伤识别方法是振动损伤识别领域中应用最为广泛的方法.利用模态参数灵敏度构建结构损伤方程组,对其进行求解可以识别结构损伤位置和程度.由于实际工程中模态参数不完备性和噪声的影响,结构损伤方程易出现病态问题,直接求解可能产生错误的结果.为了解决这一问题,可以引入正则化方法进行求解.然而,各类正则化方法的基本...     

二维热弹性力学边界元法中几乎奇异积分的正则化

针对二维热弹性力学边界元法中近边界点的几乎强奇异和超奇异积分,采用一种通用算法,将其实施正则化．该方法适用于线性单元,与近边界点邻近的单元上的积分采用正则化积分公式计算,远处单元的积分仍保持常规高斯积分．算例证明了该法的有效性和精确性．     

Damage detection using a new regularization method with variable parameter

In this research, a sensitivity approach to finite element model updating is used to determine stiffness reduction factors from measured structural response. The used method causes a set of nonlinear ill-conditioned equations that need to be linearized and regularized in order to find the solution. A new approach to solve the problem is presented using variable regularization parameter. Utilization of variable regularization parameter eliminates dependency on the number of iterations and prevents the loss of regularization effect due to iterations. A new stopping criteria is used which is based on the difference between mean and variance of last iterations. Furthermore the results show that using wavelet transform to update the model yields better results than modal parameters. Expedient performance of the proposed method is shown through a numerical simulation.     

基于拓扑与形状优化的小型化金属天线设计方法

天线小型化设计需要基于先进的设计方法,基于拓扑优化的设计往往存在灰度单元,因此设计结果无法直接应用,需要进一步规整设计。而对于电磁金属结构,粗糙的规整方法会引起结构性能的很大变化以致偏离最优结果。提出一种拓扑优化和形状优化相结合的方法,用于金属天线结构的小型化设计。该方法通过拓扑优化获得金属天线结构的概念构型,进而利用形状优化对概念构型进行边界规整和精细化设计。形状优化方法采用多控制点贝塞尔曲线描述拓扑概念构型,通过贝塞尔曲线控制点的移动实现天线构型的调控。给出了贝塞尔曲线控制点的设置原则,基于拓扑优化得到场量分布结果,利用较少的贝塞尔曲线控制点实现天线拓扑构型结构特征的有效调控。该方法可以获得无灰度单元残留的拓扑结果,同时可有效避免密度阈值规整方法中天线性能改变的问题,并且获得的拓扑构型边界光滑。数值算例表明拓扑优化和形状优化相结合方法的有效性。此外,该方法可拓展到其他类型电磁器件的优化设计中。