过载性损伤与防护生物力学

损伤与防护生物力学(injury and protection biomechanics)是研究生物组织或器官损伤机理及其防护方法的一门交叉性学科,属于现代生物力学的重要分支.其研究目标是降低载荷环境下组织或器官的损伤程度,主要内容包括载荷造成生物组织和器官的损伤机制、损伤耐受极限以及损伤过程中的生物力学动态响应、如何改善组织和器官所处的力学环境降低其损伤程度、有效的防护装备优化设计思路.高过载性载荷由于其作用短时性和爆发性具有较高致命性, 因此,人在过载环境下的抗损伤能力已越来越成为航空器研制、汽车性能提升、运动员竞技能力提升与充分发挥的瓶颈;尤其是更快、更灵活新型飞机的出现,超音速弹射救生、大过载高增长率的机动飞行防护等问题向损伤与防护生物力学研究提出了新的挑战,同时也为损伤与防护生物力学的发展提供了新发展机遇.随着科技不断进步,航空航天、交通事故、体育运动乃至日常生活中老年人跌倒等过程中人体冲击过载性损伤越来越呈现发生率高、防护效率低等问题,一方面由于人体耐限实验会造成损伤而难以获得真实数据,另一方面生物组织具有复杂非线性及黏弹性、可再生和重建特性,涉及到如何精准描述生物组织或器官的本构关系、组织解剖学特征与其力学特性之间相关性,建立不同尺度的组织或器官损伤机理与耐受极限、防护方法及防护装备设计准则.为此,本文将主要总结过载性损伤与防护生物力学的主要研究内容和研究方法,并在此基础上针对人体在复杂过载环境下的损伤类型、损伤机制(包括生物力学和力学生物学响应)、损伤耐限及防护方法进行回顾,包括近年来该领域国内外的主要进展, 并提出该领域发展趋势.过载性损伤与防护生物力学研究对于保障和提高复杂过载环境下人体安全性具有重要意义,可为解决航空航天、交通、体育运动中广泛涉及的骨肌多轴向损伤评价方法与标准制定提供科学依据,对指导防护装备优化设计具有重要理论价值,同时该方面研究在工程仿生材料和防护装备方面具有潜在实用价值和广阔应用前景.      

半圆柱壳挠曲问题解析解精确性的研究

壳体力学已于上世纪由多位专家发展成熟,其中简支柱壳挠曲问题采用改进莱维解法的三角级数法解出,但是其解法复杂,手算难以完成．为讨论其结果的精确性,通过编写运行基于MATLAB的运算程序导出实例化解析解,与基于力学基本理论的推想假设对比,再引入有限元计算结果进行比较研究．最终发现,理论解析解应力和位移具有分布形式大致准确性,但仍存在不容忽视的细节与局部性问题．研究表明,理论解法工程意义有限,结果尚需改进．     

骨陷窝-骨细胞形状和方向对骨单元内液体流动行为的影响

骨组织内的流体流动不仅为骨细胞的生存提供了充足营养供应及代谢物排放途径,也在骨重建过程中起到关键作用. 为了更精确地阐明骨内液体流动的具体形式,这项研究利用骨陷窝-骨细胞的密度,形态和方向等参数来计算骨单元内液体的流动行为. 首先,计算出不同形状和方向的骨陷窝周围骨小管的数量及分布情况,其次利用算出的参数以及骨组织其他微结构数据来估计骨组织的渗透率和孔隙率等参数,最后根据计算所得的参数建立骨单元的多孔弹性力学有限元模型,并分析了在轴向位移载荷作用下骨陷窝形状和方向对骨单元内液体渗流行为的影响. 结果表明,在所研究的参数范围内不同骨单元模型的相同区域上,骨陷窝形状影响下的骨单元最大压力和流速比最小的分别增加了86%和18%;骨陷窝方向影响下的最大压力和流速比最小的分别增加了125%和56%. 伸长形骨陷窝对单个骨单元局部压力的影响远大于扁平形和圆形骨陷窝. 骨陷窝从0°绕$x$轴旋转到90°过程中压力是逐渐降低的,且30°,45°和60°的模型对骨单元内局部流速有显著影响. 该模型表示骨陷窝的形状和方向以及骨小管的三维分布对骨单元内液体压力和流速幅值及沿不同方向的流动差异有显著的影响. 这项研究将有助于精确量化描述骨内液体的流体行为.      

闽江福州河段环境工程地质问题的评价研究

