高速船舱壁加筋板流固耦合振动分析

将结构有限元法和流体有限元法结合起来,根据Garlerkin法和Hamilton变分原理分别推导出离散的流体和结构的运动方程,用湿模态法求得流体附加质量矩阵,对高速船舱壁加筋板流固耦合系统进行了模态分析。并进行模态试验,研究不同水深时舱壁加筋板结构流固耦合振动模态,同时与有限元计算结果进行比较,从而揭示了该类结构流固耦合振动的动力特性,对舱壁加筋板结构的动力设计具有指导意义。     

静水中反应堆吊篮结构的振动特性研究

将吊篮与压力容器视为包含粘性流体的有限长同心圆柱壳体结构。从圆柱壳的运动控制微分方程和流体的 N- S方程出发 ,在流体不可压缩粘性假设下 ,导出吊篮结构自由振动时 ,流体作用于壳体结构上的附加质量、附加阻尼矩阵。再结合固体有限元技术 ,建立吊篮结构的流固耦合振动的有限元模型。用此模型对吊篮结构在静水中的固有振动特性进行分析 ,计算结果与试验实测结果的比较表明 ,本文方法能较好地满足工程分析要求     

Explicit dynamics in impact simulation using a NURBS contact interface

In this paper, the impact problem and the subsequent wave propagation are considered. For the contact discretization an intermediate non-uniform rational B-spline (NURBS) layer is added between the contacting finite element bodies, which allows a smooth contact formulation and efficient element-based integration. The impact event is ill-posed and requires a regularization to avoid propagating stress oscillations. A nonlinear mesh-dependent penalty regularization is used, where the stiffness of the penalty regularization increases upon mesh refinement. Explicit time integration methods are well suited for wave propagation problems, but are efficient only for diagonal mass matrices. Using a spectral element discretization in combination with a NURBS contact layer the bulk part of the mass matrix is diagonal.     

Reduction of the degrees of freedom in solving dynamic problems by the finite element method

The method of lowering the order of a matrix equation which describes the dynamics of a finite element model is presented. The finite element with a ‘truncated’ mass matrix is obtained for calculating thin plate vibrations. Such an element has one vibrational degree of freedom at each nodal point. In the case of uniform systems the accuracy provided by the suggested element is no less than that provided by the non-conforming elements, which have three vibrational degrees of freedom at each nodal point, and in some cases it is greater. The finite elements with a ‘truncated’ mass matrix have essential advantages in the study of non-uniform systems.     

Direct and sparse construction of consistent inverse mass matrices: general variational formulation and application to selective mass scaling

Classical explicit finite element formulations rely on lumped mass matrices. A diagonalized mass matrix enables a trivial computation of the acceleration vector from the force vector. Recently, non‐diagonal mass matrices for explicit finite element analysis (FEA) have received attention due to the selective mass scaling (SMS) technique. SMS allows larger time step sizes without substantial loss of accuracy. However, an expensive solution for accelerations is required at each time step. In the present study, this problem is solved by directly constructing the inverse mass matrix. First, a consistent and sparse inverse mass matrix is built from the modified Hamiltons principle with independent displacement and momentum variables. Usage of biorthogonal bases for momentum allows elimination of momentum unknowns without matrix inversions and directly yields the inverse mass matrix denoted here as reciprocal mass matrix (RMM). Secondly, a variational mass scaling technique is applied to the RMM. It is based on the penalized Hamiltons principle with an additional velocity variable and a free parameter. Using element‐wise bases for velocity and a local elimination yields variationally scaled RMM. Thirdly, examples illustrating the efficiency of the proposed method for simplex elements are presented and discussed. Copyright © 2014 John Wiley & Sons, Ltd.     

发展附加质量模型应用于储液罐的动力分析

提出一种新的用于表述附加质量分布规律的公式和有限元单元模型,接入ABAQUS有限元软件对储液罐动力试验模型进行了数值模拟,分析对比了弹塑性屈曲“象足效应”和“钻石效应”的实验数据,证明所提出的附加质量公式是合理和有效的,并指出原经验公式在分析动水压力问题时的缺陷。模拟中考虑了储液罐壁动力响应加速度的方向对液体惯性作用力的影响,从而准确地模拟了储液罐壁与液体的相互作用。动力时程和频率响应的定量分析表明:地震荷载对储液罐计算模型所做的外力功是导致罐壁塑性变形、产生“象足效应”的决定性因素。     

