Adhesively-bonded Patch Repair with Composites

Adhesively-bonded composite patch repairs over cracked or corrosion-damaged metallic aircraft structures have shown great promise for extending life of ageing structures. This study presents the numerical investigation into the interface behaviour of adhesively-bonded cracked aluminum alloy substrate patched with fibre-reinforced composite material. The adhesive is modelled as an elasto-plastic bilinear material to characterise the debond behaviour, while the defective substrate is regarded as linear elastic continuum. Two typical patch shapes were selected based on information available in the literature. Geometric and material nonlinear analyses for square and octagonal patches were performed to capture peel and shear stresses developed between the substrate and the patch to examine the possibility of interface delamination/debonding. Parametric studies on adhesive thickness and patch thickness were carried out to predict their infuence on damage tolerance of repaired structures.

Effect of interaction of macrocracks on the stress intensity factor in a beam

Beams with a central edge crack, as well as other noncentral vertical and inclined edge cracks distributed symmetrically, subjected to three-point as well as four-point bending, are analysed using the finite element technique. Values of stress intensity factor K1 at the central crack tip for a crack-to-beam depth ratio Image equal to 0.5, are calculated for various cracked-beam configurations, using the compliance calibration technique as well as method of strain energy release rate. These are compared with the value of K1 for the case of a beam with central edge crack alone. Results of the present parametric study are used to specify the range of values pertaining to basic parameters such as crack-to-beam depth ratios, geometry and position with respect to central edge crack, of other macrocracks for which interaction exists. Accordingly, the macrocracks are classified as either interacting or noninteracting types. Hence for noninteracting types of cracks, analytical expressions available for the determination of K1 in the case of beam with a central edge crack alone, are applicable.

Stress intensity factor determination of radially cracked circular rings subjected to tension using photoelastic technique

A parametric study was carried out to determine the Stress Intensity Factor (SIF) in a cracked circular ring by using the photoelastic technique. The stress intensity factors for mode I deformation were determined by subjecting the specimens to the tensile loading from inner boundary and through the holes. The results of Non-Dimensional Stress Intensity Factor (NDSIF) variation with non-dimensional crack length for both methods of loading are compared with each other and with published results.

Generalized Burgers equations and Euler–Painlevé transcendents. III

It was proposed earlier [P. L. Sachdev, K. R. C. Nair, and V. G. Tikekar, J. Math. Phys. 27, 1506 (1986); P. L. Sachdev and K. R. C. Nair, ibid. 28, 977 (1987)] that the Euler–Painlevé equations  y(d2y/dη2)+a(dy/dη)2 +f(η)y(dy/dη)+g(η)y2+b(dy/dη) +c=0 represent generalized Burgers equations (GBE’s) in the same way as Painlevé equations represent the Korteweg–de Vries type of equations. The earlier studies were carried out in the context of GBE’s with damping and those with spherical and cylindrical symmetry. In the present paper, GBE’s with variable coefficients of viscosity and those with inhomogeneous terms are considered for their possible connection to Euler–Painlevé equations. It is found that the Euler–Painlevé equation, which represents the GBE ut+uβux=(δ/2)g(t)uxx, g(t)=(1+t)n, β>0, has solutions, which either decay or oscillate at η=±∞, only when −1parametric value β=βn seems to bifurcate the families of solutions, which remain bounded at η=±∞. Other GBE’s considered here are also found to be reducible to Euler–Painlevé equations.

A dragonfly inspired flapping wing actuated by electro active polymers

An energy-based variational approach is used for structural dynamic modeling of the IPMC (Ionic Polymer Metal Composites) flapping wing. Dynamic characteristics of the wing are analyzed using numerical simulations. Starting with the initial design, critical parameters which have influence on the performance of the wing are identified through parametric studies. An optimization study is performed to obtain improved flapping actuation of the IPMC wing. It is shown that the optimization algorithm leads to a flapping wing with dimensions similar to the dragonfly Aeshna Multicolor wing. An unsteady aerodynamic model based on modified strip theory is used to obtain the aerodynamic forces. It is found that the IPMC wing generates sufficient lift to support its own weight and carry a small payload. It is therefore a potential candidate for flapping wing of micro air vehicles.

Methods in general model localization