江河环境工程地质问题是环境工程地质学研究的重要课题。江河沿岸环境地质条件是影响和决定沿江城市环境质量的重要背景 ;是控制和影响沿江河岸城市的环境容量和质量的重要因素。本文以闽江下游福州河段为例 ,概要评述了闽江福州河段的水动力条件和流水地质作用类型、强度等状况 ;分析研究了该河段的流水地质作用和地貌环境变化动态及其与沿江河岸人类各种工程建设活动的互相作用的特征和规律。对促进和提高沿江河岸“人—地系统”和谐协调平衡发展水平 ,正确规划、预测和评价沿江河岸各类工程建设发展的经济和社会效益及其前景 ,都具有重要意义.     

大型柔性航天器动力学与振动控制研究进展

随着航天重大工程的逐步实施,航天器正朝着超高速、超大尺度、多功能的方向发展,其面临的发射和运行环境也更加恶劣.航天器发射过程中的振动及其主/被动控制、在轨运行中大型柔性航天器动力学建模与动态响应分析、结构振动与飞行器姿态的混合控制等问题越来越复杂且难于处理;航天器结构的大型化和柔性化(如大阵面天线和太阳翼等)也对其地面试验和半实物仿真提出了挑战.本文着重介绍大型柔性航天器涉及到的动力学与振动控制问题,包括航天器发射过程中的整星隔振,大型柔性结构动力学建模与振动响应分析,大型柔性航天器的结构振动与姿轨控耦合动力学及其混合控制等.提炼出航天动力学与控制领域中亟待解决的若干基础科学问题,包括:多刚柔体系统动力学建模与模型降阶(涉及大变形柔性体动力学建模、多求解器合作仿真、模型降阶、组合结构动力学建模的解析方法等);复杂结构状态空间模型构建方法与能控性(涉及状态空间模型构建的理论与实验方法、复杂结构振动控制系统的能观性与能控性等);航天器姿态运动与大型柔性结构振动的混合控制律设计(涉及姿态机动与结构振动的鲁棒混合控制、执行机构与压电控制器的协同控制等).      

分数导数粘弹性模型描述的土中管桩竖向振动的频域解

在分数导数粘弹性本构模型的基础上综合考虑桩周土和桩芯土的平衡方程和几何方程建立了桩周土和桩芯土的竖向运动的控制方程．在频率域内利用分离变量法和分数导数的性质求解了桩周土和桩芯土竖向振动控制方程．考虑管桩与桩周土、管桩与桩芯土的边界连续性条件以及三角函数的正交性得到了分数导数粘弹性模型描述的土中管桩的竖向振动,通过数值分析研究了管桩和土体模型参数和几何参数对管桩的桩顶复刚度的影响规律．结果显示：桩芯土本构模型的分数导数的阶数对管桩竖向振动的影响较桩周土本构模型的阶数要小,且与频率有一定关系;桩芯土与桩周土的模型参数比τ1 和τ2 对等效阻尼的影响较对刚度因子的影响要大;管桩桩周和桩芯的直径比d 越小,管桩复刚度的实部和虚部就越大;土体力学性能对管桩竖向振动的影响要比管桩桩身力学性能的影响小．     

On the methodologies of stress analysis of composite structures: Part 1: State of the art—phenomenological and physical methodologies

The rapidly increasing technological importance of composite materials and composite structures is leading to the development of new, more advanced models of their actual response to mechanical and thermal loads. This in turn results in the development of new experimental and analytical methods for determination of the mechanical and thermal responses of such structures and materials to various loads. In this respect the reliability and the predictive power of various methods and techniques of stress analysis become very important since all the analytical, experimental and numerical methods used for the determination, prediction and optimization of the actual mechanical responses of composite structures and materials are based on the concepts of strain and stress. Because of the inherently three-dimensional stress and strain states in composite materials and structures and the wide use of viscoelastic polymers as the matrix and some reinforcing fiber materials, a more rigorous type of modelling than had been common in the past is needed of all the involved physical phenomena which influence the strain and stress states at the local and global levels. Also, a more rigorous analysis of practical consequences of the physical and mathematical simplifications is required to assure reliability and accuracy of various methods of stress analysis. The influence of the above-mentioned factors on the reliability and applicability of analytical and experimental procedures is illustrated by examples of actual material responses.Part 2 of this paper presents theories and techniques of three new methods of strain/stress analysis which have been developed on the basis of comprehensive physical models of involved phenomena: the isodyne, strain gradient and thermoelastic effect methods. Presented examples illustrate the efficacy of these methods.     

