Dynamics Modeling and Parameter Identification for a Coupled-Drive Dual-Arm Nursing Robot

A dual-arm nursing robot can gently lift patients and transfer them between a bed and a wheelchair.With its lightweight design,high load-bearing capacity,and smooth surface,the coupled-drive joint is particularly well suited for these robots.However,the coupled nature of the joint disrupts the direct linear relationship between the input and output torques,posing challenges for dynamic modeling and practical applications.This study investigated the transmission mechanism of this joint and employed the Lagrangian method to construct a dynamic model of its internal dynamics.Building on this foundation,the Newton-Euler method was used to develop a dynamic model for the entire robotic arm.A continuously differentiable friction model was incorporated to reduce the vibrations caused by speed transitions to zero.An experimental method was designed to compensate for gravity,inertia,and modeling errors to identify the parameters of the friction model.This method establishes a mapping relationship between the friction force and motor current.In addition,a Fourier series-based excitation trajectory was developed to facilitate the identification of the dynamic model parameters of the robotic arm.Trajectory tracking experiments were conducted during the experimental validation phase,demonstrating the high accuracy of the dynamic model and the parameter identification method for the robotic arm.This study presents a dynamic modeling and parameter identification method for coupled-drive joint robotic arms,thereby establishing a foundation for motion control in humanoid nursing robots.

Parameter calibration of the tensile-shear interactive damage constitutive model for sandstone failure

The tensile-shear interactive damage(TSID)model is a novel and powerful constitutive model for rock-like materials.This study proposes a methodology to calibrate the TSID model parameters to simulate sandstone.The basic parameters of sandstone are determined through a series of static and dynamic tests,including uniaxial compression,Brazilian disc,triaxial compression under varying confining pressures,hydrostatic compression,and dynamic compression and tensile tests with a split Hopkinson pressure bar.Based on the sandstone test results from this study and previous research,a step-by-step procedure for parameter calibration is outlined,which accounts for the categories of the strength surface,equation of state(EOS),strain rate effect,and damage.The calibrated parameters are verified through numerical tests that correspond to the experimental loading conditions.Consistency between numerical results and experimental data indicates the precision and reliability of the calibrated parameters.The methodology presented in this study is scientifically sound,straightforward,and essential for improving the TSID model.Furthermore,it has the potential to contribute to other rock constitutive models,particularly new user-defined models.

Construction of an Arrhenius constitutive model for Mg-Y-Nd-Zr-Gd rare earth magnesium alloy based on the Zener-Hollomon parameter and objective evaluation of its accuracy in the twinning-rich intervals

According to a high-temperature compression test of rare earth magnesium alloy(WE43),a strain-compensated constitutive model of the Arrhenius equation based on Zener-Hollomon parameters was established,and the rheological behaviors were predicted.The model exhibited relatively serious prediction distortion in the low-temperature and high-strain rate parameter interval,and its accuracy was still unsatisfactory even after modification by a correction operator considering the coupling of temperature and strain rate.The microstructure characterization and statistical analysis showed that a large number of twinning occurred in the parameter intervals with prediction deviation.The occurrence of twinning complicated the local internal stress distribution by drastically changing the crystal orientation and led to significant fluctuations in the macroscopic strain-stress and hardening curves relative to the rheological processes dominated by the dislocation and softening mechanisms,making the logarithm of the strain rate and stress deviate from the linear relationship.This twinning phenomenon was greatly influenced by the temperature and strain rate.Herein,the influence mechanism on twinning behavior was analyzed from the perspective of the interaction of dislocation and twinning.

Numerical Simulation and Parameter Estimation of Fractional-Order Dynamic Epidemic Model for COVID-19

The outbreak of COVID-19 in 2019 resulted in numerous infections and deaths. In order to better study the transmission of COVID-19, this article adopts an improved fractional-order SIR model. Firstly, the properties of the model are studied, including the feasible domain and bounded solutions of the system. Secondly, the stability of the system is discussed, among other things. Then, the GMMP method is introduced to obtain numerical solutions for the COVID-19 system and combined with the improved MH-NMSS-PSO parameter estimation method to fit the real data of Delhi, India from April 1, 2020 to June 30, 2020. The results show that the fitting effect is quite ideal. Finally, long-term predictions were made on the number of infections. We accurately estimate that the peak number of infections in Delhi, India, can reach around 2.1 million. This paper also compares the fitting performance of the integer-order COVID-19 model and the fractional-order COVID-19 model using the real data from Delhi. The results indicate that the fractional-order model with different orders, as we proposed, performs the best.

