具有非局部体力矩的非局部弹性理论

本文基于非局部连续统场论的公理系统,建立了具有非局部体力矩作用的非局部弹性理论,我们证明了,在非局部弹性固体中存在着非局部体力矩,非局部体力矩引起了应力的非对称和非局部体力矩是由材料中的共价键产生的。     

非局部—非对称拟连续统理论

本文基于原子点阵模型,应用广义函数论,建立了非局部一非对称拟连续统理论.并讨论了非对称应力的物理—力学模型,即微观—宏观相结合的模型理论.     

用非局部弹性场论研究圆盘状裂纹问题

本文采用非局部弹性理论。用Love位移函数导出三维轴对称问题的非局部弹性应力的一般形式解,并求解了圆盘裂纹问题。得到了裂纹尖端区的应力是有界的,再次证实了非局部理论模型求解断裂力学问题的正确性。     

关于非局部场论的两点注记

研究了非局部场论中尚未完全解决的两个基本问题：其一为局部化体力，力矩残余之间的相关性，由此得到了一个描述两者关系的定理；其二为线性非局部弹性理论的应力边界条件的提法；文中所得到的应力边界条件不仅解释了在裂纹混合边界值问题中线性非局部弹性理论方程的解在常应力边界条件下不存在的问题，而且可以给出裂纹尖端的分子内聚力模型。     

薄宽带材局部纵向屈曲的辛弹性力学求解方法

局部纵向屈曲是普遍存在于薄宽带材生产过程的板形缺陷,是屈曲研究的难点,精确的解析求解方法对局部纵向屈曲形成机理的研究和板形质量的提高具有重要意义。本文将任意位置的局部纵向屈曲分为带材边部和内部两类,采用辛弹性力学方法直接推导得到了局部纵向屈曲区域承受不同边界约束条件时的临界屈曲应力和屈曲挠度函数,并将求解结果与有限元和相关文献结果进行了对比。结果表明：辛弹性力学方法与有限元方法相比具有相同计算精度和更高的计算效率,计算精度高于传统能量法;带材边界的约束条件对临界屈曲应力、屈曲区域几何形状和屈曲挠度函数均存在显著影响,验证了传统能量法求解的不足,有利于提高局部屈曲计算精度。     

线性非局部弹性理论的修正

本文通过考虑局部化残余力的影响对线性非局部弹性理论进行了修正，由修正后的理论所导出的应力边界条件包含了物体微观结构的长程力的作用，这个结果不仅解释了在裂纹混合边界值问题中线性非局部弹性理论方程的解在常应力边界条件下不存在的问题，而且可以自然地得到裂纹尖端的Ｂａｒｅｎｂｌａｔｔ分子内聚力模型。     

基于能量等效原理的应变局部化分析:Ⅱ.有限元解法

以非局部塑性理论为基础,应用状态空间理论,通过局部和非局部两个状态空间的塑性能量耗散率等效原理,提出了一种求解应变局部化问题的新方法,以得到与网格无关的数值解.针对二维问题的屈服函数和流动法则导出了求解非局部内变量的一般方程,并提出了在有限元环境中求解应变局部化问题的应力更新算法.为了验证所提出的方法,对1个一维拉杆和3个二维平面应变加载试件进行了有限元分析.数值结果表明,塑性应变的分布和载荷-位移曲线都随着网格的变小而稳定地收敛,应变局部化区域的尺寸只与材料内尺度有关,而对有限元网格的大小不敏感.对于一维问题,当有限元网格尺寸减小时,数值解收敛于解析解.对于二维剪切带局部化问题,数值解随着网格尺寸的减小而稳定地向唯一解收敛.当网格尺寸减小时,剪切带的宽度和方向基本上没有变化.而且得到的塑性应变分布和网格变形是平滑的.这说明,所提方法可以克服经典连续介质力学模型导致的网格相关性问题,从而获得具有物理意义的客观解.此模型只需要单元之间的位移插值函数具有C~0连续性,因而容易在现有的有限元程序中实现而无需对程序作大的修改.     

关于线性非局部弹性理论的应力边界条件