On the methodologies of stress analysis of composite structures: Part 2: New experimental approaches

The rapidly increasing technological importance of composite materials and composite structures is leading to the development of new, more advanced models of their actual response to mechanical and thermal loads. This in turn results in the development of new experimental and analytical methods for determination of the mechanical and thermal responses of such structures and materials to various loads. In this respect the reliability and the predictive power of various methods and techniques of stress analysis become very important since all the analytical, experimental and numerical methods used for the determination, prediction and optimization of the actual mechanical responses of composite structures and materials are based on the concepts of strain and stress. Because of the inherently three-dimensional stress and strain states in composite materials and structures and the wide use of viscoelastic polymers as the matrix and some reinforcing fiber materials, a more rigorous type of modelling than had been common in the past is needed of all the involved physical phenomena which influence the strain and stress states at the local and global levels. Also, a more rigorous analysis of practical consequences of the physical and mathematical simplifications is required to assure reliability and accuracy of various methods of stress analysis. The influence of the above-mentioned factors on the reliability and applicability of analytical and experimental procedures is illustrated by examples of actual material responses.Part 2 of this paper presents theories and techniques of three new methods of strain/stress analysis which have been developed on the basis of comprehensive physical models of involved phenomena: the isodyne, strain gradient and thermoelastic effect methods. Presented examples illustrate the efficacy of these methods.     

水力压裂解析模型裂缝扩展参数敏感性分析

水力压裂形成复杂裂缝网络是致密储层油气开采的重要技术,掌握水压裂缝扩展机理是控制压裂行为和优化压裂效果的关键.水压裂缝动态扩展行为涉及储层岩体、注入压裂液、压裂实施工艺等方面,其中水力压裂扩展时间、压裂液流体动力粘度系数、压裂液流体注入流速、储层岩石剪切模量成为决定裂缝扩展长度和裂缝开度的重要因素.本研究采用KGD、P...     

干燥环境下表面黏附性能调控研究进展

实现黏附可控是解决微/纳机械产品黏附失效问题和发展先进黏附材料及新型微操作/微装配技术等的有效途径. 基于该研究背景,针对干燥不带电系统,对近年来发展的几类黏附调控方法进行了综述,从表面改形、材料改性及外场控制三方面介绍了织构、涂层/薄膜和复合材料、速度、温度、外力和磁场在黏附调控方面的研究进展,分析和总结了各类方法的黏附调控机理、关键影响因素和特点,并指出了应用领域,介绍了多手段复合调控黏附方法的最新研究成果. 最后,探讨了当前黏附调控方面存在的问题和未来可能的发展方向.      

Reliability domain of the Timoshenko model in dynamics of rods and plates

Comparing the frequency spectrum of rods and plates calculated with the rotation inertia and shear strain taken into account (the Timoshenko model) with the results of computations based on the solution of the stationary dynamic problem of the theory of elasticity for canonical bodies in rectangular and cylindrical coordinates, we determine the boundaries of the reliability domain for the Timoshenko model: we obtain the admissible (in the framework of accepted refined kinematic hypotheses) values of minimum flexibilities for rods of annular and quadratic cross-sections and of maximal relative thicknesses for plates with annular and square cross-sections. We present the basic dependencies for the solution of the problems of natural frequencies and mode shapes for rods and plates mentioned above (the Timoshenko model) and briefly describe the solutions of similar problems by methods of the theory of elasticity.     

神经动力学研究进展和若干思考

神经动力学是动力学与控制学科的基础性分支,属于力学与脑科学、智能科学的国际前沿交叉学科领域,主要是通过动力学与控制的基本理论和方法,建立合理的模型来探究神经系统电生理动力学行为和脑认知功能的机理.近年来,国内外学者在神经动力学的基础研究方面取得了显著成果,包括神经元和神经元网络动力学行为的深入研究、大脑不同功能结构的建模分析以及神经疾病关联脑区的网络动力学建模与控制等.本文首先对国内外神经动力学研究领域取得的进展做了较全面的概括分析,特别是给出了建模方面的发展历程.进而,基于解析生物神经网络及其动力学的研究成果,对神经动力学未来的研究方向提出了一些思考展望,期望神经动力学的研究将助力具备较强可解释性和泛化能力的类脑智能原理和方法的突破及在重大工程中的应用.     

