含抗转能力散粒体的宏微观力学特性数值分析

颗粒间抗转因素对散粒体宏微观特性的影响十分显著。本文采用离散单元法,将蒋明镜（2005年）等提出的抗转动接触模型植入离散元分析软件PFC2D中,模拟了大量刚性边界下的砂土双轴试验,研究了压缩过程中试样的宏观力学性质及部分微观参量的变化情况,用以定性分析颗粒间抗转作用对散粒体的宏微观力学特性的影响。研究表明：颗粒间的抗转...     

颗粒离散元模拟堆积碎石土变形的参数灵敏度分析

采用颗粒离散单元方法,实现任意形状块石的模拟,基于室内试验数据标定滑坡坡体物质的堆(残、坡)积碎石土的细观参数,并考虑试样尺度效应,对影响堆积碎石土宏观变形特征的细观参数进行虚拟试验敏感度分析。研究表明,1以三轴剪切试验获得测试数据并标定颗粒离散元细观参数具有相对误差在5%以内的可靠性;2虚拟模型1(101mm×200mm)与模型2(300mm×600mm)对标定的细观参数具有尺度效应,但其相对误差控制在9%以内;3离散元颗粒间的摩擦系数与碎石土的内摩擦角、抗剪强度及残余强度为非线性正相关,摩擦系数每增加0.1,峰值主应力差平均增加118.85kPa,残余强度平均增加90.44kPa;4围压越大,材料的剪胀性越弱,围压在100kPa~500kPa时,剪胀特征值K在3~6之间变化,随着围压增加,模型破坏时的粘结力以近线性增加;5杨氏模量越大,碎石土的抗剪强度越大,但两者之间并不成线性关系,且不同杨氏模量对材料的残余强度没有显著影响。     

复杂应力条件下脆性材料的受拉破坏准则

根据脆性材料微粒对之间距离超过特定距离而破坏的微观机理,结合宏细观试验结果指出的绝大多数的破坏颗粒对的方向与σ1的方向之间的夹角都大于0°这一事实,导出了新的脆性材料张破坏准则: σ1 ασ2 βσ3≤[σ],式中系数β与材料的优势角θ(破坏颗粒对的方向与σ1的方向之间的夹角的某种平均) 和泊松比μ有关．同时还指出了经典强度理论由于没有考虑优势角从而导致了它们与宏观实验结果不符．实验结果表明,在张破坏条件下此准则比经典的准则有更高的精度．     

颗粒材料三维应力路径下的接触组构特性

颗粒材料的宏观应力变形特征与其微观接触力、组构等紧密相关.一般而言,强接触系统属于颗粒内部体系的传力结构,其对应的组构张量是影响宏观应力性质的重要因素.细观数值方法(如离散单元法)能够反映物理试验的基本规律,并且可以方便地提取宏微观数据来研究颗粒体系的应力变形机制.采用离散单元法(discrete element method,DEM)进行一系列等$p$等$b$应力路径下颗粒材料的真三轴试验,在此基础上研究了三维应力路径下颗粒材料的宏微观力学参数的演化过程、三维组构张量与应力张量多重联系以及强接触体系反映的宏观应力特征.研究表明:颗粒体系偏应力峰值状态和临界状态均存在与加载路径无关的宏微观特征;三维应力路径下组构张量与应力张量存在非共轴性,但其联合不变量演化过程表现出加载路径无关的特征;与弱接触系统的组构张量相比,强接触系统的组构张量更能反映宏观应力张量的特征;强弱接触体系的组构张量对颗粒体系宏观响应的贡献不同,其分界点存在一定取值范围,但采用平均接触力较为简单合理.      

颗粒材料三维应力路径下的接触组构特性

