# The Analysis of Stress Waves at a Junction of Beam and String

Authors: Mu Chen; Ronghui Wang; Niujing Ma; Zhuojie Zhang
Journal: Advances in Civil Engineering
 (2018)

---

## Abstract

In the bridge engineering, there are some problems about the dynamics that traditional theory cannot solve. So, the theory about stress waves is introduced to solve the related problems. This is a new attempt that the mechanic theory is applied to practical engineering. The stress wave at a junction of the structure composed of beams and strings is investigated in this paper. The structure is studied because the existence of a soft rope makes the transmission of the force in the bridge structure different from the traditional theory, and it is the basis for further research. The equilibrium equations of the displacement and the internal force are built based on the hypothesis. The fast Fourier transform (FFT) numerical algorithm is used to express an incident pulse of arbitrary shape. The analytical solutions are substantiated by comparing with the finite element programs. The conclusion that if the cross section of the string is relatively small, then the energy density of the structure is relatively large, which is disadvantageous to the structure, can be obtained from this paper.


---

## Body

## 1. Introduction

Complex structures, such as arch bridge, suspension bridge, and cable-stayed bridge, are used increasingly in the world. Different from the traditional theory, the vehicle load is regarded as the vibration load in the bridge structure. The theory about stress waves is introduced to study the influence of the vibration load. The propagation speed of longitudinal waves and transverse waves in the concrete and steel is in the range of 1000 m/s to 5000 m/s, and in those large structures, the time scale of vehicle load and stress wave in balance are in the same order of magnitude. So, different from the traditional theory of structural dynamics that ignores stress wave propagation and interaction before the static equilibrium, the spread of the stress process is what we concern. Investigation of stress waves at a junction composed of beams and strings is the most important and foundational work.Lee and Kolsky [1] investigate the generation of stress pulses at the junction of two noncollinear rods. They used the simple one-dimensional theory to describe the propagation of longitudinal waves and used the Timoshenko equation for flexural propagation. The analytic relations between the incident pulse and the four generated pulses are derived both for an incident longitudinal pulse and an incident flexural pulse. Atkins and Hunter [2] study the propagation of longitudinal elastic waves around right-angled corners in rods of square cross section. In their work, the rod is assumed to behave as a rigid body, and it was demonstrated that longitudinal elastic waves in rods of square cross section are transmitted in part around right-angled corners. Doyle [3] used the fast Fourier transform algorithm that allows the force history to be determined to establish a strain relation for the transverse impact of beams.Mace [4] studied wave reflection and transmission in beams. In his paper, it is seen that incident near fields can give rise to substantial propagation components. Doyle and Kamle [5, 6] had taken an experimental study of the reflection and transmission of flexural waves at an arbitrary T-joint. It could be obtained from their study that the dynamic response is not sensitive to particular values of the joint model, and it depends primarily on the member arrangement and only secondarily on the particular shape and mass of the joint. Horner and White [7] predicted vibrational power transmission through bends and joints in beam-like structures. They built some models which determine the wave type which carries most power in each section of the system. By establishing the wave type, it is then possible to apply the most suitable vibration control technique. Sack [8] studied transverse oscillations in travelling strings. He found that the frequency response to a simple harmonic disturbance and the envelope of the oscillations are similar to those of a clamped string, but the phase of oscillation varies from point to point so that the instantaneous configuration of the string is not a sine curve. Little and Heywoodb [9] had described the axial stress waves in an elastic-perfectly-plastic bar. Different from the classical elastic solution, the load is greater than half the yield load in their study. The paper describes and verifies the model and presents the predictions that it gives for various levels of load. Chen [10] had thoroughly reviewed the research of transverse vibrations of axially moving strings and their control. Chen presents the governing equations with large amplitude and reviews progress on discretized or direct approximate analytical analyses and numerical approaches based on the Galerkin method or the finite difference method. Guo et al. [11] investigated the elastic wave and energy propagation in angled beams. Experimental and numerical methods were taken. It could demonstrate that, after several reflections from and transmission across the bends, energy is progressively smeared throughout the entire beam and does not concentrate at any particular segment. Li et al. [12] developed a fractal damage joint model based on the fractal damage theory to investigate the transmission and reflection of stress waves across joints. The fractal geometrical characteristics of joint surfaces are investigated by using a laser profilometer to scan the joint surfaces in their study. Lancioni et al. [13] studied the scattering of seismic waves induced by the presence of underground cavities in homogeneous soils by using the numerical method. Singh and Goya [14] established the governing equations of a transversely isotropic microstretch material, which are specialized in x-z plane. Plane wave solutions of these governing equations result in a biquadratic velocity equation. They found that, for a specific material, numerical simulation in the presence as well as in the absence of microstretch shows that the coupled longitudinal displacement (CLD) wave is the fastest wave and the coupled transverse microrotational (CTM) wave is the observed slowest wave. Piccolroaza et al. [15] presented a novel analysis of Floquet–Bloch flexural waves in a periodic lattice-like structure consisting of flexural beam ligaments. A special feature of this structure is in the presence of the rotational inertia, which is commonly neglected in conventional models of the Euler–Bernoulli type. Numerical examples were taken to illustrate directional localization, negative refraction, localization at an interface, and neutrality for propagating plane waves across a structured interface for a frequency range corresponding to a Dirac cone.The previous studies focused on the distribution of the stress wave in the junction of the beams. According to the need of practical engineering, the distribution of the stress wave that was produced by transverse impact in the junction of the beam and string is the focus of the study in this paper.


