油藏数值模拟有限元并行计算方法研究

针对目前油藏数值模拟普遍采用的有限差分法计算精度低的问题,提出了兼顾计算精度、计算速度问题的有限元油藏数值模拟方法,即在建立了油藏数值模拟数学模型的基础上通过有限元数值分析方法建立有限元数值模型,但有限元在油藏数值模拟时存在单机计算困难、计算时间长的问题,为此提出了利用区域分解技术的油藏数值模拟并行计算方法,最后将该方法通过实例进行检验,取得了良好的加速比和并行效率。     

Simulink behavioral modeling of a 10-bit pipelined ADC

The increasing architecture complexity of data converters makes it necessary to use behavioral models to simulate their electrical performance and to determine their relevant data features. For this purpose, a specific data converter simulation environment has been developed which allows designers to perform time-domain behavioral simulations of pipelined analog to digital converters (ADCs). All the necessary blocks of this specific simulation environment have been implemented using the popular Matlab simulink environment. The purpose of this paper is to present the behavioral models of these blocks taking into account most of the pipelined ADC non-idealities, such as sampling jitter, noise, and operational amplifier parameters (white noise, finite DC gain, finite bandwidth, slew rate, and saturation voltages). Simulations, using a 10-bit pipelined ADC as a design example, show that in addition to the limits analysis and the electrical features extraction, designers can determine the specifications of the basic blocks in order to meet the given data converter requirements.     

嵌入多层黏弹性胶膜复合材料阻尼工字梁的多目标设计优化

针对嵌入多层黏弹性胶膜的复合材料阻尼工字梁,传统的混合单元法在进行结构动态分析与设计优化时存在很大困难的问题,采用基于离散层理论的多层梁单元建模分析复合材料阻尼工字梁.通过对正交各向异性铺层和腹板进行等效处理的方法,对工字梁凸缘嵌入单层阻尼层模型进行参数化分析;分别对嵌入多层或单层阻尼层的模型建立多目标优化模型,优化目标为模态损耗因子和固有频率最大化,设计变量为阻尼层层数（厚度）及其嵌入位置;应用多目标遗传算法进行优化求解.结果表明：基于离散层理论的阻尼梁单元计算精度好且易于优化,对于嵌入单层较厚阻尼和嵌入多层较薄阻尼的复合材料工字梁,获得的阻尼效果与动刚度损失基本相当,但对于高阻尼的方案,前者比后者的动刚度损失更大.     

Damping optimization of viscoelastic laminated sandwich composite structures

Recent developments on the optimization of passive damping for vibration reduction in sandwich structures are presented in this paper, showing the importance of appropriate finite element models associated with gradient based optimizers for computationally efficient damping maximization programs. A new finite element model for anisotropic laminated plate structures with viscoelastic core and laminated anisotropic face layers has been formulated, using a mixed layerwise approach. The complex modulus approach is used for the viscoelastic material behavior, and the dynamic problem is solved in the frequency domain. Constrained optimization is conducted for the maximization of modal loss factors, using gradient based optimization associated with the developed model, and single and multiobjective optimization based on genetic algorithms using an alternative ABAQUS finite element model. The model has been applied successfully and comparative optimal design applications in sandwich structures are presented and discussed.     

A nonlinear visco-elastic constitutive model for large rotation finite element formulations

Although all known materials have internal damping that leads to energy dissipation, most existing large deformation visco-elastic finite element formulations are based on linear constitutive models or on nonlinear constitutive models that can be used in the framework of an incremental co-rotational finite element solution procedure. In this investigation, a new nonlinear objective visco-elastic constitutive model that can be implemented in non-incremental large rotation and large deformation finite element formulations is developed. This new model is based on developing a simple linear relationship between the damping forces and the rates of deformation vector gradients. The deformation vector gradients can be defined using the decomposition of the matrix of position vector gradients. In this paper, the decomposition associated with the use of the tangent frame that is equivalent to the QR decomposition is employed to define the matrix of deformation gradients that enter into the formulation of the viso-elastic constitutive model developed in this investigation. Using the relationship between the deformation gradients and the components of the Green–Lagrange strain tensor, it is shown that the damping forces depend nonlinearly on the strains and linearly on the classical strain rates. The relationship between the damping forces and strains and their rates is used to develop a new visco-elastic model that satisfies the objectivity requirements and leads to zero strain rates under an arbitrary rigid body displacement. The linear visco-elastic Kelvin–Voigt model frequently used in the literature can be obtained as a special case of the proposed nonlinear model when only two visco-elastic coefficients are used. As demonstrated in this paper, the use of two visco-elastic coefficients only leads to viscous coupling between the deformation gradients. The model developed in this investigation can be used in the framework of large deformation and large rotation non-incremental solution procedure without the need for using existing co-rotational finite element formulations. The finite element absolute nodal coordinate formulation (ANCF) that allows for straightforward implementation of general constitutive material models is used in the validation of the proposed visco-elastic model. A comparison with the linear visco-elastic model is also made in this study. The results obtained in this investigation show that there is a good agreement between the solutions obtained using the proposed nonlinear model and the linear model in the case of small deformations.     

