Authors:Friedrich Pfeiffer; Johannes Mayet Pages: 1411 - 1426 Abstract: Abstract Chain fountains are known since long time, and many efforts have been taken to model and to explain the dynamics of such a chain fountain. A chain consists of many small elements starting from an inertial container, forms this specific arc and comes to an inertial position again after a rather long vertical distance. As the chain elements are all connected by a bearing-type structure, they all have to move with the same velocity v. In the following we shall consider only the stationary case and not the evolution from a state of rest to the fountain with velocity v. Most models known use the idea of a continuous model. In the following we shall apply multibody theory in addition by modeling each bead with its connections separately. Results confirm the approach. PubDate: 2017-09-01 DOI: 10.1007/s00419-017-1260-y Issue No:Vol. 87, No. 9 (2017)

Authors:Guoqi Zhao; Wenfeng Hao; Xixuan Sheng; Yin Luo; Guangping Guo Pages: 1427 - 1438 Abstract: Abstract In this study, the interactions of matrix crack with inclusions of different shapes were investigated using the method of caustic. First, the specimens with inclusions of different shapes were prepared, where glass was used as inclusion and epoxy was used as matrix. Then, caustic experiments were conducted, and the typical caustic spots at the crack tip with varied distances from three different shapes of inclusions were obtained. Ultimately, the stress intensity factors of the cracks shielded by different shapes of inclusions were extracted from the caustic spots, and finite element simulations were conducted to verify the experimental results using the ABAQUS software. The results show that the stress intensity factors measured by the caustic method are in good agreement with the finite element simulation results. PubDate: 2017-09-01 DOI: 10.1007/s00419-017-1261-x Issue No:Vol. 87, No. 9 (2017)

Authors:Jialin Tian; Yinglin Yang; Lin Yang Pages: 1439 - 1451 Abstract: Abstract For the nonlinear vibration of the drill string in the drilling process, vibration characteristics analysis and experimental study of the drill string are conducted, which are to analyze the drill string dynamic characteristics with wellbore random friction force on the basis of the horizontal well. Firstly, considering the wellbore random friction force, the analysis models of the drill string vibration and the drilling efficiency of the horizontal well are established. Then, the establishment method of the random wellbore friction field is also obtained. With the combination of solution expressions of each force in the vibration equation, the discrete method of the dynamic model is established. According to the experimental test, the key input parameters are determined, and then, the example analysis of the vibration model is conducted. With the comparison of the experimental test and theoretical calculation, the influences of key parameters on the dynamic characteristics of the drill string are analyzed to verify the accuracy of the analysis model. The results can provide a new insight to the researches of the drill string dynamics, especially for the complex well, such as ultra-deep well, branch well and directional well, in which the wellbore friction has a significant influence on the result of the drill string kinetics. Moreover, the results can offer an important guidance for the design and application of new downhole tools. PubDate: 2017-09-01 DOI: 10.1007/s00419-017-1262-9 Issue No:Vol. 87, No. 9 (2017)

Authors:Xuhui Liao; Shunming Li; Lianying Liao; Haodong Meng Pages: 1453 - 1463 Abstract: Abstract An important procedure in transfer path analysis (TPA) is to measure the frequency response functions (FRFs) of the decoupled passive subsystem. The classical TPA method obtains the passive subsystem’s FRFs by direct measuring when the system is disassembled. The main shortcoming of the classical method to measure the FRFs is that it is time-consuming due to the necessity to dismount the active part. In this paper, a novel method is proposed to estimate the passive subsystem’ FRF matrix without disassembling the coupled mechanical structure. The key idea of this method is that the effect of a coupled subsystem will be canceled out if the links which connect this subsystem with the other one have no deformation, since the coupled systems influence each other only through the links which can be regarded as combinations of connecting springs and dampers. Following this idea, the expression of the passive subsystem’s FRF matrix can be deduced from the entire system’s FRF matrix directly. The proposed method in this paper is called the virtual decoupling method, since the decoupling is not ‘real’ but ‘virtual’. Obviously, the actual decoupling procedure is avoided so that the shortcoming mentioned above is overcome. The method is validated by a numerical model and a finite element model. PubDate: 2017-09-01 DOI: 10.1007/s00419-017-1264-7 Issue No:Vol. 87, No. 9 (2017)