## 2. Theoretical Analysis

### 2.1. Wave Propagation through the Junction

We will here consider the reflection of a stress pulse at the boundary of beams and strings of different diameters and materials (Figure1(a)). There is an angle β between two beams and α between beam 1 and string 2. We call the flexural displacement due to the incident wave v0. We denote the displacement of the reflected longitudinal wave by u1, the lateral displacement produced by the flexural wave by v1 in beam 1, the longitudinal displacement in string 2 by u2, the lateral displacement in string 2 by v2, the longitudinal displacement in beam 3 by u3, and the lateral displacement in beam 3 by v3. The joint is modeled as a rigid body, and the internal forces are shown in Figure 1(b).Figure 1
The displacement and stress of the three-member intersection. (a) The displacement of the components. (b) The internal forces of the junction.
(a)
(b)Continuity of the displacement of beam 1 and string 2 leads to(1)u1=u2cosα−v2sinα,(2)v0−v1=u2sinα+v2cosα.Continuity of the displacement of beams 1 and 3 leads to(3)u1=u3cosβ−v3sinβ+L2sinβ∂v3∂x3,(4)v0−v1=u3sinβ+v3cosβ−L21+cosβ∂v3∂x3,(5)∂v0∂x1−∂v1∂x1=−∂v3∂x3.The balance of the forces at the junction leads to(6)E1I1∂2v0−v1∂x12−L2E1I1∂3v1∂x13=E3I3∂3v3∂x33+L2E3I3∂3v3∂x33−Ij∂v1∂x1,(7)E1A1∂u1∂x1=E2A2∂u2∂x2cosα+E3A3∂u3∂x3cosβ+E3I3∂3v3∂x33sinβ−mj∂2u1∂t2,(8)E1I1∂3v0−v1∂x13=E2A2∂u2∂x2sinα−E3A3∂u3∂x3sinβ−E3I3∂3v3∂x33cosβ−mj∂u2∂t,where E1I1 and E3I3 are bending stiffness of beams 1 and 3. E1A1, E2A2, and E3A3 are axial stiffness of beams 1 and 3 and string 2. L is a parameter characterizing the size of the joint. Ij and mj are the moment of inertia and quality of the joint.


### 2.2. Characterization of the Wave

