含有多个涂层压电纤维复合材料的平面问题研究

压电纤维在未来的复合材料结构健康监测中具有重要作用.本文基于横观各向同性压电材料位移和应力连续条件以及经典的复势函数理论,讨论了同时受到平面内机械载荷和出平面电载荷作用时含有多个带涂层压电纤维的无限大线弹性基体的平面力学问题.首先将线弹性基体、涂层和压电纤维的应力场、位移场表示成复势函数,然后通过横观各向同性压电材料和线弹性材料的位移和应力连续条件确定复势函数表达式.将得到的复势函数表达式代入线弹性基体、涂层和压电纤维的的应力场、位移场公式可确定其应力场和位移场.最后,通过定量的案例讨论了涂层的材料属性对线弹性基体应力场的影响.案例分析表明涂层的材料属性对压电复合材料的应力场有重要的影响.     

铁中Ⅰ型裂纹裂尖场与裂尖形变机制分析

应用有限元及连续介质力学计算了平面应交条件下体心立方铁中不同取向裂纹的裂尖应力场和弹性应变能密度分布,从宏微观相结合角度分析了裂尖场分布与裂尖微结构演化的相互关联.指出裂尖塑性变形的具体形式与裂尖滑移面上分切应力的大小密切相关,并从能量角度解释了裂尖相变产生的原因,最后结合相关分子动力学研究成果探讨了裂尖奇异性区域.     

不同排水条件下非饱和土中柱孔扩张问题的解析分析

现有的圆柱孔扩张理论已可为诸如石油工程中井筒稳定性鉴定、 及旁压和圆锥贯入实验分析等提供理论依据, 但在非饱和地基压力注浆, 复合地基处理等实际工程问题中却鲜有应用. 基于弹塑性理论和非饱和土力学原理, 采用统一强度理论, 对非饱和土中柱形小孔扩张问题进行了解析研究. 首先将柱孔周围土体分为弹性区和塑性区, 并考虑在弹性区遵循小应变理论, 在塑性区遵循大应变理论, 同时考虑了中间主应力及粒间吸力对非饱和土体强度的影响. 其次应用有效应力表示的统一强度准则, 在本构关系、几何方程、动量平衡方程等基本方程的基础上, 结合相应的边界条件, 最终获得了不同排水条件下柱孔扩张时周围弹塑性区域内的应力场、应变场、位移场及极限扩孔压力的解析表达式. 通过数值算例和参数分析, 在与现有的饱和及非饱和土中柱孔扩张理论进行退化验证的同时, 分析了吸力、剪胀参数、中主应力效应参数及初始径向有效应力等对弹塑性区域内的应力场、应变场及位移场的影响规律, 验证了本文理论的正确性及有效性, 以期为实际工程问题提供合理的理论依据.      

一维六方准晶中纳米尺度开裂孔洞的III型断裂力学

研究了一维六方准晶中纳米尺度开裂孔洞的Ⅲ型断裂力学问题。基于复变弹性理论和表面弹性理论获得了考虑表面效应时椭圆孔边裂纹的应力场、应力强度因子和能量释放率的解析表达；讨论了缺陷尺寸、裂纹/孔洞比、耦合系数和施加载荷对应力强度因子和能量释放率的影响。研究表明：考虑表面效应且缺陷的尺寸在纳米尺度时，声子场和相位子场的无量纲应力强度因子以及无量纲能量释放率具有明显的尺寸依赖；裂纹相对尺寸较小时，表面效应对声子场和相位子场的无量纲应力强度因子影响较小；纳米尺度时无量纲能量释放率随耦合系数的增加而增大；耦合系数一定时，无量纲能量释放率受到椭圆孔尺寸影响；随着声子场载荷的增大，无量纲能量释放率先减小后增加，最后趋于稳定；无量纲能量释放率随相位子场载荷的增大单调减小，非常小和非常大的声子场载荷(或相位子场载荷)屏蔽了相位子场载荷(或声子场载荷)的影响。     

一维六方准晶中纳米尺度开裂孔洞的III型断裂力学