Authors:S. H. Mirtalaie; M. A. Hajabasi Pages: 1465 - 1494 Abstract: Abstract The nonlinear axial-lateral-torsional free vibration of the rotating shaft is analyzed by employing the Rayleigh beam theory. The effects of lateral, axial and torsional deformations, gyroscopic forces and rotary inertia are taken into account, but the shear deformations are neglected. In the new developed dynamic model, the nonlinearities are originated from the stretching of beam centerline, nonlinear curvature and twist and inertial terms which leads to the coupling between the axial, lateral and torsional deformations. The deformed configuration of the cross section of the beam is represented by the axial and lateral deformations, also the geometry of the beam in the deformed configuration is represented by Euler angles. A system of coupled nonlinear differential equations is obtained which is examined by the method of multiple scales and the nonlinear natural frequencies are determined. The accuracy of the solutions is inspected by comparing the free vibration response of the system with the numerical integration of the governing equations. The effect of the spin speed and radius-to-length ratio of the rotating shaft on the free vibrational behavior of the system is inspected. The study demonstrates the effect of axial-lateral-torsional coupling on the nonlinear free vibrations of the rotating shaft. PubDate: 2017-09-01 DOI: 10.1007/s00419-017-1265-6 Issue No:Vol. 87, No. 9 (2017)

Authors:B. Emek Abali; Wolfgang H. Müller; Francesco dell’Isola Pages: 1495 - 1510 Abstract: Abstract In continuum mechanics, there exists a unique theory for elasticity, which includes the first gradient of displacement. The corresponding generalization of elasticity is referred to as strain gradient elasticity or higher gradient theories, where the second and higher gradients of displacement are involved. Unfortunately, there is a lack of consensus among scientists how to achieve the generalization. Various suggestions were made, in order to compare or even verify these, we need a generic computational tool. In this paper, we follow an unusual but quite convenient way of formulation based on action principles. First, in order to present its benefits, we start with the action principle leading to the well-known form of elasticity theory and present a variational formulation in order to obtain a weak form. Second, we generalize elasticity and point out, in which term the suggested formalism differs. By using the same approach, we obtain a weak form for strain gradient elasticity. The weak forms for elasticity and for strain gradient elasticity are solved numerically by using open-source packages—by using the finite element method in space and finite difference method in time. We present some applications from elasticity as well as strain gradient elasticity and simulate the so-called size effect. PubDate: 2017-09-01 DOI: 10.1007/s00419-017-1266-5 Issue No:Vol. 87, No. 9 (2017)

Authors:Arpita Mandal; Chaitali Ray; Salil Haldar Pages: 1511 - 1523 Abstract: Abstract The present paper deals with free vibration analysis of laminated composite skew plates with and without cut-outs. The experimental investigation along with numerical simulation has been presented in this paper for complete understanding of the dynamic behaviour of laminated skew plates with cut-out. Glass fibre-reinforced laminated composite plates have been prepared by resin infusion process using vacuum bagging technique in the laboratory for experimental analysis. The experimental studies have been carried out on skew composite plates with varying size of cut-out placed at the centre. The numerical analysis has been carried out by developing a computer code in MATLAB. Special attention is drawn on the formulation of mass matrix by considering effect of rotary inertia. The results obtained by the finite element formulation using nine-noded isoparametric plate-bending elements are validated by comparing the results from relevant published literature. The numerical and experimental data are then compared for experimental verification of present investigation. The consistency of mode shapes between experimental and numerical investigations is checked by using modal assurance criteria. PubDate: 2017-09-01 DOI: 10.1007/s00419-017-1267-4 Issue No:Vol. 87, No. 9 (2017)