应力边界条件的提法是线性非局部弹性理论尚未解决的一个理论问题。文中针对这一问题进行了研究，所导出的应力边界条件包含了物体微观结构的长程相互作用，这个结果不仅解释了在裂纹混合边界值问题中非线性局部弹性理论方程的解在常应力边界条件下不存在的问题，而且可以自然地得到裂纹尖端的分子内聚力模型。     

非局部模型与变形局部化数值模拟

基于非线性身体场论,建立了材料非局部连续模型、变分方程及相应的实时拖带系大变形有限元数值模型,设计了这一模型的数值卷积算法,由于广义函数弱收敛定理和卷积理论,证明所提出的非局部连续模型具备收敛性和稳定性。并阐明了材料特征尺度数学物理意义,统计加权函数的选择原则,数值结果表明非局部模型描述变形局部化问题是适当的。     

混凝土单轴受拉的非局部本构模型

混凝土受拉本构行为存在很强的局部软化现象,使得单轴受拉试验无法给出应力-应变关系,而只能给出应力-位移关系。本文根据内变量理论和等效应变假设建立了基于真实应变的混凝土单轴受力本构方程,并根据Weibull分布可以描述混凝土等脆性材料断裂过程的试验现象,建立了关于弹性应变的损伤演化规律。然后,通过假设平均应变与真实弹性应变的函数关系,在应力-平均应变的本构关系中采用平均弹性应变以描述其非局部行为,而在材料的损伤演化规律中采用真实弹性应变以描述其局部行为,由此建立了单轴受拉荷载条件下的非局部本构模型。最后,对一个单调受拉试验和一个反复受拉试验的仿真结果表明所提出的非局部本构模型可以准确地模拟试验结果。     

均布载荷作用下带边缘大波纹膜片的非线性弯曲

采用轴对称旋转壳体的简化Reissner方程，研究了在均布载荷作用下具有硬中心的带边缘大波纹膜片的非线性弯曲问题．应用积分方程方法，获得了具有夹紧固定和滑动固定两种外边界的膜片的特征关系，即荷载-中心挠度曲线．作为算例，给出了夹紧固定膜片中的应力分布．     

脆性材料中应力波衰减规律与层裂实验设计的数值模拟

对混凝土、岩石类脆性材料的层裂实验进行了有限元模拟,研究了应力波在此类材料中传播的衰减规律,包括两类机制:弹性波因大尺寸试样的几何弥散产生的小幅度线性衰减、与应变率相关的黏塑性波因本构关系导致的指数衰减。在此基础上,提出了包含常数项的指数型应力波峰值拟合公式。建议采用可以忽略应力波衰减影响的细长形试样进行层裂实验。混凝土类脆性材料层裂破坏模拟结果显示,有限元模拟得到的层裂片厚度与一维应力波理论得到的结果非常吻合,验证了按一维应力波理论确定层裂强度的实验方法的有效性。通过对比3种不同入射波形下层裂片的形状和净拉应力波形,发现不对称的入射波形状更有利于实验获得平直的层裂断面和较准确的层裂强度。     

Two-phase non-linear model for blood flow in asymmetric and axisymmetric stenosed arteries

The pulsatile flow of a two-phase model for blood flow through axisymmetric and asymmetric stenosed narrow arteries is analyzed, treating blood as a two-phase model with the suspension of all the erythrocytes in the core region as the Herschel-Bulkley material and plasma in the peripheral layer as the Newtonian fluid. The perturbation method is applied to solve the resulting non-linear implicit system of partial differential equations. The expressions for various flow quantities are obtained. It is found that the pressure drop, plug core radius, wall shear stress increase as the yield stress or stenosis height increases. It is noted that the velocity increases, longitudinal impedance decreases as the amplitude increases. For asymmetric stenosis, the wall shear stress increases non-linearly with the increase of the axial distance. The estimates of the increase in longitudinal impedance to flow of the two-phase Herschel-Bulkley material are significantly lower than those of the single-phase Herschel-Bulkley material. The results show the advantages of two-phase flow over single-phase flow in small diameter arteries with stenosis.     

