基于数据重建的三角网格模型简化优化方法

由CT切片数据重建得到的三角网格模型常存在数据量大、狭长三角形多等问题,针对这些问题,研究了一种保持特征的高质量三角网格模型的简化优化方法。该方法分为网格简化和网格优化两个阶段。首先,采用基于曲面变化的二次误差度量计算边折叠代价,并按代价值的大小进行迭代的折叠简化,可较好地保持模型表面的特征;其次,通过二阶加权伞算子对简化模型中局部存在的狭长三角形进行优化处理,改善三角网格模型的质量。实验结果表明,该方法能够较好地保持特征区域的细节信息,并可靠地生成高质量、低几何误差的简化模型。     

基于工业CT图像重构的网格模型的保特征简化方法

对基于工业CT图像重构的网格模型进行网格简化时,大多数现有网格模型简化算法会丢失特征,出现网格质量不好的问题。因此提出一种网格模型保特征简化方法,该方法用三角形折叠法对原始模型进行简化,当简化后模型的平均二面角角度误差达到允许误差后,再使用边折叠法对模型进行简化。在三角形折叠法中提出了利用被折叠三角形的法向量、各个顶点的高斯曲率及其在周边三角形上的投影确定该三角形的折叠点,利用局部体积误差与二面角角度误差的无因次化和确定折叠代价的方法;在边折叠法中提出了将二面角角度误差引入到二次误差测度(QEM)法的折叠代价中的改进QEM法。实验结果表明:与其他算法相比,该方法能够生成保特征、高质量、低几何误差的网格模型。     

保特征的自适应三角网格规范化算法

为了提高基于工业CT图像重构的三角网格质量，提出了顶点调整和特征因子相结合的网格规范化算法。引入模型特征因子和局部网格质量提高程度作为网格调整的控制条件，保留原始模型的局部细节特征；采用法矢量阈值按规律递增的方式自适应控制网格调整，实现不同曲率特征的自动识别。实验结果表明，与现有方法相比，该方法能更好地规范模型的三角网格，同时保留了原始模型的细节特征。     

Finite element analysis of the dynamic behavior of a laminated windscreen with frequency dependent viscoelastic core

This paper presents numerical and experimental validation of results obtained by a shell finite element, which has been developed for modeling of the dynamic behavior of sandwich multilayered structures with a viscoelastic core. The proposed shell finite element is very easy to implement in existing finite element solvers, since it uses only the displacements as degrees of freedom at external faces and at inter-layer interfaces. The displacement field is linearly interpolated in the thickness direction of each layer, and analytical integration is made in the thickness direction in order to avoid meshing of each sandwich layer by solid elements. Only the two dimensional mid-surface of reference is meshed, facilitating the mesh generation task. A simplified modal approach using a real modal basis is also proposed to efficiently calculate the dynamic response of the sandwich structure. The proposed method reduces the memory size and computing time and takes into account the frequency-dependence of the polymer core mechanical properties. Results obtained by the proposed element in conjunction with the simplified modal method have been numerically and experimentally validated by comparison to results obtained by commercial software codes (MSC/nastran and ESI/rayon-vtm), and to measurements done on automobile windscreens.     

Finite element prediction of wave motion in structural waveguides

A method is presented by which the wavenumbers for a one-dimensional waveguide can be predicted from a finite element (FE) model. The method involves postprocessing a conventional, but low order, FE model, the mass and stiffness matrices of which are typically found using a conventional FE package. This is in contrast to the most popular previous waveguide/FE approach, sometimes termed the spectral finite element approach, which requires new spectral element matrices to be developed. In the approach described here, a section of the waveguide is modeled using conventional FE software and the dynamic stiffness matrix formed. A periodicity condition is applied, the wavenumbers following from the eigensolution of the resulting transfer matrix. The method is described, estimation of wavenumbers, energy, and group velocity discussed, and numerical examples presented. These concern wave propagation in a beam and a simply supported plate strip, for which analytical solutions exist, and the more complex case of a viscoelastic laminate, which involves postprocessing an ANSYS FE model. The method is seen to yield accurate results for the wavenumbers and group velocities of both propagating and evanescent waves.     

Corotational Finite Element Formulation for Virtual-Reality Based Surgery Simulators