ChromAlign: A two-step algorithmic procedure for time alignment of three-dimensional LC-MS chromatographic surfaces

We present an algorithmic approach to align three-dimensional chromatographic surfaces of LC-MS data of complex mixture samples. The approach consists of two steps. In the first step, we prealign chromatographic profiles: two-dimensional projections of chromatographic surfaces. This is accomplished by correlation analysis using fast Fourier transforms. In this step, a temporal offset that maximizes the overlap and dot product between two chromatographic profiles is determined. In the second step, the algorithm generates correlation matrix elements between full mass scans of the reference and sample chromatographic surfaces. The temporal offset from the first step indicates a range of the mass scans that are possibly correlated, then the correlation matrix is calculated only for these mass scans. The correlation matrix carries information on highly correlated scans, but it does not itself determine the scan or time alignment. Alignment is determined as a path in the correlation matrix that maximizes the sum of the correlation matrix elements. The computational complexity of the optimal path generation problem is reduced by the use of dynamic programming. The program produces time-aligned surfaces. The use of the temporal offset from the first step in the second step reduces the computation time for generating the correlation matrix and speeds up the process. The algorithm has been implemented in a program, ChromAlign, developed in C++ language for the .NET2 environment in WINDOWS XP. In this work, we demonstrate the applications of ChromAlign to alignment of LC-MS surfaces of several datasets: a mixture of known proteins, samples from digests of surface proteins of T-cells, and samples prepared from digests of cerebrospinal fluid. ChromAlign accurately aligns the LC-MS surfaces we studied. In these examples, we discuss various aspects of the alignment by ChromAlign, such as constant time axis shifts and warping of chromatographic surfaces.     

Sustained-release pellets prepared by combination of wax matrices and double-layer coatings for extremely water-soluble drugs

This study was performed in order to develop a sustained-release pellet formulation containing venlafaxine hydrochloride (VEN), an extremely water-soluble drug, prepared by combination of wax matrices and double-layer coatings. The influence of both double-layer polymeric coats and wax matrices on the release of VEN from sustained-release pellets was investigated. The pellets were prepared by wet mass extrusion spheronization methods and then coated with a fluidized bed coater. For the pellets coated with Eudragit NE30D alone, a coating level of nearly 40% was required to pass the dissolution test compared with commercial product, and it was accompanied by an unacceptable lag time. The application of an alcohol-soluble polymeric subcoat, Opadry I, was added before the Eudragit NE30D coating process, which resulted in a marked delay in drug release. However, a faster release was observed for the formulation coated with a high subcoat level (10%) at the end of the dissolution test. A further delay in drug release was observed when a wax matrix, octadecanol, was added to the core pellet formulation. The kinetics of drug release changed from the Higuchi model to a zero order model and the predominant mechanism controlling drug release changed from diffusion to dissolution upon increasing the amount of octadecanol within the matrix pellets. In addition, the drug release was markedly influenced by the drug to matrix ratio. In conclusion, the 40% drug-loaded core pellets with double-layer coatings (8% Opadry I and 12% Eudragit NE30D) and 20% octadecanol matrix produced the desired profile for once-daily sustained release compared with the commercial product, and these pellets remained stable during storage.     

A mixture model for wet bones—I theory

A mixture model for wet bones is given in this paper. In developing the theory, it is assumed that a wet bone is mainly composed of electrolytes, bone salts and the solid matrix bone which has the piezoelectric property. Based on the fact that the electrical conductance and calcification in the wet bones are, in turn, due to monovalent ions (Na, Cl, HCO₃) and mass transfer evolving between the piezoelectric matrix bone and the bone salts (Ca₃, (PO₄)₂), the governing balance equations and the constitutive relations (both nonlinear and linear) of wet bones are derived and various special cases are discussed.     

Refined triangular discrete Kirchhoff plate element for thin plate bending,vibration and buckling analysis

A refined triangular discrete Kirchhoff thin plate bending element RDKT which can be used to improve the original triangular discrete Kirchhoff thin plate bending element DKT is presented. In order to improve the accuracy of the analysis a simple explicit expression of a refined constant strain matrix with an adjustable constant can be introduced into its formulation. The new element displacement function can be used to formulate a mass matrix called combined mass matrix for calculation of the natural frequency and in the same way a combined geometric stiffness matrix can be obtained for buckling analysis. Numerical examples are presented to show that the present methods indeed, can improve the accuracy of thin plate bending, vibration and buckling analysis. © 1998 John Wiley & Sons, Ltd.     

Added mass computation by the boundary integral method

