Article
Refine
Year of publication
Institute
- IfB - Institut für Bioengineering (390) (remove)
Has Fulltext
- no (390) (remove)
Language
- English (390) (remove)
Document Type
- Article (390) (remove)
Keywords
- CellDrum (3)
- Stiffness (2)
- bacterial cellulose (2)
- damage (2)
- locomotion (2)
- muscle fascicle behavior (2)
- prebiotic (2)
- ultrasound (2)
- ultrasound imaging (2)
- Achilles tendon (1)
Wearable EEG has gained popularity in recent years driven by promising uses outside of clinics and research. The ubiquitous application of continuous EEG requires unobtrusive form-factors that are easily acceptable by the end-users. In this progression, wearable EEG systems have been moving from full scalp to forehead and recently to the ear. The aim of this study is to demonstrate that emerging ear-EEG provides similar impedance and signal properties as established forehead EEG. EEG data using eyes-open and closed alpha paradigm were acquired from ten healthy subjects using generic earpieces fitted with three custom-made electrodes and a forehead electrode (at Fpx) after impedance analysis. Inter-subject variability in in-ear electrode impedance ranged from 20 kΩ to 25 kΩ at 10 Hz. Signal quality was comparable with an SNR of 6 for in-ear and 8 for forehead electrodes. Alpha attenuation was significant during the eyes-open condition in all in-ear electrodes, and it followed the structure of power spectral density plots of forehead electrodes, with the Pearson correlation coefficient of 0.92 between in-ear locations ELE (Left Ear Superior) and ERE (Right Ear Superior) and forehead locations, Fp1 and Fp2, respectively. The results indicate that in-ear EEG is an unobtrusive alternative in terms of impedance, signal properties and information content to established forehead EEG.
Exercise training effectively mitigates aging-induced health and fitness impairments. Traditional training recommendations for the elderly focus separately on relevant physiological fitness domains, such as balance, flexibility, strength and endurance. Thus, a more holistic and functional training framework is needed. The proposed agility training concept integratively tackles spatial orientation, stop and go, balance and strength. The presented protocol aims at introducing a two-armed, one-year randomized controlled trial, evaluating the effects of this concept on neuromuscular, cardiovascular, cognitive and psychosocial health outcomes in healthy older adults. Eighty-five participants were enrolled in this ongoing trial. Seventy-nine participants completed baseline testing and were block-randomized to the agility training group or the inactive control group. All participants undergo pre- and post-testing with interim assessment after six months. The intervention group currently receives supervised, group-based agility training twice a week over one year, with progressively demanding perceptual, cognitive and physical exercises. Knee extension strength, reactive balance, dual task gait speed and the Agility Challenge for the Elderly (ACE) serve as primary endpoints and neuromuscular, cognitive, cardiovascular, and psychosocial meassures serve as surrogate secondary outcomes. Our protocol promotes a comprehensive exercise training concept for older adults, that might facilitate stakeholders in health and exercise to stimulate relevant health outcomes without relying on excessively time-consuming physical activity recommendations.
Air-pulse corneal applanation signal curve parameters for the characterisation of keratoconus
(2011)
Altered gastrocnemius contractile behavior in former achilles tendon rupture patients during walking
(2022)
Achilles tendon rupture (ATR) remains associated with functional limitations years after injury. Architectural remodeling of the gastrocnemius medialis (GM) muscle is typically observed in the affected leg and may compensate force deficits caused by a longer tendon. Yet patients seem to retain functional limitations during—low-force—walking gait. To explore the potential limits imposed by the remodeled GM muscle-tendon unit (MTU) on walking gait, we examined the contractile behavior of muscle fascicles during the stance phase. In a cross-sectional design, we studied nine former patients (males; age: 45 ± 9 years; height: 180 ± 7 cm; weight: 83 ± 6 kg) with a history of complete unilateral ATR, approximately 4 years post-surgery. Using ultrasonography, GM tendon morphology, muscle architecture at rest, and fascicular behavior were assessed during walking at 1.5 m⋅s–1 on a treadmill. Walking patterns were recorded with a motion capture system. The unaffected leg served as control. Lower limbs kinematics were largely similar between legs during walking. Typical features of ATR-related MTU remodeling were observed during the stance sub-phases corresponding to series elastic element (SEE) lengthening (energy storage) and SEE shortening (energy release), with shorter GM fascicles (36 and 36%, respectively) and greater pennation angles (8° and 12°, respectively). However, relative to the optimal fascicle length for force production, fascicles operated at comparable length in both legs. Similarly, when expressed relative to optimal fascicle length, fascicle contraction velocity was not different between sides, except at the time-point of peak series elastic element (SEE) length, where it was 39 ± 49% lower in the affected leg. Concomitantly, fascicles rotation during contraction was greater in the affected leg during the whole stance-phase, and architectural gear ratios (AGR) was larger during SEE lengthening. Under the present testing conditions, former ATR patients had recovered a relatively symmetrical walking gait pattern. Differences in seen AGR seem to accommodate the profound changes in MTU architecture, limiting the required fascicle shortening velocity. Overall, the contractile behavior of the GM fascicles does not restrict length- or velocity-dependent force potentials during this locomotor task.
Altered neurovascular coupling as measured by optical imaging: a biomarker for Alzheimer’s disease
(2017)
Retinal vessels are similar to cerebral vessels in their structure and function. Moderately low oscillation frequencies of around 0.1 Hz have been reported as the driving force for paravascular drainage in gray matter in mice and are known as the frequencies of lymphatic vessels in humans. We aimed to elucidate whether retinal vessel oscillations are altered in Alzheimer's disease (AD) at the stage of dementia or mild cognitive impairment (MCI). Seventeen patients with mild-to-moderate dementia due to AD (ADD); 23 patients with MCI due to AD, and 18 cognitively healthy controls (HC) were examined using Dynamic Retinal Vessel Analyzer. Oscillatory temporal changes of retinal vessel diameters were evaluated using mathematical signal analysis. Especially at moderately low frequencies around 0.1 Hz, arterial oscillations in ADD and MCI significantly prevailed over HC oscillations and correlated with disease severity. The pronounced retinal arterial vasomotion at moderately low frequencies in the ADD and MCI groups would be compatible with the view of a compensatory upregulation of paravascular drainage in AD and strengthen the amyloid clearance hypothesis.
Edge-based and face-based smoothed finite element methods (ES-FEM and FS-FEM, respectively) are modified versions of the finite element method allowing to achieve more accurate results and to reduce sensitivity to mesh distortion, at least for linear elements. These properties make the two methods very attractive. However, their implementation in a standard finite element code is nontrivial because it requires heavy and extensive modifications to the code architecture. In this article, we present an element-based formulation of ES-FEM and FS-FEM methods allowing to implement the two methods in a standard finite element code with no modifications to its architecture. Moreover, the element-based formulation permits to easily manage any type of element, especially in 3D models where, to the best of the authors' knowledge, only tetrahedral elements are used in FS-FEM applications found in the literature. Shape functions for non-simplex 3D elements are proposed in order to apply FS-FEM to any standard finite element.
Replacement tissues, designed to fill in articular cartilage defects, should exhibit the same properties as the native material. The aim of this study is to foster the understanding of, firstly, the mechanical behavior of the material itself and, secondly, the influence of cultivation parameters on cell seeded implants as well as on cell migration into acellular implants. In this study, acellular cartilage replacement material is theoretically, numerically and experimentally investigated regarding its viscoelastic properties, where a phenomenological model for practical applications is developed. Furthermore, remodeling and cell migration are investigated.
Three-dimensional (3D) full-field measurements provide a comprehensive and accurate validation of finite element (FE) models. For the validation, the result of the model and measurements are compared based on two respective point-sets and this requires the point-sets to be registered in one coordinate system. Point-set registration is a non-convex optimization problem that has widely been solved by the ordinary iterative closest point algorithm. However, this approach necessitates a good initialization without which it easily returns a local optimum, i.e. an erroneous registration. The globally optimal iterative closest point (Go-ICP) algorithm has overcome this drawback and forms the basis for the presented open-source tool that can be used for the validation of FE models using 3D full-field measurements. The capability of the tool is demonstrated using an application example from the field of biomechanics. Methodological problems that arise in real-world data and the respective implemented solution approaches are discussed.
Analysis of the long-term effect of the MBST® nuclear magnetic resonance therapy on gonarthrosis
(2016)
Bacterial cellulose (BC) is a promising material for biomedical applications due to its unique properties such as high mechanical strength and biocompatibility. This article describes the microbiological synthesis, modification, and characterization of the obtained BC-nanocomposites originating from symbiotic consortium Medusomyces gisevii. Two BC-modifications have been obtained: BC-Ag and BC-calcium phosphate (BC-Ca3(PO4)2). Structure and physicochemical properties of the BC and its modifications were investigated by scanning electron microscopy (SEM), energy-dispersive X-ray spectroscopy (EDX), atomic force microscopy (AFM), and infrared Fourier spectroscopy as well as by measurements of mechanical and water holding/absorbing capacities. Topographic analysis of the surface revealed multicomponent thick fibrils (150–160 nm in diameter and about 15 µm in length) constituted by 50–60 nm nanofibrils weaved into a left-hand helix. Distinctive features of Ca-phosphate-modified BC samples were (a) the presence of 500–700 nm entanglements and (b) inclusions of Ca3(PO4)2 crystals. The samples impregnated with Ag nanoparticles exhibited numerous roundish inclusions, about 110 nm in diameter. The boundaries between the organic and inorganic phases were very distinct in both cases. The Ag-modified samples also showed a prominent waving pattern in the packing of nanofibrils. The obtained BC gel films possessed water-holding capacity of about 62.35 g/g. However, the dried (to a constant mass) BC-films later exhibited a low water absorption capacity (3.82 g/g). It was found that decellularized BC samples had 2.4 times larger Young’s modulus and 2.2 times greater tensile strength as compared to dehydrated native BC films. We presume that this was caused by molecular compaction of the BC structure.
Biocompatibility, flexibility and durability make polydimethylsiloxane (PDMS) membranes top candidates in biomedical applications. CellDrum technology uses large area, <10 µm thin membranes as mechanical stress sensors of thin cell layers. For this to be successful, the properties (thickness, temperature, dust, wrinkles, etc.) must be precisely controlled. The following parameters of membrane fabrication by means of the Floating-on-Water (FoW) method were investigated: (1) PDMS volume, (2) ambient temperature, (3) membrane deflection and (4) membrane mechanical compliance. Significant differences were found between all PDMS volumes and thicknesses tested (p < 0.01). They also differed from the calculated values. At room temperatures between 22 and 26 °C, significant differences in average thickness values were found, as well as a continuous decrease in thicknesses within a 4 °C temperature elevation. No correlation was found between the membrane thickness groups (between 3–4 µm) in terms of deflection and compliance. We successfully present a fabrication method for thin bio-functionalized membranes in conjunction with a four-step quality management system. The results highlight the importance of tight regulation of production parameters through quality control. The use of membranes described here could also become the basis for material testing on thin, viscous layers such as polymers, dyes and adhesives, which goes far beyond biological applications.
Background
Osteoporosis is associated with the risk of fractures near the hip. Age and comorbidities increase the perioperative risk. Due to the ageing population, fracture of the proximal femur also proves to be a socio-economic problem. Preventive surgical measures have hardly been used so far.
Methods
10 pairs of human femora from fresh cadavers were divided into control and low-volume femoroplasty groups and subjected to a Hayes fall-loading fracture test. The results of the respective localization and classification of the fracture site, the Singh index determined by computed tomography (CT) examination and the parameters in terms of fracture force, work to fracture and stiffness were evaluated statistically and with the finite element method. In addition, a finite element parametric study with different position angles and variants of the tubular geometry of the femoroplasty was performed.
Findings
Compared to the control group, the work to fracture could be increased by 33.2%. The fracture force increased by 19.9%. The used technique and instrumentation proved to be standardized and reproducible with an average poly(methyl methacrylate) volume of 10.5 ml. The parametric study showed the best results for the selected angle and geometry.
Interpretation
The cadaver studies demonstrated the biomechanical efficacy of the low-volume tubular femoroplasty. The numerical calculations confirmed the optimal choice of positioning as well as the inner and outer diameter of the tube in this setting. The standardized minimally invasive technique with the instruments developed for it could be used in further comparative studies to confirm the measured biomechanical results.
This study aims to quantify the kinematics, kinetics and muscular activity of all-out handcycling exercise and examine their alterations during the course of a 15-s sprint test. Twelve able-bodied competitive triathletes performed a 15-s all-out sprint test in a recumbent racing handcycle that was attached to an ergometer. During the sprint test, tangential crank kinetics, 3D joint kinematics and muscular activity of 10 muscles of the upper extremity and trunk were examined using a power metre, motion capturing and surface electromyography (sEMG), respectively. Parameters were compared between revolution one (R1), revolution two (R2), the average of revolution 3 to 13 (R3) and the average of the remaining revolutions (R4). Shoulder abduction and internal-rotation increased, whereas maximal shoulder retroversion decreased during the sprint. Except for the wrist angles, angular velocity increased for every joint of the upper extremity. Several muscles demonstrated an increase in muscular activation, an earlier onset of muscular activation in crank cycle and an increased range of activation. During the course of a 15-s all-out sprint test in handcycling, the shoulder muscles and the muscles associated to the push phase demonstrate indications for short-duration fatigue. These findings are helpful to prevent injuries and improve performance in all-out handcycling.
Purpose
This study aims to investigate the biomechanics of handcycling during a continuous load trial (CLT) to assess the mechanisms underlying fatigue in upper body exercise.
Methods
Twelve able-bodied triathletes performed a 30-min CLT at a power output corresponding to lactate threshold in a racing recumbent handcycle mounted on a stationary ergometer. During the CLT, ratings of perceived exertion (RPE), tangential crank kinetics, 3D joint kinematics, and muscular activity of ten muscles of the upper extremity and trunk were examined using motion capturing and surface electromyography.
Results
During the CLT, spontaneously chosen cadence and RPE increased, whereas crank torque decreased. Rotational work was higher during the pull phase. Peripheral RPE was higher compared to central RPE. Joint range of motion decreased for elbow-flexion and radial-duction. Integrated EMG (iEMG) increased in the forearm flexors, forearm extensors, and M. deltoideus (Pars spinalis). An earlier onset of activation was found for M. deltoideus (Pars clavicularis), M. pectoralis major, M. rectus abdominis, M. biceps brachii, and the forearm flexors.
Conclusion
Fatigue-related alterations seem to apply analogously in handcycling and cycling. The most distal muscles are responsible for force transmission on the cranks and might thus suffer most from neuromuscular fatigue. The findings indicate that peripheral fatigue (at similar lactate values) is higher in handcycling compared to leg cycling, at least for inexperienced participants. An increase in cadence might delay peripheral fatigue by a reduced vascular occlusion. We assume that the gap between peripheral and central fatigue can be reduced by sport-specific endurance training.
It was generally believed that coal sources are not favorable as live-in habitats for microorganisms due to their recalcitrant chemical nature and negligible decomposition. However, accumulating evidence has revealed the presence of diverse microbial groups in coal environments and their significant metabolic role in coal biogeochemical dynamics and ecosystem functioning. The high oxygen content, organic fractions, and lignin-like structures of lower-rank coals may provide effective means for microbial attack, still representing a greatly unexplored frontier in microbiology. Coal degradation/conversion technology by native bacterial and fungal species has great potential in agricultural development, chemical industry production, and environmental rehabilitation. Furthermore, native microalgal species can offer a sustainable energy source and an excellent bioremediation strategy applicable to coal spill/seam waters. Additionally, the measures of the fate of the microbial community would serve as an indicator of restoration progress on post-coal-mining sites. This review puts forward a comprehensive vision of coal biodegradation and bioprocessing by microorganisms native to coal environments for determining their biotechnological potential and possible applications.
Using the OpenSim software and verified anatomical data, a computer model for the calculation of biomechanical parameters is developed and used to determine the effect of a reattachment of the Supraspinatus muscle with a medial displacement of the muscle attachment point, which may be necessary for a rupture of the supraspinatus tendon. The results include the influence of the operation on basic biomechanical parameters such as the lever arm, as well as the calculated the muscle activations for the supraspinatus and deltoid. In addition, the influence on joint stability is examined by an analysis of the joint reaction force. The study provides a detailed description of the used model, as well as medical findings to a reattachment of the supraspinatus.
Mit der Software OpenSim und überprüften anatomischen Daten wird ein Computermodell zur Berechnung von biomechanischen Parametern entwickelt und genutzt, um den Effekt einer Refixierung des Supraspinatusmuskels mit einer medialen Verschiebung des Muskelansatzpunktes zu ermitteln, wie sie unter anderem nach einem Riss der Supraspinatussehne notwendig sein kann. Die Ergebnisse umfassen hierbei den Einfluss der Operation auf grundlegende biomechanische Parameter wie den Hebelarm sowie die berechneten Muskelaktivierungen für den Supraspinatus und Deltoideus. Zusätzlich wird der Einfluss auf die Gelenkstabilität betrachtet und durch eine Analyse der Gelenkreaktionskraft untersucht. Die Studie bietet eine detaillierte Beschreibung des genutzten Modells, sowie medizinische Erkenntnisse zu einer Refixierung des Supraspinatus.
To better understand what kinds of sports and exercise could be beneficial for the intervertebral disc (IVD), we performed a review to synthesise the literature on IVD adaptation with loading and exercise. The state of the literature did not permit a systematic review; therefore, we performed a narrative review. The majority of the available data come from cell or whole-disc loading models and animal exercise models. However, some studies have examined the impact of specific sports on IVD degeneration in humans and acute exercise on disc size. Based on the data available in the literature, loading types that are likely beneficial to the IVD are dynamic, axial, at slow to moderate movement speeds, and of a magnitude experienced in walking and jogging. Static loading, torsional loading, flexion with compression, rapid loading, high-impact loading and explosive tasks are likely detrimental for the IVD. Reduced physical activity and disuse appear to be detrimental for the IVD. We also consider the impact of genetics and the likelihood of a ‘critical period’ for the effect of exercise in IVD development. The current review summarises the literature to increase awareness amongst exercise, rehabilitation and ergonomic professionals regarding IVD health and provides recommendations on future directions in research.
Can vascular function be assessed by the interpretation of retinal vascular diameter changes?
(2011)
Aneurysmal subarachnoid hemorrhage (aSAH) is associated with early and delayed brain injury due to several underlying and interrelated processes, which include inflammation, oxidative stress, endothelial, and neuronal apoptosis. Treatment with melatonin, a cytoprotective neurohormone with anti-inflammatory, anti-oxidant and anti-apoptotic effects, has been shown to attenuate early brain injury (EBI) and to prevent delayed cerebral vasospasm in experimental aSAH models. Less is known about the role of endogenous melatonin for aSAH outcome and how its production is altered by the pathophysiological cascades initiated during EBI. In the present observational study, we analyzed changes in melatonin levels during the first three weeks after aSAH.
This study evaluates neuromechanical control and muscle-tendon interaction during energy storage and dissipation tasks in hypergravity. During parabolic flights, while 17 subjects performed drop jumps (DJs) and drop landings (DLs), electromyography (EMG) of the lower limb muscles was combined with in vivo fascicle dynamics of the gastrocnemius medialis, two-dimensional (2D) kinematics, and kinetics to measure and analyze changes in energy management. Comparisons were made between movement modalities executed in hypergravity (1.8 G) and gravity on ground (1 G). In 1.8 G, ankle dorsiflexion, knee joint flexion, and vertical center of mass (COM) displacement are lower in DJs than in DLs; within each movement modality, joint flexion amplitudes and COM displacement demonstrate higher values in 1.8 G than in 1 G. Concomitantly, negative peak ankle joint power, vertical ground reaction forces, and leg stiffness are similar between both movement modalities (1.8 G). In DJs, EMG activity in 1.8 G is lower during the COM deceleration phase than in 1 G, thus impairing quasi-isometric fascicle behavior. In DLs, EMG activity before and during the COM deceleration phase is higher, and fascicles are stretched less in 1.8 G than in 1 G. Compared with the situation in 1 G, highly task-specific neuromuscular activity is diminished in 1.8 G, resulting in fascicle lengthening in both movement modalities. Specifically, in DJs, a high magnitude of neuromuscular activity is impaired, resulting in altered energy storage. In contrast, in DLs, linear stiffening of the system due to higher neuromuscular activity combined with lower fascicle stretch enhances the buffering function of the tendon, and thus the capacity to safely dissipate energy.
Objective
Hemodialysis patients show an approximately threefold higher prevalence of cognitive impairment compared to the age-matched general population. Impaired microcirculatory function is one of the assumed causes. Dynamic retinal vessel analysis is a quantitative method for measuring neurovascular coupling and microvascular endothelial function. We hypothesize that cognitive impairment is associated with altered microcirculation of retinal vessels.
Methods
152 chronic hemodialysis patients underwent cognitive testing using the Montreal Cognitive Assessment. Retinal microcirculation was assessed by Dynamic Retinal Vessel Analysis, which carries out an examination recording retinal vessels' reaction to a flicker light stimulus under standardized conditions.
Results
In unadjusted as well as in adjusted linear regression analyses a significant association between the visuospatial executive function domain score of the Montreal Cognitive Assessment and the maximum arteriolar dilation as response of retinal arterioles to the flicker light stimulation was obtained.
Conclusion
This is the first study determining retinal microvascular function as surrogate for cerebral microvascular function and cognition in hemodialysis patients. The relationship between impairment in executive function and reduced arteriolar reaction to flicker light stimulation supports the involvement of cerebral small vessel disease as contributing factor for the development of cognitive impairment in this patient population and might be a target for noninvasive disease monitoring and therapeutic intervention.
Combining physiological relevance and throughput for in vitro cardiac contractility measurement
(2020)
Despite increasing acceptance of human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) in safety pharmacology, controversy remains about the physiological relevance of existing in vitro models for their mechanical testing. We hypothesize that existing signs of immaturity of the cell models result from an improper mechanical environment. We cultured hiPSC-CMs in a 96-well format on hyperelastic silicone membranes imitating their native mechanical environment, resulting in physiological responses to compound stimuli.We validated cell responses on the FLEXcyte 96, with a set of reference compounds covering a broad range of cellular targets, including ion channel modulators, adrenergic receptor modulators and kinase inhibitors. Acute (10 - 30 min) and chronic (up to 7 days) effects were investigated. Furthermore, the measurements were complemented with electromechanical models based on electrophysiological recordings of the used cell types.hiPSC-CMs were cultured on freely-swinging, ultra-thin and hyperelastic silicone membranes. The weight of the cell culture medium deflects the membranes downwards. Rhythmic contraction of the hiPSC-CMs resulted in dynamic deflection changes which were quantified by capacitive distance sensing. The cells were cultured for 7 days prior to compound addition. Acute measurements were conducted 10-30 minutes after compound addition in standard culture medium. For chronic treatment, compound-containing medium was replaced daily for up to 7 days. Electrophysiological properties of the employed cell types were recorded by automated patch-clamp (Patchliner) and the results were integrated into the electromechanical model of the system.Calcium channel agonist S Bay K8644 and beta-adrenergic stimulator isoproterenol induced significant positive inotropic responses without additional external stimulation. Kinase inhibitors displayed cardiotoxic effects on a functional level at low concentrations. The system-integrated analysis detected alterations in beating shape as well as frequency and arrhythmic events and we provide a quantitative measure of these.
Comparison of different training algorithms for the leg extension training with an industrial robot
(2018)
In the past, different training scenarios have been developed and implemented on robotic research platforms, but no systematic analysis and comparison have been done so far. This paper deals with the comparison of an isokinematic (motion with constant velocity) and an isotonic (motion against constant weight) training algorithm. Both algorithms are designed for a robotic research platform consisting of a 3D force plate and a high payload industrial robot, which allows leg extension training with arbitrary six-dimensional motion trajectories. In the isokinematic as well as the isotonic training algorithm, individual paths are defined i n C artesian s pace by sufficient s upport p oses. I n t he i sotonic t raining s cenario, the trajectory is adapted to the measured force as the robot should only move along the trajectory as long as the force applied by the user exceeds a minimum threshold. In the isotonic training scenario however, the robot’s acceleration is a function of the force applied by the user. To validate these findings, a simulative experiment with a simple linear trajectory is performed. For this purpose, the same force path is applied in both training scenarios. The results illustrate that the algorithms differ in the force dependent trajectory adaption.
This paper presents a novel numerical procedure for computing limit and shakedown loads of structures using a node-based smoothed FEM in combination with a primal–dual algorithm. An associated primal–dual form based on the von Mises yield criterion is adopted. The primal-dual algorithm together with a Newton-like iteration are then used to solve this associated primal–dual form to determine simultaneously both approximate upper and quasi-lower bounds of the plastic collapse limit and the shakedown limit. The present formulation uses only linear approximations and its implementation into finite element programs is quite simple. Several numerical examples are given to show the reliability, accuracy, and generality of the present formulation compared with other available methods.