Surgical simulation provides a means for trainees to develop surgical competence that encompasses requisite knowledge, technical and cognitive skills and decision-making ability. Considering virtual-reality based surgery simulators, the key requirement is sufficiently accurate and numerically efficient computation of deformation behavior of soft tissues, which is highly nonlinear. The paper offers a simplified geometrically nonlinear corotational finite element formulation to meet the imposed requirements. The approach is used in combination with a rather simple type of finite element and an appropriate solver is chosen for fast computation of dynamical behavior. The finite element formulation is enriched with a coupled-mesh technique to enable modelling of complex geometries by relatively simple computational models. A few examples of models of internal organs are provided to discuss the aspects of the developed tools.     

Finite element analysis of quasistatic crack propagation in brittle media with voids or inclusions

A three-level finite element scheme is proposed for simulation of crack propagation in heterogeneous media including randomly distributed voids or inclusions. To reduce total degrees of freedom in the view of mesh gradation, the entire domain is categorized into three regions of different-level meshes: a region of coarse-level mesh, a region of intermediate-level mesh, and a region of fine-level mesh. The region of coarse-level mesh is chosen to be far from the crack to treat the material inhomogeneities in the sense of coarse-graining through homogenization, while the region near the crack is composed of the intermediate-level mesh to model the presence of inhomogeneities in detail. Furthermore, the region very near the crack tip is refined into the fine-level mesh to capture a steep gradient of elastic field due to the crack tip singularity. Variable-node finite elements are employed to satisfy the nodal connectivity and compatibility between the neighboring different-level meshes. Local remeshing is needed for readjustment of mesh near the crack tip in accordance with crack growth, and this is automatically made according to preset values of parameters determining the propagation step size of crack, and so the entire process is fully automatic. The effectiveness of the proposed scheme is demonstrated through several numerical examples. Meanwhile, the effect of voids and inclusions on the crack propagation is discussed in terms of T-stresses, with the aid of three-level adaptive scheme.     

自适应间断有限元方法求解三维欧拉方程

将龙格库塔间断有限元方法(RDDG)与自适应方法相结合,求解三维欧拉方程.区域剖分采用非结构四面体网格,依据数值解的变化采用自适应技术对网格进行局部加密或粗化,减少总体网格数目,提高计算效率.给出四种自适应策略并分析不同自适应策略的优缺点.数值算例表明方法的有效性.     

某型内转子光电雷达稳定平台的模态分析与优化

为了更好地设计某型内转子光电雷达稳定平台,运用Pro-E软件对其进行了三维实体建模,合理简化后,进入Ansys软件建立有限元模型.应用有限元理论和Ansys软件对该平台进行了模态分析,得到稳定平台前八阶模态的固有频率和振型.其中对比了不同网格、不同材料产生的不同模态,对同类分析提供了参考,也为此稳定平台的设计优化提供了...     

Finite element modeling of heating phenomena of cracks excited by high-intensity ultrasonic pulses

A three-dimensional thermo-mechanical coupled finite element model is built up to simulate the phenomena of dynamical contact and frictional heating of crack faces when the plate containing the crack is excited by high-intensity ultrasonic pulses.In the finite element model,the high-power ultrasonic transducer is modeled by using a piezoelectric thermal-analogy method,and the dynamical interaction between both crack faces is modeled using a contact-impact theory.In the simulations,the frictional heating taking place at the crack faces is quantitatively calculated by using finite element thermal-structural coupling analysis,especially,the influences of acoustic chaos to plate vibration and crack heating are calculated and analysed in detail.Meanwhile,the related ultrasonic infrared images are also obtained experimentally,and the theoretical simulation results are in agreement with that of the experiments.The results show that,by using the theoretical method,a good simulation of dynamic interaction and friction heating process of the crack faces under non-chaotic or chaotic sound excitation can be obtained.     

基于质心Voronoi图重构的UDSM边折叠简化

城市数字表面模型网格(UDSM)的相邻网格常常出现曲率剧变,而这些位置是UDSM的细节部分,简化过程中应当尽量保持。针对该情况,引入了质心Voronoi图重划分网格,将曲率较小的表面的点云密度大大降低。重划分的网格表面细节与周围的平滑表面的三角网格尺寸悬殊,在该基础上使用二次误差矩阵边折叠进行LOD构建时网格发生明显变化,范围大大减少。算法在时间性能与网格误差与直接边折叠相近的前提下,更多地保存简化后的网格细节。     

