NUMERICAL SIMULATION OF FIBROBLAST DUROTAXIS MIGRATION INDUCED BY STIFFNESS GRADIENT OF SUBSTRATE1)

建立了一种细胞趋硬性迁移的理论模型和有限元分析框架,为连续变刚度人工基质的试验设计提供理论依据。考虑了细胞体的黏弹性属性,以及细胞与基质间的配受体动态反应过程,并以配受体合成时间为时间步长,将细胞运动方程化为静力学形式进行求解。对有限元过程提出一种动约束,便于消除其结构矩阵的奇异性。结果表明,模型能够模拟黏着斑内部力的快速波动现象,细胞的运动速度与观测数据一致,可有效模拟20,h以上的长时程问题。     

Residual‐based stabilization of the finite element approximation to the acoustic perturbation equations for low Mach number aeroacoustics

The acoustic perturbation equations (APE) are suitable to predict aerodynamic noise in the presence of a non‐uniform mean flow. As for any hybrid computational aeroacoustics approach, a first computational fluid dynamics simulation is carried out from which the mean flow characteristics and acoustic sources are obtained. In a second step, the APE are solved to get the acoustic pressure and particle velocity fields. However, resorting to the finite element method (FEM) for that purpose is not straightforward. Whereas mixed finite elements satisfying an appropriate inf–sup compatibility condition can be built in the case of no mean flow, that is, for the standard wave equation in mixed form, these are difficult to implement and their good performance is yet to be checked for more complex wave operators. As a consequence, strong simplifying assumptions are usually considered when solving the APE with FEM. It is possible to avoid them by resorting to stabilized formulations. In this work, a residual‐based stabilized FEM is presented for the APE at low Mach numbers, which allows one to deal with the APE convective and reaction terms in its full extent. The key of the approach resides in the design of the matrix of stabilization parameters. The performance of the formulation and the contributions of the different terms in the equations are tested for an acoustic pulse propagating in sheared‐solenoidal mean flow, and for the aeolian tone generated by flow past a two‐dimensional cylinder. Copyright © 2016 John Wiley & Sons, Ltd.     

Selectivity of seismic electric signal (SES) of the 2000 Izu earthquake swarm: a 3D FEM numerical simulation model

Although seismic electric signal (SES) has been used for short-term prediction of earthquakes, selectivity of SES still remains as one of the mysterious features. As a case study, we made a numerical simulation based on a 3D finite element method (FEM) on the selectivity of SES observed in the case of the 2000 Izu earthquake swarm. Our numerical results indicated that the existence of conductive channel under Niijima island could explain the reported SES selectivity.     

A refined numerical model for determining inflation-burst behavior of composite membrane structures

The airship structures made of multi-layer composite fabrics or membranes can offer the platform for earth observations, wireless communications and space research due to light weight and good mechanical performance. The structural safety and serviceability strongly depend on material properties and working conditions. Available studies are limited within service stress limits or are lack of suitable biaxial tensile constitutive models for understanding structural behavior. This paper thus focuses on a refined numerical model for determining inflation-burst behavior of composite airship structures considering new biaxial constitutive equations, novel failure criteria and manufacture factors.The differences between ideal and real forms of airship structures, e.g. volume difference, demonstrate the necessity for incorporating cutting-pattern effects in the initial numerical model. For structural analysis, stress distributions on real structural forms are different from those on ideal forms because of welding parts that can enhance local stiffness. The ultimate pressures are 56.7 kPa and 59.5 kPa for ideal and real structural forms. Structural breaking initiated at the maximum diameter of ideal structural forms propagates fast while welding parts can prevent breaking propagation for real structural forms. Therefore, the refined numerical model can reveal basic structural behavior and safety performance of airship structures in the inflation-burst processes.     

Snap-through path in a bistable dielectric elastomer actuator

The dielectric elastomer (DE) has attracted significant attention due to its desired features, including large deformation, fast response, and high energy density. However, for a DE actuator (DEA) utilizing a snap-through deformation mode, most existing theoretical models fail to predict its deformation path. This paper develops a new finite element method (FEM) based on the three-parameter Gent-Gent model suitable for capturing strain-stiffening behaviors. The simulation results are verified by experiments, indicating that the FEM can accurately characterize the snap-through path of a DE. The method proposed in this paper provides theoretical guidance and inspiration for designing and applying DEs and bistable electroactive actuators.     

物理参数变化对短脉冲激光激励温度场的影响

 为研究多物理参数（耦合系数、电子热导率、电子热容、晶格热容）同时随温度变化对短脉冲激光辐照金属材料产生温度场分布的影响,基于双温耦合理论,建立了短脉冲激光辐照金属材料金的加热过程的有限元求解模型。在同时考虑脉冲激光的空间、时间分布和多参数同时随温度变化的情况下,得到短脉冲激光辐照金属材料金激励产生的温度场二维瞬态分布,并进一步比较了多物理参数同时随温度变化和采用室温物理参数两种情况下温度场分布的区别。数值结果表明：多物理参数同时随温度变化使电子温度和晶格温度的上升变快,最大值变大,而且使得材料中激光穿透直接辐照到的区域温度变高。     

An alternative alpha finite element method with discrete shear gap technique for analysis of laminated composite plates

This paper presents an alternative alpha finite element method using triangular meshes (AαFEM) for static, free vibration and buckling analyses of laminated composite plates. In the AαFEM, an assumed strain field is carefully constructed by combining compatible strains and additional strains with an adjustable parameter α which can produce an effectively softer stiffness formulation compared to the linear triangular element. The stiffness matrices are obtained based on the strain smoothing technique over the smoothing domains and the constant strains on triangular sub-domains associated with the nodes of the elements. The discrete shear gap (DSG) method is incorporated into the AαFEM to eliminate transverse shear locking and an improved triangular element termed as AαDSG3 is proposed. Several numerical examples are then given to demonstrate the effectiveness of the AαDSG3.     

A posteriori error estimates for optimal distributed control governed by the evolution equations

We describe a technique for a posteriori error estimates suitable to the optimal control problem governed by the evolution equations solved by the method of lines. It is applied to the control problem governed by the parabolic equation, convection-diffusion equation and hyperbolic equation. The error is measured with the aid of the L²-norm in the space-time cylinder combined with a special time weighted energy norm.