In the paper by Doyle and Kamle [5, 6], the fast Fourier transform (FFT) numerical algorithm was used to express an incident pulse of arbitrary shape. The FFT is defined as(9)Xk=∑n=0N−1xne−2πikn/N,n,k=0,1,2,…,N−1,where x(n) is the value of the signal at time (n − 1)Δt and X(k) is the FFT at frequency k/NΔt. Different from the previous study, the Timoshenko beams are used in this paper. The displacement can be written as(10)v=Ae−ik1x+Be−k2xeiwt,where(11)k1=k1+14wkcs4+12wkcs21/2,k=ρAEI1/4w,k2=k1+14wkcs4−12wkics21/2,cs=Gρ.The displacements and strains can be written as(12)vx,t=∑n=0NAne−ik1x+Bnek2xeiwt,ex,t=−h2∂2v∂x2=h2∑n=0NAnk12e−ik1x−Bnk22ek2xeiwt,where An, Bn, and k1, k2 are understood to be frequency dependent. If the FFT is taken of the strain history, then it can be approximated with a finite set of frequency components, that is,(13)ex,t=∑n=0Ne^x,wneiwnt,where e^x,wn is the FFT component at frequency wn of the strain history at x. If the strain is sampled at two locations, then An and Bn at discrete frequency can be solved as follows:(14)An=1k2n2e^x1,wne−k2nx1−e^x2,wne−k2nx2Q,Q=k1n2k2n2e−ik1nx1−k2nx1−e−ik1nx2−k2nx2,Bn=1k1n2e^x1,wneik1nx1−e^x2,wneik1nx2L,L=k2n2k1n2ek2nx2+ik1nx1−ek2nx1+ik1nx1.This is different from that reported by Doyle and Kamle [5, 6]. They assumed that the accuracy of the simpler Bernoulli–Euler model was up to the requirements of the study. This is generally true when the wavelength is larger than the size of the beam section. However, if the wavelength and the size of the beam section are approximately equal, then the Timoshenko model can get more accurate results.Finally, the incident and transmitted displacements are expanded in the forms:(15)u1=∑n=0ND1neiwt−x1/c0,(16)v1=∑n=0ND2ne−ik11x1+D3ne−k21x1eiwt,(17)u2=∑n=0ND4neiwt−x2/c2,(18)v2=∑n=0ND5neiwt−x2/c′2,(19)u3=∑n=0ND6neiwt−x3/c3,(20)v3=∑n=0ND7ne−ik13x3+D8ne−k23x3eiwt,where c2 and c′2 are velocities of the longitudinal and flexural waves in string 2. The hypotheses that v2 is small and perpendicular to the x-axis at any time are taken in this paper. So, the tension T of the string does not change over time. c2 and c′2 can be expressed in the forms:(21)c2=E2ρ2,c′2=Tρ2.c
3 in (19) is the velocity of the longitudinal wave and it can be found using the following equation:(22)c3=E3ρ3.For convenience, the joint is located atx = 0. From (1)–(8), the eight unknowns D1n,D2n,…,D8n can be determined in terms of the coefficients An and Bn which are given by the initial conditions. Due to the limitation of the space, the expression of the coefficients will not be given in this article.


## 2.1. Wave Propagation through the Junction

We will here consider the reflection of a stress pulse at the boundary of beams and strings of different diameters and materials (Figure1(a)). There is an angle β between two beams and α between beam 1 and string 2. We call the flexural displacement due to the incident wave v0. We denote the displacement of the reflected longitudinal wave by u1, the lateral displacement produced by the flexural wave by v1 in beam 1, the longitudinal displacement in string 2 by u2, the lateral displacement in string 2 by v2, the longitudinal displacement in beam 3 by u3, and the lateral displacement in beam 3 by v3. The joint is modeled as a rigid body, and the internal forces are shown in Figure 1(b).Figure 1
The displacement and stress of the three-member intersection. (a) The displacement of the components. (b) The internal forces of the junction.
(a)
(b)Continuity of the displacement of beam 1 and string 2 leads to(1)u1=u2cosα−v2sinα,(2)v0−v1=u2sinα+v2cosα.Continuity of the displacement of beams 1 and 3 leads to(3)u1=u3cosβ−v3sinβ+L2sinβ∂v3∂x3,(4)v0−v1=u3sinβ+v3cosβ−L21+cosβ∂v3∂x3,(5)∂v0∂x1−∂v1∂x1=−∂v3∂x3.The balance of the forces at the junction leads to(6)E1I1∂2v0−v1∂x12−L2E1I1∂3v1∂x13=E3I3∂3v3∂x33+L2E3I3∂3v3∂x33−Ij∂v1∂x1,(7)E1A1∂u1∂x1=E2A2∂u2∂x2cosα+E3A3∂u3∂x3cosβ+E3I3∂3v3∂x33sinβ−mj∂2u1∂t2,(8)E1I1∂3v0−v1∂x13=E2A2∂u2∂x2sinα−E3A3∂u3∂x3sinβ−E3I3∂3v3∂x33cosβ−mj∂u2∂t,where E1I1 and E3I3 are bending stiffness of beams 1 and 3. E1A1, E2A2, and E3A3 are axial stiffness of beams 1 and 3 and string 2. L is a parameter characterizing the size of the joint. Ij and mj are the moment of inertia and quality of the joint.


