This paper presents an approach to the Euler-Bernoulli beam theory (EBBT) using the finite difference method (FDM). The EBBT covers the case of small deflections, and shear deformations are not considered. The FDM is an approximate method for solving problems described with differential equations. The FDM does not involve solving differential equations; equations are formulated with values at selected points of the structure. Generally, the finite difference approximations are derived based on fourth-order polynomial hypothesis (FOPH) and second-order polynomial hypothesis (SOPH) for the deflection curve; the FOPH is made for the fourth and third derivative of the deflection curve while the SOPH is made for its second and first derivative. In addition, the boundary conditions and not the governing equations are applied at the beam’s ends. In this paper, the FOPH was made for all of the derivatives of the deflection curve, and additional points were introduced at the beam’s ends and positions of discontinuity (concentrated loads or moments, supports, hinges, springs, etc.). The introduction of additional points allowed us to apply the governing equations at the beam’s ends and to satisfy the boundary and continuity conditions. Moreover, grid points with variable spacing were also considered, the grid being uniform within beam segments. First-order analysis, second-order analysis, and vibration analysis of structures were conducted with this model. Furthermore, tapered beams were analyzed (element stiffness matrix, second-order analysis). Finally, a direct time integration method (DTIM) was presented. The FDM-based DTIM enabled the analysis of forced vibration of structures, with damping taken into account. The results obtained in this paper showed good agreement with those of other studies, and the accuracy was increased through a grid refinement. Especially in the first-order analysis of uniform beams, the results were exact for uniformly distributed and concentrated loads regardless of the grid. Further research will be needed to investigate polynomial refinements (higher-order polynomials such as fifth-order, sixth-order…) of the deflection curve; the polynomial refinements aimed to increase the accuracy, whereby non-centered finite difference approximations at beam’s ends and positions of discontinuity would be used.

The rehabilitation of the road embankments alongside the irrigation canal, which is constructed on a Bangkok soft clay foundation, has resulted from land use growth for several decades. As a result, the embankment is always confronted with the problem of instability. Numerous articles have demonstrated that rapid drawdown of the water level in the canal could be one of the most significant reasons, whereas vehicle cyclic loading may potentially contribute to some effect of the embankment failure. To explore and assess the effects of the influence factors, the finite element method (FEM) and conventional and dynamic triaxial approaches were employed. The study also investigated the consequences of varying thicknesses of soft clay and cyclic loading with different magnitudes and frequencies. Consequently, cyclic loading can lead to an augmentation in strain distribution and excess pore water pressure. Likewise, more significant deformation and excess pore water pressure were induced by a lower frequency and deeper thickness of soft clay compared to a higher frequency and thinner thickness. In brief, vehicle repetitive loading is one of the primary factors contributing to embankment failure, with the threshold stress being approximately three-fourths of the undrained shear strength. In contrast, threshold strain is approximately one to two-fold of the monotonic strain condition.

Abstract: This paper outlines the methodology and results for a two species finite volume scalar computational drug transport model developed for simulating the mass transport of Poly(lactic-co-glycolic acid (PLGA) ) from a half-embedded single strut implanted in a coronary arterial vessel wall. The mathematical drug transport model incorporates the convection-diffusion equation in scalar form (dimensionless) with a two species (free-drug and bound-drug) mass transport setup including reversible equilibrium reaction source terms for the free and bound-drug states to account for the pharmaco-kinetic reactions in the arterial wall. The relative reaction rates of the added source terms control the interconversion of drug between the free and bound-drug states. The model is solved by a 2-D finite-volume method for discretizing and solving the free and bound drug transport equations with anisotropic vascular drug diffusivities. This model is an improvement over our previously developed model using the finite-difference method. A dimensionless characteristic scaling pre-analysis was conducted a priori to evaluate the significance of implementing the reaction source terms in the transport equations. This paper reports the findings of an investigation of the interstitial flow profile into the arterial wall and the free and bound drug diffusion profiles with a parametric study of varying polymer drug concentration (low and high), tortuosity, porosity, and Peclet and DamKöhler numbers over the course of 400 hours (16.67 days). The results also reveal how a single species drug delivery model that neglects both a reversible binding reaction source term and the porosity and tortuosity of the arterial wall cannot accurately predict the distribution of both the free and bound drug.

