• Coupled systems mechanics
• /
• v.12 no.5
• /
• pp.419-428
• /
• 2023

We have developed a model for estimating the parameters of viscous materials from indirect tensile tests for asphalt. This is a simple Burger nonlinear rheological two-cell model or standard model. At the same time, we begin to develop a more versatile and complex multi-cell model. The simple model is validated using experimental load-displacement results from laboratory tests: The recorded displacements are used as input values and the measured force data are simulated with the model. The optimal model parameters are estimated using the Levenberg-Marquardt method and a very good agreement between the experimental results and the model calculations is shown. However, not all parts of the model are active in the loading phase of the experiment, so we extended the validation of the model to the simulation of the relaxation behaviour. In this stage, the other model parameters are activated and the simulation results are consistent with the literature. At this stage, we have estimated the parameters only for the two-cell uniaxial model, but further work will include results for the multi-cell model.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.3
• /
• pp.243-259
• /
• 2016

The Richards equation for water movement in unsaturated soil is highly nonlinear partial differential equations which are not solvable analytically unless unrealistic and oversimplifying assumptions are made regarding the attributes, dynamics, and properties of the physical systems. Therefore, conventionally, numerical solutions are the only feasible procedures to model flow in partially saturated porous media. The standard Finite element numerical technique is usually coupled with an Euler time discretizations scheme. Except for the fully explicit forward method, any other Euler time-marching algorithm generates nonlinear algebraic equations which should be solved using iterative procedures such as Newton and Picard iterations. In this study, lumped mass and distributed mass in the frame of Picard and Newton iterative techniques were evaluated to determine the most efficient method to solve the Richards equation with finite element model. The accuracy and computational efficiency of the scheme and of the Picard and Newton models are assessed for three test problems simulating one-dimensional flow processes in unsaturated porous media. Results demonstrated that, the conventional mass distributed finite element method suffers from numerical oscillations at the wetting front, especially for very dry initial conditions. Even though small mesh sizes are applied for all the test problems, it is shown that the traditional mass-distributed scheme can still generate an incorrect response due to the highly nonlinear properties of water flow in unsaturated soil and cause numerical oscillation. On the other hand, non oscillatory solutions are obtained and non-physics solutions for these problems are evaded by using the mass-lumped finite element method.

• Structural Engineering and Mechanics
• /
• v.43 no.5
• /
• pp.561-582
• /
• 2012

Despite popularity of FEM in analysis of static and dynamic structural problems and the routine applicability of FE softwares, analytical methods based on simple mathematical relations is still largely sought by many researchers and practicing engineers around the world. Development of such analytical methods for analysis of free vibration of non-prismatic beams is also of primary concern. In this paper a new and simple method is proposed for determination of vibration frequencies of non-prismatic beams under variable axial forces. The governing differential equation is first obtained and, according to a harmonic vibration, is converted into a single variable equation in terms of location. Through repetitive integrations, integral equation for the weak form of governing equation is derived. The integration constants are determined using the boundary conditions applied to the problem. The mode shape functions are approximated by a power series. Substitution of the power series into the integral equation transforms it into a system of linear algebraic equations. Natural frequencies are determined using a non-trivial solution for system of equations. Presented method is formulated for beams having various end conditions and is extended for determination of the buckling load of non-prismatic beams. The efficiency and convergence rate of the current approach are investigated through comparison of the numerical results obtained to those obtained using available finite element software.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.37 no.9
• /
• pp.1069-1075
• /
• 2013

An effective method is necessary for the real-time analysis of an unmanned military robot. To achieve highly efficient simulations, a subsystem synthesis method has been developed. The subsystem synthesis method separately generates equations of motion for the base body and for the subsystem. The equations of motion are expressed by DAE, which consist of differential equations and algebraic equations. To increase the accuracy and efficiency of solutions, DAE solvers such as the Direct, CS (Constraint Stabilization), and GCP (Generalized Coordinate Partitioning) method are employed. In this study, the subsystem synthesis method is applied for effective multi-body dynamics analysis of an unmanned military robot, and a comparative study of three different DAE solvers is carried out.

• Smart Structures and Systems
• /
• v.18 no.6
• /
• pp.1125-1143
• /
• 2016

Forced vibration analysis of a simple supported viscoelastic nanobeam is studied based on modified couple stress theory (MCST). The nanobeam is excited by a transverse triangular force impulse modulated by a harmonic motion. The elastic medium is considered as Winkler-Pasternak elastic foundation.The damping effect is considered by using the Kelvin-Voigt viscoelastic model. The inclusion of an additional material parameter enables the new beam model to capture the size effect. The new non-classical beam model reduces to the classical beam model when the length scale parameter is set to zero. The considered problem is investigated within the Timoshenko beam theory by using finite element method. The effects of the transverse shear deformation and rotary inertia are included according to the Timoshenko beam theory. The obtained system of differential equations is reduced to a linear algebraic equation system and solved in the time domain by using Newmark average acceleration method. Numerical results are presented to investigate the influences the material length scale parameter, the parameter of the elastic medium and aspect ratio on the dynamic response of the nanobeam. Also, the difference between the classical beam theory (CBT) and modified couple stress theory is investigated for forced vibration responses of nanobeams.

