A new quasi-continuum constitutive model for crack growth in an isotropic solid

A new quasi-continuum constitutive model is established based on the randomized cohesive bonds model proposed by Gao and Klein (1998). This model bridges the microscopic discrete constitution characters and the macroscopic mechanical properties of material. In the presented constitutive model, both the bond stretch energy potential and the rotation energy potential are considered, which makes the presented constitutive model applicable to different Poisson-ratio and Young's modulus materials. By establishing a phenomenological bond stiffness function according to the complete stress–strain relationship of uniaxial tension test, the fracture criterion is directly incorporated into the constitutive model. The method requires no external fracture criterion when simulating fracture initiation and propagation, which brings convenience in the numerical simulation. At last, the presented constitutive model is applied to an example of crack growth in an isotropic solid.     

A lattice-based model of the kinetics of twin boundary motion

In this paper we derive a macroscopic kinetic law for twin boundary motion from a lattice dynamical model. The model is developed for compound and type-1 twins and it is explicitly illustrated for a Cu-Al-Ni shape memory alloy. The governing multiple-well energy is calculated using an effective interatomic potential; a Frenkel-Kontorowa type model is developed for the dynamics at the lattice scale; and a quasi-continuum approximation is used to determine the continuum-scale kinetics. The model predicts that compound twins in the Cu-Al-Ni shape memory alloy are an order of magnitude more mobile than type-1 twins which is consistent with experimental observations.     

A homogenization analysis of the field theoretic approach to the quasi-continuum method

Using the orbital-free density functional theory as a model theory, we present an analysis of the field theoretic approach to quasi-continuum method. In particular, by perturbation method and multiple scale analysis, we provide a formal justification for the validity of numerical coarse-graining of various fields in the quasi-continuum reduction of field theories by taking the homogenization limit. Further, we derive the homogenized equations that govern the behavior of electronic fields in regions of smooth deformations. Using Fourier analysis, we determine the far-field solutions for these fields in the presence of local defects, and subsequently estimate cell-size effects in computed defect energies.     

A field theoretical approach to the quasi-continuum method

The quasi-continuum method has provided many insights into the behavior of lattice defects in the past decade. However, recent numerical analysis suggests that the approximations introduced in various formulations of the quasi-continuum method lead to inconsistencies—namely, appearance of ghost forces or residual forces, non-conservative nature of approximate forces, etc.—which affect the numerical accuracy and stability of the method. In this work, we identify the source of these errors to be the incompatibility of using quadrature rules, which is a local notion, on a non-local representation of energy. We eliminate these errors by first reformulating the extended interatomic interactions into a local variational problem that describes the energy of a system via potential fields. We subsequently introduce the quasi-continuum reduction of these potential fields using an adaptive finite-element discretization of the formulation. We demonstrate that the present formulation resolves the inconsistencies present in previous formulations of the quasi-continuum method, and show using numerical examples the remarkable improvement in the accuracy of solutions. Further, this field theoretic formulation of quasi-continuum method makes mathematical analysis of the method more amenable using functional analysis and homogenization theories.     

Fluid–structure interaction modelling of nonlinear aeroelastic structures using the finite element corotational theory

The application of the finite element corotational theory to model geometric nonlinear structures within a fluid–structure interaction procedure is proposed. A dynamic corotational approximately-energy-conserving algorithm is used to solve the nonlinear structural response and it is shown that this algorithm's application with a four-node flat finite element is more stable than the nonlinear implicit Newmark method. This structural dynamic algorithm is coupled with the unsteady vortex-ring method using a staggered technique. These procedures were used to obtain aeroelastic results of a nonlinear plate-type wing subjected to low speed airflow. It is shown that stable and accurate numerical solutions are obtained using the proposed fluid–structure interaction algorithm. Furthermore, it is illustrated that geometric nonlinearities lead to limit cycle oscillations.     

Amplitude dependence of filler-reinforced rubber: Experiments,constitutive modelling and FEM – Implementation

