A computational framework for fluid–porous structure interaction with large structural deformation

We study the effect of poroelasticity on fluid–structure interaction. More precisely, we analyze the role of fluid flow through a deformable porous matrix in the energy dissipation behavior of a poroelastic structure. For this purpose, we develop and use a nonlinear poroelastic computational model and apply it to the fluid–structure interaction simulations. We discretize the problem by means of the finite element method for the spatial approximation and using finite differences in time. The numerical discretization leads to a system of non-linear equations that are solved by Newton’s method. We adopt a moving mesh algorithm, based on the Arbitrary Lagrangian–Eulerian method to handle large deformations of the structure. To reduce the computational cost, the coupled problem of free fluid, porous media flow and solid mechanics is split among its components and solved using a partitioned approach. Numerical results show that the flow through the porous matrix is responsible for generating a hysteresis loop in the stress versus displacement diagrams of the poroelastic structure. The sensitivity of this effect with respect to the parameters of the problem is also analyzed.

     

Fluid pressure response in poroelastic materials subjected to cyclic loading

When cyclic loading is applied to poroelastic materials, a transient stage of interstitial fluid pressure occurs, preceding a steady state. In each stage, the fluid pressure exhibits a characteristic mechanical behavior. In this study, an analytical solution for fluid pressure in two-dimensional poroelastic materials, which is assumed to be isotropic, under cyclic axial and bending loading is presented, based on poroelasticity. The obtained analytical solution contains transient and steady-state responses. Both of these depend on three dimensionless parameters: the dimensionless stress coefficient; the dimensionless frequency; and, the axial-bending loading ratio. We focus particularly on the transient behavior of interstitial fluid pressure with changes in the dimensionless frequency and the axial-bending loading ratio. The transient properties, such as half-value period and contribution factor, depend largely on the dimensionless frequency and have peak values when its value is about 10. This suggests that, under these conditions, the transient response can significantly affect the mechanical behavior of poroelastic materials.     

An analytical poroelastic model for laboratorial mechanical testing of the articular cartilage (AC)

The articular cartilage (AC) can be seen as a biphasic poroelastic material. The cartilage deformation under compression mainly leads to an interstitial fluid flow in the porous solid phase. In this paper, an analytical poroelastic model for the AC under laboratorial mechanical testing is developed. The solutions of interstitial fluid pressure and velocity are obtained. The results show the following facts. (i) Both the pressure and fluid velocity amplitudes are proportional to the strain loading amplitude. (ii) Both the amplitudes of pore fluid pressure and velocity in the AC depend more on the loading amplitude than on the frequency. Thus, in order to obtain the considerable fluid stimulus for the AC cell responses, the most effective way is to increase the loading amplitude rather than the frequency. (iii) Both the interstitial fluid pressure and velocity are strongly affected by permeability variations. This model can be used in experimental tests of the parameters of AC or other poroelastic materials, and in research of mechanotransduction and injury mechanism involved interstitial fluid flow.     

Transient response of fluid pressure in a poroelastic material under uniaxial cyclic loading

Poroelasticity is a theory that quantifies the time-dependent mechanical behavior of a fluid-saturated porous medium induced by the interaction between matrix deformation and interstitial fluid flow. Based on this theory, we present an analytical solution of interstitial fluid pressure in poroelastic materials under uniaxial cyclic loading. The solution contains transient and steady-state responses. Both responses depend on two dimensionless parameters: the dimensionless frequency Ω that stands for the ratio of the characteristic time of the fluid pressure relaxation to that of applied forces, and the dimensionless stress coefficient H governing the solid-fluid coupling behavior in poroelastic materials. When the phase shift between the applied cyclic loading and the corresponding fluid pressure evolution in steady-state is pronounced, the transient response is comparable in magnitude to the steady-state one and an increase in the rate of change of fluid pressure is observed immediately after loading. The transient response of fluid pressure may have a significant effect on the mechanical behavior of poroelastic materials in various fields.     

