基于最小加速度原理的刚塑性动力问题

应用有限变形的最小加速度原理，推导出了小变形或有限变形构件受到刚性飞射物撞击时的刚塑性运动方程，指出该方法具有直接、简捷和可靠的特点。     

有限变形弹塑性理论及一致性算法

通过对对数应变的引用，成功地将小变形弹塑性理论移植到有限变形之中，形成了一套完整的有限变形弹塑性理论，进一步发展了具有一阶精度的一致性弹塑性算法。最后，给出几个算例以证明此理论的正确性及算法的有效性。     

考虑晶界效应的多晶体有限变形分析

将晶界及其影响区综合考虑，建立了考虑晶界效应的力学模型，结合晶体塑性理论，利用有限变形有限元对多晶体进行数值模拟，数值结果显示了细观层次下晶粒变形场的特点，理论计算同实验定性一致。     

超静定梁变形计算的有限差分法

推导了超静定梁变形计算的有限差分方程，研究了边界条件，编制了计算程序，计算 了超静定梁的变形. 文中工作扩大了有限差分法的应用范围.     

等截面梁有限变形的传递函数增量算法

本文介绍了一种计算等截面梁有限变形的新方法－传递函数增量算法，它是一种半解析数值计算方法。此算法充分利用增量失空法Ｇａｕｓｓ求积公式计算非线性有限变形的特点，并将这些特点与传递函数方法，有效地结合起来，既避免了数值方法计算量大的困难，又使得求解高阶非线性微分方程的解析解成为可能，算例分析表明，这是一种易编程，计算量小，收敛快，求解精度高的行之有效的计算方法。     

有限空心截圆锥扭转变形问题解析解

对有限空心截圆锥扭转变形问题进行了力学分析，确定其外锥面的内柱面上的载荷分布规律，应用三维弹性理论求解了有限空心截圆锥扭转变形问题，得到了应力解析式。文中给出了实际算例。     

框架结构P-△效应分析的微分求积单元法

采用一种新的数值方法——微分求积单元法分析框架结构的P-△效应。微分求积单元法采用微分求积法直接求解微分方程的技术，并结合有限分割技术而形成。首先建立考虑剪切变形和轴力二阶效应的框架结构单元平衡微分方程，通过微分求积离散而得到梁单元的一般弹性刚度方程；同时考虑变形后节点的平衡条件和变形协调条件，导出框架结构整体二阶分析的微分求积单元法力学模型。由于该分析模型中包括了单元及结构的所有离散形式的控制方程，因此采用该模型进行结构分析可得出较为精确的解。数值算例的分析比较，表明了该法用于框架结构P-△效应分析的正确性和有效性。本文导出的框架结构二阶分析的微分求积单元法力学模型可用于框架结构剪切变形与几何非线性的耦合效应分析。     

软化梁的动力行为

由于材料的应变软化，或是由于薄壁梁在挠曲变形时的截面畸变，梁的弯矩一曲率关系可能呈现先强化后软化的性态。近几年来北京大学一个研究组的工作证实，具有软化特性的半无限长梁在自由端受常速度冲击时，软化可能发展为一个有限长的区段，并在梁内移动，同时，英国ＵＭＩＳＴ的一个研究组以核电站中的管道甩动问题为工程背景，集中研究了具有软化特性的有限长梁的动力行为。提出了计及弹性、强化、软化和梁的几何大变形的理论模型，并发展了相应的有限差分程序。这一理论模型能够预报软化区的萌生、发展和局部化，同管道甩动实验观测极为符合。     

三维有限变形线性光弹性材料研究

制成了１５％应变量的三维有限变形线性光弹性材料，可供对橡胶类零件的应力分析之用；在该材料中，以尼龙６作内含体，可以制作初应力很小的三维有限变形复合光弹性模型。     

冻土的有限变形本构关系的实验研究

根据有限变形理论,给出了冻土三轴蠕变实验数据处理时所需Green应变和Kirchhoff应力的计算公式,并根据冻土三轴蠕变实验结果给出冻土的有限变形本构关系及其蠕变参数。通过对冻土的实验数据对比分析可知:在相同试验条件下,应变较小时,小应变ε1 和Green应变Ez的数值几乎相等,随着应变的增加,两种表征方式的计算结果差别越来越大,可见对于冻土如果不考虑试件变形前后长度的变化,计算得出的应变偏离实际变形情况较大,因此冻土材料的本构关系采用有限变形表征更确切。根据冻土的有限变形本构关系计算得到的冻结壁内侧最大径向位移更接近实测结果,因此冻结壁设计计算依据冻土的有限变形本构关系更为合理。对今后冻结工程的设计具有参考价值。     