颗粒材料的宏观应力变形特征与其微观接触力、组构等紧密相关.一般而言,强接触系统属于颗粒内部体系的传力结构,其对应的组构张量是影响宏观应力性质的重要因素.细观数值方法 (如离散单元法)能够反映物理试验的基本规律,并且可以方便地提取宏微观数据来研究颗粒体系的应力变形机制.采用离散单元法(discrete element method, DEM)进行一系列等p等b应力路径下颗粒材料的真三轴试验,在此基础上研究了三维应力路径下颗粒材料的宏微观力学参数的演化过程、三维组构张量与应力张量多重联系以及强接触体系反映的宏观应力特征.研究表明:颗粒体系偏应力峰值状态和临界状态均存在与加载路径无关的宏微观特征;三维应力路径下组构张量与应力张量存在非共轴性,但其联合不变量演化过程表现出加载路径无关的特征;与弱接触系统的组构张量相比,强接触系统的组构张量更能反映宏观应力张量的特征;强弱接触体系的组构张量对颗粒体系宏观响应的贡献不同,其分界点存在一定取值范围,但采用平均接触力较为简单合理.     

基于三维离散元堆积碎石土细-宏观力学参数相关性研究

以室内三轴压缩试验为基础,运用PFC3D离散元软件生成含有不同形状碎石块的1∶1堆积碎石土数值模型,标定得到相应的细观参数(相对误差控制在10%以内)。改变细观参数,定量探讨其与宏观力学参数(粘聚力和内摩擦角)关系,得到细观参数与粘聚力、内摩擦角的细-宏观参数关系式,并对其可靠性进行了验证。研究发现,(1)粘聚力与接触粘结力、颗粒摩擦系数呈线性正相关;同时受粘结力比值M=σc/τc的影响,当M1时,粘聚力随着M的增加而线性增加;M≥4时,粘聚力趋于稳定。(2)内摩擦角与颗粒摩擦系数呈对数关系,与接触粘结力呈抛物线关系。(3)细-宏观参数关系式是在含水率为7%室内试验结果基础上研究得出,将关系式运用到含水率为9%和11%的室内试验结果进行标定验证:依据含水率为9%和11%室内试验得到的c和φ值,运用细-宏观参数关系式计算得到细观参数,经验证其数值-室内试验结果相对误差在10%以内,说明通过细-宏观力学关系式得到的细观参数符合精度要求,关系式可靠。为后续土石混合体数值试验研究者进行细观参数标定工作提供了参考依据。     

被动桩中土拱效应特征与影响参数研究

采用平面有限元方法,对被动桩在粘性、无粘性土体条件下的土拱形成机理进行了分析。从土拱效应的4个主要方面:土拱形状、桩周土体塑性(拉裂)区的分布、竖向位移等值线及桩后土体残余荷载分担比,探讨了改变桩、土参数,如粘聚力、内摩擦角、弹性模量、泊松比、剪胀角、桩土接触面特性时,土拱效应的变化规律。研究表明,土性不同时,土拱形成过程中土体的受力和破坏形式不同;泊松比、剪胀角、桩土接触特征对土拱效应的影响最明显;强度高、剪胀角大、桩土接触面粗糙或泊松比小的土体,易形成土拱。根据研究结果,指出Tom io Ito排桩塑性绕流土压力理论计算公式尚存不足。     

颗粒破碎的微观力学模拟

基于颗粒尺度建立了考虑颗粒破碎过程中能量耗散的微观力学模型。为了模拟颗粒体的破碎过程和侧向应力的影响,把胶结在一起的三个等粒径的颗粒体当作颗粒破碎体。基于颗粒破碎准则和热力学原理,考虑了颗粒破碎过程中颗粒接触点之间的塑性滑移和颗粒转动以及侧向应力作用,通过颗粒破碎体内部颗粒之间的接触破碎来模拟颗粒体的破碎。计算结果表明:本文所提出的微观力学模型可以模拟胶结在一起的大颗粒的破碎过程;颗粒之间的微观力学参数对颗粒破碎过程中的力、塑性滑移、颗粒的转动角有不同程度的影响;力-位移关系曲线在破碎前是线性的,而颗粒破碎时竖向力迅速降低到较小值,并伴有颗粒的转动。     

考虑颗粒转矩的接触网络诱发各向异性分析