Behavior of cell aggregates under force-controlled compression

In this paper we study the mechanical behavior of multicellular aggregates under a cycle of compressive loads and releases. Some analytical properties of the solution are discussed and numerical results are presented for a compressive test under constant force imposed on a cylindrical specimen and for a cycle of compressions and releases. We show that a steady loaded configuration is achieved. The analytical determination of the steady state value allows to obtain mechanical parameters of the cellular structure that cannot be obtained on the basis of creep tests at constant stress.     

不可压饱和多孔弹性梁、杆动力响应的数学模型

基于多孔介质理论,首先建立了饱和多孔弹性杆件弯曲与轴向变形时动力响应的数学模型.其次,基于多孔弹性梁弯曲变形的数学模型,利用Laplace变换,分析了两端可渗透的饱和多孔弹性悬臂梁在自由端受阶梯载荷作用下的动静力响应,给出了梁弯曲时挠度、弯矩以及孔隙流体压力等效力偶等物理量随时间的响应曲线.发现不可压多孔弹性梁的拟静态响应亦存在Mandel-Cryer现象,多孔弹性梁的挠度具有与粘弹性梁挠度类似的蠕变特征,然而,其应力响应不同于粘弹性梁,随着时间的增加,梁拟静态响应的弯矩逐渐增加,并达到一个稳态值.这些结果有助于揭示植物根茎等力学行为的机理.     

Identification of the elastic stiffness of composites using the virtual fields method and digital image correlation

This paper presents an effective methodology for characterizing the mechanical parameters of composites using digital image correlation combined with the virtual fields method.By using a three-point bending test configuration,this method can identify all mechanical parameters of the material with merely a single test.Successful results verified that this method is especially effective for characterizing composite materials.In this study,the method is applied to measure the orthotropic elastic parameters of fiber-reinforced polymer-matrix composites before and after the hygrothermal aging process.The results indicate that the hygrothermal aging environment significantly influences the mechanical property of a composite.The components of the parameters in the direction of the fiber bundle decreased significantly.From the accuracy analysis,we found that the actual measurement accuracy is sensitive to a shift of the horizontal edges and rotation of the vertical edges.     

Instability analysis of non-homogeneous materials under biaxial loading

In this paper, we formulate the instability problem of non-homogeneous materials under biaxial loading, in the sense that the mechanical parameters, and more specifically the strain hardening or softening and the strain-rate sensitivity, are spatially dependent functions. We first study the time behavior of the process and present the evolution of strain non-uniformities before the critical time, to show that it depends on the interplay between the spatial distribution of mechanical parameters and the initial non-uniformities. We next study the instability modes of a material exhibiting inhomogeneities along one direction by an effective instability analysis [Dudzinski, D., Molinari, A., 1991. Perturbation analysis of thermoviscoplastic instabilities in biaxial loading. International Journal of Solids and Structures 27 (5), 601–628] adapted to the non-homogeneous case. This method is based on a suitable non-uniform reference solution and allows to select the localization modes activated at different deformation levels for different “deformation paths” of the material zones.     

可重构力学超材料的设计与波动特性研究

力学超材料中的弯曲梁双稳态结构由于其主动调控性强且调控精度高等优点近年来受到广泛关注.文章利用中心受压弯曲梁的不稳定性设计了六角型双稳态结构,首先建立了等效弯曲梁模型,基于梁变形微分方程及能量最低原理探明了结构双稳态特性的产生基理,之后利用有限元数值计算研究了结构几何参数对其整体力学性能的影响,分别得到了具备自恢复及双稳态性能的结构几何参数范围,绘制了几何参数与力学性能之间的相图.同时,可重构结构的可控变形能力有助于调整整体的色散特性,利用数值仿真研究了具备双稳态特性的结构在拉伸和压缩两种构型下的色散关系,对比分析了不同结构几何参数及构型转变对结构产生的带隙位置及范围的影响,之后对由不同构型单胞组成的周期性结构进行了频响分析来验证带隙计算的准确性.通过六角型可重构结构的力学特性、色散特性研究及频响分析表明可以通过结构几何参数的设计实现对结构整体性能的主动调控,为可逆向设计的弹性波超材料结构研究分析提供了一条可靠路径.     