Numerical Models and Methods of Atmospheric Parameters Originating in the Formation of the Earth’s Climatic Cycle

Atmospheric models are physical equations based on the ideal gas law. Applied to the atmosphere, this law yields equations for water, vapor (gas), ice, air, humidity, dryness, fire, and heat, thus defining the model of key atmospheric parameters. The distribution of these parameters across the entire planet Earth is the origin of the formation of the climatic cycle, which is a normal climatic variation. To do this, the Earth is divided into eight (8) parts according to the number of key parameters to be defined in a physical representation of the model. Following this distribution, numerical models calculate the constants for the formation of water, vapor, ice, dryness, thermal energy (fire), heat, air, and humidity. These models vary in complexity depending on the indirect trigonometric direction and simplicity in the sum of neighboring models. Note that the constants obtained from the equations yield 275.156˚K (2.006˚C) for water, 273.1596˚K (0.00963˚C) for vapor, 273.1633˚K (0.0133˚C) for ice, 0.00365 in/s for atmospheric dryness, 1.996 in2/s for humidity, 2.993 in2/s for air, 1 J for thermal energy of fire, and 0.9963 J for heat. In summary, this study aims to define the main parameters and natural phenomena contributing to the modification of planetary climate. .

Runoff Modeling in Ungauged Catchments Using Machine Learning Algorithm-Based Model Parameters Regionalization Methodology

Model parameters estimation is a pivotal issue for runoff modeling in ungauged catchments.The nonlinear relationship between model parameters and catchment descriptors is a major obstacle for parameter regionalization,which is the most widely used approach.Runoff modeling was studied in 38 catchments located in the Yellow–Huai–Hai River Basin(YHHRB).The values of the Nash–Sutcliffe efficiency coefficient(NSE),coefficient of determination(R2),and percent bias(PBIAS)indicated the acceptable performance of the soil and water assessment tool(SWAT)model in the YHHRB.Nine descriptors belonging to the categories of climate,soil,vegetation,and topography were used to express the catchment characteristics related to the hydrological processes.The quantitative relationships between the parameters of the SWAT model and the catchment descriptors were analyzed by six regression-based models,including linear regression(LR)equations,support vector regression(SVR),random forest(RF),k-nearest neighbor(kNN),decision tree(DT),and radial basis function(RBF).Each of the 38 catchments was assumed to be an ungauged catchment in turn.Then,the parameters in each target catchment were estimated by the constructed regression models based on the remaining 37 donor catchments.Furthermore,the similaritybased regionalization scheme was used for comparison with the regression-based approach.The results indicated that the runoff with the highest accuracy was modeled by the SVR-based scheme in ungauged catchments.Compared with the traditional LR-based approach,the accuracy of the runoff modeling in ungauged catchments was improved by the machine learning algorithms because of the outstanding capability to deal with nonlinear relationships.The performances of different approaches were similar in humid regions,while the advantages of the machine learning techniques were more evident in arid regions.When the study area contained nested catchments,the best result was calculated with the similarity-based parameter regionalization scheme because of the high catchment density and short spatial distance.The new findings could improve flood forecasting and water resources planning in regions that lack observed data.

Modeling of Computer Virus Propagation with Fuzzy Parameters

Typically,a computer has infectivity as soon as it is infected.It is a reality that no antivirus programming can identify and eliminate all kinds of viruses,suggesting that infections would persevere on the Internet.To understand the dynamics of the virus propagation in a better way,a computer virus spread model with fuzzy parameters is presented in this work.It is assumed that all infected computers do not have the same contribution to the virus transmission process and each computer has a different degree of infectivity,which depends on the quantity of virus.Considering this,the parametersβandγbeing functions of the computer virus load,are considered fuzzy numbers.Using fuzzy theory helps us understand the spread of computer viruses more realistically as these parameters have fixed values in classical models.The essential features of the model,like reproduction number and equilibrium analysis,are discussed in fuzzy senses.Moreover,with fuzziness,two numerical methods,the forward Euler technique,and a nonstandard finite difference(NSFD)scheme,respectively,are developed and analyzed.In the evidence of the numerical simulations,the proposed NSFD method preserves the main features of the dynamic system.It can be considered a reliable tool to predict such types of solutions.

