具有损伤耦合效应的弹塑性蠕变问题结构分析的变分原理

本文基于连续介质损伤力学中有效应力的概念,研究弹塑性蠕变问题中的损伤耦合,并应用由最优控制理论基本思想发展起来的参变量变分原理建立起用于弹塑性蠕变损伤问题结构分析的变分原理.文中给出了原理的证明.该原理的物理意义明确,表达式简单且规范,容易为数值手段实现.     

参变量变分原理的提出、发展与应用

参变量变分原理及其参数二次规划算法是由钟万勰院士1985年针对弹性接触边界非线性问题首次提出来的,经过将近40年的不断发展,目前参变量变分原理已经成功应用于各个领域,其中包括弹塑性分析、接触问题、润滑力学、岩土力学、变刚度杆系结构、先进材料性能分析、材料的蠕变与损伤、柔性结构力学和LQ最优控制等各个工程领域。本文首先回顾了参变量变分原理的起源,介绍了参变量变分原理的基本概念,然后以弹塑性分析问题为例,阐明建立参变量变分原理的理论模型以及实现数值参数二次规划求解原理,最后详细回顾了参变量变分原理的基本理论与相应数值算法在各个领域的发展及其工程应用,展示了参变量变分原理在求解各类非线性问题的特色与优势。     

考虑损伤累积的热弹塑性问题变分原理及其有限元方法

本文基于连续介质损伤力学中有效应力的概念，研究建筑结构抗火分析中遇到的热弹塑性问题与损伤的耦合，并运用参变量变分原理，建立起用于热弹塑性损伤问题结构分析的变分原理。本文给出了原理应用的有限元列式。具有明确的物理意义，表达形式规范，便于数值手段实现。     

弹性接触问题的变分原理及参数二次规划求解

本文给出了平面与空间接触问题的带参变量的变分极值原理。接触是考虑库伦摩擦的。参变量二次规划可以用于精确求解,并且通过数例说明了计算方法。     

理想粘塑性流体润滑问题的参变量变分原理

本文给出了理想粘塑性流体润滑问题的参变量变分原理.在膜厚方向压力为常数的假设下,塑性剪切滑移面将发生在固液交界面上,因而可以选择边界速度滑移量为参变量(控制变量).文中讨论了采用有限元求解时的实施过程,原问题最后可化为求解带约束条件的参数二次规划问题.该方法简单可靠,具有良好的工程应用前景.     

一类铁电模型的参变量变分原理及其应用

将参变量变分原理引入铁电问题。对一类借用了经典弹塑性理论中的概念和方法的多轴铁电模型建立基于Helmholtz自由能的参变量变分原理,可以有效处理传统变分原理中由非关联流动法则或屈服面不考虑材料系数变化所引起的切线模量非对称困难。相应于参变量变分原理,引入参数二次规划算法,可获得具有可靠数值稳定性的一套铁电算法。将该算法应用于一个具体的铁电模型,数值计算结果表明本文方法的有效性。     

岩土力学中的参变量最小余能原理

本文对于岩土材料的流动法则与塑性势不相关连而引起的Drucker公设不成立,因而塑性力学中的经典的变分原理也不再适用的问题,提出了参变量余能变分原理。如同弹性接触问题及极限分析中那样,该原理可望用于有限元数值分析中。     

岩土中弹塑性渗透固结问题的参变量变分原理

本文给出渗透固结过程的参变量变分原理,可以用来处理岩土弹塑性渗透固结问题,并不受 Drucker 假设的限制,对于关连或非关连塑性流动情形皆可。     

弹塑性岩土介质中渗透固结问题的参数二次规划法

岩土类介质在承受荷载时,不仅产生弹塑性变形,还伴随着渗流和固结,是一个关于时间的动态过程。本文建立起处理该问题的参变量变分原理以及相应的有限元方法,这样将原问题化为求解带约束条件(本构状态方程)的二次规划问题.文中讨论了单元的选取形式及具体的实施过程,还给出了一个实例.     

解三维摩擦接触问题的一个二阶锥线性互补法

针对三维摩擦接触问题的求解,给出了一种基于参变量变分原理的二阶锥线性互补法. 首先,基于三维Coulomb摩擦锥在数学表述上属于二阶锥的事实,利用二阶锥规划对偶理论,建立了三维Coulomb摩擦接触条件的参变量二阶锥线性互补模型,它是二维Coulomb摩擦接触条件参变量线性互补模型在三维情形下的自然推广;随后,利用参变量变分原理与有限元方法,建立了求解三维摩擦接触问题的二阶锥线性互补法. 较之于将三维Coulomb摩擦锥进行显式线性化的线性互补法,该方法无需对三维Coulomb摩擦锥进行线性化,因而在保证精度的前提下所解问题的规模要小很多. 最后通过算例展示了该方法的特点.      