## 2.2. Characterization of the Wave

In the paper by Doyle and Kamle [5, 6], the fast Fourier transform (FFT) numerical algorithm was used to express an incident pulse of arbitrary shape. The FFT is defined as(9)Xk=∑n=0N−1xne−2πikn/N,n,k=0,1,2,…,N−1,where x(n) is the value of the signal at time (n − 1)Δt and X(k) is the FFT at frequency k/NΔt. Different from the previous study, the Timoshenko beams are used in this paper. The displacement can be written as(10)v=Ae−ik1x+Be−k2xeiwt,where(11)k1=k1+14wkcs4+12wkcs21/2,k=ρAEI1/4w,k2=k1+14wkcs4−12wkics21/2,cs=Gρ.The displacements and strains can be written as(12)vx,t=∑n=0NAne−ik1x+Bnek2xeiwt,ex,t=−h2∂2v∂x2=h2∑n=0NAnk12e−ik1x−Bnk22ek2xeiwt,where An, Bn, and k1, k2 are understood to be frequency dependent. If the FFT is taken of the strain history, then it can be approximated with a finite set of frequency components, that is,(13)ex,t=∑n=0Ne^x,wneiwnt,where e^x,wn is the FFT component at frequency wn of the strain history at x. If the strain is sampled at two locations, then An and Bn at discrete frequency can be solved as follows:(14)An=1k2n2e^x1,wne−k2nx1−e^x2,wne−k2nx2Q,Q=k1n2k2n2e−ik1nx1−k2nx1−e−ik1nx2−k2nx2,Bn=1k1n2e^x1,wneik1nx1−e^x2,wneik1nx2L,L=k2n2k1n2ek2nx2+ik1nx1−ek2nx1+ik1nx1.This is different from that reported by Doyle and Kamle [5, 6]. They assumed that the accuracy of the simpler Bernoulli–Euler model was up to the requirements of the study. This is generally true when the wavelength is larger than the size of the beam section. However, if the wavelength and the size of the beam section are approximately equal, then the Timoshenko model can get more accurate results.Finally, the incident and transmitted displacements are expanded in the forms:(15)u1=∑n=0ND1neiwt−x1/c0,(16)v1=∑n=0ND2ne−ik11x1+D3ne−k21x1eiwt,(17)u2=∑n=0ND4neiwt−x2/c2,(18)v2=∑n=0ND5neiwt−x2/c′2,(19)u3=∑n=0ND6neiwt−x3/c3,(20)v3=∑n=0ND7ne−ik13x3+D8ne−k23x3eiwt,where c2 and c′2 are velocities of the longitudinal and flexural waves in string 2. The hypotheses that v2 is small and perpendicular to the x-axis at any time are taken in this paper. So, the tension T of the string does not change over time. c2 and c′2 can be expressed in the forms:(21)c2=E2ρ2,c′2=Tρ2.c
3 in (19) is the velocity of the longitudinal wave and it can be found using the following equation:(22)c3=E3ρ3.For convenience, the joint is located atx = 0. From (1)–(8), the eight unknowns D1n,D2n,…,D8n can be determined in terms of the coefficients An and Bn which are given by the initial conditions. Due to the limitation of the space, the expression of the coefficients will not be given in this article.


## 3. Numerical Validation of the Cosine Wave

