高性能混凝土框架结构滞回全过程分析

基于非线性有限元原理，编制了高性能混凝土梁基于本构关系的滞回全过程分析程序，程序设计中考虑了材料非线性、几何非线性、预应力作用、混凝土裂面效应、轴力二次矩、材料双切线模量场等因素的影响，并考虑了高性能混凝土的特点，程序计算值与试验结果吻合较好。在此基础上，编制了高性能混凝土框架结构基于恢复力模型的滞回全过程分析程序，对该程序进行了试验验证，计算表明高性能混凝土框架结构具有良好的抗震性能。     

预应力混凝土桥梁结构非线性仿真研究

针对现代预应力桥梁结构受力复杂、非线性效应明显的问题,从理论上对其非线性仿真方法进行了探讨。通过在实体等参元中引入梁和板壳等理论,提出了用于复杂混凝土桥梁结构仿真建模的实体退化系列单元;从预应力筋的实际作用机理出发,提出了预应力效应分析的二阶段法;建立了混凝土徐变的三参数粘弹性模型,并推导了三维徐变效应分析的递推公式;从几何非线性、材料非线性及求解方法等三个方面对混凝土开裂及极限承载能力分析的具体实现方式进行了讨论。用三个算例分别对预应力、徐变及极限承载状态等的分析方法的有效性进行了验证。结果表明,本文方法能更真实有效地模拟实际结构,客观分析预应力及徐变的三维效应,可以准确预测结构的开裂及极限承载状态,具有很强的适用性。     

预应力等效荷载的应用及非线性阶段计算方法

等效荷载概念在预应力结构的研究和设计中存在广泛应用．依据等效荷载的经典定义,对线性变换后预应力等效荷载是否变化进行论述,指出采用等效荷载法计算预应力混凝土受弯构件在使用阶段的预加力反拱值时,应取用净截面刚度．针对工程界提出的一些观点,用解析方法证明直线、抛物线组合线型力筋的预应力作用仍可用等效荷载法计算,以及等效荷载法并不存在与预应力筋线型相关的局限性．为回答无粘结、有粘结预应力筋应力增长后等效荷载是否变化的问题,提出实时等效荷载的概念,建立了预应力混凝土梁实时等效荷载的计算方法．推导基于实时等效荷载的承载力极限状态平衡方程,并证明其与将预应力筋完全作为抗力材料的平衡方程具有等价性．最后,介绍了实时等效荷载运用于预应力结构非线性阶段变形计算时所体现的优越性．     

混凝土率型内时损伤本构模型

混凝土是一种典型的率敏感材料,为了更好地描述混凝土结构在动力、冲击荷载作用下的强度和变形特征,本文结合内时理论和损伤理论建立了一种考虑混凝土率效应的内时损伤本构模型。该模型的特点:将混凝土材料的受力软化效应分解为密实状态的塑性效应和由微裂缝扩展引起的刚度退化效应。前者由内时理论来描述,这使该模型摆脱了一般弹塑性模型中屈服面的概念,从而更符合混凝土的变形特性,并且简化了非线性计算过程;后者由损伤理论来描述,根据混凝土的动力试验结果建立了增量型的损伤演变方程,从而使该模型能够较好地反映混凝土的动力特性。最后,应用本文建议的模型对一钢筋混凝土简支梁进行了非线性分析,结果表明:当结构承受快速荷载作用时,应变率对结构的受力性能影响较大,在进行结构分析时必须予以考虑。     

一个改进的混凝土非线性双剪切强度准则

通过对混凝土、岩石等材料强度理论的分析研究,基于普遍形式强度理论及双剪切强度理论,考虑混凝土在多轴荷载作用下材料强度性质在π平面和拉伸压缩子午线上不同的变化规律,将混凝土的强度准则定义为两个较大主剪切应力平方和的函数;并考虑其对应剪切应力作用面上的正应力和静水压力的影响,推导了一个考虑中间主应力效应和复杂极限面情况的非线性双剪切强度准则;通过实验数据获得了相关参数。结果表明该强度准则可以较好地描述混凝土材料的破坏特性。     

桥梁空间分析中预应力效应分析方法研究

预应力技术在桥梁中的广泛使用,使得桥梁结构受力状况复杂,对预应力混凝土桥梁进行空间分析是必要的;但由于预应力筋数量庞大且形状复杂,对预应力空间效应的模拟是预应力混凝土桥梁空间分析的难点。从预应力筋的实际作用机理出发,提出在预应力张拉阶段用等效节点荷载来模拟预应力对结构的作用;而在预应力管道灌浆后,将预应力筋作为结构的一部分,用实体退化单元进行结构分析,并用两算例对本文的方法进行了验证。结果表明,该方法能真正模拟预应力的空间效应,可用于预应力混凝土桥的空间分析。     

预应力混凝土平面杆系结构的有限元方法

建立了基于有限元方法的考虑材料和几何非线性的任意截面预应力混凝土平面杆系结构的数值分析模型,可用于模拟预应力混凝土大跨度梁、单向偏压细长柱等的非线性全过程结构响应。引入修正的Rodriguez截面模型确定截面切线刚度,其中混凝土的贡献通过截面边界顶点定义的梯形单元来实现;在此基础上利用传统的平面非线性杆单元导出了标准有限元公式。通过两个算例验证了该模型的可靠性和适用性。     

模拟酸雨腐蚀下预应力混凝土梁的抗弯性能研究

