• Journal of applied mathematics & informatics
• /
• v.15 no.1_2
• /
• pp.1-28
• /
• 2004

We first apply a first order splitting to a semilinear reaction-diffusion equation and then discretize the resulting system by an $H^1$-Galerkin mixed finite element method in space. This semidiscrete method yields a system of differential algebraic equations (DAEs) of index one. A priori error estimates for semidiscrete scheme are derived for both differ-ential as well as algebraic components. For fully discretization, an implicit Runge-Kutta (IRK) methods is applied to the temporal direction and the error estimates are discussed for both components. Finally, we conclude the paper with a numerical example.

• Coupled systems mechanics
• /
• v.10 no.1
• /
• pp.1-19
• /
• 2021

Composite nature of the masonry structures in general causes complex and non-linear behaviour, especially in intense vibration conditions. The presence of different types and forms of structural elements and different materials is a major problem for the analysis of these type of structures. For this reason, the analysis of the behaviour of masonry structures requires a combination of experimental tests and non-linear mathematical modelling. The famous UNESCO Heritage Old Bridge in Mostar was selected as an example for the analysis of the global behaviour of reinforced stone arch masonry bridges. As part of the experimental research, a model of the Old Bridge was constructed in a scale of 1:9 and tested on a shaking table platform for different levels of seismic excitation. Non-linear mathematical modelling was performed using a combined finite-discrete element method (FDEM), including the effect of connection elements. The paper presents the horizontal displacement of the top of the arch and the failure mechanism of the Old Bridge model for the experimental and the numerical phase, as well as the comparison of the results. This research provided a clearer insight into the global behaviour of stone arch masonry structures reinforced with steel clamps and steel dowels, which is significant for the structures classified as world cultural heritage.

• Computers and Concrete
• /
• v.3 no.2_3
• /
• pp.79-90
• /
• 2006

This paper describes the results of an investigation into high mass-low velocity impact behaviour of reinforced concrete beams. Tests have been conducted on fifteen 2.7 m or 1.5 m span beams under drop-weight loads. A high-speed video camera has been used at rates of up to 4,500 frames per second in order to record the crack formation, propagation, particle spallation and scabbing. In some tests the strain in the reinforcement has been recorded using "Durham" strain gauged bars, a technique developed by Scott and Marchand (2000) in which the strain gauges are embedded in the bars, so that the strains in the reinforcement can be recorded without affecting the bond between the concrete and the reinforcement. The impact force acting on the beams has been measured using a load cell placed within the impactor. A high-speed data logging system has been used to record the impact load, strains, accelerations, etc., so that time histories can be obtained. This research has led to the development of computational techniques based on combined continuum/discontinuum methods (finite/discrete element methods) to permit the simulation of impact loaded reinforced concrete beams. The implementation has been within the software package ELFEN (2004). Beams, similar to those tested, have been analysed using ELFEN a good agreement has been obtained for both the load-time histories and the crack patterns.

• Coupled systems mechanics
• /
• v.11 no.2
• /
• pp.167-198
• /
• 2022

In this paper we deal with classical instability problems of heterogeneous Euler beam under conservative loading. It is chosen as the model problem to systematically present several possible solution methods from simplest deterministic to more complex stochastic approach, both of which that can handle more complex engineering problems. We first present classical analytic solution along with rigorous definition of the classical Euler buckling problem starting from homogeneous beam with either simplified linearized theory or the most general geometrically exact beam theory. We then present the numerical solution to this problem by using reduced model constructed by discrete approximation based upon the weak form of the instability problem featuring von Karman (virtual) strain combined with the finite element method. We explain how such numerical approach can easily be adapted to solving instability problems much more complex than classical Euler's beam and in particular for heterogeneous beam, where analytic solution is not readily available. We finally present the stochastic approach making use of the Duffing oscillator, as the corresponding reduced model for heterogeneous Euler's beam within the dynamics framework. We show that such an approach allows computing probability density function quantifying all possible solutions to this instability problem. We conclude that increased computational cost of the stochastic framework is more than compensated by its ability to take into account beam material heterogeneities described in terms of fast oscillating stochastic process, which is typical of time evolution of internal variables describing plasticity and damage.