Authors:Shaowei Hu; Jiang Yu; Congjie Wei; Zhaoguang Zhang Pages: 1525 - 1539 Abstract: Abstract Due to non-uniform distribution of structural deformation along its transverse width direction, shear lag behavior widely exists in composite structure with multi-box and leads to structural instability and destruction. To in-depth explore its mechanical mechanism, a type of steel–concrete composite structure with double-box (the DBSCCS model) is proposed, and its longitudinal warping shape functions are set up. Based on the minimum potential energy principle, governing differential equations of the DBSCCS model and its boundary conditions are deduced by means of the variational method. And then, its strain functions and shear lag coefficients are also obtained under concentrated loading and symmetrical loading, respectively. What is more, experimental verification and its related parametric sensitivity analysis are launched based on deduced longitudinal strain functions and shear lag coefficients. Through this analysis, it shows that this method can be used to illustrate and predict shear lag characteristics for this type of the DBSCCS model. That further suggests that it provides a more reference value for engineering design and structure optimization in some extent for the composite structures. PubDate: 2017-09-01 DOI: 10.1007/s00419-017-1268-3 Issue No:Vol. 87, No. 9 (2017)

Authors:Pei Zhu; Xingmin Ren; Weiyang Qin; Yongfeng Yang; Zhiyong Zhou Pages: 1541 - 1554 Abstract: Abstract In this paper, we study the characteristics of a tri-stable energy harvester (TEH) that is realized by the effect of magnetic attractive forces. The electromechanical model is established, and the corresponding coupling equations are derived by Euler–Lagrange equation. The potential energy indicates that the TEH’s potential well depths are the determinant factors for performance and can be designed such that the snap-through is easy to be elicited. We find that the TEH exhibits the best performance when the three potential well’s depths are nearly identical. To highlight the advantage of the TEH in harvesting energy, the comparisons between the tri-stable energy harvester and the bi-stable energy harvester (BEH) are carried out in simulations and experiments. The results prove that the TEH is preferable to the BEH in energy harvesting. The validation experiments show that the TEH owns a wide range of frequency of snap-through and high output voltage. PubDate: 2017-09-01 DOI: 10.1007/s00419-017-1270-9 Issue No:Vol. 87, No. 9 (2017)

Authors:Abdolrasoul Ranjbaran; Mohammad Ranjbaran Pages: 1555 - 1565 Abstract: Abstract Beam, column, plate, and any other structure, under full or partial compressive loading, are prone to failure by the buckling phenomenon. At the instant of failure, the structure may be in unpredictable elastic, elastic–plastic, full plastic, cracked, or other forms of deterioration state. Therefore, in spite of so much study, there is no definite solution to the problem. In this paper a unified, simple, and exact theory is proposed where buckling is considered as the change of state of structure between intact and collapsed states, and then the buckling capacity is innovatively expressed via states and phenomena functions, which are explicitly defined as functions of state variable. The state variable is determined by calibration of the structure slenderness ratio. The efficacy of the work is verified via concise mathematical logics, and comparison of the results with those of the others via seven examples. PubDate: 2017-09-01 DOI: 10.1007/s00419-017-1273-6 Issue No:Vol. 87, No. 9 (2017)

Authors:C. B. Silbermann; J. Ihlemann Abstract: Abstract Continuum dislocation theory (CDT) allows the consideration of dislocation ensembles by introducing the dislocation density tensor. Though the kinematics of geometrically linear CDT are well established, the closure of governing field equations is not finished yet. The present study now brings together different principles for such a closure: It is shown how the field equations for the CDT can be obtained from potential energy minimization and from the phase field approach. These two energetic methods are integrated into a generic thermodynamic framework with twofold benefit: First, the rigorous thermodynamic treatment allows clarifying physical consequences of the energetic methods, among them the proof of thermodynamic consistency. Second, the framework provides a basis for consistent extensions of CDT. In this way, a new dynamic formulation of CDT is presented, which enables the analysis of the evolution of dislocation structures during plastic deformation. Moreover, a variety of possible dissipative phenomena is considered and the mechanical balance laws are deduced. For two special cases, the field equations are derived in the strong form and the stability of the solution is analyzed. Next, a flexible numerical solution algorithm is presented using the finite difference method. Solutions of various initial boundary value problems are presented for the case of plane deformations. Therefore, some of the dissipative phenomena are further investigated and two distinct sources of the Bauschinger effect are identified. Special attention is also given to different boundary conditions and their effect on the solution. For the case of uniaxial compression, the numerical results are confronted with experimental data. Thus, the simulations are validated and a new consistent interpretation of the experimental results is achieved. PubDate: 2017-09-20 DOI: 10.1007/s00419-017-1296-z