Computational techniques for the treatment of fluid-structure interaction effects by discrete boundary integral methods are examinde. Attention is focused on the computation of the added mass matrix by finite element methods for a structure submerged in an infinite, inviscid, incompressible fluid. A general computational procedure is presented that is based upon a variational approach involving the assumption of constant source strength over each surface element. This is followed by an analysis of the discretization error for a spherical body that is then used to develop a hierarchy of computational schemes. These schemes are than evaluated numerically in terms of ‘fluid boundary modes’ for a submerged spherical surface. One scheme has been found to be surprisingly accurate in relation to its computational demands.     

A new boundary element formulation for wave load analysis

A new boundary element (BEM) formulation is proposed for wave load analysis of submerged or floating bodies. The presented formulation, through establishing an impedance relation, permits the evaluation of the hydrodynamic coefficients (added mass and damping coefficients) and the coefficients of wave exciting forces systematically in terms of system matrices of BEM without solving any special problem, such as, unit velocity or unit excitation problem. It also eliminates the need for scattering analysis in the evaluation of wave exciting forces. The imaginary and real parts of impedance matrix give, respectively, added mass and damping matrices whose elements describe the fluid resistance against the motion of the body. The formulation is explained through the use of a simple fluid-solid system under wave excitations, which involves a uniform fluid layer containing a solid cylindrical body. In the formulation, the solid body is taken first as deformable, then, it is specialized when it is rigid. The validity of the proposed method is verified by comparing its result with those available in literature for rigid submerged or floating bodies.     

The iterative solution of Taylor—Galerkin augmented mass matrix equations

This paper investigates the convergence properties of iterative schemes for the solution of finite element mass matrix equations that arise through the application of a Taylor-Galerkin algorithm to solve instationary Navier-Stokes equations. This is a time-stepping algorithm that involves Galerkin mass matrix equations at fractional stages within each time-step. Plane Poiseuille flow and shear-driven cavity flow are selected as benchmark problems on which to investigate the effects of various choice of scheme and time-step dependency. The iterative convergence of each mass matrix equation for a single fractional stage is studied, both at the element and the system matrix level. The underlying theory is confirmed and it is shown how optimal iterative convergence rates may be achieved for a Jacobi scheme by employing an appropriate acceleration factor. Moreover, this factor is trivial to compute. The consequential effects on the convergence of the time-stepping procedure to reach steady-state are also considered where non-linear effects are present.     

Sustained-Release Pellets Prepared by Combination of Wax Matrices and Double-Layer Coatings for Extremely Water-Soluble Drugs

This study was performed in order to develop a sustained-release pellet formulation containing venlafaxine hydrochloride (VEN), an extremely water-soluble drug, prepared by combination of wax matrices and double-layer coatings. The influence of both double-layer polymeric coats and wax matrices on the release of VEN from sustained-release pellets was investigated. The pellets were prepared by wet mass extrusion spheronization methods and then coated with a fluidized bed coater. For the pellets coated with Eudragit® NE30D alone, a coating level of nearly 40% was required to pass the dissolution test compared with commercial product, and it was accompanied by an unacceptable lag time. The application of an alcohol-soluble polymeric subcoat, Opadry® I, was added before the Eudragit® NE30D coating process, which resulted in a marked delay in drug release. However, a faster release was observed for the formulation coated with a high subcoat level (10%) at the end of the dissolution test. A further delay in drug release was observed when a wax matrix, octadecanol, was added to the core pellet formulation. The kinetics of drug release changed from the Higuchi model to a zero order model and the predominant mechanism controlling drug release changed from diffusion to dissolution upon increasing the amount of octadecanol within the matrix pellets. In addition, the drug release was markedly influenced by the drug to matrix ratio. In conclusion, the 40% drug-loaded core pellets with double-layer coatings (8% Opadry® I and 12% Eudragit® NE30D) and 20% octadecanol matrix produced the desired profile for once-daily sustained release compared with the commercial product, and these pellets remained stable during storage.     

ITERATIVE METHOD USING CONSISTENT MASS MATRIX IN AXISYMMETRICAL FINITE ELEMENT ANALYSIS OF HYPERVELOCITY IMPACT

We present in this paper an iterative method using consistent mass matrix in axisymmetrical finite element analysis of hypervelocity impact. To retain the advantage of integration on an element-by-element basis which is at the heart of modern hydrocodes, we suggest that the first step should be to solve for accelerations at an advanced time step by using the lumped mass approach, then iterate using a consistent mass matrix to improve the estimate. Examples are given to show the improved resolution with the new method.     