研究了一维六方准晶中纳米尺度开裂孔洞的Ⅲ型断裂力学问题。基于复变弹性理论和表面弹性理论获得了考虑表面效应时椭圆孔边裂纹的应力场、应力强度因子和能量释放率的解析表达；讨论了缺陷尺寸、裂纹/孔洞比、耦合系数和施加载荷对应力强度因子和能量释放率的影响。研究表明：考虑表面效应且缺陷的尺寸在纳米尺度时，声子场和相位子场的无量纲应力强度因子以及无量纲能量释放率具有明显的尺寸依赖；裂纹相对尺寸较小时，表面效应对声子场和相位子场的无量纲应力强度因子影响较小；纳米尺度时无量纲能量释放率随耦合系数的增加而增大；耦合系数一定时，无量纲能量释放率受到椭圆孔尺寸影响；随着声子场载荷的增大，无量纲能量释放率先减小后增加，最后趋于稳定；无量纲能量释放率随相位子场载荷的增大单调减小，非常小和非常大的声子场载荷(或相位子场载荷)屏蔽了相位子场载荷(或声子场载荷)的影响。     

橡胶圆柱体的平面外剪切变形分析

利用Gao Y C给出的一类不可压缩的应变能函数,分析了不均匀橡胶圆柱体受平面外剪切变形问题.通过平面外剪切变形问题的变形模式,并结合应力和位移的边界条件以及界面连续条件得到了位移场和应力场在瞬时构型中的解析解,并讨论了材料的不均匀性对位移场及应力场分布的影响.结果表明:不均匀模型可以方便地退化成均匀模型;橡胶圆柱体的不均匀性对位移场和应力场的分布有着重要影响.     

应力诱发界面扩散下夹杂演化的有限元模拟

应力作用下内连导线中的夹杂等缺陷会迁移和变形,从而影响电子器件的可靠性.论文基于应力诱发表面扩散机制下的弱解描述,推导了固-固界面在界面扩散机制下的有限单元控制方程,并数值分析了应力诱发界面扩散下金属内连导线中夹杂的形貌演化.研究结果表明:在拉压应力下,夹杂会发生分节或者圆形化;在双向等值拉应力下,夹杂存在着沿长轴生长...     

理想正交各向异性弹塑性材料平面应力条件下Ⅰ型静止裂纹尖端场

在本文中,以 Hill 的塑性理论为基础,详细地讨论了理想正交各向异性弹塑性材料,平面应力条件下Ⅰ型静止裂纹尖端场解。裂纹尖端应力场不包含应力间断线,但包含弹性区。分析的结果表明(i)对于平面应力静止裂纹问题,应力场解不是唯一的,场解中的自由参数必须由远场条件来确定。(ii)裂纹尖端的应力、应变的奇异性,无论是各向异性材料还是各向同性材料,都是相同的。但在各向异性材料中,各向异性参数影响着应力、应变的幅度和分布。     

理想正交各向异性弹塑性材料平面应力条件下Ⅰ型静止裂纹尖端场

在本文中,以 Hill 的塑性理论为基础,详细地讨论了理想正交各向异性弹塑性材料,平面应力条件下Ⅰ型静止裂纹尖端场解。裂纹尖端应力场不包含应力间断线,但包含弹性区。分析的结果表明(i)对于平面应力静止裂纹问题,应力场解不是唯一的,场解中的自由参数必须由远场条件来确定。(ii)裂纹尖端的应力、应变的奇异性,无论是各向异性材料还是各向同性材料,都是相同的。但在各向异性材料中,各向异性参数影响着应力、应变的幅度和分布。     

单晶体塑性滑移有限变形下的应力计算

为了探讨和发展单晶金属材料的非弹性有限变形分析方法,提出一种单晶体各向异性弹塑性分析的计算格式.该方法是一种以初始构形为变形计算参考构形的描述方法,它对单晶体塑性构形的演化用增量计算以跟随加载路径,而在应力计算时在卸载构形的基础上用Hencky对数弹性应变来计算总量的应力以保证计算的稳定和收敛;通过求解满足瞬时屈服条件和应力与弹性应变关系的广义胡克定律的非线性方程组来搜索激活滑移系.     

粗糙面在梯度表面层上滑动接触的应力分布

对粗糙面在梯度表面层上的滑动过程进行应力分布研究,以模拟实际摩擦过程中,考虑塑性变形情况下,梯度覆层体中的应力分布规律,同时与均质体及单覆层体进行比较研究,分析了在表面载荷相同时滑动接触的应力分布。结果表明覆层体出现塑性变形后,在接触表面上的压力分布与弹性变形时有很大变化,在界面处梯度层的应力分布比单层膜更为理想,其应变梯度也较小;受相同表面载荷作用下产生塑性变形时,梯度层膜在基体产生塑性变形较小     