• Structural Engineering and Mechanics
• /
• v.45 no.1
• /
• pp.95-110
• /
• 2013

In the recent decade, practical of wavelet technique is being utilized in various domain of science. Particularly, engineers are interested to the wavelet solution method in the time series analysis. Fundamentally, seismic responses of structures against time history loading such as an earthquake, illustrates optimum capability of systems. In this paper, a procedure using particularly discrete Haar wavelet basis functions is introduced, to solve dynamic equation of motion. In the proposed approach, a straightforward formulation in a fluent manner is derived from the approximation of the displacements. For this purpose, Haar operational matrix is derived and applied in the dynamic analysis. It's free-scaled matrix converts differential equation of motion to the algebraic equations. It is shown that accuracy of dynamic responses relies on, access of load in the first step, before piecewise analysis added to the technique of equation solver in the last step for large scale of wavelet. To demonstrate the effectiveness of this scheme, improved formulations are extended to the linear and nonlinear structural dynamic analysis. The validity and effectiveness of the developed method is verified with three examples. The results were compared with those from the numerical methods such as Duhamel integration, Runge-Kutta and Wilson-${\theta}$ method.

• Proceedings of the KIEE Conference
• /
• 2000.07b
• /
• pp.635-637
• /
• 2000

To analyze the transient state of an induction motor, there have been studies for using the magnetic equivalent circuit method (MECM) instead of the time differential finite-element method, MECM which analyzes magnetic equivalent circuits after converting each part of an electric machine into the magnetic circuit elements. has the merits of short calculation-time and comparatively accurate results. To analyze an electric machine with MECM, we have to replace stator and rotor with the magnetic elements and express the air gap, where electromechanical energy conversion takes place, with the permeance. So in this paper, to analyze an Induction Motor with MECM, we express the magnetic equivalent circuit as algebraic equations and then as the matrix for solving easily them.

• Structural Engineering and Mechanics
• /
• v.52 no.5
• /
• pp.997-1017
• /
• 2014

According to the structure characteristics of the non-uniform beam bridge, a practical model for calculating the vibration equation of the non-uniform beam bridge is given and the application scope of the model includes not only the beam bridge structure but also the non-uniform beam with added masses and elastic supports. Based on the Bernoulli-Euler beam theory, extending the application of the modal perturbation method and establishment of a semi-analytical method for solving the vibration equation of the non-uniform beam with added masses and elastic supports based is able to be made. In the modal subspace of the uniform beam with the elastic supports, the variable coefficient differential equation that describes the dynamic behavior of the non-uniform beam is converted to nonlinear algebraic equations. Extending the application of the modal perturbation method is suitable for solving the vibration equation of the simply supported and continuous non-uniform beam with its arbitrary added masses and elastic supports. The examples, that are analyzed, demonstrate the high precision and fast convergence speed of the method. Further study of the timesaving method for the dynamic characteristics of symmetrical beam and the symmetry of mode shape should be developed. Eventually, the effects of elastic supports and added masses on dynamic characteristics of the three-span non-uniform beam bridge are reported.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.4
• /
• pp.583-594
• /
• 1989

The strong nonliner dynamic behavior of mechanical systems in the presence of clearances are studied. The nonlinearity is induced from the assumed symmetric piecewise-linear characteristics for stiffness and damping by the contact and uncontact. Based on Stoker's assertion concering the reasoning beyond the occurrence of subharmonics, the nonlinear differential equation is converted to four nonlinear algebraic equations form the boundary conditions at the contact points. For a single contact per half exciting period, under the assumption of symmetric response, the steady-state solutions obtained are in agreement with those of numerical integration. Also a nondimen-sionalized formulation is made for the purpose of parametric studies.

• Proceedings of the KIEE Conference
• /
• 2000.11b
• /
• pp.330-332
• /
• 2000

To analyze the transient state of an induction motor, there have been studies for using the magnetic equivalent circuit method(MECD) instead of the time differential finite-element method. MECD which analyzes magnetic equivalent circuits after converting each part of an electric machine into the magnetic circuit elements, has the merits of short calculation-time and comparatively accurate results. To analyze an electric machine with MECM, we have to replace stator and rotor with the magnetic elements and express the air gap, where electromechanical energy conversion takes place, with the permeance. So in this study, to analyze an Induction Motor with HECM, we express the magnetic equivalent circuit as algebraic equations and then as the matrix for solving easily them. In particular, all relations are formed with matrixes to solve Mathematically them in the programming process later. As a result, this theory will be the basis on the static and dynamic analysis of an Induction Motor.