High dimensional model representation method for fuzzy structural dynamics

Uncertainty propagation in multi-parameter complex structures possess significant computational challenges. This paper investigates the possibility of using the High Dimensional Model Representation (HDMR) approach when uncertain system parameters are modeled using fuzzy variables. In particular, the application of HDMR is proposed for fuzzy finite element analysis of linear dynamical systems. The HDMR expansion is an efficient formulation for high-dimensional mapping in complex systems if the higher order variable correlations are weak, thereby permitting the input-output relationship behavior to be captured by the terms of low-order. The computational effort to determine the expansion functions using the α-cut method scales polynomically with the number of variables rather than exponentially. This logic is based on the fundamental assumption underlying the HDMR representation that only low-order correlations among the input variables are likely to have significant impacts upon the outputs for most high-dimensional complex systems. The proposed method is first illustrated for multi-parameter nonlinear mathematical test functions with fuzzy variables. The method is then integrated with a commercial finite element software (ADINA). Modal analysis of a simplified aircraft wing with fuzzy parameters has been used to illustrate the generality of the proposed approach. In the numerical examples, triangular membership functions have been used and the results have been validated against direct Monte Carlo simulations. It is shown that using the proposed HDMR approach, the number of finite element function calls can be reduced without significantly compromising the accuracy.     

Accelerated simulations of aromatic polymers: application to polyether ether ketone (PEEK)

For aromatic polymers, the out-of-plane oscillations of aromatic groups limit the maximum accessible time step in a molecular dynamics simulation. We present a systematic approach to removing such high-frequency oscillations from planar groups along aromatic polymer backbones, while preserving the dynamical properties of the system. We consider, as an example, the industrially important polymer, polyether ether ketone (PEEK), and show that this coarse graining technique maintains excellent agreement with the fully flexible all-atom and all-atom rigid bond models whilst allowing the time step to increase fivefold to 5 fs.     

弧形调强放射治疗对剂量计算方法的要求