The focus of the present paper is the experimental investigation, the constitutive representation and the numerical simulation of the amplitude dependence of filler-reinforced elastomers. A standard way to investigate the dynamic properties of viscoelastic materials is via the dynamic modulus which is obtained from stress signals due to harmonic strain excitations. Based on comprehensive experimental data, an amplitude-dependent constitutive model of finite viscoelasticity is developed. The model is based on a modified Maxwell chain with process-dependent viscosities which depend on additional internal state variables. The evaluation of this thermodynamically consistent model is possible in both the time domain, via stress-time signals, and in the frequency domain, via the dynamic modulus. This property is very profitable for the parameter identification process. The implementation of the constitutive model into the commercial finite element code ANSYS with the user-programmable feature (UPF) USERMAT for large deformations in updated Lagrange formulation is presented. This implementation allows simulating the time-dependent behaviour of rubber components under arbitrary transient loading histories. Due to physical and geometrical nonlinearities, these simulations are not possible in the frequency domain. But, transient FEM computations of large loading histories are sometimes not possible in an acceptable time. In the context of the parameter identification the fundamental ideas are presented, how this problem has been solved. Transient FEM simulations of real rubber components are also shown to visualize the properties of the model in the context of the transient material behaviour.     

SIZE EFFECTS OF ELASTIC MODULUS OF FCC METALS BASED ON THE CAUCHY-BORN RULE AND NANOPLATE MODELS

In the present research, a simple quasi-continuum model, the Cauchy-Born rule model, is used to investigate the size effects of elastic modulus for fcc metals. By considering a nanoplate model and calculating the strain energy for the nano-sized plate under tension and bending, the relationship between the elastic modulus and the plate thickness is found. Size effects of the elastic modulus are displayed by the relative differences of the elastic modulus between the nano-sized plate sample and the bulk sample. By comparing the present results with those of others, the effectiveness of the Cauchy-Born rule model in studying the size effects of material properties are shown.     

Propagation,Dispersion, and Localization of Elastic Waves in Low-Symmetric Anisotropic Waveguides

The results of a theoretical investigation into ultrasonic linear normal wavefields in anisotropic three-dimensional bodies with mechanical orthorhombic symmetry are presented. A detailed study is made of the dispersion dependencies for higher normal wave modes in a single-layered orthorhombic plate-type waveguide. Moreover, the distribution of complex branches corresponding to the edge standing wave modes is studied and their role in the transformation of the entire spectrum upon a change in the travel direction in the waveguide plane is analyzed. A new type of localization of the higher modes of high-frequency short-range normal waves in a crystal layer is described. An efficient method is developed for studying the spectrum of ultrasonic normal waves in a circular cylindrical waveguide made of an orthorhombic monocrystal     

Changes in the natural frequency of a ferromagnetic rod in a magnetic field due to magnetoelastic interaction

Based on the magnetoelastic generalized variational principle and Hamilton＇s principle, a dynamic theoretical model characterizing the magnetoelastic interaction of a soft ferromagnetic medium in an applied magnetic field is developed in this paper. From the variational manipulation of magnetic scale potential and elastic displacement, all the fundamental equations for the magnetic field and mechanical deformation, as well as the magnetic body force and magnetic traction for describing magnetoelastic interaction are derived. The theoretical model is applied to a ferromagnetic rod vibrating in an applied magnetic field using a perturbation technique and the Galerkin method. The results show that the magnetic field will change the natural frequencies of the ferromagnetic rod by causing a decrease with the bending motion along the applied magnetic field where the magnetoelastic buckling will take place, and by causing an increase when the bending motion of the rod is perpendicular to the field. The prediction by the mode presented in this paper qualitatively agrees with the natural frequency changes of the ferromagnetic rod observed in the experiment.     

微重力下成一定夹角平板间的表面张力驱动流动的研究

空间微重力环境中,由于重力基本消失,表面张力等次级力发挥主要作用,流体行为与地面迥异,因此有必要深入探究微重力环境中的流体行为规律和特征.板式贮箱利用板式组件在微重力环境中对流体进行管理,从而为推力器提供不夹气的推进剂,这对航天器精确进行姿态控制、轨道调整具有重要意义.板式组件中常包含成一定夹角的平板结构,比如蓄液叶片...     