Finite Element Guidelines for Simulation of Delamination Dominated Failures in Composite Materials Validated by Case Studies

The focus of this paper is on the computational modelling of progressive damage in composite structures of fibre reinforced laminae. A general review of modelling approaches to failure in the context of the finite element method is first presented, with an emphasis on models based on continuum damage mechanics. The way in which delamination and matrix splitting (that may or may not interact with fibre-tension damage) should be addressed in the framework of a commercial finite element code is considered next. An important feature of the analysis is it does not rely on customized user-subroutines but solely on the analysis capabilities of the general purpose software Abaqus, thus ensuring that the numerical results can be universally reproduced. It is shown that the finite element simulations can accurately represent the physical mechanisms controlling damage development and progression and reproduce a number of phenomena including delamination, laminate in-plane failure and behaviour at notches. The paper ends giving guidelines for the generalized modelling methodology using Abaqus without user-subroutines.     

Design of an SOI-MEMS high resolution capacitive type single axis accelerometer

In recent years, high resolution micro accelerometers are increasingly finding various applications in different segments of life. The current work deals with the design of a high resolution single axis accelerometer based on SOI-MEMS technology. Accordingly two different approaches for designing comb type, capacitive, SOI-MEMS accelerometer is presented. Initially a system level approach using a simulation platform from SABER is carried out to obtain the basic design. Later, a device level design is carried out by building three dimensional (3D) geometric models using finite element (FE) simulations through CoventorWare software. Different design parameters like mechanical and electrical sensitivity, capacitance values, resonant frequency, etc. are obtained in either of the cases and compared. The design is optimized based on the overall sensitivity and the system noise level both electrical and mechanical, respectively. The complete design is worked out in accordance with the silicon on insulator based multiuser MEMS fabrication processes (MUMPs) technology from MEMSCAP foundry.     

A Digital Input Shaper for Stable and Transparent Haptic Interaction

This paper presents frequency-dependent digital damping referred to a digital input shaper (DIS). Compared to the previously proposed analog input shaper, the DIS makes the implementation of damping much more flexible and significantly reduces the electrical noise. Moreover, the DIS incorporates a position control method at the initial contact to improve the perceived contact hardness. Then, through a series of simulations and benchmark virtual wall experiments, we show that the DIS can significantly increase the impedance range that a haptic interface can stably display and that the initial contact hardness that a user perceives can be noticeably improved by using the proposed position control method. Finally, the DIS is applied to scaled teleoperation using an atomic force microscope and through contact experiments it is shown that it can indeed increase the scaled impedance of a real environment that can be stably displayed. ©     

Design sensitivity analysis of acoustical damping and its application to design of musical bells

Damping characteristics of a musical bell plays an important role in characterizing the musical sound. The total damping consists of acoustical damping and internal damping. Acoustical damping depends upon resonating frequencies and vibration patterns while internal damping is a material property. The acoustical damping of a vibrating structure is formulated via boundary element method and finite element method using eigenmode decomposition. The design sensitivity of acoustical damping is derived using an adjoint variable method of the eigenvalue problem. Design optimization of a musical bell is then performed in terms of acoustical parameters. The goal of the optimization problem is to design a harmonically tuned bell with given acoustical damping values. The proposed automated design process integrates finite element analysis, boundary element analysis, design sensitivity analysis, mode-tracking algorithm and optimization module, seamlessly. It is demonstrated by numerical examples to show practical applications.     

Shape optimization of a vehicle hat-shelf: Improving acoustic properties for different load cases by maximizing first eigenfrequency

This paper presents design optimization of the geometry of a vehicle hat-shelf. At first two existing finite element discretizations are investigated for two different element types. The structural model is then parameterized. Only four design variables have been chosen to control the shape modification of the hat-shelf. The aim of this paper is to decrease the vehicle interior noise due to three different excitations for two cases of fluid damping. With respect to the support conditions of the hat-shelf these three load cases and the two cases of different damping are considered simultaneously by maximizing the lowest eigenfrequency of the structural model. Although remarkable differences in the natural frequencies are discovered for the four discretizations, a similar dependence of the objective function in terms of the design variables is observed. Thus, a multigrid strategy can be applied. The coarsest mesh is used to obtain suitable initial sets of optimization variables, one of the finer meshes serves for a pre-optimization and the finest mesh is optimized to find the final set of parameters. While the lowest eigenfrequency of the original model is found at about 31 Hz, the corresponding value in the optimized variant exceeds 100 Hz being the upper bound of the frequency range under consideration. Evaluation of the noise transfer function proves that this strategy decreases its average between 4.4 and 13.9 dB.     