This paper presents an approach to the two-dimensional analysis of elastic isotropic deep beams using the finite difference method (FDM). Deep beams are subjected to in-plane loading and present a shear span to height ratio of less than 2.50; consequently, Euler-Bernoulli beam theory and Timoshenko beam theory do not apply. Deep beams analysis is generally conducted using numerical methods such as the finite element method and to a lesser extent the FDM; the strut-and-tie model and the stress field method are also widely utilized. Analytical approaches usually make use of the Airy stress function, where stresses are formulated in terms of the stress function; however, the exact solution of this function satisfying all of the boundary conditions can hardly be found, even for simple cases. In this paper, deep beams were analyzed using the FDM. The FDM is an approximate method for solving problems described with differential equations. The FDM does not involve solving differential equations; equations are formulated with values at selected nodes of the structure. Therefore, the deep beam was discretized with a two-dimensional grid, and additional nodes were introduced at the boundaries and at positions of discontinuity (openings, brutal change of material properties, non-uniform grid spacing), the number of additional nodes corresponding to the number of boundary conditions at the node of interest. The introduction of additional nodes allowed us to apply the governing equations at boundary nodes and satisfy the boundary and continuity conditions. An Airy stress function approach and a displacement potential function approach were considered in this study whereby strong formulations of equations (equilibrium, kinematic, and constitutive) were set. Stress and stability analyses were carried out with this model; furthermore, deep beams of varying stiffness, layered beams, and beams having openings were analyzed. For slender beams, the results obtained with the Airy stress function approach showed good agreement with those of the Euler-Bernoulli beam theory, and for deep beams, the shapes of stress distributions were in good agreement with a proper understanding of the behavior of structures. On the other hand, the displacement potential function approach delivered unsatisfactory results, probably due to the use of an inefficient equation solver; a more powerful tool will be needed in future research for this purpose.

The combination of interstitial hyperthermia treatment (IHT) with high dose rate brachytherapy (HDR-BT) for prostate cancer treatment and has the potential to improve clinical outcome, since it highly enhances the efficiency of cell kill, especially when applied simultaneously. Therefore, we have developed the ThermoBrachy applicators. To effectively apply optimal targeted IHT, treatment planning is considered essential. However, treatment planning in IHT is rarely applied since it is regarded difficult to accurately calculate the deposited energy in the tissue in a short enough time for clinical practice. In this study, we investigated various time-efficient methods for fast computation of the electromagnetic (EM) energy deposition resulting from the ThermoBrachy applicators. Initially, we investigated the use of an electro-quasistatic solver. Next, we extended our investigation to the application of geometric simplifications. Furthermore, we investigated the validity of the superpositioning principle, which can enable adaptive treatment plan optimization without the need for continuous recomputation of the EM field. Finally, we evaluated the accuracy of the methods by comparing them to the golden standard Finite-Difference Time-Domain calculation method using gamma-index analysis. The simplifications considerably reduced the computation time needed, improving from >12 h to a few seconds. All investigated methods showed excellent agreement with the golden standard by showing a >99% passing rate with 1%/0.5 mm Dose Difference and Distance-to-Agreement criteria. These results allow the proposed electromagnetic simulation method to be used for fast and accurate adaptive treatment planning.

This paper presents an approach to the Kirchhoff-Love plate theory (KLPT) using the finite difference method (FDM). The KLPT covers the case of small deflections, and shear deformations are not considered. The FDM is an approximate method for solving problems described with differential equations. The FDM does not involve solving differential equations; equations are formulated with values at selected points of the structure. Generally in the case of KLPT, the finite difference approximations are derived based on the fourth-order polynomial hypothesis (FOPH) and second-order polynomial hypothesis (SOPH) for the deflection surface. The FOPH is made for the fourth and third derivative of the deflection surface while the SOPH is made for its second and first derivative; this leads to a 13-point stencil for the governing equation. In addition, the boundary conditions and not the governing equations are applied at the plate edges. In this paper, the FOPH was made for all of the derivatives of the deflection surface; this led to a 25-point stencil for the governing equation. Furthermore, additional nodes were introduced at plate edges and at positions of discontinuity (continuous supports/hinges, incorporated beams, stiffeners, brutal change of stiffness, etc.), the number of additional nodes corresponding to the number of boundary conditions at the node of interest. The introduction of additional nodes allowed us to apply the governing equations at the plate edges and to satisfy the boundary and continuity conditions. First-order analysis, second-order analysis, buckling analysis, and vibration analysis of plates were conducted with this model. Moreover, plates of varying thickness and plates with stiffeners were analyzed. Finally, a direct time integration method (DTIM) was presented. The FDM-based DTIM enabled the analysis of forced vibration of structures, with damping taken into account. In first-order, second-order, buckling, and vibration analyses of rectangular plates, the results obtained in this paper were in good agreement with those of well-established methods, and the accuracy was increased through a grid refinement.