Characterization of piezo-film sensors for direct vibration and impact measurements

Piezo-film sensors were employed in determining the dynamic response of [(0 deg/90 deg)₄]_s s-glass/epoxy laminates and 2024 aluminum specimens. Simple beam- and plate-type sensor equations were derived based on classical plate theory incorporating the linear piezoelastic constitutive relationship. A series of vibration and impact tests were conducted for the determination of structural dynamic response. Piezo-film sensors, with a thickness and area of 110 μm and 1×1 cm², respectively, were connected directly to a voltage measurement device in these tests. The first three bending frequencies of the glass-fiber-reinforced plastic (GFRP) cantilever specimen were examined. Experimental results and those simulated by the MARC finite-element code were found to be in good agreement, with the difference between the two being less than five percent. At frequencies above 3 kHz, piezo-film transducers are capable of closely detecting structural dynamic response in the absence of charge amplification. At frequencies lower than approximately 3 kHz, however, the voltage measured directly from a piezo-film sensor underestimates structural response. A modified piezo-film sensor equation is thus proposed for lower frequency measurements. Effect of frequency and piezo-film's size on lower frequency attenuation is explicitly formulated based on a simple RC circuit analogy. Drop tests were also performed on clamped [(0 deg/90 deg)₄]_s GFRP laminates and aluminum targets, with nine piezo-film sensors being glued to the specimen's distal surface in order to determine the low-speed impact response. Specimen transient central deflection subjected to impact loading was identified based on test findings and the plate-type piezo-film sensor equation. The results were found to be in good agreement with the numerical solution obtained from the MARC finite-element code.     

Dynamic Systems Analysis and Control Based on a Configuration Manifold Model

This paper presents a configuration manifold model for the analysis of dynamic systems and the development of control algorithms from both geometrical and topological points of view. The fundamental theory of surfaces and differential manifolds endowed with Riemannian metrics is overviewed. The concepts of configuration manifolds (C-manifolds) and their immersions and embeddings are then introduced and applied to dynamic systems modeling. An explicit form of the smooth embedding for a given dynamic system with its C-manifold is derived. In an open serial-chain robotic system, a topological equivalence, i.e. a homeomorphism, is found and shown to be useful for dynamic model reduction. With topology being viewed as the structure of geometry, we discover and prove that the kinematics of a dynamic system determines its topology so that the kinematics is virtually a structure of the system's dynamics. This key point of view is further extended to the development of an adaptive control strategy. A computer simulation study is finally performed to verify the proposed model and adaptive control scheme.     

Finite and impulsive motion of constrained mechanical systems via Jourdain's principle: discrete and hybrid parameter models

The main purpose of this paper is to present a unified analytical dynamics framework for the analysis of finite and impulsive motion of mechanical systems using Jourdain's principle. Emphasis is given to the general case when a mechanical system is described by a hybrid (discrete-distributed) parameter model. A large group of finite and impulsive, generally non-holonomic, constraints are analysed in detail and a so-called extended Appellian classification is presented for these constrained motion problems. The fundamental dynamic equation of constrained systems is developed in terms of velocity variations (Jourdain's principle). Based on this equation and the constraints, the methods of quasivelocities and Lagrangian multipliers are adopted and interpreted for the finite motion of hybrid parameter models of mechanical systems; and the methods of independent quasivelocity variations and Lagrangian multipliers are introduced for the analysis of impulsive motion of such models. To illustrate the proposed material, an example of a one-link flexible arm intercepting and capturing a moving target is considered.     

纳米力学进展

概述在固态下纳米力学的若干研究内容.首先对纳米力学及其范畴 进行界定，然后介绍纳米力学方法，包括属于纳观计算力学范畴的大规 模分子动力学算法、连续介质/分子动力学交叠层算法、准连续介质算法 和LMPM方法；及属于纳观实验力学范畴的纳米云纹法和纳米压痕法.随即 阐述纳米力学的新兴研究领域：包括纳米晶体的超塑性变形、纳观断裂 力学、纳米管力学和纳米压痕力学.     