Effect of Milling Parameters on DEM Modeling of a Planetary Ball Mill

The effects of the milling parameters involving shape of powder particles, rotation speed, and ball-to-powder diameter (BPDR) on DEM modeling in the planetary ball mill were investigated. BPDR was varied from 1 to 10. The results revealed that the size and shape of the powder particles do not give a significant change in simulation results when BPDR attains maximum value of 10. The increasing of BPDR leads to the increase of simulation time and size. Hence, the effect of change of the powder particle shape on the calculated data size is not significant. The results also revealed that the increasing rotation speed increases impact energy between powder particles.

Transient elastography with controlled attenuation parameter for the diagnosis of colorectal polyps in patients with nonalcoholic fatty liver disease

BACKGROUND The severity of nonalcoholic fatty liver disease(NAFLD)and lipid metabolism are related to the occurrence of colorectal polyps.Liver-controlled attenuation parameters(liver-CAPs)have been established to predict the prognosis of hepatic steatosis patients.AIM To explore the risk factors associated with colorectal polyps in patients with NAFLD by analyzing liver-CAPs and establishing a diagnostic model.METHODS Patients who were diagnosed with colorectal polyps in the Department of Gastroenterology of our hospital between June 2021 and April 2022 composed the case group,and those with no important abnormalities composed the control group.The area under the receiver operating characteristic curve was used to predict the diagnostic efficiency.Differences were considered statistically significant when P<0.05.RESULTS The median triglyceride(TG)and liver-CAP in the case group were significantly greater than those in the control group(mmol/L,1.74 vs 1.05;dB/m,282 vs 254,P<0.05).TG and liver-CAP were found to be independent risk factors for colorectal polyps,with ORs of 2.338(95%CI:1.154–4.733)and 1.019(95%CI:1.006–1.033),respectively(P<0.05).And there was no difference in the diagnostic efficacy between liver-CAP and TG combined with liver-CAP(TG+CAP)(P>0.05).When the liver-CAP was greater than 291 dB/m,colorectal polyps were more likely to occur.CONCLUSION The levels of TG and liver-CAP in patients with colorectal polyps are significantly greater than those patients without polyps.Liver-CAP alone can be used to diagnose NAFLD with colorectal polyps.

Total ionizing dose effect modeling method for CMOS digital-integrated circuit

Simulating the total ionizing dose(TID)of an electrical system using transistor-level models can be difficult and expensive,particularly for digital-integrated circuits(ICs).In this study,a method for modeling TID effects in complementary metaloxide semiconductor(CMOS)digital ICs based on the input/output buffer information specification(IBIS)was proposed.The digital IC was first divided into three parts based on its internal structure:the input buffer,output buffer,and functional area.Each of these three parts was separately modeled.Using the IBIS model,the transistor V-I characteristic curves of the buffers were processed,and the physical parameters were extracted and modeled using VHDL-AMS.In the functional area,logic functions were modeled in VHDL according to the data sheet.A golden digital IC model was developed by combining the input buffer,output buffer,and functional area models.Furthermore,the golden ratio was reconstructed based on TID experimental data,enabling the assessment of TID effects on the threshold voltage,carrier mobility,and time series of the digital IC.TID experiments were conducted using a CMOS non-inverting multiplexer,NC7SZ157,and the results were compared with the simulation results,which showed that the relative errors were less than 2%at each dose point.This confirms the practicality and accuracy of the proposed modeling method.The TID effect model for digital ICs developed using this modeling technique includes both the logical function of the IC and changes in electrical properties and functional degradation impacted by TID,which has potential applications in the design of radiation-hardening tolerance in digital ICs.

A Novel On-Site-Real-Time Method for Identifying Characteristic Parameters Using Ultrasonic Echo Groups and Neural Network