THE PARAMETRIC VARATIONAL PRINCIPLE AND NON-LINEAR FINITE ELEMENT METHOD FOR ANALYSIS OF ASTROMESH ANTENNA STRUCTURES

针对大型周边桁架式索网天线由拉索拉压模量不同引起的本构非线性和结构大变形引起的几何非线性问题,给出了基于参变量变分原理的几何非线性有限元方法. 首先针对含预应力索单元拉压模量不同分段描述的本构关系,通过引入参变量,导出了基于参变量及其互补方程的统一描述形式,避免了传统算法需要根据当前变形对索单元张紧/松弛状态的预测,提高了算法收敛性. 然后利用拉格朗日应变描述索网天线结构大变形问题,结合几何非线性有限元法,建立了基于参变量的非线性平衡方程和线性互补方程;并给出了牛顿-拉斐逊迭代法与莱姆算法相结合的求解算法. 数值算例验证了本文提出的算法比传统算法具有更稳定的收敛性和更高的求解精度,特别适合于大型索网天线结构的高精度变形分析和预测.     

索网天线的参变量变分及非线性有限元方法

针对大型周边桁架式索网天线由拉索拉压模量不同引起的本构非线性和结构大变形引起的几何非线性问题,给出了基于参变量变分原理的几何非线性有限元方法. 首先针对含预应力索单元拉压模量不同分段描述的本构关系,通过引入参变量,导出了基于参变量及其互补方程的统一描述形式,避免了传统算法需要根据当前变形对索单元张紧/松弛状态的预测,提高了算法收敛性. 然后利用拉格朗日应变描述索网天线结构大变形问题,结合几何非线性有限元法,建立了基于参变量的非线性平衡方程和线性互补方程;并给出了牛顿-拉斐逊迭代法与莱姆算法相结合的求解算法. 数值算例验证了本文提出的算法比传统算法具有更稳定的收敛性和更高的求解精度,特别适合于大型索网天线结构的高精度变形分析和预测.      

分段线性强化问题的全量无迭代解

本文应用[1]提出的塑性全量理论中的控制变量(也称参变量)变分原理,对各级分段线性强化规律建立了相应的本构状态方程以及有限元求解公式,可使分段线性强化全量问题的数值解不需迭代.对于一般的(?)—(?)曲线,适当选择分段数目可达到足够的精度。本文给出了有一定代表性的手工及数值算例。     

纤维增强复合材料层压板的失效过程分析——变参量变分方法

本文建立了关于具有逐层失效不连续本构关系的纤维增强复合材料层压板的参变量变分原理,所构造的受参变量控制的能量泛函的驻值问题可转化为一个拟二次规划问题.本文给出了一个求解此规则问题的有效算法.本文方法的优点在于避免了分析层压板强度的冗长迭代过程,可在一个加载步内求出对应各单层失效的一系列强度比值.     

压气机超速制造时弹塑性接触的数值模拟

“叶轮超速工艺“是叶轮制造中的一种重要方法.以某柴油机涡轮增压器的压气机为例,采用有限元参数二次规划法结合多重子结构技术分析模拟超速预加载制造时叶轮、轴套与轴之间接触应力的相应分布规律.针对不同的超速转速进行了大量计算,获得了超速转速与接触应力;超速转速与压气机内孔最大残余变形量;超速转速与轴套壁厚之间的关系.结果表明本文所采用的分析方法是可行的,弹塑性接触分析结果可作为叶轮设计及加工工艺的参考.     

A second strain gradient elasticity theory with second velocity gradient inertia – Part I: Constitutive equations and quasi-static behavior