Error Analysis and Stochastic Modeling of Low-cost MEMS Accelerometer

This paper presents the error analysis and stochastic modeling of commercial low-cost MEMS Accelerometer. Although Micro Electro Mechanical Systems (MEMS) based sensors have been utilized for the development of low-cost integrated navigation systems on the benefits of low inherent cost, small size, low power consumption, and solid reliability, it is significantly important to characterize the error behaviors of MEMS-based sensors and to construct more sophisticated mathematical modeling methods. The errors of MEMS-based accelerometer have been identified into deterministic and stochastic error sources and the stochastic error part was the focus to be discussed in this paper using discrete parameter models of stationary random process. Appropriate Autoregressive (AR) models have been analyzed which can be used to help the development of appropriate optimal algorithm for multiple sensor integration.     

Model updating of multistory shear buildings for simultaneous identification of mass,stiffness and damping matrices using two different soft-computing methods

In this research, two novel methods for simultaneous identification of mass–damping–stiffness of shear buildings are proposed. The first method presents a procedure to estimate the natural frequencies, modal damping ratios, and modal shapes of shear buildings from their forced vibration responses. To estimate the coefficient matrices of a state-space model, an auto-regressive exogenous excitation (ARX) model cooperating with a neural network concept is employed. The modal parameters of the structure are then evaluated from the eigenparameters of the coefficient matrix of the model. Finally, modal parameters are used to identify the physical/structural (i.e., mass, damping, and stiffness) matrices of the structure. In the second method, a direct strategy of physical/structural identification is developed from the dynamic responses of the structure without any eigenvalue analysis or optimization processes that are usually necessary in inverse problems. This method modifies the governing equations of motion based on relative responses of consecutive stories such that the new set of equations can be implemented in a cluster of artificial neural networks. The number of neural networks is equal to the number of degree-of-freedom of the structure. It is shown the noise effects may partially be eliminated by using high-order finite impulse response (FIR) filters in both methods. Finally, the feasibility and accuracy of the presented model updating methods are examined through numerical studies on multistory shear buildings using the simulated records with various noise levels. The excellent agreement of the obtained results with those of the finite element models shows the feasibility of the proposed methods.     

Simulation of Thin-Film Damping and Thermal Mechanical Noise Spectra for Advanced Micromachined Microphone Structures

In many micromachined sensors the thin (2-10 thick) air film between a compliant diaphragm and backplate electrode plays a dominant role in shaping both the dynamic and thermal noise characteristics of the device. Silicon microphone structures used in grating-based optical-interference microphones have recently been introduced that employ backplates with minimal area to achieve low damping and low thermal noise levels. Finite-element based modeling procedures based on 2-D discretization of the governing Reynolds equation are ideally suited for studying thin-film dynamics in such structures which utilize relatively complex backplate geometries. In this paper, the dynamic properties of both the diaphragm and thin air film are studied using a modal projection procedure in a commonly used finite element software and the results are used to simulate the dynamic frequency response of the coupled structure to internally generated electrostatic actuation pressure. The model is also extended to simulate thermal mechanical noise spectra of these advanced sensing structures. In all cases simulations are compared with measured data and show excellent agreement - demonstrating 0.8 and 1.8 thermal force and thermal pressure noise levels, respectively, for the 1.5 mm diameter structures under study which have a fundamental diaphragm resonance-limited bandwidth near 20 kHz.     

矢量水听器的芯片级减震结构的设计

在MEMS矢量水听器现有结构的基础上,设计出一种新型弹簧减震结构,期望利用该结构提高水听器抗噪能力。根据Ansys有限元仿真,确定出该弹簧的几何尺寸、劲度系数(k≈5 N/m)及弹簧的所在位置——弹簧设置为内外圈相互垂直的两对。仿真结果表明,该减震结构抑制流噪声能力可提高3倍以上,同时几乎不影响水听器的接收灵敏度。该芯...     

Neurocomputing method based on structural finite element analysis of discrete model

