Acoustoelastic effect simulation by time–space finite element formulation based on quadratic interpolation of the acceleration

Acoustoelastic effect describes the change of ultrasound velocity due to the initial stress. Its simulation involves a numerical analysis of nonlinear elastodynamics and requires high accuracy in the time domain. A time–space finite element formulation, derived from the quadratic interpolation of the acceleration within a time segment, is proposed for an accurate simulation of the acoustoelastic effect in the present study. Ten different integration schemes are generated based on this formulation and nine of them are found to be conditionally stable. Among the nine stable schemes, one is found to obtain a spectral radius of one when the normalized step ratio is less than 5.477, indicating no numerical dissipation or numerical divergence. Compared with integration schemes from previous studies, this integration scheme demonstrates better performance in calculation accuracy and energy conservation. A two-stage approach, namely the static stage and the dynamic stage, has been employed in the simulation of the acoustoelastic effect. The former stage is adopted to obtain the initial stress and the latter stage, where the proposed integration scheme is implemented, is adopted to simulate the ultrasound propagation in an initial stress state. The simulation results of the dynamic stage show that the ultrasound velocity increases in a compression stress state and decreases in a tension stress state for aluminum alloy, which is in good agreement with previous experimental studies. Together with the simulation result of the static stage, it is conjectured that the acoustoelastic effect results from the stress-dependent elastic modulus.     

Aeroelastic analysis of cantilever plates using Peters’ aerodynamic model,and the influence of choosing beam or plate theories as the structural model

In this paper, the aeroelastic analyses of a rectangular cantilever plate of varying aspect ratio is presented. The classical plate theory has been selected as the structural model. The main point that distinguishes this study from previously reported research is employing Peters’ theory to model aerodynamic effect which is not straightforward. The Peters’ aerodynamic model was originally developed to provide lift and moment, which is only applicable to the structural model based on the beam theories. In this study, using the basic concept of the Peters’ aerodynamic model in addition to utilizing the Fourier series, the pressure distribution is derived, which makes Peters’ model applicable to structural models based on plate theory. This combination provides a much simpler state–space aeroelastic model for plates in comparison to the prevalent panel methods, which could lead to a significant reduction in computational time. In addition, the aeroelastic response of the plate with respect to changes in the structural model from the beam theory to the plate theory is evaluated. By using data from an experiment carried out at Duke University, the theoretical results are evaluated. Furthermore, the differences in structural models obtained from the plate and beam theories can be divided into two distinct parts, which are responsible for differences in bending and torsional behaviors of the structure, separately. This approach enables us to measure the effects of differences of each behavior separately, which could provide with a new insight into the problem. It has been determined that the flutter speeds obtained from the beam and plate aeroelastic models are little affected by the difference in bending behavior, but rather is mainly caused by the difference in torsional frequencies.     

SYSTEM OF VECTOR QUASI-EQUILIBRIUM PROBLEMS AND ITS APPLICATIONS

A new system of vector quasi-equilibrium problems is introduced and its existence of solution is proved. As applications, some existence results of weak Pareto equilibrium for both constrained multicriteria games and multicriteria games without constrained correspondences are also shown.     

Dynamic Analysis of an Infinite Cylindrical Hole in a Saturated Poroelastic Medium

In this paper, the dynamic response of an infinite cylindrical hole embedded in a porous medium and subjected to an axisymmetric ring load is investigated. Two scalar potentials and two vector potentials are introduced to decouple the governing equations of Biot’s theory. By taking a Fourier transform with respect to time and the axial coordinate, we derive general solutions for the potentials, displacements, stresses and pore pressures in the frequency-wave-number domain. Using the general solutions and a set of boundary conditions applied at the hole surface, the frequency-wave-number domain solutions for the proposed problem are determined. Numerical inversion of the Fourier transform with respect to the axial wave number yields the frequency domain solutions, while a double inverse Fourier transform with respect to frequency as well as the axial wave number generates the time-space domain solution. The numerical results of this paper indicate that the dynamic response of a porous medium surrounding an infinite hole is dependant upon many factors including the parameters of the porous media, the location of receivers, the boundary conditions along the hole surface as well as the load characteristics.     

Series solutions of unsteady MHD flows above a rotating disk