The aim of this study was to evaluate and test methods which could improve local estimates of a general model fitted to a large area. In the first three studies, the intention was to divide the study area into sub-areas that were as homogeneous as possible according to the residuals of the general model, and in the fourth study, the localization was based on the local neighbourhood. According to spatial autocorrelation (SA), points closer together in space are more likely to be similar than those that are farther apart. Local indicators of SA (LISAs) test the similarity of data clusters. A LISA was calculated for every observation in the dataset, and together with the spatial position and residual of the global model, the data were segmented using two different methods: classification and regression trees (CART) and the multiresolution segmentation algorithm (MS) of the eCognition software. The general model was then re-fitted (localized) to the formed sub-areas. In kriging, the SA is modelled with a variogram, and the spatial correlation is a function of the distance (and direction) between the observation and the point of calculation. A general trend is corrected with the residual information of the neighbourhood, whose size is controlled by the number of the nearest neighbours. Nearness is measured as Euclidian distance. With all methods, the root mean square errors (RMSEs) were lower, but with the methods that segmented the study area, the deviance in single localized RMSEs was wide. Therefore, an element capable of controlling the division or localization should be included in the segmentation-localization process. Kriging, on the other hand, provided stable estimates when the number of neighbours was sufficient (over 30), thus offering the best potential for further studies. Even CART could be combined with kriging or non-parametric methods, such as most similar neighbours (MSN).

Laminar compressible boundary layer flow at a three-dimensional stagnation point with vectored mass transfer

The effect of surface mass transfer velocities having normal, principal and transverse direction components (�vectored� suction and injection) on the steady, laminar, compressible boundary layer at a three-dimensional stagnation point has been investigated both for nodal and saddle points of attachment. The similarity solutions of the boundary layer equations were obtained numerically by the method of parametric differentiation. The principal and transverse direction surface mass transfer velocities significantly affect the skin friction (both in the principal and transverse directions) and the heat transfer. Also the inadequacy of assuming a linear viscosity-temperature relation at low-wall temperatures is shown.

Lahopuumäärän ennustaminen ja kartoitus lentokonelaserkeilauksen avulla

Lahopuun määrästä ja sijoittumisesta ollaan kiinnostuneita paitsi elinympäristöjen monimuotoisuuden, myös ilmakehän hiilen varastoinnin kannalta. Tutkimuksen tavoitteena oli kehittää aluepohjainen laserkeilausdataa hyödyntävä malli lahopuukohteiden paikantamiseksi ja lahopuun määrän estimoimiseksi. Samalla tutkittiin mallin selityskyvyn muuttumista mallinnettavan ruudun kokoa suurennettaessa. Tutkimusalue sijaitsi Itä-Suomessa Sonkajärvellä ja koostui pääasiassa nuorista hoidetuista talousmetsistä. Tutkimuksessa käytettiin harvapulssista laserkeilausdataa sekä kaistoittain mitattua maastodataa kuolleesta puuaineksesta. Aineisto jaettiin siten, että neljäsosa datasta oli käytössä mallinnusta varten ja loput varattiin valmiiden mallien testaamiseen. Lahopuun mallintamisessa käytettiin sekä parametrista että ei-parametrista mallinnusmenetelmää. Logistisen regression avulla erikokoisille (0,04, 0,20, 0,32, 0,52 ja 1,00 ha) ruuduille ennustettiin todennäköisyys lahopuun esiintymiselle. Muodostettujen mallien selittävät muuttujat valittiin 80 laserpiirteen ja näiden muunnoksien joukosta. Mallien selittävät muuttujat valittiin kolmessa vaiheessa. Aluksi muuttujia tarkasteltiin visuaalisesti kuvaamalla ne lahopuumäärän suhteen. Ensimmäisessä vaiheessa sopivimmiksi arvioitujen muuttujien selityskykyä testattiin mallinnuksen toisessa vaiheessa yhden muuttujan mallien avulla. Lopullisessa usean muuttujan mallissa selittävien muuttujien kriteerinä oli tilastollinen merkitsevyys 5 % riskitasolla. 0,20 hehtaarin ruutukoolle luotu malli parametrisoitiin muun kokoisille ruuduille. Logistisella regressiolla toteutetun parametrisen mallintamisen lisäksi, 0,04 ja 1,0 hehtaarin ruutukokojen aineistot luokiteltiin ei-parametrisen CART-mallinnuksen (Classification and Regression Trees) avulla. CARTmenetelmällä etsittiin aineistosta vaikeasti havaittavia epälineaarisia riippuvuuksia laserpiirteiden ja lahopuumäärän välillä. CART-luokittelu tehtiin sekä lahopuustoisuuden että lahopuutilavuuden suhteen. CART-luokituksella päästiin logistista regressiota parempiin tuloksiin ruutujen luokituksessa lahopuustoisuuden suhteen. Logistisella mallilla tehty luokitus parani ruutukoon suurentuessa 0,04 ha:sta(kappa 0,19) 0,32 ha:iin asti (kappa 0,38). 0,52 ha:n ruutukoolla luokituksen kappa-arvo kääntyi laskuun (kappa 0,32) ja laski edelleen hehtaarin ruutukokoon saakka (kappa 0,26). CART-luokitus parani ruutukoon kasvaessa. Luokitustulokset olivat logistista mallinnusta parempia sekä 0,04 ha:n (kappa 0,24) että 1,0 ha:n (kappa 0,52) ruutukoolla. CART-malleilla määritettyjen ruutukohtaisten lahopuutilavuuksien suhteellinen RMSE pieneni ruutukoon kasvaessa. 0,04 hehtaarin ruutukoolla koko aineiston lahopuumäärän suhteellinen RMSE oli 197,1 %, kun hehtaarin ruutukoolla vastaava luku oli 120,3 %. Tämän tutkimuksen tulosten perusteella voidaan todeta, että maastossa mitatun lahopuumäärän ja tutkimuksessa käytettyjen laserpiirteiden yhteys on pienellä ruutukoolla hyvin heikko, mutta vahvistuu hieman ruutukoon kasvaessa. Kun mallinnuksessa käytetty ruutukoko kasvaa, pienialaisten lahopuukeskittymien havaitseminen kuitenkin vaikeutuu. Tutkimuksessa kohteen lahopuustoisuus pystyttiin kartoittamaan kohtuullisesti suurella ruutukoolla, mutta pienialaisten kohteiden kartoittaminen ei onnistunut käytetyillä menetelmillä. Pienialaisten kohteiden paikantaminen laserkeilauksen avulla edellyttää jatkotutkimusta erityisesti tiheäpulssisen laserdatan käytöstä lahopuuinventoinneissa.