Based on structural finite element analysis of discrete models, a neurocomputing strategy is developed in this paper. Dynamic iterative equations are constructed in terms of neural networks of discrete models. Determination of the iterative step size, which is important for convergence, is investigated based on the positive definiteness of the finite element stiffness matrix. Consequently, a method of choosing the step size of dynamic equations is proposed and the computational formula of the best step size is derived. The analysis of the computing model shows that the solution of finite element system equations can be obtained by the method of neural network computation efficiently. The proposed method can be used for parallel computation of structural finite element in a large-scale integrated circuit (LSI).     

Anatomical and dynamic volume spline model applied to facial soft tissue

Biomechanical modeling of soft tissue is a complex problem for achieving realistic surgical simulations, surgical planning, and scientific analysis. In the literature, three categories of biomechanical models: spline based models, spring models, and finite element models (FEMs) are mainly used for dealing with this problem. Among these, spline based models offer relatively fast and realistic soft tissue simulations by utilizing both the spring and FEMs. In this paper, a new dynamic volume spline model for human face skin is proposed and the performance of our model is discussed by estimating the results of facial surgery of three different patients. Face models of the patients are obtained from 3D CT/MR scans by segmenting the skull, muscle, and skin layers. In these face models, the skull and the muscle layers are considered as the rigid boundary for the skin layer and the skin layer is modeled by our dynamic volume spline. The control points of the dynamic volume spline are localized masses with viscoelastic material properties (stiffness, damping, and mass). These parameters are computed from the skin material properties that were published in the literature. Once the face models are generated, facial surgery plannings are simulated. Infact, the pre‐surgery face models are modified according to the surgical plans and the estimated post‐surgery face models are compared with the actual post‐surgery face models. Moreover, in order to discuss the performance of our dynamic volume spline model, the same analyses are performed on the post‐surgery estimations of a conventional tool. Copyright © 2011 John Wiley & Sons, Ltd.     

Finite element-based analysis of shunted piezoelectric structures for vibration damping

Piezoelectric patches shunted with passive electrical networks can be attached to a host structure for reduction of structural vibrations. This approach is frequently called “shunted piezo damping” and has the advantage of guaranteed stability and low complexity in implementation. For numerical treatment of such structures, a finite element modelling methodology is presented that incorporates both the piezoelectric coupling effects of the patches and the electrical dynamics of the connected passive electrical circuits. It allows direct computation of the achieved modal damping ratios as a major design criterion of interest. The damping ratios are determined from the eigenvalue problem corresponding to the coupled model containing piezoelectric structure and passive electrical circuit. The model includes local stiffening and mass effects as a result of the attached patches and, therefore, enables accurate prediction of the natural frequencies and corresponding modal damping ratios. This becomes crucial for choosing the patch thickness to achieve optimal modal damping for a given host structure. Additionally, structures with complex geometry or spatially varying material properties can easily be handled. Furthermore, the use of a finite element formulation for the coupled model of piezoelectric patches and a host structure facilitates design modifications and systematic investigations of parameter dependencies. In this paper, the impact of parameters of the passive electrical network on modal damping ratios as well as the variation of the patch thickness are studied. An application of this modelling method is realized by commercial software packages by importing fully coupled ANSYS^© – models in MATLAB^©. Afterwards, modal truncation is applied, the dynamic equations of the passive electrical network are integrated into the piezoelectric model and eigenvalue problems are solved to extract the increase in modal damping ratios. The numerical results are verified by experiments.     

Robust method for identifying material parameters based on virtual fields in elastodynamics

We develop an inverse method with the purpose of extracting elastic properties of materials in the framework of transient dynamics. To this end, we create a small linear system based on a set of well-chosen time-dependent virtual fields (VF) and measurement data. The parameters are the solutions of this system and can be quickly extracted. We compare this new method with the classical finite element model updating (FEMU) method for different case studies. In our study, the measurements are synthetic, i.e, they are calculated using a fine finite element (FE) model. Uniform white noise is added to model measurement uncertainties. Results, based on Monte Carlo simulations, show that our method is more robust and accurate than the FEMU method for an acceptable noise level. Our new method appears well-adapted to linear elasticity in transient dynamics.     

Hexahedral mesh generation for biomedical models in SCIRun

Biomedical simulations are often dependent on numerical approximation methods, including finite element, finite difference, and finite volume methods, to model the varied phenomena of interest. An important requirement of the numerical approximation methods above is the need to create a discrete decomposition of the model geometry into a ‘mesh’. Historically, the generation of these meshes has been a critical bottleneck in efforts to efficiently generate biomedical simulations which can be utilized in understanding, planning, and diagnosing biomedical conditions. In this paper we discuss a methodology for generating hexahedral meshes for biomedical models using an algorithm implemented in the SCIRun Problem Solving Environment. The method is flexible and can be utilized to build up conformal hexahedral meshes ranging from models defined by single isosurfaces to more complex geometries with multi-surface boundaries.