On-site and real-time non-destructive measurement of elastic constants for materials of a component in a in-service structure is a challenge due to structural complexities,such as ambiguous boundary,variable thickness,nonuniform material properties.This work develops for the first time a method that uses ultrasound echo groups and artificial neural network(ANN)for reliable on-site real-time identification of material parameters.The use of echo groups allows the use of lower frequencies,and hence more accommodative to structural complexity.To train the ANNs,a numerical model is established that is capable of computing the waveform of ultrasonic echo groups for any given set of material properties of a given structure.The waveform of an ultrasonic echo groups at an interest location on the surface the structure with material parameters varying in a predefined range are then computed using the numerical model.This results in a set of dataset for training the ANN model.Once the ANN is trained,the material parameters can be identified simultaneously using the actual measured echo waveform as input to the ANN.Intensive tests have been conducted both numerically and experimentally to evaluate the effectiveness and accuracy of the currently proposed method.The results show that the maximum identification error of numerical example is less than 2%,and the maximum identification error of experimental test is less than 7%.Compared with currently prevailing methods and equipment,the proposefy the density and thickness,in addition to the elastic constants.Moreover,the reliability and accuracy of inverse prediction is significantly improved.Thus,it has broad applications and enables real-time field measurements,which has not been fulfilled by any other available methods or equipment.

Analysis and comparison of retinal vascular parameters under different glucose metabolic status based on deep learning

AIM:To develop a deep learning-based model for automatic retinal vascular segmentation,analyzing and comparing parameters under diverse glucose metabolic status(normal,prediabetes,diabetes)and to assess the potential of artificial intelligence(AI)in image segmentation and retinal vascular parameters for predicting prediabetes and diabetes.METHODS:Retinal fundus photos from 200 normal individuals,200 prediabetic patients,and 200 diabetic patients(600 eyes in total)were used.The U-Net network served as the foundational architecture for retinal arteryvein segmentation.An automatic segmentation and evaluation system for retinal vascular parameters was trained,encompassing 26 parameters.RESULTS:Significant differences were found in retinal vascular parameters across normal,prediabetes,and diabetes groups,including artery diameter(P=0.008),fractal dimension(P=0.000),vein curvature(P=0.003),C-zone artery branching vessel count(P=0.049),C-zone vein branching vessel count(P=0.041),artery branching angle(P=0.005),vein branching angle(P=0.001),artery angle asymmetry degree(P=0.003),vessel length density(P=0.000),and vessel area density(P=0.000),totaling 10 parameters.CONCLUSION:The deep learning-based model facilitates retinal vascular parameter identification and quantification,revealing significant differences.These parameters exhibit potential as biomarkers for prediabetes and diabetes.

Optimizing near-carbon-free nuclear energy systems:advances in reactor operation digital twin through hybrid machine learning algorithms for parameter identification and state estimation

Accurate and efficient online parameter identification and state estimation are crucial for leveraging digital twin simulations to optimize the operation of near-carbon-free nuclear energy systems.In previous studies,we developed a reactor operation digital twin(RODT).However,non-differentiabilities and discontinuities arise when employing machine learning-based surrogate forward models,challenging traditional gradient-based inverse methods and their variants.This study investigated deterministic and metaheuristic algorithms and developed hybrid algorithms to address these issues.An efficient modular RODT software framework that incorporates these methods into its post-evaluation module is presented for comprehensive comparison.The methods were rigorously assessed based on convergence profiles,stability with respect to noise,and computational performance.The numerical results show that the hybrid KNNLHS algorithm excels in real-time online applications,balancing accuracy and efficiency with a prediction error rate of only 1%and processing times of less than 0.1 s.Contrastingly,algorithms such as FSA,DE,and ADE,although slightly slower(approximately 1 s),demonstrated higher accuracy with a 0.3%relative L_2 error,which advances RODT methodologies to harness machine learning and system modeling for improved reactor monitoring,systematic diagnosis of off-normal events,and lifetime management strategies.The developed modular software and novel optimization methods presented offer pathways to realize the full potential of RODT for transforming energy engineering practices.

Artificial Intelligence Based Meteorological Parameter Forecasting for Optimizing Response of Nuclear Emergency Decision Support System