Non-elliptic deformation field near the tip of a mixed-mode crack in a compressible hyperelastic material

This paper gives an asymptotic analysis of the deformation field near the tip of an arbitrary mixed-mode crack in a compressible hyperelastic harmonic material which loses ellipticity at sufficiently large deformations. It is found that the near-tip deformation field is characterized by a localized non-elliptic deformation band issuing from the crack-tip and bounded by two curves of discontinuous deformation gradient. Explicit expression for the near-tip deformation field is obtained both inside and outside the localized deformation band. In particular, a simple relation is derived that determines the orientation of the deformation band in terms of two complex governing parameters of the near-tip fields inside and outside the deformation band, respectively.     

软材料接触力学问题的实验研究

针对硫化硅橡胶在刚性压头作用下的大变形接触问题,研究了橡胶类软材料接触区区域变形场的实验测量技术。为了研究硫化硅橡胶与刚性体在试件内部接触的大变形问题,本文提出了一种体内数字栅云纹实验技术,给出了该技术的原理和实现方法。然后运用数字图像处理方法得到了橡胶材料大变形位移u场、v场,并结合橡胶材料在压头作用下的理论分区模型,分析了软材料接触区域的大变形场,讨论了接触区域的变形分区规律,并初步给出了接触区域应变计算方法。     

Effects of unsteady deformation of flapping wing on its aerodynamic forces

Effects of unsteady deformation of a flapping model insect wing on its aerodynamic force production are studied by solving the Navier-Stokes equations on a dynamically deforming grid.Aerodynamic forces on the flapping wing are not much affected by considerable twist,but affected by camber deformation.The effect of combined camber and twist deformation is similar to that of camber deformation.With a deformation of 6% camber and 20°twist(typical values observed for wings of many insects),lift is increased bv 10%～20%and lift-to-drag ratio by around 10%compared with the case of a rigid flat-plate wing.As a result.the deformation can increase the maximum lift coefficient of an insect.and reduce its power requirement for flight.For example,for a hovering bumblebee with dynamically deforming wings(6?mber and 20°twist),aerodynamic power required is reduced by about 16%compared with the case of rigid wings.     

Numerical simulation of effect of convection-diffusion on oxygen transport in microcirculation

The entire process of oxygen transport in microcirculation by developing a 3D porous media model is calculated numerically with coupled solid deformation-fluid seepage-convection and diffusion . The principal novelty of the model is that it takes into account volumetric deformation of both capillary and tissues resulting from capillary fluctuation. How solid deformation, fluid seepage, and convection-diffusion combine to affect oxygen transport is examined quantitatively: (1) Solid deformation is more significant in the middle of capillary, where the maximum value of volumetric deformation reaches about 0.5%. (2) Solid deformation has positive influence on the tissue fluid so that it flows more uniformly and causes oxygen to be transported more uniformly, and eventually impacts oxygen concentration by 0.1%-0.5%. (3) Convection-diffusion coupled deformation and seepage has a maximum (16%) and average (3%) increase in oxygen concentration, compared with pure molecular diffusion. Its more significant role is to allow oxygen to be transported more evenly.     

Uniaxial deformation behavior of different polypropylene cast films at temperatures near the melting point

A new apparatus for stretching polypropylenes at elevated temperatures below the melting point at high deformation speeds (up to 750 mm/s) is described. In the temperature range of 140-160 °C the tensile behavior of polypropylene undergoes a shift from the ductile to the quasi-rubber-like deformation behavior. Furthermore, the deformation behavior is strongly affected by the strain rate. The homogeneity of sample deformation increases with the deformation rate. Furthermore, the influence of the cast film morphology on the uniaxial deformation behavior is investigated and discussed in terms of the degree of crystallinity and crystallite size. The degree of the residual crystallinity plays a decisive role for the deformation behavior. At similar degree of crystallinity, the yield stress is significantly influenced by the crystallite size in a way that the larger the crystallite size the higher the yield stress.     

相近韦伯数条件下激波波后液滴初期变形的影响机制

以实验结合数值模拟与理论分析的方法,研究韦伯数在2 100~2 700区间内,不同组合流动参数对液滴破碎初期变形的影响与作用机制。实验中通过高速摄影捕捉到一系列具有明显差异的液滴变形模态,表明在相近韦伯数下液滴的初期变形仍受到气流速度、密度等具体流动参数的显著影响。以刚性球体替代液滴进行外流数值模拟,利用球体表面气动力分布推算出的液滴表面变形趋势与实际变形形态吻合,表明液滴的初期变形特征与外流流动分离和涡特征具有一致性。对流场和理论变形数据的分析显示,流动分离发展阶段和稳定阶段对液滴作用力以及它们所诱导的液滴变形特征存在很大差异;分离发展与液滴变形过程的特征时间之比可由气液密度比的平方根表示,它决定了液滴早期变形的基本形态。分离发展阶段所占时间比例越高,即实验中气液密度比越高,则液滴更倾向于发展出单个显著的环形突起,反之则趋于形成多个相对均衡的突起。     