A meshfree method for transverse vibrations of anisotropic plates

A meshfree approach, called displacement boundary method, for anisotropic Kirchhoff plate dynamic analysis is presented. This method is deduced from a variational principle, which uses a modified hybrid functional involving the generalized displacements and generalized tractions on the boundary and the lateral deflection in the domain as independent variables. The discretization process is based on the employment of the fundamental solutions of the static problem operator for the expression of the variables involved in the functional. The stiffness and mass matrices obtained for the dynamic model are frequency-independent, symmetric and positive definite and their computation involves boundary integrals of regular kernels only. Due to its features, the final resolving system can be solved with the classical approaches by using standard numerical procedures. To assess the formulation, the free vibrations of some anisotropic plates were calculated and the results compared with those obtained using other solution techniques. The present results are in good agreement with those found in the literature showing the accuracy and effectiveness of the proposed approach.     

Dynamics and stability of turbocharger rotors

The paper discusses the bifurcation and stability behavior of (automotive) turbochargers with full-floating ring bearings. Turbocharger rotors exhibit a highly nonlinear behavior due to the nonlinearities introduced by the floating ring bearings. A flexible multibody model of the rotor/bearing system is presented. Numerical run-up simulations are compared with corresponding test rig measurements. The nonlinear oscillation effects are thoroughly investigated by means of simulated and measured rotor vibrations. The influence of various system parameters on the bifurcation behavior of the rotor/bearing system is analyzed. The article examines rotors supported in full-floating ring bearings with plain circular bearing geometry in the inner and outer oil gap. By recapitulating the well-known oil whirl and oil whip phenomena for single and double oil film bearings, the paper gives an overview on the fundamental dynamic effects occurring in turbocharger systems.     

A nonlinear threshold model applied to spallation analysis

Spallation in heterogeneous media is a complex, dynamic process. Generally speaking, the spallation process is relevant to multiple scales and the diversity and coupling of physics at different scales present two fundamental difficulties for spallation modeling and simulation. More importantly, these difficulties can be greatly enhanced by the disordered heterogeneity on multi-scales. In this paper, a driven nonlinear threshold model for damage evolution in heterogeneous materials is presented and a trans-scale formulation of damage evolution is obtained. The damage evolution in spallation is analyzed with the formulation. Scaling of the formulation reveals that some dimensionless numbers govern the whole process of deformation and damage evolution. The effects of heterogeneity in terms of Weibull modulus on damage evolution in spallation process are also investigated.     

多尺度模拟与计算研究进展

简要介绍了多尺问题与研究方法.重点论述了两类常见多尺度问题的模拟计算方法与研究进展,分析了各自的优缺点和使用范围.对现有研究的局限性和存在的问题进行分析,指出了进一步研究多尺度模拟与计算的必要性.介绍了求解含有孤立缺陷问题的非局部准连续体法,MAAD等方法以及求解基于微观模型本构模拟问题的局部连续体法、HMM等方法.文章对多尺度模拟与计算的前景进行展望,提出了一些亟待解决的问题.     

Qualitative criteria in non-linear dynamic buckling and stability of autonomous dissipative systems

A general qualitative approach for dynamic buckling and stability of autonomous dissipative structural systems is comprehensively presented. Attention is focused on systems which under the same statically applied loading exhibit a limit point instability or an unstable branching point instability with a non-linear fundamental path. Using the total energy equation, the theory of point and periodic attractors of the basin of attraction of a stable equilibrium point, of local and global bifurcations, of the inset and outset manifolds of a saddle and of the geometry of the channel of motion, the stability of the fundamental equilibrium path and the mechanism of dynamic buckling are thoroughly discussed. This allows us to establish useful qualitative criteria leading to exact, approximate and upper/lower bound buckling estimates without integrating the highly non-linear initial-value problem. The individual and coupling effect of geometric and material non-linearities of damping and mass distribution on the dynamic buckling load are also examined. A comparison of the results of the above qualitative analysis with those obtained via numerical simulation is performed on several two- and three-degree-of-freedom models of engineering importance.