多粗糙峰弹塑性接触的有限元分析

采用具有一定数目圆形粗糙峰的刚性表面对弹塑性半元限体进行压下的模型来模拟多粗糙峰接触,并用有限元法对该模型进行了弹生分析,揭示了我地接触区应力状态的影响规律,发现中心接触区的变莆主要受到一定数目邻近粗糙峰的影响,而处于较远自找影响较小,同弹性接触相比,在弹塑性接触过程中有更多的邻近粗糙对中心接触区发生作用,改变粗糙的间距、曲率半径和压下深度都会对其产生影响。     

等截面梁有限变形的传递函数解法

本文应用传递函数方法对等截面梁的有限变形进行了分析,对于等截面梁的有限变形问题,该方法从变分方程出发把问题表述为状态空间的形式,然后利用Gauss积分对轴力进行加速迭代求解,不需要进行增量迭代即可取得具有良好计算精度的数值结果,对简单受力的等截面梁情况该解可以看作是所讨论问题的精确解。对于受力比较复杂或者阶梯变截面梁情况,为减少运算量,可以和有限元法类似,采用多个单元进行拼接,从而得到问题的解。数值算例表明,本方法具有半解析、精度高、收敛快等特点。     

Nonlinear dynamic analysis of circular plates with varying thickness

The BEM is developed for nonlinear free and forced vibrations of circular plates with variable thickness undergoing large deflections. General boundary conditions are considered, which may be also nonlinear. The problem is formulated in terms of displacements. The solution is based on the concept of the analog equation, according to which the two coupled nonlinear differential equations with variable coefficients pertaining to the in-plane radial and transverse deformation are converted to two uncoupled linear ones of a substitute beam with unit axial and unit bending stiffness, respectively, under fictitious quasi-static load distributions. Numerical examples are presented which illustrate the method and demonstrate its accuracy.     

A theory of large-strain isotropic thermoplasticity based on metric transformation tensors

Summary A formulation of isotropic thermoplasticity for arbitrary large elastic and plastic strains is presented. The underlying concept is the introduction of a metric transformation tensor which maps a locally defined six-dimensional plastic metric onto the metric of the current configuration. This mixed-variant tensor field provides a basis for the definition of a local isotropic hyperelastic stress response in the thermoplastic solid. Following this fundamental assumption, we derive a consistent internal variable formulation of thermoplasticity in a Lagrangian as well as a Eulerian geometric setting. On the numerical side, we discuss in detail an objective integration algorithm for the mixed-variant plastic flow rule. The special feature here is a new representation of the stress return and the algorithmic elastoplastic moduli in the eigenvalue space of the Eulerian plastic metric for plane problems. Furthermore, an algorithm for the solution of the coupled problem is formulated based on an operator split of the global field equations of thermoplasticity. The paper concludes with two representative numerical simulations of thermoplastic deformation processes.     

Elastica of a variable-arc-length circular curved beam subjected to an end follower force

This paper presents postbuckling behaviors of a variable-arc-length (VAL) circular curved beam subjected to an end follower force. One end of the VAL circular curved beam is hinged while the other end is supported by a frictionless slot, which is fixed horizontally and vertically but is allowed to rotate corresponding to loading direction. When the VAL circular curved beam is deformed, the total arc-length of the circular curved beam varies. Two approaches have been applied for the solution of this problem. The first approach is an elliptic integrals method based on elastica theory, which yields the exact closed-form solution in terms of the first and second kinds of elliptic integrals. For validation of the results, the shooting method is employed for a numerical solution by developing the set of nonlinear governing differential equations together with boundary conditions, and then integrating them by using the fourth-order Runge–Kutta algorithm. The results from both approaches are in very good agreement. From the results, it is found that the VAL circular curved beam subjected to an end follower force can be deformed in many mode shapes. For the first and third modes, the beam exhibits both stable and unstable configurations, whereas for the second mode only an unstable configuration exists. The influences of initial curvature on the critical load and the deformed configurations are highlighted.     

Research note on evaluation of the plastic deformation of beams under transverse impact

The permanent deformation of a rigid perfectly plastic beam subjected to a time-dependent line load is calculated. The beam has finite length, and the load is applied transversely in the middle of the span. At the initial time, the load is zero. Thereafter, the load has a triangular time dependence. For sufficiently large values of the peak load, static and moving plastic hinges appear, and the beam is subjected to plastic bending. Previous studies of this problem have provided a numerical solution of a system of approximate equations. In this article, a numerical solution of the system of exact equations is provided. The numerical results represent a significant improvement over those that were previously available.     