According to the theory of Fourier transform, the arbitrary wave produced by vehicle vibration load can be explained as a superposition of the different frequency and amplitude cosine waves. In order to simplify calculation, consider a beam-string system subjected to a cosine wave. The length of the string and beam is 1 m. A square whose side length is 25 mm is used as a beam cross section, and a circle whose radius is 1 mm is used as a string cross section. The top of the string and the right of the beam are fixed. The load is applied in the left side of the beam. In order to simplify the process of calculation, a cycle of cosine wave is calculated in this article, as shown in Figure2(a). The circular frequency and the amplitude of the wave are 2000π and 3 mm. For convenience of the calculation, two beams adopt the same material and cross section. The beam is made of the steel with Young’s modulus E = 207,800 N/mm2, Poisson’s ratio = 0.3, and density = 7.8e−9 t/mm3. The string is made of the steel with Young’s modulus E = 195,000 N/mm2, Poisson’s ratio = 0.3, and density = 7.8e−9 t/mm3.Figure 2
A cycle of cosine load (a) and the finite element of the model (b).
(a)
(b)The angleα changed with the different mechanical system of the bridge. To investigate the angle α between the beam and string on the distribution of the force in the joint, three angles corresponding to π/2, π/3, and π/6 are investigated, respectively. The angle β is kept as zero. Figure 2(b) shows the finite element of the model when α = π/2. The right end of the beam and the left end of the cable are consolidated. The load is applied to the left of the beam. Because the section size of the string is much smaller and the string is always tensioned under the cosine load, a two-node linear Timoshenko beam element (B31) was used and each bent beam and the string were discretized into 100 elements of uniform size.Figures3(a)–3(c) show the difference between the two methods (theoretical analysis and finite element analysis by Abaqus programs) in three locations (the left and right of the joint and the string). We choose the half cycle when the string bears the pulling force to discuss. The phase is ignored because it is not our concern, and the time domain is from 0 to 0.5 ms. It can be seen from Figure 3 that the amplitudes of Q1, N2, and Q3 are 23940 N, 405 N, and 23220 N by theoretical analysis and 20836 N, 328 N, and 20128 N by finite element analysis. It is noticed that the amplitude of the internal force by finite element analysis is slightly smaller than that by theoretical analysis. This is not surprising because of the loss in the process of the wave propagation and the error of the algorithm. It can be seen from the results of the comparison that the analytical solutions are in good agreement with finite element analysis prediction.Figure 3
The internal force by two methods (α = π/2). (a) Reflected wave. (b) The string. (c) The transmission wave.
(a)
(b)
(c)Figure4 illustrates the effect of α on the internal force of the joint. It is shown in Figures 4(a), 4(b), 4(e), and 4(f) that the effect of α on the shearing forces on both sides of the joint is minor compared to that on the pulling force of the string whatever the analytical method is. It means that the pulling force N2 exhibits a great relevance to the angle α. When α = π/2, π/3, and π/6, N2 = 405, 352, and 212 N by theoretical analysis and 328, 268, and 148 N by finite element analysis. It can be concluded that the pulling force has a nonlinear relationship with the angle α. It can also be observed that the rate of change of the pulling force N2 by theoretical analysis is close to that by finite element analysis. This can also prove the accuracy of the theoretical method on the other side.Figure 4
Effect ofα on the internal force of the joint (β = 0). (a) Q1 by theoretical analysis. (b) Q1 by finite element analysis. (c) N2 by theoretical analysis. (d) N2 by finite element analysis. (e) Q3 by theoretical analysis. (f) Q3 by finite element analysis.
(a)
(b)
(c)
(d)
(e)
(f)The energy density is another important research object in this paper.α = π/2, β = 0, and the cosine wave are taken to study the influence of the radius of the string on the energy density. The energy density can be calculated from mechanical energy as follows:(23)Wk1=12ξΔx∂v0−v1∂t2=−12ξΔxw2A−D2−D3eiwt2,(24)Wp1=12c′0ξΔx∂v0−v1∂x12=12c′0ξΔxAiwc′0+D2ik11+D3k212eiwt2,(25)w1=Wp1+Wk1ΔV=−12ρ1A−D2−D3w2−c′0Aiwc′0+D2ik11+D3k212eiwt2,(26)Wk2=12ξΔxπ∂u2∂t2=−12ρΔxD42w2eiwt2,(27)Wp2=12c22ξΔx∂u2∂x22=−12ρ2ΔxD42w2eiwt2,(28)w2=Wk1+Wp1ΔV=−ρ2D42w2eiwt2,(29)Wk3=12ξΔx∂v3∂t2=−12ξΔxw2D7+D82eiwt2,(30)Wp3=12c32ξΔx∂v3∂x32=12c32ξΔxD7ik13+D8k232eiwt2,(31)w3=Wk3+Wp3ΔV=12ρ3−D7+D82w2−c32D7ik31+D8k322eiwt2,where Wk, Wp, and w are the kinetic energy, potential energy, and energy density. ξ is the quality of the unit length. ΔV is the unit volume. The subscript numbers 1, 2, and 3 represent the locations on the left side of the joint, the string, and the right of the joint. From (25), (28), and (31), the energy densities w1,w2, and w3 can be obtained by the coefficients A, D2, D3, D4, D7, and D8.Figure5 illustrates the influence of the radius of the string on the energy density. It is shown in Figure 5 that the energy density w3 is much larger than that in locations 1 and 2. It can be explained that the reflected wave v1 has an unloading effect on the incident wave v0. It also can be observed that the energy density w3 varies little with the radius of the string from Figure 5(c), and that the energy densities w1 and w2 decrease with the increase of the radius of the string from Figures 5(a) and 5(b). Such a discrepancy is attributed to the sensitivity to the tensile strength of the string. As a matter of fact, with the increase of the radius, the string should take on more energy. This can also explain the phenomenon that the energy density w1 decreases. But the growth rate of the tensile strength of the string is much greater than that of the energy they sustain. This indicates that increasing the rigidity of the string properly can make the structure more reasonable. In practical engineering, the choice of the string will also be limited by other factors. Therefore, meeting the various conditions of the diameter of the string is the best choice.Figure 5
The influence of the radius of the string on the energy density (α = π/2, β = 0). (a) The energy density caused by the flexural wave v0 − v1. (b) The energy density caused by the longitudinal wave u2. (c) The energy density caused by the flexural wave v3.
(a)
(b)
(c)Referring to (25), the expressions of the energy density w4 caused by the longitudinal wave u1, w5 caused by the flexural wave v2, and w6 caused by the longitudinal wave u3 are as follows:(32)w4=−ρ1D12w2eiwt2,w5=−ρ2D52w2eiwt2,w6=−ρ1D62w2eiwt2.Figure6 plots the influence of α on the energy density when the radius of the string r = 1 mm. It can be seen from Figure 6(a) that the energy density w4 by the flexural wave v0 − v1 decreases with the increase of the angle α. The energy density w2 by the longitudinal wave u2 is opposite to the change in energy density w1. The energy density w3 by the flexural wave v3 is not influenced by angular change at all. It can be observed from Figures 6(d) and 6(f) that the energy densities w4 and w5 by the longitudinal wave reach maximum when the angle α is π/4. The reason is that the mechanical energy reaches maximum in that angle. The effect of the longitudinal wave on the whole energy is very small because the energy density caused by the longitudinal wave is two orders of magnitude smaller than that caused by the flexural wave. From Figure 6(e), we can obtain that the energy density w5 caused by the flexural wave v2 decreases with the increase of the angle. On the contrary, the energy density w1 increases with the increase of the angle, as shown in Figure 6(b). So, we can come to a conclusion that the flexural wave is the main wave pattern in the string when the angle is relatively small; otherwise, it is the longitudinal wave.Figure 6
The influence ofα on the energy density (r = 1 mm). (a) The energy density by the flexural wave v0 − v1. (b) The energy density by the longitudinal wave u2. (c) The energy density by the longitudinal wave v3. (d) The energy density by the longitudinal wave u1. (e) The energy density by the flexural wave v2. (f) The energy density by the longitudinal wave u3.
(a)
(b)
(c)
(d)
(e)
(f)


## 4. Conclusions

The stress waves at a joint of a structure composed of the string and beam are investigated in this paper. The governing equations are established by displacement balance and internal force balance. The fast Fourier transform (FFT) numerical algorithm is used to deal with the incident wave of arbitrary wave type. The unknown coefficients of the reflected wave and transmitted wave are obtained by governing equations. The analytical solutions are compared with the finite element analysis by Abaqus program, and some parameters are analyzed in this paper. We can yield the following conclusions: first, the results of theoretical analysis are slightly larger than those of the finite element analysis, and the error is about 10 percent. Because the wave is dissipative in the finite element analysis, the result is reasonable, and the theoretical analysis is correct. Second, through the analysis of the influence of the angleα on the internal force and energy density, we can come to a conclusion that the flexural wave is the main wave pattern in the string when the angle is relatively small; otherwise, it is the longitudinal wave. The last conclusion is that the tensile strength of the string has a great influence on the stress of the structure. If the tensile strength is relatively small, then the energy density of the string and beam will be relatively large, which is not conducive to the long-term use of the structure. This is important for practical engineering and further study.



---

*Source: 4863503-2018-03-07.xml*