Ratkaiseeko raha? : Tulospalkkiojärjestelmään ja työympäristöön liittyvien kokemusten yhteys esimiesten työtyytyväisyyteen ja työmotivaatioon

The aim of this study was to investigate the relationship between merit pay system and work environment and foremen´s work satisfaction and work motivation. There has been a lot of investigation on rewarding. Less research has been done on previous surveys among the merit pay systems and motivation investigations. According to former surveys, rewarding systems cannot be released from its context. Therefore this survey expanded to deal with work environment. It was also essential to investigate different dimensions of extrinsic and intrinsic motivation and equity of rewarding. Investigation or work motivation and work satisfaction was challenging because both of these concepts have been investigated under quite traditional frame of reference of work motivation theories. In some surveys, the concepts have not been even separated or they have been used even as synonyms. The data were collected with the 193 foremen working in the profit centers of the different chains of the company in the field of retail trade. The questions were: Are the experiences of merit pay system and work environment related to foremen´s work satisfaction and work motivation? Are the backround variables related to foremen´s work satisfaction and work motivation? The data collection was carried out by an electronic inquiry during May 2010. 137 replied from foremen working under merit pay system. The research material was analyzed with PASW-software. Various analyzing methods were used: factor analyses, regression analyses and group of different parametric and non-parametric analyses. In contrast to theoretical framework in the factor analyses work satisfaction and work motivation clustered into the same dimension. As a main result the atmosphere, possibilities to influence and the atmosphere of leading were strongly positively related to foremen´s work satisfaction and work motivation. According to regression analyses these factors were able to explain 55 % of the foremen´s work satisfaction and work motivation. The best explanatory variable was atmosphere. Instead, the backround variables (age, sex, working years, group of profession, education) were not associated with work satisfaction and work motivation.

Heat transfer correlations for PCM-based heat sinks with plate fins

Characterization of melting process in a Phase Change Material (PCM)-based heat sink with plate fin type thermal conductivity enhancers (TCEs) is numerically studied in this paper. Detailed parametric investigations are performed to find the effect of aspect ratio of enclosure and the applied heat flux on the thermal performance of the heat sinks. Various non-dimensional numbers, such as Nusselt number (Nu), Rayleigh number (Ra), Stefan number (Ste) and Fourier number (Fo) based on a characteristic length scale, are identified as important parameters. The half fin thickness and the fin height are varied to obtain a wide range of aspect ratios of an enclosure. It is found that a single correlation of Nu with Ra is not applicable for all aspect ratios of enclosure with melt convection taken into account. To find appropriate length scales, enclosures with different aspect ratios are divided into three categories, viz. (a) shallow enclosure, (b) rectangular enclosure and (c) tall enclosure. Accordingly, an appropriate characteristic length scale is identified for each type of enclosure and correlation of Nu with Ra based on that characteristic length scale is developed. (C) 2010 Elsevier Ltd. All rights reserved.

Seismic Stability Analysis of a Himalayan Rock Slope

The seismic slope stability analysis of the right abutment of a railway bridge proposed at about 350 m above the ground level, crossing a river and connecting two huge hillocks in the Himalayas, India, is presented in this paper. The rock slopes are composed of highly jointed rock mass and the joint spacing and orientation are varying at different locations. Seismic slope stability analysis of the slope under consideration is carried out using both pseudo-static approach and time response approach as the site is located in seismic zone V as per the earth quake zonation maps of India. Stability of the slope is studied numerically using program FLAC. The results obtained from the pseudo-static analysis are presented in the form of Factor of Safety (FOS) and the results obtained from the time response analysis of the slope are presented in terms of horizontal and vertical displacements along the slope. The results obtained from both the analyses confirmed the global stability of the slope as the FOS in case of pseudo-static analysis is above 1.0 and the displacements observed in case of time response analysis are within the permissible limits. This paper also presents the results obtained from the parametric analysis performed in the case of time response analysis in order to understand the effect of individual parameters on the overall stability of the slope.