颗粒材料的宏观力学行为与接触网络的组构各向异性密切相关, 根据接触点的滑动与否、转动与否和强弱力情况, 可以将颗粒间的接触系统分为不同的子接触网络. 一般而言, 不同的子接触网络在颗粒体系中的传力机制不同, 对宏观力学响应的贡献也有不同. 采用离散单元法(discrete element method, DEM)模拟了不同抗转动系数$\mu_r$下颗粒材料三轴剪切试验, 分析了剪切过程中不同子接触网络的组构张量的演变规律, 并探究了颗粒抗转动效应对子接触网络各向异性指标演变规律的影响. 研究发现: 剪切过程中转动、非转动接触的组构张量变化不是独立的, 受到颗粒间滑动与否的影响; 非滑动、强接触网络是颗粒间的主要传力结构, 非滑动接触网络的接触法向和法向接触力各向异性均随$\mu_r$的增大而增大, 其对宏观应力的贡献程度随$\mu_r$的增大而减小;强接触网络的接触法向各向异性随$\mu_r$的增大而增大, 但法向接触力各向异性随$\mu_r$的增大无明显变化, 强接触网络对宏观应力的贡献程度在不同$\mu_r$情况下均相同.      

岩土材料细观、宏观强度参数的关系研究

岩土材料细观与宏观概念的界定有个相对尺度的区别.当前对于岩土体宏观与细观参数间规律的探索主要集中在本构关系和计算分析两个方面.论文利用颗粒流,设计了无摩擦作用下3个颗粒间黏结强度增加、弱黏结和强黏结作用下19个内摩擦角在0-90度单调增加等共计41组单轴压缩数值试验.分析发现:(1)颗粒间黏结强度的增加,使得粒间较难破裂,微裂纹数量急剧减少.(2)黏结强度一定时,45度前后破坏模式截然不同,与剪切破坏角有关.(3)强黏结峰后曲线跌落明显,脆性破坏特征明显.(4)颗粒间强度参数的变化,不改变能量转移模式,黏结强度的改变使得应变能占输入能量的比例提高,即材料破坏的塑性渐变减小.弱黏结时,随着内摩擦角的增加应变能占输入能量比例逐步减小.强黏结时,应变能占输入能量比例近似,减小不明显.     

大学物理《力学》中的自由度分析

对《力学》中的物体自由度进行多方面分析,以深化教学、提高学生正 确分析物理问题的能力.使用实际教学分析的研究方法,在《力学》范围内讨论自由度与坐标、 自由与约束的关系并得以下结论: (1) 同一物体的自由度随其所在的``空间'不同而不同, 不因坐标系的选取不同而 异, 在同类参考系中不因参考系的动静而有别;(2)自由度遵循叠加原理. 讨论了质点系的总自由度及相关计算问题,并指出研究《力学》中自由度的意义.     

两自由度振动系统的斜碰撞分析

研究斜碰撞振动系统动力学的一个关键问题是对系统在碰撞前后的状态进行合理描述和正确计算．针对两弹性体斜碰撞问题，基于瞬间碰撞假设，提出了采用步进冲量来分析和求解斜碰撞前后的状态关系；并以弹簧摆和振子组成的两自由度斜碰撞振动系统为例，具体介绍了该算法如何实现．用解析方法讨论了该系统在斜碰撞过程中可能出现的各种力学现象，将冲量步进算法得到的数值解与解析结果进行对比，取得了完全一致的结果．该数值方法能适应多种斜碰撞问题的计算．     

The distribution of zeroes of solutions of neutral equations

1.IntroductionThestudyofoscillator}'behaviorofneutraldifferentialequationsisarelativealnewfieldandisveryinterestinginapplications.Inrecentyears.theoscillatorytheoryfortheseequationshasbeenextensivelydeveloped.WerefertotherecentbookbyGyoriandLadaSI'I,Baino…     

Application of the theory of branching of solutions of nonlinear equations in the estimation of pole perturbance of linear closed-loop systems

When the state and input matrices of a multivariable linear time-invariant system are perturbed, the problem of the estimation of pole perturbance of closed-loop system is considered by making use of the theory of branching of solutions for nonlinear equations.     