全柔性空间机器人运动振动一体化输入受限重复学习控制

探究基座、臂、关节全柔性影响下空间机器人动力学模拟、运动控制及基座、臂、关节三重柔性振动主动抑制的问题, 设计了不基于系统模型信息的运动振动一体化输入受限重复学习控制算法. 将柔性基座与关节等效为线性弹簧与扭转弹簧, 柔性臂视为欧拉-伯努利梁模型, 利用拉格朗日方程与假设模态法建立动力学模型, 然后, 用奇异摄动理论将模型分解为包含刚性变量与臂柔性振动的慢变子系统, 包含基座、关节柔性振动的快变子系统, 并分别设计相应的子控制器, 构成了带关节柔性补偿的一体化控制算法. 针对慢变子系统, 提出输入受限重复学习控制算法, 由双曲正切函数, 饱和函数与重复学习项构成, 双曲正切函数与饱和函数实现输入力矩受限要求, 重复学习项补偿周期性系统误差, 以完成对基座姿态、关节铰周期轨迹的渐进稳定追踪. 然而, 为了同时抑制慢变子系统臂的柔性振动, 运用虚拟力的概念, 构造同时反映臂柔性振动与系统刚性运动的混合轨迹, 提出了基于虚拟力概念的输入受限重复学习控制器, 保证基座、关节轨迹精确追踪的同时, 对臂的柔性振动主动抑制. 针对快变子系统, 采用线性二次最优控制算法抑制基座与关节的柔性振动. 仿真结果表明: 控制器适用于一般柔性非线性系统, 满足输入力矩受限要求, 实现对周期信号的高精度追踪, 有效抑制基座、臂、关节的柔性振动, 证实算法的可行性.      

基于模态柔度曲率的砖石古塔损伤识别研究

损伤识别是结构性能评估的有效途径之一,砖石古塔的模态柔度曲率是进行损伤识别的有效参数之一。为研究古塔结构损伤机制与模态柔度曲率的关系,依据兴教寺窥基塔结构模型振动台试验及数值模拟结果,分析周期性荷载作用下塔体动力响应及破坏特征,建立了塔体损伤程度与模态参数的关系模型。结果表明,古塔结构开裂破坏后,刚度减小,模态柔度曲率增大;数值计算与试验测试所得加速度及频域响应规律一致,塔体模型结构试件的开裂损伤随加载工况增多而持续累积,加载前期损伤参数变化缓慢,而后保持稳定,试验及模拟所得塔体的破坏过程一致。因此,塔体模态柔度曲率变化与古塔结构损伤累积具有一致性,可用于有效可靠识别砖石古塔结构损伤识别。     

基于分数导数本构的粘弹性土层—隧洞衬砌系统的稳态动力响应

将混凝土衬砌材料视为具有分数导数本构的粘弹性体,在频率域研究了深埋圆形隧洞粘弹性土的稳态动力响应。利用衬砌的内边界条件和混凝土衬砌与土体界面处的连续性条件,得到了粘弹性土体和衬砌稳态振动的应力和位移解析解。进行了隧洞衬砌厚度、土体阻尼比、本构阶数、材料参数比等对土体的位移和应力幅值影响的算例分析,结果表明:材料参数比对系统动力特性有较大影响,随着材料参数比的增加,位移和应力幅值减小;随着土体的阻尼比、衬砌厚度的增加,位移和应力幅值减小;当材料参数比Tσ/Tε=3时,随着分数导数阶数增大,响应幅值减小。     

含气泡油滴撞击油膜壁面时气泡的变形与破裂

油--气润滑过程中润滑油液滴受高速气流扰动易形成含气泡油滴,微气泡将对油滴撞击壁面时的运动过程以及壁面油膜 层的形成质量产生重要影响. 基于耦合的水平集--体积分数 方法,对含气泡油滴撞击油膜壁面行为进行数值模拟研究, 考察含气泡油滴撞击油膜壁面时气泡的变形运动过程,探讨气泡破裂的动力学机制,分析气泡大小、碰撞速度和液体黏度等因素对含气 泡油滴撞壁过程中气泡变形特征参数的影响规律. 研究表明:含气泡油滴撞击油膜壁面后气泡会发生变形,并破裂形成膜液滴;气泡随同 液滴运动过程中,气泡内外压力和速度梯度变化是使气泡发生破裂的主要诱因. 气泡大小对气泡破裂方式影响较大,气泡较小时发生单 点破裂,而气泡较大时更容易发生多处破裂. 不同大小气泡受力差异较大,气泡大小与破裂发生时刻没有明显相关性. 碰撞速度和液体 黏度对气泡的变形、破裂和破裂发生时刻都具有一定的影响. 碰撞速度越大,油滴动能越大,更容易产生气泡变形和破裂现象. 液体黏 度增大,在油滴撞壁运动前期促进气泡变形,而在运动后期可以阻延气泡破裂行为发生.      