为了探讨酸雨腐蚀环境对预应力混凝土梁抗弯性能的影响,制作了8根预应力混凝土梁,从混凝土强度等级、预应力度和钢绞线腐蚀率等方面,进行8根预应力混凝土梁的三分点静载试验。试验结果及分析表明,混凝土强度对屈服后试验梁的抗弯性能影响较大,对弹性阶段的刚度影响较小;随着预应力度的增加,试验梁的开裂荷载逐渐增大,进入裂缝阶段后,刚度下降速度逐渐加快;钢绞线腐蚀率越大,试验梁的极限抗弯承载力越低;腐蚀率在一定范围内(3.22%),随着腐蚀率的增大,预应力混凝土梁试件的开裂弯矩、屈服弯矩和极限弯矩明显降低。根据试验结果,提出了模拟酸雨环境下预应力混凝土梁试件抗弯承载力的计算方法,同时采用ANSYS软件对预应力混凝土梁进行非线性有限元分析,将试验梁的计算值和模拟值与试验值进行比较,均吻合较好。     

模拟酸雨腐蚀下预应力混凝土梁的抗弯性能研究

为了探讨酸雨腐蚀环境对预应力混凝土梁抗弯性能的影响，制作了8根预应力混凝土梁，从混凝土强度等级、预应力度和钢绞线腐蚀率等方面，进行8根预应力混凝土梁的三分点静载试验。试验结果及分析表明，混凝土强度对屈服后试验梁的抗弯性能影响较大，对弹性阶段的刚度影响较小；随着预应力度的增加，试验梁的开裂荷载逐渐增大，进入裂缝阶段后，刚度下降速度逐渐加快；钢绞线腐蚀率越大，试验梁的极限抗弯承载力越低；腐蚀率在一定范围内（>3.22%），随着腐蚀率的增大，预应力混凝土梁试件的开裂弯矩、屈服弯矩和极限弯矩明显降低。根据试验结果，提出了模拟酸雨环境下预应力混凝土梁试件抗弯承载力的计算方法，同时采用ANSYS软件对预应力混凝土梁进行非线性有限元分析，将试验梁的计算值和模拟值与试验值进行比较，均吻合较好。     

部分预应力混凝土梁抗冲击性能试验研究

通过在跨中施加一落锤冲击荷载的方法,研究了3根无粘结部分预应力混凝土梁的抗冲击性能,给出了典型的实测波形曲线,测得了无粘结部分预应力混凝土试验梁在落锤冲击荷载作用下的动态位移延性和极限承载能力,分析了落锤重力势能对预应力混凝土梁动力性能的影响,并根据相似模拟理论,得出了与试验梁相似的12 m跨度的实际无粘结部分预应力混凝土梁的极限抗冲击承载能力。肯定了无粘结部分预应力混凝土试验梁具有很好的抗冲击性能,不会发生脆性破坏,完全可用于大跨度地下结构。     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

Effet de l'anisotropie élastique cristalline sur la distribution des facteurs de Schmid à la surface des polycristauxInfluence of crystalline elasticity anisotropy on Schmid factor distribution at the free surface of polycrystals

Experimental studies of the plasticity mechanisms of polycrystals are usually based on the Schmid factor distribution supposing crystalline elasticity isotropy. A numerical evaluation of the effect of crystalline elasticity anisotropy on the apparent Schmid factor distribution at the free surface of polycrystals is presented. Cubic elasticity is considered. Order II stresses (averaged on all grains with the same crystallographic orientation) as well as variations between averages computed on grains with the same crystallographic orientation but with different neighbour grains are computed. The Finite Element Method is used. Commonly studied metals presenting an increasing anisotropy degree are considered (aluminium, nickel, austenite, copper). Concerning order II stresses in strongly anisotropic metals, the apparent Schmid factor distribution is drifted towards small Schmid factor values (the maximum Schmid factor is equal to 0.43 instead of 0.5) and the slip activation order between characteristic orientations of the crystallographic standard triangle is modified. The computed square deviations of the stresses averaged on grains with the same crystallographic orientation but with different neighbour grains are a bit higher than the second order ones (inter-orientation scatter). Our numerical evaluations agree quantitatively with several observations and measures of the literature concerning stress and strain distribution in copper and austenite polycrystals submitted to low amplitude loadings. Hopefully, the given apparent Schmid factor distributions could help to better understand the observations of the plasticity mechanisms taking place at the free surface of polycrystals. To cite this article: M. Sauzay, C. R. Mecanique 334 (2006).     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).     

Dynamics and Breakup of Viscoelastic Liquids (A Review)

The phenomena of hydrodynamic breakup of liquid jets, drops, films, bridges, and filaments are reviewed for liquids with viscoelastic properties. The reasons for breakup are capillary instabilities, collisions with rigid obstacles, and other forms of dynamic action. The relationship between the properties of the liquids and the features of the breakup process is discussed.     

Flow stability analysis and excitation using pulsating jets

Classical flow stability applied to transition from laminar to turbulent flow may also describe the behavior of vorticity fluctuations created by a pulsating jet placed along a solid boundary. A numerical laminar flow experiment involving a pulsating jet placed along the surface of a duct with flow separation downstream, resulted in eliminating most part of the separated flow region. Applying the same approach to a turbulent flow, it was possible to develop a turbulent stability flow formulation and apply successfully turbulent pulsating jet flow separation control. To cite this article: D. Skamnakis, K. Papailiou, C. R. Mecanique 333 (2005).