基于微凸体接触压力分布的粗糙表面法向接触建模

提出一种考虑微凸体弹塑性接触变形影响的粗糙表面法向接触力学模型.采用有限元模拟微凸体弹塑性接触过程,分析不同塑性屈服条件对微凸体接触载荷和实际接触面积的影响.再根据微凸体接触面上压力分布的变化规律,将微凸体的接触状态分为完全弹性接触阶段、弹塑性接触阶段、完全塑性接触阶段.分析接触面压力变化规律对微凸体法向接触载荷-变形的影响,再利用GW模型的数理统计分析的方法得到粗糙表面的法向接触载荷.将论文提出的模型与完全弹性模型、CEB模型、ZMC模型、KE模型、JG模型进行对比,并且研究了塑性指数对粗糙面接触载荷-平均高度距离的影响.结果表明,论文提出的模型能够更好地描述微凸体法向接触载荷与接触变形的变化趋势,模型预测粗糙表面法向载荷与ZMC、KE模型具有较好的一致性;粗糙面接触载荷随着平均接触距离增加而减少,随着塑性指数的增加,不同模型预测的法向接触载荷差异逐渐增大.     

Stress analysis of a layered elastic solid in contact with a rough surface exhibiting fractal behavior

A contact stress analysis is presented for a layered elastic half-space in contact with a rough surface exhibiting self-affine (fractal) behavior. Relationships for the mean contact pressure versus representative strain and the real half-contact width versus elastic properties of the layer and the substrate, asperity radius, layer thickness, and truncated half-contact width were derived from finite element simulations of a layered medium compressed elastically by a rigid cylindrical asperity. These relationships were incorporated in a numerical algorithm that was used to obtain the contact pressure distributions and stresses generated by the asperity contacts formed at the interface of the layered medium and the fractal surface. Analytical solutions illustrate the significance of the elastic material properties, layer thickness, and surface topography (roughness) on global parameters such as normal load and real contact area. Results for the contact pressure distribution and the surface and subsurface stresses provide insight into the initiation of yielding and the tendency for cracking in the layered medium. It is shown that cracking at the surface and the layer/substrate interface is more likely to occur in the case of a stiff layer, whereas surface cracking is more prominent for a relatively compliant layer.     

固定接触界面切向静弹性刚度问题研究

根据两球体单峰同时受法向、切向载荷时微滑切向应力的分布以及MB模型,给出了界面的总切向接触静弹性条件刚度、总条件法向载荷的解析解。将切向接触静弹性条件刚度的解析解嵌入到有限元软件中,获得了整机的理论模态。以一款八四七厂华中工学院XHK5140型自动换刀计算机数控立式镗铣床上的结合部为研究对象,通过实验对解析解进行了定量验证。研究结果表明：在理论振型与实验振型一致的条件下,界面模型的相对误差在-19．2％～16．8％之间。     

Contact analysis in the presence of an ellipsoidal inhomogeneity within a half space

Many materials contain inhomogeneities or inclusions that may greatly affect their mechanical properties. Such inhomogeneities are for example encountered in the case of composite materials or materials containing precipitates. This paper presents an analysis of contact pressure and subsurface stress field for contact problems in the presence of anisotropic elastic inhomogeneities of ellipsoidal shape. Accounting for any orientation and material properties of the inhomogeneities are the major novelties of this work. The semi-analytical method proposed to solve the contact problem is based on Eshelby’s formalism and uses 2D and 3D Fast Fourier Transforms to speed up the computation. The time and memory necessary are greatly reduced in comparison with the classical finite element method. The model can be seen as an enrichment technique where the enrichment fields from the heterogeneous solution are superimposed to the homogeneous problem. The definition of complex geometries made by combination of inclusions can easily be achieved. A parametric analysis on the effect of elastic properties and geometrical features of the inhomogeneity (size, depth and orientation) is proposed. The model allows to obtain the contact pressure distribution – disturbed by the presence of inhomogeneities – as well as subsurface and matrix/inhomogeneity interface stresses. It is shown that the presence of an inclusion below the contact surface affects significantly the contact pressure and subsurfaces stress distributions when located at a depth lower than 0.7 times the contact radius. The anisotropy directions and material data are also key elements that strongly affect the elastic contact solution. In the case of normal contact between a spherical indenter and an elastic half space containing a single inhomogeneity whose center is located straight below the contact center, the normal stress at the inhomogeneity/matrix interface is mostly compressive. Finally when the axes of the ellipsoidal inclusion do not coincide with the contact problem axes, the pressure distribution is not symmetrical.     