在弧形调强放射治疗的治疗计划设计中, 由于包含有很多照射方向, 调强最优化的射束元矩阵计算需要很大的计算量和存储量, 为提高计算效率常使用简化剂量计算模型计算射束元矩阵, 因此有必要研究简化模型对治疗计划质量产生影响。 对一个模拟例子和一个临床实例, 使用没考虑散射效应的原射线剂量计算模型计算射束元矩阵, 由此进行最优化计算。 在得到最优化强度分布后, 通过比较原射线剂量计算模型和微分卷积剂量计算模型得到的剂量分布, 研究了不同射束数目条件下, 使用简化剂量计算模型计算射束元剂量矩阵对最终的剂量分布质量的影响。 结果表明, 在射线束很多的情况下(对应弧形调强照射）, 用简化的剂量计算模型, 即不考虑散射来计算射束元剂量矩阵, 会导致靶区剂量分布的质量大大低于预期的剂量分布质量, 因此, 弧形调强放射治疗的最优化计算中, 有效考虑散射的影响是必要的。 In the treatment planning for arc intensity modulated radiation therapy, because many irradiation directions are involved, the computing time and storage space needed for calculating beamlet dose matrices in optimization is quite heavy. In order to improve the computation efficiency, the simplified dose calculation is often used for the calculation of the dose matrices. Thus, it is deserved to study how this simplification could influence the quality of the treatment plan. In this paper, a simulation and a clinical case are adopted. Using the primary dose calculation model without taking into account the scattering effect to generate the dose matrices of beamlets, the optimization for beam intensity profile are firstly carried out. Then, based on the obtained intensity profile, the dose distributions are recalculated by using the primary dose calculation model and the differential convolution superposition dose calculation model which is more accurate but more time consuming. By comparing dose distributions obtained by this two models, the influence of using simplified model for dose matrix calculation on beam profile optimization is studied. The results demonstrate that when the beam number is large(corresponding to the arc modulated radiation), using the simplified model for the calculation of dose matrix of beamlets will reduce the quality of dose distribution greatly comparing with the expected dose distribution quality. Thus it is very necessary to correctly take into account the scattering effect in beam profile optimization for the arc intensity modulated radiation therapy.     

Two-δ- Function Generalized Plasma Pole Model for First Principle Calculations of Quasiparticle Energies in Many-Electron Systems

To evaluate the dynamical effects of the screened interaction in the calculations of quasiparticle energies in many-electron systems a two-δ-function generalized plasma pole model (GPP) is introduced to simulate the dynamical dielectric function. The usual single δ-function GPP model has the drawback of over simplifications and for the crystals without.the center of symmetry is inappropriate to describe the finite frequency behavior for dielectric function matrices. The discrete frequency summation method requires too much computation to achieve converged results since ab initio calculations of dielectric function matrices are to be carried out for many different frequencies. The two-δ GPP model is an op timization of the two approaches. We analyze the two-δ GPP model and propose a method to determine from the first principle calculations the amplitudes and effective frequencies of these δ functions. Analytical solutions are found for the second order equations for the parameter matrices entering the model. This enables realistic applications of the method to the first principle quasiparticle calculations and makes the calculations truly adjustable parameter free.     

A finite element model using multi-layered shell element in laser forming

Due to its enormously high flexibility laser forming has been gaining importance in recent years. This rapidness and flexibility demand very precise controlling strategies especially when simulating the process of large plates and challenging the limited computation power of the current workstation. A simple, robust and accurate modeling method of laser forming has been demonstrated to solve this problem. The simplified model is meshed by multi-layered shell element, simulated with a more real scanning method and fewer parameters. The intelligent meshing strategies have reduced the number of elements dramatically. Thus the simulation efficiency has been improved significantly. By comparing the simulation results under the simplified model with the results under the traditional model for laser forming, the applicability of proposed method has been proven. The method of these simplified models is also suitable to simulate complex finite element models, which take much time to simulate. It would throw some light on the thermal mechanically coupled-field simulation of large sheet.     

Semi-analytical finite element analysis of elastic waveguides subjected to axial loads

Predicting the influence of axial loads on the wave propagation in structures such as rails requires numerical analysis. Conventional three-dimensional finite element analysis has previously been applied to this problem. The process is tedious as it requires that a number of different length models be solved and that the user identify the computed modes of propagation. In this paper, the more specialised semi-analytical finite element method is extended to account for the effect of axial load. The semi-analytical finite element method includes the wave propagation as a complex exponential in the element formulation and therefore only a two-dimensional mesh of the cross-section of the waveguide is required. It was found that the stiffness matrix required to describe the effect of axial load is proportional to the mass matrix, which makes the extension to existing software trivial. The method was verified by application to an aluminium rod, where after phase and group velocities of propagating waves in a rail were computed to demonstrate the method.     

Torsional vibration of crankshaft in an engine–propeller nonlinear dynamical system

Theoretical and experimental studies on torsional vibration of an aircraft engine–propeller system are presented in this paper. Two system models—a rigid body model and a flexible body model, are developed for predicting torsional vibrations of the crankshaft under different engine powers and propeller pitch settings. In the flexible body model, the distributed torsional flexibility and mass moment of inertia of the crankshaft are considered using the finite element method. The nonlinear autonomous equations of motion for the engine–propeller dynamical system are established using the augmented Lagrange equations, and solved using the Runge–Kutta method after a degrees of freedom reduction scheme is applied. Experiments are carried out on a three-cylinder four-stroke engine. Both theoretical and experimental studies reveal that the crankshaft flexibility has significant influence on the system dynamical behavior.     

BGK方法在非结构网格上的应用

采用旋转局部坐标的方法,发展了一种针对非结构网格的BGK计算方法.该方法属于有限体积法,大致分为两个步骤:①空间离散:②通量计算及时间推进.在第①步中,采用基于最小二乘法的高阶ENO格式来获得宏观物理量的高阶导数;在第②步中,采用旋转坐标轴的方法来计算非结构网格单元各边的通量.并得出了后台阶绕流(Backward Facing Step)及翼型绕流(Flow Over an Airfoil)两个算例的计算结果.     

An adaptive high-order minimum action method

In this work, we present an adaptive high-order minimum action method for dynamical systems perturbed by small noise. We use the hp finite element method to approximate the minimal action path and nonlinear conjugate gradient method to solve the optimization problem given by the Freidlin–Wentzell least action principle. The gradient of the discrete action functional is obtained through the functional derivative and the moving mesh technique is employed to enhance the approximation accuracy. Numerical examples are given to demonstrate the efficiency and accuracy of the proposed numerical method.