How to deal with sparse macroseismic data

This study discusses the scope of historical earthquake analysis in low-seismicity regions. Examples of non-damaging earthquake reports are given from the Eastern Baltic (Fennoscandian) Shield in north-eastern Europe from the 16th to the 19th centuries. The information available for past earthquakes in the region is typically sparse and cannot be increased through a careful search of the archives. This study applies recommended rigorous methodologies of historical seismology developed using ample data to the sparse reports from the Eastern Baltic Shield. Attention is paid to the context of reporting, the identity and role of the authors, the circumstances of the reporting, and the opportunity to verify the available information by collating the sources. We evaluate the reliability of oral earthquake recollections and develop criteria for cases when a historical earthquake is attested to by a single source. We propose parametric earthquake scenarios as a way to deal with sparse macroseismic reports and as an improvement to existing databases.

MPD arcjet cathode surface current density

The radial current density distribution on the cathode longitudinal surface of magnetoplasmadynamic arcjets for axisymmetric geometries has been obtained by simultaneous solution of the electromagnetic equations for a given uniform gas dynamic field. The problem formulation permits a parametric study of the effects of the Hall parameter and the magnetic Reynolds number. The solution for the current density distribution displays current concentrations at two locations, that is, at the upstream and downstream ends of the cathode. This result is in conformity with known experimental data. The parameters responsible for these current concentrations are identified. It is shown that the effect of the magnetic Reynolds number on the current density distribution is different depending on whether or not the Hall effect is included. This result is also found to be consistent with experimental data.

A Smooth Finite Element Method Based on Reproducing Kernel DMS-Splines

The element-based piecewise smooth functional approximation in the conventional finite element method (FEM) results in discontinuous first and higher order derivatives across element boundaries Despite the significant advantages of the FEM in modelling complicated geometries, a motivation in developing mesh-free methods has been the ease with which higher order globally smooth shape functions can be derived via the reproduction of polynomials There is thus a case for combining these advantages in a so-called hybrid scheme or a `smooth FEM' that, whilst retaining the popular mesh-based discretization, obtains shape functions with uniform C-p (p >= 1) continuity One such recent attempt, a NURBS based parametric bridging method (Shaw et al 2008b), uses polynomial reproducing, tensor-product non-uniform rational B-splines (NURBS) over a typical FE mesh and relies upon a (possibly piecewise) bijective geometric map between the physical domain and a rectangular (cuboidal) parametric domain The present work aims at a significant extension and improvement of this concept by replacing NURBS with DMS-splines (say, of degree n > 0) that are defined over triangles and provide Cn-1 continuity across the triangle edges This relieves the need for a geometric map that could precipitate ill-conditioning of the discretized equations Delaunay triangulation is used to discretize the physical domain and shape functions are constructed via the polynomial reproduction condition, which quite remarkably relieves the solution of its sensitive dependence on the selected knotsets Derivatives of shape functions are also constructed based on the principle of reproduction of derivatives of polynomials (Shaw and Roy 2008a) Within the present scheme, the triangles also serve as background integration cells in weak formulations thereby overcoming non-conformability issues Numerical examples involving the evaluation of derivatives of targeted functions up to the fourth order and applications of the method to a few boundary value problems of general interest in solid mechanics over (non-simply connected) bounded domains in 2D are presented towards the end of the paper

A geodesic-based triangulation technique for damage location in metallic and composite plates

Lamb-wave-based damage detection methods using the triangulation technique are not suitable for handling structures with complex shapes and discontinuities as the parametric/analytical representation of these structures is very difficult. The geodesic concept is used along with the triangulation technique to overcome the above problem. The present work is based on the fundamental fact that a wave takes the minimum energy path to travel between two points on any multiply connected surface and this reduces to the shortest distance path or geodesic. The geodesics are computed on the meshed surface of the structure using the fast marching method. The wave response matrix of the given sensor configuration for the healthy and the damaged structure is obtained experimentally. The healthy and damage response matrices are compared and their difference gives the time information about the reflection of waves from the damage. A wavelet transform is used to extract the arrival time information of the wave scattered by the damage from the acquired Lamb wave signals. The computed geodesics and time information are used in the ellipse algorithm of triangulation formulation to locate the loci of possible damage location points for each actuator-sensor pair. The results obtained for all actuator-sensor pairs are combined and the intersection of multiple loci gives the damage location result. Experiments were conducted in aluminum and composite plate specimens to validate this method.