This paper presents a novel artificial intelligence (AI) based approach to predict crucial meteorological parameters such as temperature,pressure,and wind speed,typically calculated from computationally intensive weather research and forecasting (WRF) model.Accurate meteorological data is indispensable for simulating the release of radioactive effluents,especially in dispersion modeling for nuclear emergency decision support systems.Simulation of meteorological conditions during nuclear emergencies using the conventional WRF model is very complex and time-consuming.Therefore,a new artificial neural network (ANN) based technique was proposed as a viable alternative for meteorological prediction.A multi-input multi-output neural network was trained using historical site-specific meteorological data to forecast the meteorological parameters.Comprehensive evaluation of this technique was conducted to test its performance in forecasting various parameters including atmospheric pressure,temperature,and wind speed components in both East-West and North-South directions.The performance of developed network was evaluated on an unknown dataset,and acquired results are within the acceptable range for all meteorological parameters.Results show that ANNs possess the capability to forecast meteorological parameters,such as temperature and pressure,at multiple spatial locations within a grid with high accuracy,utilizing input data from a single station.However,accuracy is slightly compromised when predicting wind speed components.Root mean square error (RMSE) was utilized to report the accuracy of predicted results,with values of 1.453℃for temperature,77 Pa for predicted pressure,1.058 m/s for the wind speed of U-component and 0.959 m/s for the wind speed of V-component.In conclusion,this approach offers a precise,efficient,and wellinformed method for administrative decision-making during nuclear emergencies.

Improved cat swarm optimization for parameter estimation of mixed additive and multiplicative random error model

To estimate the parameters of the mixed additive and multiplicative(MAM)random error model using the weighted least squares iterative algorithm that requires derivation of the complex weight array,we introduce a derivative-free cat swarm optimization for parameter estimation.We embed the Powell method,which uses conjugate direction acceleration and does not need to derive the objective function,into the original cat swarm optimization to accelerate its convergence speed and search accuracy.We use the ordinary least squares,weighted least squares,original cat swarm optimization,particle swarm algorithm and improved cat swarm optimization to estimate the parameters of the straight-line fitting MAM model with lower nonlinearity and the DEM MAM model with higher nonlinearity,respectively.The experimental results show that the improved cat swarm optimization has faster convergence speed,higher search accuracy,and better stability than the original cat swarm optimization and the particle swarm algorithm.At the same time,the improved cat swarm optimization can obtain results consistent with the weighted least squares method based on the objective function only while avoiding multiple complex weight array derivations.The method in this paper provides a new idea for theoretical research on parameter estimation of MAM error models.

Effects of gravel on the water absorption characteristics and hydraulic parameters of stony soil

The eastern foothills of the Helan Mountains in China are a typical mountainous region of soil and gravel,where gravel could affect the water movement process in the soil.This study focused on the effects of different gravel contents on the water absorption characteristics and hydraulic parameters of stony soil.The stony soil samples were collected from the eastern foothills of the Helan Mountains in April 2023 and used as the experimental materials to conduct a one-dimensional horizontal soil column absorption experiment.Six experimental groups with gravel contents of 0%,10%,20%,30%,40%,and 50%were established to determine the saturated hydraulic conductivity(K_(s)),saturated water content(θ_(s)),initial water content(θ_(i)),and retention water content(θ_(r)),and explore the changes in the wetting front depth and cumulative absorption volume during the absorption experiment.The Philip model was used to fit the soil absorption process and determine the soil water absorption rate.Then the length of the characteristic wetting front depth,shape coefficient,empirical parameter,inverse intake suction and soil water suction were derived from the van Genuchten model.Finally,the hydraulic parameters mentioned above were used to fit the soil water characteristic curves,unsaturated hydraulic conductivity(K_(θ))and specific water capacity(C(h)).The results showed that the wetting front depth and cumulative absorption volume of each treatment gradually decreased with increasing gravel content.Compared with control check treatment with gravel content of 0%,soil water absorption rates in the treatments with gravel contents of 10%,20%,30%,40%,and 50%decreased by 11.47%,17.97%,25.24%,29.83%,and 42.45%,respectively.As the gravel content increased,inverse intake suction gradually increased,and shape coefficient,K_(s),θ_(s),andθ_(r)gradually decreased.For the same soil water content,soil water suction and K_(θ)gradually decreased with increasing gravel content.At the same soil water suction,C(h)decreased with increasing gravel content,and the water use efficiency worsened.Overall,the water holding capacity,hydraulic conductivity,and water use efficiency of stony soil in the eastern foothills of the Helan Mountains decreased with increasing gravel content.This study could provide data support for improving soil water use efficiency in the eastern foothills of the Helan Mountains and other similar rocky mountainous areas.

Analysis of Modeling the Influence of Electromagnetic Fields Radiated by Industrial Static Converters and Impacts on Operators Using Maxwell’s Equations