Monochromatic surface waves on impermeable boundaries in two-component poroelastic media

The dispersion relation for surface waves on an impermeable boundary of a fully saturated poroelastic medium is investigated numerically over the whole range of applicable frequencies. To this aim a linear simplified model of a two-component poroelastic medium is used. Similarly to the classical Biot’s model, it is a continuum mechanical model but it is much simpler due to the lack of coupling of stresses. However, results for bulk waves following for these two models agree very well indeed which motivates the application of the simplified model in the analysis of surface waves. In the whole range of frequencies there exist two modes of surface waves corresponding to the classical Rayleigh and Stoneley waves. The numerical results for velocities and attenuations of these waves are shown for different values of the bulk permeability coefficient in different ranges of frequencies. In particular, we expose the low and high frequency limits, and demonstrate the existence of the Stoneley wave in the whole range of frequencies as well as the leaky character of the Rayleigh wave.     

参数灵敏度分析的神经网络方法及其工程应用

在系统分析中，参数灵敏度分析不仅为判断各系统参数的重要性大小提供了依据，量化的灵敏度指标也是后续参数估计的前提。然而，在多效实际系统中，系统参数与系统状态间的显式函数关系不易得到，导致一阶灵敏度指标无法直接求取。简化的单因素分析方法亦存在模型粗糙、精度不高的缺点。本文研究采用人工神经网络的高精度泛化映射，通过少量样本的训练，建立复杂系统中多个系统参数与系统状态间的近似映射关系，继而推导得到统一的灵敏度计算列式。简单结构的神经网络方法和解析方法的对比计算显示了方法的有效性和可靠性。最后，应用该法对某斜拉桥结构的荷载参数和刚度参数进行了考查，得到一般性结论。     

折纸超材料折展稳态特性研究

Origami metamaterials can regulate and control macroscopic deformation by continuous deformation of their internal microstructure. Hence, mechanical properties such as Poisson's ratio, stiffness, and modulus of the metamaterials can be adjusted and designed. This work theoretically studied the mechanical behavior of origami metamaterials under the synergy of complex internal configuration and folding motion, with the methods of torsion spring equivalent and energy principle. Here, we established a mechanical model to describe the folding deformation of the metamaterials, and analyzed the influence of geometric parameters on external loading. Through parameter analysis, it is found that the external loading exhibits monotonicity during the unfolding process and does not have stability. The external loading can show three situations such as monotonicity, mono-stability and bi-stability during the folding process, which are closely related to the parameters. This work provides important guidance to origami metamaterials for improving the design of the configuration and the regulation of their performance.      

三维编织复合材料模量的双尺度有限元计算

针对三维编织复合材料的力学性能进行了双尺度有限元(TSA)数值计算,给出了计算模型和算法过程,并将数值结果与文献中的实验数据进行了比较,验证了算法的物理准确性。编织复合材料的力学性能不仅依赖于材料的基本组份,也与细观构造相关。双尺度有限元计算可以数值模拟出三维编织复合材料的整体力学性能,从而为材料的研发提供指导。本文的双尺度有限元三维数值计算方法可以推广到其他增强／孔隙等多相复合材料的数值模拟。     

DYNAMIC RESPONSE OF SATURATED POROUS ELASTIC FOUNDATION UNDER POROSITY ANISOTROPY$^{\bf 1)}$