Authors:Fabio C. Figueiredo; Lavinia A. Borges Abstract: Abstract The aim of this paper is to propose a limit analysis formulation concerning prescription of non-homogeneous velocities and unilateral conditions with friction at structures’ contact interfaces. This formulation is especially suitable for determining the limit state conditions in structures in which the external action is defined by prescribed velocities on boundaries, particularly if the contact interface is not planar and the force distribution is not known a priori. The requirement of body’s non-penetrability is attended by applying the unilateral conditions at normal direction and a sliding rule based on Coulomb friction law at tangential direction. Under limit state, if there is sliding between the contact surfaces, the external collapse power is consumed by plastic and friction dissipation. As applications, the influence of friction coefficient at tool–specimen interface at scratch test problem and the lateral resistance of a soil due to lateral movement of a partially embedded pipe are investigated. PubDate: 2017-09-20 DOI: 10.1007/s00419-017-1304-3

Authors:Helal Chowdhury; Konstantin Naumenko; Holm Altenbach; Manja Krüger Abstract: Abstract Determining critical stresses for different slip systems is one of the most important parts in crystal plasticity modeling of anisotropy. However, the task of finding individual critical resolved shear stress (CRSS) for every single slip system, if not impossible, is formidable and a delicate one especially if the microstructure is very complex. Slip family-based, mechanism-based and morphology-based (e.g., phase interface) slip systems classification and hence determining CRSS consistent with experimental measurements are often used in crystal plasticity. In this work, a novel approach to determining CRSS at high homologous temperature has been proposed by crystal plasticity modeling of rate-dependent anisotropy. Two-internal-variable-based phenomenological crystal viscoplasticity model is adopted for simulating isothermal, two-phase, single-crystal-like Al-rich lamellar Ti–61.8at.%Al binary alloy at high-temperature compression state ( \(1050\,^\circ \hbox {C}\) ) by employing finite strain and finite rotation framework. To the best of authors’ knowledge, this is the first micromechanical modeling attempt with long-period superstructures. Conventional approaches related to CRSS estimation are also compared with the proposed one. Our material parameters are based on calibrating three different sets of compressive stain rate-controlled plasticity data taken from the loading of two different lamellar directions. It is revealed that the proposed approach works fine for rate-dependent anisotropy modeling, while other conventional approaches highly under- or overestimate available anisotropic experimental behavior of this alloy. PubDate: 2017-09-20 DOI: 10.1007/s00419-017-1291-4

Authors:Roman Bogacz; Włodzimierz Kurnik Abstract: Abstract The paper contains description of dynamical problems connected with the self-excitation and kinematic excitation of rail vehicle and railway track. Some phenomena are presented which may create high loads resulting in track degradation and fatigue of wheelset axles. An alternative approach to rail vehicle hunting is proposed. Examples of experimental investigations are given which show that the dynamical load acting on the track can be much higher than the static load. PubDate: 2017-09-19 DOI: 10.1007/s00419-017-1298-x

Authors:İsa Çömez Abstract: Abstract In this study, frictional moving contact problem for a rigid cylindrical punch and an elastic layer is considered. The punch is subjected to concentrated normal and tangential force, and moves steadily with a constant subsonic velocity on the boundary. The problem is reduced to a singular integral equation of the second kind, in which the contact stress and the contact area are the unknowns, and it is treated using Fourier transforms and the boundary conditions for the problem. The numerical solution of the singular integral equation is obtained by using the Gauss–Jacobi integration formulas. Numerical results for the contact stress and the contact area are given. The results show that with increasing values of relative moving velocity, contact width between the moving punch and the layer increases, whereas contact stress decreases. PubDate: 2017-09-16 DOI: 10.1007/s00419-017-1306-1