The study of Electromagnetic Compatibility is essential to ensure the harmonious operation of electronic equipment in a shared environment. The basic principles of Electromagnetic Compatibility focus on the ability of devices to withstand electromagnetic disturbances and not produce disturbances that could affect other systems. Imperceptible in most work situations, electromagnetic fields can, beyond certain thresholds, have effects on human health. The objective of the present article is focused on the modeling analysis of the influence of geometric parameters of industrial static converters radiated electromagnetic fields using Maxwell’s equations. To do this we used the analytical formalism for calculating the electromagnetic field emitted by a filiform conductor, to model the electromagnetic radiation of this device in the spatio-temporal domain. The interactions of electromagnetic waves with human bodies are complex and depend on several factors linked to the characteristics of the incident wave. To model these interactions, we implemented the physical laws of electromagnetic wave propagation based on Maxwell’s and bio-heat equations to obtain consistent results. These obtained models allowed us to evaluate the spatial profile of induced current and temperature of biological tissue during exposure to electromagnetic waves generated by this system. The simulation 2D results obtained from computer tools show that the temperature variation and current induced by the electromagnetic field can have a very significant influence on the life of biological tissue. The paper provides a comprehensive analysis using advanced mathematical models to evaluate the influence of electromagnetic fields. The findings have direct implications for workplace safety, potentially influencing standards and regulations concerning electromagnetic exposure in industrial settings.

Parameter sensitivity analysis for a biochemically-based photosynthesis model

A challenge for the development of Land Surface Models(LSMs) is improving transpiration of water exchange and photosynthesis of carbon exchange between terrestrial plants and the atmosphere, both of which are governed by stoma in leaves. In the photosynthesis module of these LSMs, variations of parameters arising from diversity in plant functional types(PFTs) and climate remain unclear. Identifying sensitive parameters among all photosynthetic parameters before parameter estimation can not only reduce operation cost, but also improve the usability of photosynthesis models worldwide. Here, we analyzed 13 parameters of a biochemically-based photosynthesis model(FvCB), implemented in many LSMs, using two sensitivity analysis(SA) methods(i.e., the Sobol’ method and the Morris method) for setting up the parameter ensemble. Three different model performance metrics, i.e.,Root Mean Squared Error(RMSE), Nash Sutcliffe efficiency(NSE), and Standard Deviation(STDEV) were introduced for model assessment and sensitive parameters identification. The results showed that among all photosynthetic parameters only a small portion of parameters were sensitive, and the sensitive parameters were different across plant functional types: maximum rate of Rubisco activity(Vcmax25), maximum electron transport rate(Jmax25), triose phosphate use rate(TPU) and dark respiration in light(Rd) were sensitive in broad leafevergreen trees(BET), broad leaf-deciduous trees(BDT) and needle leaf-evergreen trees(NET), while only Vcmax25and TPU are sensitive in short vegetation(SV), dwarf trees and shrubs(DTS), and agriculture and grassland(AG). The two sensitivity analysis methods suggested a strong SA coherence;in contrast, different model performance metrics led to different SA results. This misfit suggests that more accurate values of sensitive parameters, specifically, species specific and seasonal variable parameters, are required to improve the performance of the FvCB model.

Constructing Confidence Regions for Autoregressive-Model Parameters

We discuss formulas and techniques for finding maximum-likelihood estimators of parameters of autoregressive (with particular emphasis on Markov and Yule) models, computing their asymptotic variance-covariance matrix and displaying the resulting confidence regions;Monte Carlo simulation is then used to establish the accuracy of the corresponding level of confidence. The results indicate that a direct application of the Central Limit Theorem yields errors too large to be acceptable;instead, we recommend using a technique based directly on the natural logarithm of the likelihood function, verifying its substantially higher accuracy. Our study is then extended to the case of estimating only a subset of a model’s parameters, when the remaining ones (called nuisance) are of no interest to us.

Introduction to Mesh Based Generated Lumped Parameter Models for Electromagnetic Problems using Triangular Elements

This paper is an introduction to mesh based generated reluctance network modeling using triangular elements.Many contributions on mesh based generated reluctance networks using rectangular shaped elements have been published,but very few on those generated from a mesh using triangular elements.The use of triangular elements is aimed at extending the application of the approach to any shape of modeled devices.Basic concepts of the approach are presented in the case of electromagnetic devices.The procedure for coding the approach in the case of a flat linear permanent magnet machine is presented.Codes developed under MATLAB environment are also included.