Efficient coupling of direct forcing immersed boundary-lattice Boltzmann method and finite element method to simulate fluid-structure interactions

Fluid-structure interactions (FSI) of rigid and flexible bodies are simulated in this article. For the fluid flow, multidirect forcing immersed boundary method (IBM) is adopted to capture the moving boundary, and lattice Boltzmann method (LBM) is used to evolve the flow field. Compared with our previous no-penetration IBM, less iterations are required in this work. In addition, larger velocity in lattice units can be used and the nonphysical force oscillations are suppressed due to the C³ 6-point kernel. Multi-relaxation-time collision operator and local grid refinement are also adopted in LBM to enhance the numerical stability and efficiency. To improve the efficiency of the FSI coupling algorithm, the mesh of the deformable structure can be coarser than the Lagrangian mesh using Newton-Cotes formulas to integrate the traction on the structure surface. A variety of benchmarks, including flow around a circular cylinder with Reynold number ranging from 20 to 200, forced oscillation of a circular cylinder, vortex-induced vibration (VIV) of an elastically mounted circular cylinder in two degrees of freedom, and VIV of an elastic cantilever beam attached to a circular cylinder, are carried out to evaluate the accuracy and stability of the present coupling algorithm. For the benchmark FSI problem considered in this article, a reduction of 54% of the calculation time is achieved using coarser structure mesh. As an application of the FSI coupling algorithm, the mechanism of an elastic beam in the wake of a circular cylinder is discussed.     

固支圆板安定性的实验研究

本文利用机动法求解固支圆板在循环集中力作用下的安定问题,得到了安定载荷的理论解;利用激光全息的实验方法研究固支圆板的塑性累积现象,测得了安定载荷的上下限,并对实验结果进行了分析研究.     

Elastoplastic Deformation of Flexible Cylindrical Shells With Two Circular Holes under Axial Tension

The elastoplastic state of thin cylindrical shells with two circular holes under axial tension is analyzed considering finite deflections. The distributions of stresses along the contours of the holes and in the zone of their concentration are studied by solving doubly nonlinear boundary-value problems. The solution obtained is compared with the solutions that account for either physical nonlinearity (plastic deformations) and geometrical nonlinearity (finite deflections) alone and with a numerical solution of the linearly elastic problem. The stress-strain state near the two holes is analyzed depending on the distance between the holes and the nonlinear factors accounted for__________Translated from Prikladnaya Mekhanika, Vol. 41, No. 5, pp. 52–57, May 2005.     

Experimental investigation of circular cantilever beams subjected to impulsive loading

This paper reports the experimental results of quadrantal circular cantilever beam carrying a mass at its tip which is subjected to radial impact. It was observed that after impact a travelling plastic segment moved from the tip towards the root in a very short time, and when specimen reached the largest deformation, the beam whipped out evidently in the reverse direction. The curvature change of the curved beam reaches a peak value near the middle region of the beam and the maximum is at the root. The experimental results are compared with the solutions for the rigid-perfectly plastic model given in [12–13].The authors would like to thank Professor Yang Gui-tong for guidance and help.     

端部受冲击作用的悬臂曲梁的实验研究

本文报导了具有不同端部质量的90°圆弧悬臂梁端部受冲击载荷的实验结果.实验中观察到梁在撞击后程短时间内有塑性区段从端部附近移向根部,梁变形达最大后有明显的反向甩动.变形后梁上各点曲率改变量在梁中部附近有极大值,在根部达最大值.实验结果与理想刚塑性模型的计算结果进行了比较.     

Inelastic Deformation of Flexible Spherical Shells with Two Circular Openings

Elastoplastic analysis of thin-walled spherical shells with two identical circular openings is carried out with allowance for finite deflections. The shells are made of an isotropic homogeneous material and subjected to internal pressure of known intensity. The distributions of stresses (strains or displacements) along the contours of the openings and in the zone of their concentration are studied by solving doubly nonlinear boundary-value problems. The solution obtained is compared with the solutions that account for only physical nonlinearity (plastic deformations) and only geometrical nonlinearity (finite deflections) and with a numerical solution of the linearly elastic problem. The stress–strain state near the two openings is analyzed depending on the distance between the openings and the nonlinear factors accounted for     