A multi-cell homogenization procedure with four geometrically different groups of cell elements (respectively for the bulk, the boundary surface, the edge lines and the corner points of a body) is envisioned, which is able not only to extract the effective constitutive properties of a material, but also to assess the “surface effects” produced by the boundary surface on the near bulk material. Applied to an unbounded material in combination with the thermodynamics energy balance principles, this procedure leads to an equivalent continuum constitutively characterized by (ordinary, double and triple) generalized stresses and momenta. Also, applying this procedure to a (finite) body suitably modelled as a simple material cell system, in association with the principle of the virtual power (PVP) for quasi-static actions, an equivalent structural system is derived, featured by a (macro-scale) PVP having the typical format as for a second strain gradient material model. Due to the surface effects, the latter model does work as a combination of two subsystems, i.e. the bulk material behaving as a Cauchy continuum, and the boundary surface operating as a membrane-like boundary layer, each subsystem being in (local and global) equilibrium by its own. Further, the applied (ordinary) boundary traction splits into two (response-dependent) parts, i.e. the “Cauchy traction” transmitted to the bulk material and the “Gurtin–Murdoch traction” acting, together with all other boundary tractions, upon the boundary layer. The role of the boundary layer as a two-dimensional manifold enclosing a Cauchy continuum is elucidated, also with the aid of a discrete model. A strain gradient elasticity theory is proposed which includes a minimum total potential energy principle featuring the relevant boundary-value problem for quasi-static loads and its (unique) solution. A simple application is presented. Two appendices are included, one reports the proof of the global equilibrium of the boundary layer, the other is concerned with double and triple stresses. The paper is complemented by a companion Part II one on dynamics. Previous findings by the author [Polizzotto, C., 2012. A gradient elasticity theory for second-grade materials and higher order inertia. Int. J. Solids Struct. 49, 2121–2137] are improved and extended.     

聚合物粘结剂PVP、PAM和PVA结构与性能的分子动力学研究

在粉末打印骨支架的工艺中,粘结剂的性质是影响骨支架质量的关键因素。采用分子动力学的模拟方法对三种常用高分子粘结剂的体系进行了构建和模拟,从微观分子层面研究了聚合物粘结剂PVP,PAM和PVA的部分性质,比如密度、内聚能及力学性能,并对三种粘结剂的性能参数进行了比较,对其内在关系进行了揭示。此外,还通过建立粘结剂与羟基磷灰石的界面相互作用模型,对三种粘结剂与羟基磷灰石的界面结合能进行了计算和比较,分析了影响高聚物粘结特性的根本原因。这一工作不仅对常用粘结剂的基本性质进行了预估,而且对骨支架粉末粘结工艺中粘结剂的选择提供了理论依据。     

Rigid finite element and limit analysis

According to the lower bound theorem of limit analysis the Rigid Finite Element Method (RFEM) is applied to structural limit analysis and the linear programmings for limit analysis are deduced in this paper. Moreover, the Thermo-Parameter Method (TPM) and Parametric Variational principles (PVP) are used to reduce the computational effort while maintaining the accuracy of solutions. A better solution is also obtained in this paper. The project supported by National Natural Science Foundation of China     

Controlled preparation and characterization of nano-sized hexagonal Mg（OH）2 flame retardant

Nano-sized hexagonal magnesium hydroxide （Mg（OH）2） with good dispersibility was synthesized by a double injection-hydrothermal method, utilizing polyvinylpyrrolidone （PVP） as an additive and with optimized processing parameters. SEM and BET analysis showed that the mean particle size and specific surface area of the Mg（OH）2 particles were 174 nm and 50.77 m^2/g, respectively. The FT-IR spectra and the XRD patterns showed that PVP was adsorbed on the surface of the Mg（OH）2 crystal, thus effectively limiting particle agglomeration and hindering crystal growth along the （1 01 ） plane. TGA showed a decrease in the decomposition temperature and an increase in the weight loss of the Mg（OH）2 particles due to addition of PV/.     

Application of solvent evaporation technique for pure drug crystal spheres preparation

An innovative application of the solvent evaporation technique was suggested. Solvent evaporation technique is a technique for drug encapsulation and nanosphere preparation. The widely used technique is also facing the problem of low actual drug entrapment percent, which is not economic from the industrial view. The goal of this work is trying to use the advantage of this technique concerning the product sphericity and the ability to control particle size, to prepare a drug as pure crystals spheres. Ibuprofen is selected as a model drug. The spheres are formed by using Polyvinyl pyrrolidone (PVP) or Polyethylene glycol (PEG) as an anti-aggregating agent but not formed on using tween or span. Particle size and actual drug content depend on the concentrations the anti-aggregating agent used. Surfaces of the drug crystal spheres are porous with empty sphere internal structure on using PVP but spongy and rough on using PEG. The drug has its identity chemical form in the drug crystal spheres. IR scan of spheres prepared on using PEG showed a characteristic ether peak. DSC showed melting endothermic peak of PEG, but X-ray showed minor change in the drug crystal patterns. Drug release profiles from crystal spheres prepared with the same anti-aggregating agent are close to each other. The drug release profiles from drug crystal spheres prepared by using PEG are more controlled than that prepared by using PVP. The drug release mechanism is diffusion. It was concluded that, the same technique could be suggested for preparation of other biomedical material in pure crystals spheres with controlled particle size. These properties may encourage to prepare very small particles with spherical shape for inhalation or injection as an innovative particle technology application for the widely used technique.