Stress and electric displacement analyses in piezoelectric media with an elliptic hole and a small crack

In this paper, the interactions between an elliptic hole and an arbitrary distributed small crack in plane piezoelectric medium, which are often happened in engineering problems, are discussed. The Green’s functions in a piezoelectric plate with an elliptic hole for a generalized line dislocation and a generalized line force are presented. The small crack is represented by unknown continuous distributed dislocations. By considering traction free conditions on the surface of the small crack, the problem is then reduced to a group of singular integral equations which are solved by using a special numerical technique. Accuracy of the present method is confirmed by comparing the numerical results with those in literatures for PZT-4 when the elliptic hole is degenerated into a crack. The generalized stress intensity factors of cracks and the generalized stress on the edge of the elliptic hole are shown graphically. It is shown that the small crack may have shielding or amplifying effects on the main elliptic hole or crack, which depends on the location and orientation of the small crack. The hole near a crack can significantly reduce the stress intensity factor of the crack. The direction of the electric field is important to shielding effect.     

On the higher order asymptotic analysis of a non-uniformly propagating dynamic crack along an arbitrary path

Transient mixed-mode elastodynamic crack growth along arbitrary smoothly varying paths is considered. Asymptotically, the crack tip stress field is square root singular with the angular variation of the singular term depending weakly on the instantaneous values of the crack tip speed and on the mode-I and mode-II stress intensity factors. However, for a material particle at a small distance away from the moving crack tip, the local stress field will depend not only on the instantaneous values of the crack tip speed and stress intensity factors, but also on the past history of these time dependent quantities. In addition, for cracks propagating along curved paths the stress field is also expected to depend on the nature of the curved crack path. Here, a representation of the crack tip fields in the form of an expansion about the crack tip is obtained in powers of radial distance from the tip. The higher order coefficients of this expansion are found to depend on the time derivative of crack tip speed, the time derivatives of the two stress intensity factors as well as on the instantaneous value of the local curvature of the crack path. It is also demonstrated that even if cracks follow a curved path dictated by the criterion K ₁₁ ^d =0, the stress field may still retain higher order asymmetric components related to non-zero local curvature of the crack path.     

Simulation of unsteady laminar flow in models of terminal aneurysm of the basilar artery

Blood flow dynamics play an important role in the pathogenesis and treatment of intracranial aneurysms. The evaluation of the velocity field in the aneurysm dome and neck is important for the correct placement of endovascular coils, in addition the temporal and spatial variations of wall shear stress in the aneurysm are correlated with its growth and rupture. The present numerical investigation describes the hemodynamics in two models of terminal aneurysm of the basilar artery. Aneurysm models with an aspect ratio of 1.0 and 1.67 were studied. Each model was subject to a steady, sinusoidal and physiologically representative waveform of inflow for a mean Reynolds number of 560. Symmetric and asymmetric outflow conditions in the branches were also studied.

The three-dimensional continuity and the Navier-Stokes equations for incompressible, unsteady laminar flow with Newtonian properties were solved with a commercial software using non structured grids with 61334 and 65961 cells for models 1 and 2, respectively. The grids were primarily composed of tetrahedral elements.

The intra-aneurysmal flow was unsteady for all input conditions and in both models, the flow always showed a complex vortex structure. The inflow and outflow zones in the aneurysm neck were determined. The wall shear stress on the aneurysm showed large temporal and spatial variations. The asymmetric outflow increased the wall shear stress in both models.     

多孔介质——自由流界面应力跳跃条件下流动特性解析解