On kinematic,thermodynamic, and kinetic coupling of a damage theory for polycrystalline material

The paper proposes a new consistent formulation of polycrystalline finite-strain elasto-plasticity coupling kinematics and thermodynamics with damage using an extended multiplicative decomposition of the deformation gradient that accounts for temperature effects. The macroscopic deformation gradient comprises four terms: thermal deformation associated with the thermal expansion, the deviatoric plastic deformation attributed to the history of dislocation glide/movement, the volumetric deformation gradient associated with dissipative volume change of the material, and the elastic or recoverable deformation associated with the lattice rotation/stretch. Such a macroscopic decomposition of the deformation gradient is physically motivated by the mechanisms underlying lattice deformation, plastic flow, and evolution of damage in polycrystalline materials. It is shown that prescribing plasticity and damage evolution equations in their physical intermediate configurations leads to physically justified evolution equations in the current configuration. In the past, these equations have been modified in order to represent experimentally observed behavior with regard to damage evolution, whereas in this paper, these modifications appear naturally through mappings by the multiplicative decomposition of the deformation gradient. The prescribed kinematics captures precisely the damage deformation (of any rank) and does not require introducing a fictitious undamaged configuration or mechanically equivalent of the real damaged configuration as used in the past.     

Strain tensor for large three-dimensional surface deformation of sheet metal from an object grating

Grating techniques are used to determine the three-dimensional deformation and the tangential strain of sheet metal. A grating is fixed on the surface and taken by stereo CCD cameras in different deformation states. By suitable line-following software, the grating coordinates in the images are determined with subpixel accuracy. Using photogrammetric methods, the three-dimensional coordinates are calculated from the image coordinates. The strain usually is determined by means of a deformation gradient, which is calculated from every deformed triangle. In this paper, the gradient is determined in the center of four neighboring meshes using a polynomial approximation of the displacement function in a reference position. The influence of the nontangential deformation is considered. By simulation, a flat sheet metal is deformed to a rotational symmetric surface. The difference of the known exact strain is compared with the numerically derived strain with respect to different grating pitches. The proposed method yields good results even in the case of large spatial deformation. It is applied to the deformation of a hatlike test specimen.     

Thermal effect on the deformation of a flexible beam with large kinematical driven overall motions

The research investigates the transient longitudinal and transverse deformation of a planar flexible beam with large overall motions in a temperature field. With the increase of temperature, longitudinal deformation is caused by the thermal expansion in axial direction. Due to the coupling of longitudinal and transverse deformation, the transverse deformation is induced, which is significant in cases where temperature increases rapidly in a very short period. Furthermore, the transverse temperature gradient, which is caused by the temperature variation in the transverse direction, may lead to transverse deformation. Considering the thermal strain, equations of motion of a flexible beam with arbitrary large overall motion are derived based on virtual work principle. The high order terms of the strain tensor are taken into account, such that the geometric nonlinear deformation terms are included in the dynamic equations. Simulation results of a rotating beam are shown to reveal the thermal effect and nonlinear effect on the dynamic performance of the beam.     

考虑曲率纵向变形效应的大变形柔性梁刚柔耦合动力学建模与仿真

对在平面内做大范围转动的中心刚体柔性梁系统的动力学进行了研究,建立了考虑大变形效应的系统刚柔耦合动力学模型,并进行了动力学仿真.该动力学模型不但考虑了柔性梁横向弯曲变形和纵向变形(包含轴向拉伸变形和横向弯曲变形而引起的纵向缩短项),还考虑了纵向变形对曲率的影响,称为曲率纵向变形效应.在以往的研究中,柔性梁的横向弯曲变形能往往直接使用柔性梁横向弯曲变形来表达,并没有考虑纵向变形的影响.为了考虑柔性梁纵向变形对横向弯曲变形能的影响,在浮动坐标系下使用柔性梁参数方程形式的精确曲率公式来计算柔性梁的弯曲变形能.在此基础上建立了基于浮动坐标系的考虑曲率纵向变形效应的刚耦合动力学模型.论文给出了数值仿真算例,验证了本文所建的动力学模型既能适用于柔性梁的小变形问题,又能适用于大变形问题,且较现有高次刚柔耦合动力学模型更加适用于大变形问题的处理.论文还通过与能处理柔性梁大变形问题的绝对节点坐标法的比较,验证了模型的正确性.