变截面Timoshenko梁的单元刚度矩阵

变截面构件在工程中应用广泛,在对变截面梁进行数值计算时,需要建立变截面梁单元的刚度矩阵。该文采用势能驻值原理,考虑了轴力引起的几何非线性和剪切变形的影响,将梁截面刚度的变化率作为小量,得到了近似到二阶的单元刚度矩阵。在构造位移模式时,从梁的微分平衡方程出发,得到同样近似到二阶、分别以三次和五次多项式表示的剪切和弯曲位移模式。该文还证明了单元刚度矩阵的奇异性,给出了轴压刚度的表达式,定量论证了与某些精确解的误差,表明在一定范围内,该文的结果具有足够的精度。最后以一个计算实例说明该文的单元刚度矩阵具有较快的收敛性。     

DEM speedup： Stiffness effects on behavior of bulk material

A number of techniques exist for minimizing the computational cost of discrete element simulations （DEMs）. One such method is a reduction of particle stiffness, which allows for bigger time steps and therefore fewer iterations in a simulation. However, the limits and drawbacks of this approach are still unclear, and may lead to invalid results. This paper investigates the effect of a stiffness reduction on bulk behavior by examining three case studies. Two cases demonstrate that particle stiffness can be reduced without affecting the bulk material behavior, whereas the third test shows that a stiffness reduction influences the bulk behavior.     

金属丝网橡胶压缩力学性能研究

对金属丝网橡胶进行了静态压缩试验。利用控制变量法研究了压缩量、相对密度、金属丝丝径、丝材和承压面积对金属丝网橡胶压缩力学性能的影响,并对平均刚度和能量耗散系数随压缩量和相对密度的变化关系进行了研究。试验结果表明:随着压缩量的增加,金属丝网橡胶非线性力学特性逐渐增强;相对密度越大,金属丝网橡胶承压能力越强;金属丝的丝径和丝材主要影响金属丝网橡胶非线性阶段的力学特性,丝径越大,丝材越硬,承压能力越强;承压面积越大,金属丝网橡胶的承压和耗能性能越好;随着压缩量的增加,平均刚度增大,承压能力增强,能量耗散系数减小,减震性能降低;随着相对密度的增加,平均刚度和能量耗散系数均增大,承压能力和减震性能均增强。     

饱和粘弹性土层中悬浮桩的纵向振动

基于不可压饱和粘弹性多孔介质理论和一维杆件模型,利用虚拟土桩的概念,研究了不可缮压饱和粘弹性土层中悬浮桩纵向振动的动力刚度.通过势函数解耦和三角级数解法,得到了悬浮桩纵向谐调振动的半解析解和桩头动力刚度.在此基础上通过数值计算,研究了物性参数、几何参数等对刚度的影响,结果表明:悬浮桩纵向振动的动力刚度受桩长、土的软硬程度、激振频率等的影响,对比了细长桩和短粗桩的响应差异,这些结果可为悬浮桩-土相互作用的数值分析和工程设计提供参考.     

任意圈层径向非均质土中桩的纵向振动特性

研究三维轴对称条件下径向非均质土中桩的纵向振动. 首先将桩周土体沿径向分为任意圈层来考虑土体的径向非均质性,每个圈层土体均质; 然后结合边界条件和相邻圈层土体之间接触面上位移和应力连续条件,对任意圈层土体动力平衡方程由外而内逐层求解,进而利用桩-土完全耦合条件求解桩动力平衡方程,得到桩顶的频域响应解析解和时域响应半解析解; 最后通过对土体主要控制参数的研究,得出了土体径向非均匀性对桩-土动力响应的影响.      

单搭拉剪胶焊接头的应力和刚度分析

胶焊复合连接技术兼具胶接和点焊的优点,它因提高结构强度和刚度而使车辆轻量化成为可能.以胶接理论为基础,考虑了被粘体剪切应变,将焊点视为大剪切弹性模量胶粘剂,胶层和焊点沿其厚度方向的剪应力不变,建立了胶焊单搭拉剪接头的线弹性应力解析模型.在应力模型基础上,将接头各组成部分看作是独立承载的拉伸和剪切弹簧单元,得到了其刚度解析模型.应力模型中的正应力和剪应力与有限元解吻合得较好,证明理论模型正确,参数研究中确定了影响胶焊单搭接头应力分布特征的关键耦合参数.     

混凝土的弹塑性损伤双面本构模型

针对混凝土材料拉压应变空间下损伤机制的不同,结合连续损伤力学和塑性理论建立了一个全新的本构模型。该模型中损伤和塑性变形的演变由应变空间的同一个非弹性曲面来控制,但对拉压应变空间中非弹性曲面的演变分别采用了随动强化法则和各向同性演化规律。计算结果表明,该模型能较好地描述混凝土材料在单轴及多轴单调加载和低周反复荷载下的典型非线性特征。     

人体手臂阻抗控制研究

阻抗控制对于人体运动稳定性,尤其在进行不稳定任务时有着不可替代的贡献.本文由屈伸肌活性在实验测量和理论计算的误差得到启发,提出了即使在手臂平面自由运动中,肌肉的共同收缩对阻抗也有贡献的观点,运用神经肌肉骨骼模型及最优控制方法,计算出手臂运动轨迹,运动部分代人原有模型得到相似结论,把公式中的相互作用项解释为阻抗贡献,并综合多篇相关文献数据,推算出共同收缩阻抗,修正了旧有模型刚度椭圆单调变化的不足,通过仿真得到与实验基本一致的刚度椭圆图,并附加得出前后自由运动时手臂的共同收缩比左右运动时的幅度要大,对总体刚度椭圆的贡献更多".     

Characterization of cuttlebone for a biomimetic design of cellular structures

Cuttlebone is a natural material possessing the multifunctional properties of high porosity, high flexural stiffness and compressive strength, making it a fine example of design optimization of cellular structures created by nature. Examination of cuttlebone using scanning electron micros- copy （SEM） reveals an approximately periodic microstruc- ture, appropriate for computational characterization using direct homogenization techniques. In this paper, volume fractions and stiffness tensors were determined based on two different unit cell models that were extracted from two different cuttlefish samples. These characterized results were then used as the target values in an inverse homogenization procedure aiming to re-generate microstructures with the same properties as cuttlebone. Unit cells with similar topologies to the original cuttlebone unit cells were achieved, attaining the same volume fraction （i.e. bulk density） and the same （or very close） stiffness tensor. In addition, a range of alternate unit cell topologies were achieved also attaining the target properties, revealing the non-unique nature of this inverse homogenization problem.     

GENERALIZED REVERBERATION MATRIX FORMULATION FOR WAVE PROPAGATION IN MULTILAYERED MEDIUM

<正>Generalized reverberation matrix(GRM)formulation is presented to investigate elastic wave propagation in a complex multilayered solid by the combination of reverberation-ray matrix(RRM)method and stiffness matrix(SM)method.RRM method formulates a reverberation matrix,which reflects the reflection or refraction of the elastic waves in the multilayered solid.However,the dimension of RRM increases as the sublayer number increases,which may result in lower calculation efficiency of the generalized rays.SM formulation yields a system matrix of the constant dimension to promise higher calculation efficiency,but it is difficult to identify the generalized rays.In order to calculate the generalized rays in the complex multi-layered solid efficiently,the RRM formulation is applied to the interested sublayer for the evaluation of the generalized rays and SM formulation to the other sublayers,to construct a generalized reverberation matrix of the constant dimension,which is independent of the sublayer number.Numerical examples show that GRM formulation has higher calculation efficiency for the generalized rays in the complex multilayered-solid configuration compared with RRM formulation.