斜坡地质灾害的空间预测问题

结合典型实例阐明斜坡地质灾害的空间预测是灾害防治工作的前提条件与成败关键。继而讨论了“斜坡地质灾害空间预测”的内涵与工作程序,以及与斜坡稳定性评价预测相关的斜坡破坏型式,斜坡结构类型,斜坡破坏的环境条件和人类活动的影响等问题。在斜坡地质灾害空间预测的定义中,笔者将预测时段仅限于“工程年代”之内,且前提是保持现有斜坡结构和环境条件,以及可预料的人类活动影响;预测程序是先进行斜坡稳定性预测,然后再对稳定性差或较差的坡段进行灾害危险性(度)预测。关于斜坡破坏型式的划分,一是从它的影响范围和“量”与“质”的关系看都需重视其单次规模;二是要与斜坡的成因结构类型相结合,才能对现场地质调查与研究具有指导意义。为了展示斜坡的成因结构类型对稳定性的控制意义,本文提供了长江三峡工程库区和西部大开发中总结出来的碎屑岩和副变质岩斜坡的新类型系统,并对三类土质斜坡的典型结构特征和稳定性进行了简要比较。关于环境条件对斜坡稳定性影响的量化研究,笔者更看重多个非自体(非本坡段)天然模型观测和滑坡反算的成果。同时建议在全面准确的现场调研的基础上建立确定性地质模型并应用基础学科的相关原理进行单要素的量化分析和单元分析,以便进一步揭示斜坡变形破坏机理。     

The electrohydrodynamics of superimposed fluids subjected to a nonuniform transverse electric field

This study aims to investigate electrohydrodynamics of two superimposed fluids that are confined between a pair of two-dimensional flat plates and are exposed to a sinusoidal electric field in zero gravity. The goal is to identify the parameters that affect the flow structure and interface deformation using a simple closed form solution. The governing electrohydrodynamic equations are solved analytically for Newtonian and immiscible fluids in the framework of leaky-dielectric theory and in the limit of small electric field and fluid inertia. A detailed analysis of the electric and flow fields is presented and it is shown that the electric field induces sinusoidal electrical stresses at the interface, which lead to periodic convection cells. The parameters affecting the sense of flow circulation and strength are investigated and it is shown that the former depends on the relative magnitude of the electric permittivity and conductivity ratios while the latter is controlled by the relative thicknesses of the fluid layers and the ratio of the electric conductivities and viscosities of the fluids. The maximum flow strength is achieved at a relative thickness that is set by the competition between the electric and hydrodynamic effects. For small deformation, the distortion of the interface is examined using a normal stress balance at the interface, and it is shown that the degree of interface deformation scales with the square of the amplitude of the electric potential nonuniformity, while its wavenumber is twice that of the imposed potential nonuniformity. Furthermore, a zero-deformation curve is found, which delineates the region in the permittivity-conductivity space according to the sense of interface deformation. The results show that for certain ranges of fluid layer thicknesses and permittivity ratios, the interface will remain flat, despite the action of the nonuiform field.     

Tractive performance of a driven rigid wheel on soft ground based on the analysis of soil-wheel interaction

A new analytical method has been presented to predict the tractive performance of a rigid wheel running on soft ground. The resultant stress of the normal stress and the shear resistance applied around the peripherical contact part of the rigid wheel should be calculated by use of the dynamic pressure-sinkage curve measured from the plate loading and unloading test, considering the rolling locus of the wheel in the direction of the external resultant force of the effective driving force and the axle load. The effective driving force could be calculated as the difference of the driving force, i.e. the integration of shear resistance and the locomotion resistance calculated from the total amount of sinkage. As a result, the analytical relations between the driving force, the effective driving force and the slip ratio, the amount of sinkage and the slip ratio, the amount of eccentricity of resultant force and the slip ratio, and the entry angle, the exit angle and the slip ratio could be verified experimentally.