综合管廊早龄期混凝土温度应力仿真模拟

以汉中市综合管廊为工程背景,通过有限元软件,首先模拟分析综合管廊早龄期温度场,随后将温度场结果附加到管廊结构中求解温度应力,从而分析综合管廊早龄期混凝土温度应力与浇筑长度的关系。模拟分析发现:综合管廊在混凝土浇筑早期,温度场沿长度方向变化不显著,且随着浇筑长度的增加,温度场不变;在墙体温度上升阶段,墙体内部沿长度方向的应力主要为压应力,当墙体混凝土内部温度达到最高,混凝土开始降温,结构内变形由膨胀逐渐变为收缩,内力由压应力逐渐变为拉应力,随着时间的变化,墙体内外温差趋于稳定,墙体温度应力也逐渐平稳,此时温度应力主要受外界环境变化影响;当综合管廊浇筑长度在60m范围内,浇筑长度越大,温度应力越大,但当浇筑长度超过60m,温度应力不再增加。     

SWT-120风洞稳定段的性能测量

在普通超音速风洞中,由于受到噪声干扰,很难进行有效的边界层转捩特性试验研究,针对超声速流动特点发展较低噪声风洞十分必要,而稳定段设计的好坏直接影响到下游试验段噪声水平。本文介绍一座低噪声风洞稳定段的结构和性能测量结果。先对稳定段的结构设计做了简单介绍,然后对实验结果进行分析,实验结果表明在大角度扩散段内装置孔锥,稳定段安装消音夹层和阻尼网组等部件后,气流的速度脉动和压力脉动明显降低,其中压力脉动降低一个量级,速度脉动为1%。进一步优化设计和改进工艺,速度脉动还可进一步降低。测量结果表明SWT-120稳定段的设计是成功的,对我们以后发展更高性能的静风洞有借鉴和参考价值。     

CHARACTERISTICS OF SUBDIFFERENTIALS OF FUNCTIONS

ＣＨＡＲＡＣＴＥＲＩＳＴＩＣＳＯＦＳＵＢＤＩＦＦＥＲＥＮＴＩＡＬＳＯＦＦＵＮＣＴＩＯＮＳ（郭兴明）ＣＨＡＲＡＣＴＥＲＩＳＴＩＣＳＯＦＳＵＢＤＩＦＦＥＲＥＮＴＩＡＬＳＯＦＦＵＮＣＴＩＯＮＳ￥ＧｕｏＸｉｎｇｍｉｎｇ（ＲｅｃｅｉｖｅｄＪｕｎｅ１６，１９９５...     

Evolution of viscoelastic properties of polyurethane in the course of curing

Simultaneous measurements of the components of dynamic modulus at the set of frequencies were made by Fourier-Transform Mechanical Spectroscopy (FTMS) to investigate the curing process of a segmented polyurethane.The whole process of curing is treated as consisting of three stages. In the first stage, the material can flow and this stage is complete at the gel-point. In the second stage, the curing proceeds in the network and ends with phase separation of aggregating segments (blocks) of polymeric chains. In the third stage, the curing continues in the two-phase system. The object of investigation was the transition from the second to the third stage. It was shown that the position of the maximum of loss tangent depends on frequency and does not coincide with different characteristic points marking the time dependencies of G and G. This is due to the formation of micro-particles of the newly formed phase. The FTMS method is used for separation of the modulus changes into parts related to formation of chemical and physical networks.Special means of treating calorimetric measurements identified a weak maximum heat output at the moment corresponding to the transition point.     

One type of solutions of conduction of nervous pulses and the stability of action potential

This paper reports a type of laws which governs action poten-tial of nervous impulses,and it is discussed by general form-nonlinear dispersive process.We find that the nervous wave is aslowly varying amplitude solitary wave in the small dispersivecase.we prove that the solitary wave is not generated in the or-dinary dispersion,but a travelling wave with varying amplitudesmay be obtained.The stability of various possible action poten-tials and bifurcation in overdamped case are also discussed inthis paper.