The electrical motor design requires a lot of work and resources. In the field, motor designers need to decide on competitive motor sizing quickly to compete with other motor suppliers. Generally, it requires expensive FEM S/W to get accurate motor sizing results, but buying and maintaining many licenses is costly. Moreover, it needs some patience to get the result from FEM analysis. In this paper, the motor sizing technique using an existing production motor is presented to respond to urgent requests from customers without using FEM analysis. In order to enhance the fidelity of the result, magnetic flux loss by magnet coating and non-linear characteristics are considered during the motor sizing process. The proposed sizing method was verified by comparison to the FEM analysis and test result and showed reasonable confidence in them. Also, this method has the substantial advantage of outputting accurate results instantly without taking any time. It helps with the modification design of electrical motors without using expensive FEM S/W and a certain level of skill or experience.

Large values and gradients of stresses and deformations, triggering concentrations of stresses and deformations, arise in the corner areas of a structure. The action of forced deformations, leading to the finite rupture of the contact between the elements of a structure, also triggers the concentration of stresses, while the rupture reaches an irregular point, a line on the area boundary. The theoretical analysis of the stress-strain state (SSS) of areas with angular cutouts in the boundary under the action of discontinuous forced deformations is reduced to the study of singular solutions to the homogeneous problem of the elasticity theory that has power-related features. The calculation of stress concentration coefficients in the domain of a singular solution to the elastic problem makes no sense. It is experimentally proven that the zone, that is close to the vertex of the angular cutout in the area boundary, has substantial deformations, rotations, and it corresponds to rising values of the first and second derivatives of displacements along the radius in cases of sufficiently small radii in the neighbourhood of the irregular point of the boundary. For such areas, it is necessary to consider the plane problem of the elasticity theory, taking into account the geometric nonlinearity under the action of forced deformations. This will allow analyzing the effect of relations between orders of values of deformations, rotations, and forced deformations on the form of the equation of equilibrium. The purpose of this work is to analyze the effect of relations of deformation orders, rotations, forced deformations on the form of the equilibrium equation in the polar coordinate system for a V-shaped area under the action of forced temperature-induced deformations with regard for the geometrical non-linearity and physical linearity.

Robots in space are necessarily extremely light and lack structural stiffness resulting in natural frequencies of resonance so low as to reside inside the attitude controller’s bandwidth. A variety of input trajectories can be used to drive a controller’s attempt to ameliorate the control-structural interactions where feedback is provided by low-quality, noisy sensors. Traditionally, step functions are used as the ideal input trajectory. However, step functions are not ideal in many applications, as they are discontinuous. Alternative input trajectories are explored in this manuscript and applied to an example system that includes a flexible appendage attached to a rigid main body. The main body is controlled by a reaction wheel. The equations of motion of the flexible appendage, rigid body, and reaction wheel are derived. A feedback controller is developed to account for the rigid body modes. Additional filters are added to compensate for the system’s flexible modes. Sinusoidal trajectories are autonomously generated to feed the controller. Whiplash compensation is additionally implemented for comparison. The control method without random errors with the smallest error is the sinusoidal trajectory method, which showed a 97.39% improvement when compared to the baseline response when step trajectories were commanded, while the sinusoidal method was inferior to traditional step trajectories when sensor noise and random errors were present.

Vibration suppression in a field of frequencies and the creation of band gaps which do not allow the propagation of the waves is studied in this paper by means of microstructures designed through topology optimization. Topology optimization is formulated in frequency domain. Band gap design is based on Floquet-Bloch theory and genetic optimization. The result is appealing in view of modern 3D printing techniques.

This study addresses the enhancement of material efficiency and reduction of brittleness in timber-to-concrete adhesive connections for beam-type timber and timber-concrete composite panels. The research explores the potential benefits of adding longitudinal timber ribs to Cross-Laminated Timber (CLT) beam-type panels. Three groups of flexure-tested specimens were analysed: (1) timber panels (1400 mm × 400 mm) with two 100 mm thick CLT panels and two 60 mm thick CLT panels reinforced with 150x80 mm timber ribs; (2) eight specimens (600 mm × 100 mm × 150 mm) with CLT members (600 mm × 100 mm × 100 mm) connected to a 50 mm concrete layer using granite chips and Sikadur-31 (AB) epoxy adhesive; (3) six CLT panels (1400 mm × 400 mm × 50 mm) bonded to a 50 mm concrete layer, with two panels containing polypropylene microfibers and two panels incorporating polyethene dowels for mechanical connection. Specimens were subjected to three-point bending tests and analysed using the transformed section method, γ-method, and finite element method with ANSYS 2023R2 software. Results indicated a 53% increase in load-carrying capacity for ribbed CLT panels with no additional material consumption, a 24.8%-41.1% increase for CLT panels strengthened with a concrete layer, and improved ductility and prevention of disintegration in timber-concrete composites with polypropylene microfibers.