Authors:Sonia Parvanova; Georgi Vasilev; Petia Dineva Abstract: Abstract The paper deals with numerical evaluation of the scattered wave and dynamic stress concentration fields in a finite anisotropic solid containing multiple nano-cavities. 2D plane-strain state and in-plane wave motion are assumed. The proposed mechanical model combines classical elastodynamic theory for the bulk general anisotropic solid and the Gurtin–Murdoch theory of surface elasticity assuming localized constitutive equation for the infinitely thin interface between the cavity and the matrix. The developed computational methodology is based on the following: (a) displacement boundary integral equations along existing boundaries using the analytically derived through Radon transform fundamental solution of the equation of motion of the bulk anisotropic solid; (b) non-classical boundary conditions of the Gurtin–Murdoch model along the interface between the matrix and cavities taking into consideration a jump in the stresses as one moves from the bulk material to the cavity due to the presence of surface elasticity; and (c) elastic-viscoelastic correspondence principle. The accuracy of the developed software is proven by comparisons of the obtained results solved by boundary element method and finite element method. A detailed parametric study reveals the sensitivity of the wave field to different key factors such as size, number and configuration of the cavities, surface and bulk material properties. PubDate: 2017-09-13 DOI: 10.1007/s00419-017-1303-4

Authors:R. Bagheri Abstract: Abstract A piezoelectric half-plane weakened by several horizontal cracks is investigated under anti-plane mechanical and in-plane electrical impacts. The distributed dislocation and integral transform techniques are employed to construct integral equations of the multiple dynamic cracks embedded in the piezoelectric half-plane. At first, the stress and the electric fields in the piezoelectric half-plane are calculated by using pattern. Then, by determining distributed dislocation density on the crack surface, a system of singular integral equations with Cauchy-type singularity is derived. The dynamic field stress intensity factors are determined by using the numerical Laplace inversion and dislocation densities. Finally, several examples are solved and the effects of the geometrical parameters and cracks configuration are graphically obtained upon the dynamic field intensity factors. PubDate: 2017-09-12 DOI: 10.1007/s00419-017-1305-2

Authors:Włodzimierz Kurnik; Piotr M. Przybyłowicz; Roman Bogacz Abstract: Abstract The paper is inspired by recent experiment with a two-member discrete column subjected to dry friction force of interaction between the moving column and a moving plane. The experiment was presented in the form of a YouTube film and recommended as an experimental evidence for flutter in the Ziegler column. We show that the tested mechanism cannot be identified with the original Ziegler column. PubDate: 2017-08-31 DOI: 10.1007/s00419-017-1294-1

Authors:Airong Liu; Mark Andrew Bradford; Yong-Lin Pi Abstract: Abstract Nonlinear in-plane multiple equilibria and buckling of pinned–fixed shallow circular arches under an arbitrary radial concentrated load are investigated. Analytical solutions for the multiple nonlinear equilibria, buckling and limit points are derived. New findings are: (1) pinned–fixed shallow arches under the arbitrary concentrated load have multiple stable and unstable equilibria; (2) the position of the arbitrary concentrated load and the modified slenderness of the arch influence the number of multiple equilibria and limit points as well as the first buckling load significantly; (3) a pinned–fixed arch under the arbitrary concentrated load can buckle in a limit point instability mode, but not in a bifurcation mode; (4) when the load is located between the crown and the pinned end, the buckling load is lower than that when the load is located between the crown and the fixed end, and (5) in addition to limit points, the nonlinear equilibria of pinned–fixed arches under the arbitrary concentrated load have inflexion points, which corresponds to specific modified slenderness switching the number of equilibria and limit points or switching buckling and no buckling behaviour. The analytical solutions for inflexion points and specific modified slenderness, and for the corresponding load, axial force and displacement are also derived for the first time in the literature. Comparisons with the finite element results have shown that the analytical solutions can accurately predict the multiple equilibria, limit points, inflexion points, and buckling load of shallow pinned–fixed arches under the arbitrary concentrated load. PubDate: 2017-08-31 DOI: 10.1007/s00419-017-1300-7

Authors:Sebastian Korczak Abstract: Abstract This contribution presents some nonlinear behaviors of an underactuated mechanical system under the trajectory tracking task. Presented hovercraft model is fully controlled by the computed torque algorithm with the pseudoinverse operation and proportional-derivative feedback. General form of errors dynamic equation gives possibility to analyze their behavior. These errors present irregular behaviors because of the input force limitations. Positive values of highest Lyapunov exponent and Fourier spectrum shape prove chaotic behavior of the system. PubDate: 2017-08-30 DOI: 10.1007/s00419-017-1297-y