由于不同的沉积条件和应力状态, 天然土体通常表现出一定的各向异性特征. 文章研究地基上表面受温度载荷和机械载荷时, 孔隙率各向异性参数变化对饱和多孔弹性地基热-水-力耦合动力响应问题的影响. 基于Lord-Shulman广义热弹性理论, 结合孔隙率各向异性基本假设, 建立了孔隙率各向异性饱和多孔弹性地基热-水-力耦合动力响应模型, 利用正则模态法推导出无量纲竖向位移、超孔隙水压力、竖向应力和温度分布的解析表达式并加以图示. 正则模态法是一种利用加权残差求得解析解的方法, 相较于其他方法能快速求解偏微分方程. 当孔隙率各向异性参数为1时, 可将该各向异性耦合动力响应模型退化为热-水-力耦合动力响应模型验证该地基模型的合理性. 着重分析了孔隙率各向异性参数变化对不同物理量的影响. 结果表明: 孔隙率各向异性参数变化对物理量均有一定影响. 在地基上表面受温度载荷作用时, 对超孔隙水压力和竖向应力影响最为明显; 在地基上表面受机械载荷作用时, 对超孔隙水压力和温度影响明显. 整体而言, 无论地基上表面受何种载荷, 随着各向异性参数增大, 峰值逐渐减小, 在地基深度增加方向峰值所在位置向靠近地基上表面方向移动.     

孔隙率各向异性下饱和多孔弹性地基动力响应

由于不同的沉积条件和应力状态, 天然土体通常表现出一定的各向异性特征. 文章研究地基上表面受温度载荷和机械载荷时, 孔隙率各向异性参数变化对饱和多孔弹性地基热-水-力耦合动力响应问题的影响. 基于Lord-Shulman广义热弹性理论, 结合孔隙率各向异性基本假设, 建立了孔隙率各向异性饱和多孔弹性地基热-水-力耦合动力响应模型, 利用正则模态法推导出无量纲竖向位移、超孔隙水压力、竖向应力和温度分布的解析表达式并加以图示. 正则模态法是一种利用加权残差求得解析解的方法, 相较于其他方法能快速求解偏微分方程. 当孔隙率各向异性参数为1时, 可将该各向异性耦合动力响应模型退化为热-水-力耦合动力响应模型验证该地基模型的合理性. 着重分析了孔隙率各向异性参数变化对不同物理量的影响. 结果表明: 孔隙率各向异性参数变化对物理量均有一定影响. 在地基上表面受温度载荷作用时, 对超孔隙水压力和竖向应力影响最为明显; 在地基上表面受机械载荷作用时, 对超孔隙水压力和温度影响明显. 整体而言, 无论地基上表面受何种载荷, 随着各向异性参数增大, 峰值逐渐减小, 在地基深度增加方向峰值所在位置向靠近地基上表面方向移动.      

Three-dimensional steady-state closed form solution for multilayered fluid-saturated anisotropic finite media due to surface/internal point source

A three-dimensional(3 D)steady-state solution of fluid saturated anisotropic finite media is presented.The eigenequation method and the pseudo-Stroh formalism are used to obtain the exact solution for homogeneous saturated finite media.The propagator matrix method is introduced to deal with the corresponding multilayered poroelastic media.The poroelastic solutions due to surface or internal point fluid source are obtained.The comparison of the results of the saturated isotropic media in a half space and those obtained by the finite element method is given to illustrate the accuracy of the solution in a finite domain.Numerical solutions of a sandwich poroelastic medium are presented to analyze its hydromechanical behaviors.Two ratios of the horizontal permeability to vertical permeability and different source positions are investigated.The results show that the fluid parameters and source positions have great influence on the hydromechanical behaviors of the layered media.     

Poroelastic Response of Orthotropic Fractured Porous Media