对非对称多孔介质--自由流复合通道内多孔介质内部及多孔介质与自由流体界面处复杂质量、动量输运特性进行研究. 在多孔介质区采用Brinkman-extended Darcy模型并结合速度连续,剪切应力跳跃的界面条件对此复合通道内流体的传递现象进行求解,提出了考虑界面应力跳跃时非对称复合通道各区域流体运动速度及摩擦系数的解析式,分析了界面应力跳跃系数,达西数及无量纲多孔层偏心厚度对流体速度及摩擦系数的影响. 结果表明:改变界面性质可在一定条件下明显控制各区域流体速度分布;在达西数、多孔层偏心厚度一定情况下,界面应力系数的增大会使界面流速减小,而使流体摩擦系数增大,特别是界面应力系数小于0的情况下变化更明显,此时若不考虑界面应力系数则会造成较大误差. 当界面应力系数及多孔层偏心厚度均为较小负数值时,改变多孔层偏心厚度对界面速度的影响要大于改变界面应力系数的情况;而当界面应力系数及多孔层偏心厚度为较大正数值时,情况则相反. 较大达西数下,界面应力系数及多孔层偏心厚度对流体摩擦系数的影响均较大,继续减小达西数至一定程度时,界面应力系数对流体摩擦系数的影响可忽略不计而认为只与多孔层偏心厚度相关,且对较大多孔层偏心厚度更敏感.      

A domain perturbation study of steady flow in a cone-and-plate rheometer of non-ideal geometry

The method of domain perturbation developed by Joseph is used to calculate velocity and stress profiles in a slightly misaligned cone-and-plate rheometer where the cone is spinning and the plate is stationary. Results for a Newtonian fluid, a Criminale-Ericksen-Filbey fluid, an upper-convected Maxwell fluid, and a White-Metzner fluid are presented and compared with earlier results in which the cone is stationary and the plate is spinning (Dudgeon and Wedgewood, 1993). Streamlines calculated for the Newtonian fluid show a very small recirculation region near the stationary plate. Velocity and stress contours are symmetric around the plane of largest gap width. For the elastic fluids studied, streamlines are asymmetric. The fluid response lags where the fluid is dominated by memory effects. Much larger recirculation regions are calculated for fluids dominated by shear thinning. These recirculation regions contain a large fraction of the fluid in the apparatus and have the effect of changing the shape of the flow domain for the remaining fluid that rotates around the cone's axis. Elasticity also has a pronounced effect on the stress profile, indicating that the accuracy of the cone and plate may be compromised even for small mis-alignments.     

A comparison between different asymmetric abdominal aortic aneurysm morphologies employing computational fluid–structure interaction analysis

Considering representative asymmetric aneurysms in the abdominal aorta, the transient 3-D blood flow and pressure distributions as well as aneurysm wall stresses were numerically analyzed. To obtain more realistic and accurate results for blood flow fields and wall stress distributions, a coupled fluid-flow and solid–structure solver was employed. Geometric abdominal aortic aneurysm (AAA) variations studied included the degree of asymmetry, neck angle and bifurcation angle, and hence their impacts on the hemodynamics and biomechanics. The simulation results indicated that the assumption of symmetric AAA geometry may underestimate AAA-wall stress considerably. The neck angle influences the blood flow field substantially. A large neck angle, resulting in strong wall curvatures near the proximal neck, can produce aggravating blood flow patterns and elevated wall stresses (Von Mises). The iliac bifurcation angle affects blood flow patterns insignificantly but plays an important role in wall-stress concentrations. The wall stress of lateral asymmetric AAAs is higher than for the anterior-posterior asymmetric types. The maximum wall stress-site is located near the anterior distal side for the anterior-posterior asymmetric AAA and the distal side towards the asymmetric bulge in the lateral asymmetric AAA.     

Transient response of laminated plates with arbitrary laminations and boundary conditions under general dynamic loadings

The dynamic analysis of laminated plates with various loading and boundary conditions is presented employing generalized differential quadrature (GDQ) method. The first-order shear deformation theory is considered to model the transient response of the plate. The GDQ technique together with Newmark integration scheme is employed to solve the system of transient equations governing dynamics of the plate. Different symmetric and asymmetric lamination sequences together with various combinations of clamped, simply supported, and free boundary conditions are considered. Particular interest of this study regards to asymmetric orthotropic plates having free edge and mixed boundary conditions. It is shown that the method provides reasonably accurate results with relatively small number of grid points. Comparison of the results with those of other methods demonstrates a very good agreement. It is also revealed that the present method offers similar order of accuracy for all variables including displacements and stress resultants.