High temperature (HT) geothermal wells can provide green power 24h a day 7 days a week. Harsh environmental and operational conditions, long-term durability requirements of such wells require special cementitious composites for well construction. This paper reports a comprehensive assessment of geothermal cement composites in cyclic pressure function laboratory tests, field exposures in an HT geothermal well (300-350oC) as well as a numerical model to complement the experimental results. Performances of calcium-aluminate-cement (CAC)-based composites and calcium free cement were compared against reference Ordinary Portland Cement (OPC)/silica blend. The stability and degradation of the tested materials was characterized by crystalline composition, thermo-gravimetric, and elemental analyses, morphological studies, water-fillable porosity, and mechanical properties measurements. All CAC-based formulations outperformed the reference blend both in the function and exposure tests. The reference OPC/silica lost its mechanical properties during the 9-month well exposure through extensive HT carbonation, while properties of the CAC-based blends improved over that period. The Modified Cam-Clay (MCC) plasticity parameters of several HT cement formulations were extracted from triaxial and Brazilian tests and verified against the experimental results of function cyclic tests. These parameters can be used in well integrity models to predict field-scale behavior of the cement sheath under geothermal well conditions.

Let $f\left(1\right)=1$ , and let $f(n+1)=2^{2^{f\left(n\right)}}$ for every positive integer n. We consider the following hypothesis: if a system S ⊆ {x_i · x_j = x_k : i, j, k ∈ {1, . . . , n}} ∪ {x_i + 1 = x_k : i, k ∈{1, . . . , n}} has only finitely many solutions in non-negative integers x₁, . . . , x_n, then each such solution (x₁, . . . , x_n) satisfies x₁, . . . , x_n ≤ f (2n). We prove: (1) the hypothesisimplies that there exists an algorithm which takes as input a Diophantine equation, returns an integer, and this integer is greater than the heights of integer (non-negative integer, positive integer, rational) solutions, if the solution set is finite; (2) the hypothesis implies that there exists an algorithm for listing the Diophantine equations with infinitely many solutions in non-negative integers; (3) the hypothesis implies that the question whether or not a given Diophantine equation has only finitely many rational solutions is decidable by a single query to an oracle that decides whether or not a given Diophantine equation has a rational solution; (4) the hypothesis implies that the question whether or not a given Diophantine equation has only finitely many integer solutions is decidable by a single query to an oracle that decides whether or not a given Diophantine equation has an integer solution; (5) the hypothesis implies that if a set $\mathfrak{M}$ ⊆ has a finite-fold Diophantine representation, then $\mathfrak{M}$ is computable.

In this paper, we extend the KdVSKR equations in (1+1)-dimension and (2+1)-dimension to fractional KdVSKR equations with modified Riemann-Liouville derivative. The (2+1)-dimensional KdVSKR equation, which is a recent extension of (1+1)-dimensional KdVSKR equation, can model the resonances of solitons in shallow water. By means of the Hirota bilinear method, finite symmetry group method and consistent Riccati expansion method, many new interaction solutions have been derived. Soliton-cnoidal interaction solution of the (1+1)-dimensional fractional KdVSKR equation has been derived for the first time. For the (2+1)-dimensional frac-tional KdVSKR equation, two-wave interaction solutions and three-wave interaction solutions, including dark-soliton-sine interaction solution, bright-soliton-elliptic interaction solution, and lump-hyperbolic-sine interaction solution, have been derived. Impact of the fractional order on the shapes of these solutions has been illustrated by figures. Three-wave interaction solutions of fractional systems have not been reported in the existing references. The research idea in this paper can be applied to other fractional differential equations.

Recently, more wind turbine systems are being installed in deep waters far from the coast. Several concepts of floating wind turbine systems (FWTSs) are developed, among which the semisubmersible platform, due to its applicability in different water depths, good hydrodynamic performance, and facility in the installation process, constitutes the most explored technology compared to the others. However, a significant obstacle to the industrialization of this technology is a design of a cost-effective FWTS, which can be achieved by optimizing the geometry, size and weight of the floating platform, along with the mooring system. This is only possible by selecting a method capable of accurately analyzing the FWTS coupled hydro-aero-structural dynamics at each design stage. Accordingly, this paper aims at providing a detailed overview of the most common coupled numerical and physical methods, including their basic assumptions, formulations, limitations, and costs, used for analyzing the dynamics of FWTSs, mainly those supported by a semisubmersible, to assist the choice of the most suitable method at each design phase of FWTSs. Finally, the article discusses possible future research to address challenges in modeling FWTSs dynamics that persist to date.