全长粘结式锚杆的非局部摩擦分析

本文对Oden等提出的非局部摩擦模型进行了修正,得到了修正后的非局部摩擦模型.并应用于全长粘结式锚杆问题.文中利用Mindlin问题的位移解导出的全长粘结式锚杆沿杆体表面所受的剪应力的弹性解,对全长粘结式锚杆进行非局部摩擦分析,得到了在修正后的非局部摩擦模型下的全长粘结式锚杆的侧剪应力.文中修正的非局部摩擦模型下的全长粘结式锚杆的剪应力分布规律与试验得到的结果以及局部摩擦模型(库仑模型)下的计算结果进行了对比,结果表明修正的非局部摩擦模型下的计算结果与实验更符合.因此在描述此锚杆剪应力时,非局部摩擦模型比局部摩擦模型(库仑模型)更接近的实际.     

Equilibrium of a Prestressed Strip Reinforced with Elastic Plates

The contact problem for a prestressed elastic strip reinforced with equally spaced elastic plates is considered. The Fourier integral transform is used to construct an influence function of a unit concentrated force acting on the infinite elastic strip with one edge constrained. The transmission of forces from the thin elastic plates to the prestressed strip is analyzed. On the assumption that the beam bending model and the uniaxial stress model are valid for an elastic plate subjected to both vertical and horizontal forces, the problem is mathematically formulated as a system of integro-differential equations for unknown contact stresses. This system is reduced to an infinite system of algebraic equations solved by the reduction method. The effect of the initial stresses on the distribution of contact forces in the strip under tension and compression is studied     

A Prestressed Elastic Strip with Elastic Reinforcements

A contact problem is studied for a prestressed elastic strip with an elastic reinforcement. The integral Fourier transform is used to construct an influence function for an infinite strip with one face fixed. A unit concentrated force is applied to the strip at an arbitrary angle. The contact problem on force transfer from a thin infinite stringer to the prestressed strip is solved. The problem is mathematically formulated as a system of integro-differential equations for the unknown contact stresses on the assumption that the beam bending model and the uniaxial stress model are valid for the stringer, which is subjected to both vertical and horizontal forces. This system is solved in a closed form using the integral Fourier transform. The contact stresses are expressed in terms of Fourier integrals in a quite simple form. The influence of the initial stresses on the contact stress distribution is analyzed, and effects of concentrated load are revealed     

表面纹理磁盘滑动接触的温度和应力及退磁临界条件的研究

采用应力场和温度场耦合的有限元方法,计算磁头/磁盘滑动接触下铝基磁盘磁层内瞬态温度场和应力场及退磁临界条件,分析接触压力、滑动速度、摩擦系数以及磁盘表面纹理对磁层内最大摩擦温升值和最大应力值的影响.结果表明:波形纹理表面瞬间滑动接触所产生的温度分布呈波形特征,表面纹理越尖锐,磁层内的温度和应力越大;滑动速度对磁层内温度的影响大于对应力的影响;当磁层最大应力小于1.2 GPa时,所对应的速度和压力为安全工况;当温升大于180 K时所对应的工况将导致磁盘退磁.     

Agricultural tire deformation in the 2D case by finite element methods

The mechanical characteristics of the rubber tire and the interaction between a tire and a rigid surface were investigated by a two-dimensional (2D) finite element (FE) model. The FE model consists of a rigid rim and a rigid contact surface which interact with the elastic tire. Four distinct sets of elastic parameters are used to represent beads, sidewall, tread and lugs. Several sets of tire loads and inflation pressures were applied to the FE model as boundary conditions, together with various displacements and friction conditions. The deformation of the tire profile, the tire displacements in the vertical and lateral directions, the normal contact pressures, the frictional forces and the stress distribution of the tire components were investigated by the 2D FE model under the above boundary conditions. The calculated tire deflections were compared with the measured data. The results show a good fit between calculated and measured data, especially at high load and inflation pressure. The comparison shows that the FE analysis is suitable to predict aspects of the tire performance like its deflection and interactions with the contact surface. Compared with the experimental methods, the FE methods show many advantages in the prediction of tire deformation, contact pressure and stress distribution.