An algorithm is presented for inverting either laboratory or field poroelastic data for all the drained constants of an anisotropic (specifically orthotropic) fractured poroelastic system. While fractures normally weaken the system by increasing the mechanical compliance, any liquids present in these fractures are expected to increase the stiffness somewhat, thus negating to some extent the mechanical weakening influence of the fractures themselves. The analysis presented in this article quantifies these effects and shows that the key physical variable needed to account for the pore-fluid effects is a factor of (1 − B), where B is Skempton’s second coefficient and satisfies 0 ≤ B < 1. This scalar factor uniformly reduces the increase in compliance due to the presence of communicating fractures, thereby stiffening the fractured composite medium by a predictable amount. One further aim of the discussion is to determine the number of the poroelastic constants that needs to be known by other means to determine the rest from remote measurements, such as seismic wave propagation data in the field. Quantitative examples arising in the analysis show that, if the fracture aspect ratio a_f @ 0.1{a_f \simeq 0.1} and the pore fluid is liquid water, then for several cases considered, Skempton’s B @ 0.9{B \simeq 0.9}, and so the stiffening effect of the pore-liquid reduces the change in compliance due to the fractures by a factor 1 - B @ 0.1{1 - B \simeq 0.1}, in these examples. The results do, however, depend on the actual moduli of the unfractured elastic material, as well as on the pore-liquid bulk modulus, so these quantitative predictions are just examples, and should not be treated as universal results. Attention is also given to two previously unremarked poroelastic identities, both being useful variants of Gassmann’s equations for homogeneous—but anisotropic—poroelasticity. Relationships to Skempton’s analysis of saturated soils are also noted. The article concludes with a discussion of alternative methods of analyzing and quantifying fluid-substitution behavior in poroelastic systems, especially for those systems having heterogeneous constitution.     

纳米增强预浸料缠绕立管参数敏感性分析

复合材料立管结构参数与立管力学性能存在紧密的关联性,掌握设计参数对缠绕立管力学性能的影响规律非常必要．考虑缠绕工艺引起的孔隙缺陷,基于复合材料弹性理论及孔隙微观理论模型,构建了纳米增强预浸料缠绕立管形变模型;其次,采用局部敏感度分析法,分别讨论了缠绕立管应力与位移对孔隙含量、孔隙形状、纳米颗粒含量、纤维体积含量、缠绕层数、缠绕角度及内压的局部敏感性;最终,数值结果表明：孔隙含量对轴向应力的影响局部敏感性较强,在孔隙含量为2.8 %~3.2 %之间时,轴向应力的局部敏感性数值快速从负数增长到正值;而孔隙的引入,使得其它设计参数对立管应力与位移的局部敏感性数值的变化趋势无明显影响,但大大影响了敏感性的数值大小;此外,径向位移对孔隙、纳米颗粒含量等所有参数明显不敏感．     

Configuration Design Sensitivity Analysis of Kinematic Responses of Mechanical Systems∗

Abstract

A continuum-based configuration design sensitivity analysis method is presented for kinematically driven mechanical systems. Configuration design variable for mechanical systems are defined. The 3-1-3 Euler angle is used as the orientation design variable. The reassembly process for a mechanical system after a configuration design change is eliminated by imposing kinematic admissibility conditions of the velocity field. The direct differentiation method is used to derive design sensitivity equations. Numerical examples are presented to demonstrate the validity and effectiveness of the proposed method.     

Concurrent factors determine toughening in the hydraulic fracture of poroelastic composites

Brittle materials fail catastrophically. In consequence of their limited flaw-tolerance, failure occurs by localized fracture and is typically a dynamic process. Recently, experiments on epithelial cell monolayers have revealed that this scenario can be significantly modified when the material susceptible to cracking is adhered to a hydrogel substrate. Thanks to the hydraulic coupling between the brittle layer and the poroelastic substrate, such a composite can develop a toughening mechanism that relies on the simultaneous growth of multiple cracks. Here, we study this remarkable behaviour by means of a detailed model, and explore how the material and loading parameters concur in determining the macroscopic toughness of the system. By extending a previous study, our results show that rapid loading conveys material toughness by promoting distributed cracking. Moreover, our theoretical findings may suggest innovative architectures of flaw-insensitive materials with higher toughness.