自由梁受集中质量两点撞击的刚塑性动力响应

对矩形截面自由梁在两端同时受到完全相同的集中质量横向撞击问题进行了理论上的研究 ,通过采用刚塑性的材料模型得到了其动力响应完全解。结合数值方法给出了梁的瞬态变形 ,并讨论了输入能量、质量比等参数对梁的最终变形、能量耗散的影响。针对典型算例将完全解的结果与MSC/Dytran的计算结果进行了比较 ,两者具有合理的近似 ,但理论预测的结果略高估计了梁的最终变形。     

Meshless methods for the inverse problem related to the determination of elastoplastic properties from the torsional experiment

The problem of determining the elastoplastic properties of a prismatic bar from the given experimental relation between the torsional moment M and the angle of twist per unit length of the rod’s length θ is investigated as an inverse problem. The proposed method to solve the inverse problem is based on the solution of some sequences of the direct problem by applying the Levenberg-Marquardt iteration method. In the direct problem, these properties are known, and the torsional moment is calculated as a function of the angle of twist from the solution of a non-linear boundary value problem. This non-linear problem results from the Saint-Venant displacement assumption, the Ramberg–Osgood constitutive equation, and the deformation theory of plasticity for the stress–strain relation. To solve the direct problem in each iteration step, the Kansa method is used for the circular cross section of the rod, or the method of fundamental solutions (MFS) and the method of particular solutions (MPS) are used for the prismatic cross section of the rod. The non-linear torsion problem in the plastic region is solved using the Picard iteration.     

随机结构有限元分析的递推求解方法

提出了一种新的谱随机有限元分析方法——递推求解方法。该方法将随机结构的随机响应表示成非正交多项式展式,建立了和摄动法类似的一系列确定的递推方程,并通过确定性有限元方法对这些递推方程进行静力问题求解。算例表明,当随机量出现较大涨落时,计算结果相对于传统摄动法有不小的改进。     

多圆孔圆板问题的数值解

提出了一种求解多圆孔圆板问题的新方法。首先引入了基本特解，它由主要部分和附加部分组成。主要部分为带奇点无限平板的一个特殊弹性力学解，奇点取在内圆孔的中心处。附加部分为实心圆饭的一个特解。整个基本特解满足外圆周界为自由条件。文中把待求解取为特解系的形式，其中待定系数可用变分原理得出。最后给出了算例。     

Application of a Special Yield Stress Strain Rate Law to Structural Impulse Problems

ABSTRACT

This paper presents approximate solutions to the dynamic response of three impulsively loaded structures: a wire with an impulsively loaded end mass, an impulsively loaded circular ring, and a cantilever beam with a tip mass subjected to an impulsive load at its tip. The material is assumed to be rigid, perfectly plastic with strain rate sensitivity. A proposed power law form of yield stress strain rate relationship is used to simplify the theoretical development. Numerical solutions are presented for mild steel and are compared with previously published results. Elastic effects and wave propogations are ignored.     

梁中复合应力波的传播

采用有限差分法讨论了梁中复合应力波的传播．给出了粘塑性悬臂梁当自由端受突加弯矩载荷作用时梁内复合应力波传播的基本图象。指出，在冲击早期响应阶段．截面横向转动惯性效应起着重要作用，是不可忽视的。标志弹塑性边界的塑性效，一开始由自由端向固定端运动，但在反射卸载波的迎面作用下，会出现回退现象。在波动早期阶段，固定端主要处于弹性变形状态。此外，还对弹塑性梁中复合应力波的控制方程进行了必要的讨论。     

基于渐近均匀化方法的准周期梁结构等效刚度数值实现方法研究

相比周期梁结构,准周期梁结构沿轴向梯度变化,具有更大的设计自由度,能够获得更好的结构性能。由于其非均质性,一般将其均匀化为具有等效性质的均质梁结构,但现有工作很少涉及准周期梁结构等效性质的计算。本文针对由周期梁结构映射而成的准周期梁结构,通过引入雅可比矩阵,基于渐近均匀化方法推导的单胞方程及其等效性质计算列式,并建立了其单胞方程及等效刚度的有限元求解列式。该方法可以处理沿轴向变形的任意微单胞构型,数值算例验证了其正确性和有效性。