The series solutions of unsteady flows of a viscous incompressible electrically conducting fluid caused by an impulsively rotating infinite disk are given by means of an analytic technique, namely the homotopy analysis method. Using a set of new similarity transformations, we transfer the Navier–Stokes equations into a pair of nonlinear partial differential equations. The convergent series solutions are obtained, which are uniformly valid for all dimensionless time 0 ≤ τ < ∞ in the whole spatial region 0 ≤ η < ∞. To the best of our knowledge, such kind of series solutions have never been reported. The effect of magnetic number on the velocity is investigated.     

Parameter Determination for Reductive Dechlorination of Chlorinated Solvents

The contamination of subsurface due to the chlorinated solvents such as tetrachlorethylene (PCE) and trichlorethylene (TCE) is one of the most difficult environmental problems to treat. Bioremediation has been shown by many researchers to be a remedial alternative for this type of contamination. Chlorinated solvents are not directly mineralized but rather are transformed by microorganisms into one or more intermediate compounds before converting into a final compound. These sequential reactions, termed “reductive dehalogenation”, consist of replacing a chlorine atom by a hydrogen atom. The pathway of degradation of PCE can be expressed by the following scheme PCE → TCE → DCE → VC → ETH, where dichloroethene (DCE), vinyl chloride (VC) and finally ETH is ethylene. Since the biotransformation rate coefficients of each intermediate compound are different, they have to be determined very precisely to establish an effective treatment operation. The sequential decay can be described by Michaelis–Menten’s kinetics, which constitutes a highly nonlinear system of ordinary differential equations (ODEs). This is very sensitive to the changes of biotransformation rate coefficients. In this study we introduce a methodology how to numerically estimate the rate coefficients for Michaelis–Menten’s equations from the knowledge of the concentrations of PCE, TCE, DCE, VC and ETH. The efficiency of the proposed method is demonstrated on some examples. Estimated biotransformation coefficients are employed to predict the concentrations of chlorinated solvents. Computations and measurements show a very good agreement.     

基于微面有效应力矢量的各向异性屈服准则

基于微面模型,定义损伤变量为微面上有效承载面积的减少. 将Kachanov的一维有效 应力概念推广到三维,提出微面有效应力矢量的概念. 根据微面的有效应力矢量,将无损材 料的宏观应力张量及不变量与微面应力矢量的积分关系拓展到有损材料,得到了有损材料的 宏观有效应力张量及其不变量与宏观名义应力张量、微面面积损伤组构张量之间的关系. 将 无损材料的以应力张量不变量表示的Drucker-Prager准则推广到有损材料,建立了含缺陷 材料的各向异性屈服准则. 对有损材料,宏观有效应力张量与Murakami的有效应力张量具 有相同的形式,各向异性强度准则与Liu等提出的扩展Hill准则有相同的形式,当不考虑 静水应力对屈服的影响时,它与Hill准则具有相同的形式.     

基于非线性梁理论的有限质点法

有限质点法是以向量式力学为基础,用有限数量的质点来模拟结构的变形行为,质点的运动由牛顿运动定律来计算。在有限质点法中,质点通过构件相连,构件约束着质点的运动,并且其内力由质点的运动变量来描述。基于向量式力学的基本思想和非线性梁理论,提出了一种新的有限质点法,该方法在共旋单元坐标系中描述梁的非线性变形。以空间梁系结构为例,推导了计算构件内力的非线性公式,并考虑了弯扭耦合变形。通过两个连续欧拉角的变换公式得到共旋坐标系的旋转矩阵。与传统的有限质点法相比,本文提出的方法避免了刚体虚转动分析。通过四个结构的数值求解,验证了本文方法在计算结构大变形响应时具有较高的精度。     

基于Hopf-Cole变换的Burgers方程初边值问题的Sinc-Galerkin法

利用Sinc-Galerkin法数值求解Burgers方程的初边值问题。首先,用Hopf-Cole变换将二阶非线性的Burgers方程变换为二阶线性方程,同时把第一类边界条件变为第二类边界条件。时间上的导数采用θ加权格式离散,空间导数采用Sinc-Galerkin法离散,端点处分别引入权函数处理变换后的第二类边界条件。最后,通过数值算例验证了Sinc-Galerkin法的指数收敛性,与精确解相比,本文构造的数值格式精度高,能够有效捕捉激波等物理现象。