Vibration analysis of fluid–solid systems using a finite element displacement formulation

This report presents a finite element solution for the vibration interaction between an inviscid fluid and a solid. The equation of motion governing the inviscid fluid is expressed in terms of the displacements. This ensures that compatibility and equilibrium will be satisfied automatically along the interface of the coupled systems. To suppress circulation modes with non-zero energy, reduced integration is used when computing the element stiffness matrix contributed by the fluid. In addition, a projection is used on the element mass matrix in order to remove the spurious modes which result from the use of reduced integration. Numerical examples for both fluid and coupled fluid–solid systems are performed and the results are shown.     

Non‐linear exact geometry 12‐node solid‐shell element with three translational degrees of freedom per node

This paper presents the finite rotation exact geometry (EG) 12‐node solid‐shell element with 36 displacement degrees of freedom. The term ‘EG’ reflects the fact that coefficients of the first and second fundamental forms of the reference surface and Christoffel symbols are taken exactly at each element node. The finite element formulation developed is based on the 9‐parameter shell model by employing a new concept of sampling surfaces (S‐surfaces) inside the shell body. We introduce three S‐surfaces, namely, bottom, middle and top, and choose nine displacements of these surfaces as fundamental shell unknowns. Such choice allows one to represent the finite rotation higher order EG solid‐shell element formulation in a very compact form and to derive the strain–displacement relationships, which are objective, that is, invariant under arbitrarily large rigid‐body shell motions in convected curvilinear coordinates. The tangent stiffness matrix is evaluated by using 3D analytical integration and the explicit presentation of this matrix is given. The latter is unusual for the non‐linear EG shell element formulation. Copyright © 2011 John Wiley & Sons, Ltd.     

Inverse mass matrix via the method of localized lagrange multipliers

An efficient method for generating the mass matrix inverse of structural dynamic problems is presented, which can be tailored to improve the accuracy of target frequency ranges and/or wave contents. The present method bypasses the use of biorthogonal construction of a kernel inverse mass matrix that requires special procedures for boundary conditions and free edges or surfaces and constructs the free‐free inverse mass matrix using the standard FEM procedure. The various boundary conditions are realized by the the method of localized Lagrange multipliers. In particular, the present paper constructs the kernel inverse matrix by using the standard FEM elemental mass matrices. It is shown that the accuracy of the present inverse mass matrix is almost identical to that of a conventional consistent mass matrix or a combination of lumped and consistent mass matrices. Numerical experiments with the proposed inverse mass matrix are conducted to validate its effectiveness when applied to vibration analysis of bars, beams, and plain stress problems.     

Added mass computations by integral equation methods

This paper presents boundary integral equation procedures for calculating the added mass matrix required to determine hydroelastic vibrational modes of tanks or immersed structures. The liquid is assumed inviscid and incompressible and no surface waves are admitted. Two symmetric variational formulations of the boundary integral equation are derived, including the free surface condition and possible conditions of symmetry. These are approximated by a curved finite element method and numerical results are presented for an axisymmetric tank, immersed cantilever plates and a three-dimensional fuel storage container.     

Two-step matrix application technique to improve ionization efficiency for matrix-assisted laser desorption/ionization in imaging mass spectrometry

A novel matrix application protocol for direct tissue mass spectrometry is presented. Matrix-assisted laser desorption/ionization is a popular ionization procedure for direct tissue analysis and imaging mass spectrometry. Usually, matrixes are applied by dispensing droplets through either pipettes or automated dispensing machines, or by airbrushing. These techniques are very simple, but it was difficult to obtain uniform matrix crystals on the tissue surface, and nonuniform crystals degrade the spectrum qualities. Here we report a new matrix application protocol, which is a combination of spraying and dispensing droplets, and we have succeeded in overcoming these problems in conventional matrix applications on tissue surfaces. We call our new technique the "spray-droplet method". In this technique, tiny matrix crystals formed by spraying act as seeds for crystal growth. Our technique leads to matrix spots that are filled homogeneously with minute crystals. Such matrix crystals dramatically improve peak intensity and signal-to-noise ratio. In an example on a rat brain section, the number of detectable peaks was increased and signal intensity of m/z 5440 in our method was approximately 30.6 times higher than that in conventional methods. We used this spray-droplet method with a chemical ink-jet technology for matrix deposition to succeed in MALDI imaging of signals, which were undetectable from the conventional matrix applications.