Dynamics of Interaction of a Single Fiber with a Headset of a Sampling Drum ()

Yuldashev Jamshid^{}, Umarov Akbarjon^{}, Sarimsakov Olimjon^{}

Namangan Institute of Engineering and Technology, Namangan, Uzbekistan.

**DOI: **10.4236/eng.2020.126027
PDF HTML XML
120
Downloads
280
Views
Citations

Namangan Institute of Engineering and Technology, Namangan, Uzbekistan.

The article presents an analysis of the dynamics of interaction of a single fiber with the sampling drum of a headset during pneumomechanical spinning. The regularity is investigated and the equation of the movement of fibers on the surface of the tooth of the sampling drum headset is obtained. It was found that the dependences of the movement of fibers along the tooth at different angles of inclination, the distribution of fiber tension along the arc, and the distribution of fiber speeds in the region of the sampling arc have an increase property. As a result of the analysis of the tension and speed of the fibers, it was found that with the increase in the sampling zone, the tension force and speed of the fiber increase.

Keywords

Fibers, Sampling, Angle, Friction Force, Friction, Speed, Radius, Fiber Complex

Share and Cite:

Jamshid, Y. , Akbarjon, U. and Olimjon, S. (2020) Dynamics of Interaction of a Single Fiber with a Headset of a Sampling Drum. *Engineering*, **12**, 347-355. doi: 10.4236/eng.2020.126027.

1. Introduction: Formulation of the Problem

In the discretization zone, the headsets of the sampling drum enter the beard of the clamped fibers. The amount of immersion of the teeth of the headset depends on the shape of the supporting surface to which the fibers are pressed by the headset. The teeth of the headset act on fibers that are within their reach, and those fibers whose relationship with the mass of the tape is less than the impact of the headset are pulled out of the feed product—the tape.

In this case, the interaction of the tooth of the headset with a beard occurs according to two schemes. In the first scheme, the tooth acts on the fiber with its front face, drawing it into motion. According to the second scheme, the lateral face of the tooth creates friction forces, which when the headset is immersed in the beard, cause a certain force in the fiber, which leads to partial straightening of the curved fibers.

The fibers held by the rear ends in the clamp of the feed device are under the influence of friction forces and resistance to pulling fibers, which are in active interaction with the teeth of the headset.

In the process of interaction, one or several fibers are pulled out from the beard, which subsequently come into contact with each other and form a single complex. The separation of the fibers occurs during the implementation of mainly the second interaction scheme, where it occurs at a high frequency of exposure of individual teeth of the headset throughout the sampling area.

In this case, the headset fibers are transported from the region of high density beards in the region of low density. The orientation of individual fibers in the headset depends on their location in the beard before interacting with the teeth of the headset. If the fibers are at an angle to the direction of the headset, they can be elongated with the front edge of the tooth. Such fibers practically do not parallelize, tend to form knots. Taking into account the above mentioned mechanism of fiber movement in the sampling zone, to describe the process of the dynamic interaction of individual fibers with the teeth of the headset, we accept the following calculation schemes.

2. The Main Part

We model the fiber capture pattern with the front face of the tooth, along which it makes a movement, and thus affects the process of pulling it out of the pulp (tape). A number of studies of domestic and foreign scientists are devoted to the modeling of such processes [1] - [9].

The arrangement of a fiber of mass m on the tooth surface is shown in Figure 1, Where α denotes the angle between the front face of the tooth and the direction of the radius—the OB vector. We denote the position of the point M along the tooth by the distance r = BM.

In Figure 1. angles β, δu and γ will be equal

$\beta =\omega \cdot t-{\phi}_{0}$ ; $\gamma =\delta =\alpha +\omega \cdot t-{\phi}_{0}$

The coordinates of the material point of position (mass) M(x, y) are calculated by the formulas

$\u0445=\u041e\u0410+\u0410\u0421=R\mathrm{cos}\beta +r\mathrm{cos}\gamma =R\mathrm{cos}\left(\omega t-{\phi}_{0}\right)+r\mathrm{cos}\left(\alpha +\omega t-{\phi}_{0}\right)$ (1)

$y=\u041e{A}_{1}+{A}_{1}{C}_{1}=R\mathrm{sin}\beta +r\mathrm{sin}\gamma =R\mathrm{sin}\left(\omega t-{\phi}_{0}\right)+r\mathrm{sin}\left(\alpha +\omega t-{\phi}_{0}\right)$ (2)

We choose the distance r = BM of the generalized coordinate and find the kinetic energy of a material point of mass m by the formula:

$\begin{array}{c}\u0422=\frac{m}{2}\left({\stackrel{\dot{}}{x}}^{2}+{\stackrel{\dot{}}{y}}^{2}\right)\\ =\frac{m}{2}\{{\left[-\omega R\mathrm{sin}\left(\omega t-{\phi}_{0}\right)-r\omega \mathrm{sin}\left(\alpha +\omega t-{\phi}_{0}\right)+\stackrel{\dot{}}{r}\mathrm{cos}\left(\alpha +\omega t-{\phi}_{0}\right)\right]}^{2}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{\left[\omega R\mathrm{cos}\left(\omega t-{\phi}_{0}\right)+r\omega \mathrm{cos}\left(\alpha +\omega t-{\phi}_{0}\right)+\stackrel{\dot{}}{r}\mathrm{sin}\left(\alpha +\omega t-{\phi}_{0}\right)\right]}^{2}\}\\ =\frac{m}{2}\left({\omega}^{2}{R}^{2}+2{\omega}^{2}rR\mathrm{cos}\alpha +2\omega R\stackrel{\dot{}}{r}\mathrm{sin}\alpha +{r}^{2}{\omega}^{2}+{\stackrel{\dot{}}{r}}^{2}\right)\end{array}$ (3)

Figure 1. Layout of a fiber of mass m on the tooth surface.

The mass is affected by the gravity of the material point and the friction force created by the force of gravity and the Carioles’ force, the projections along the MV line will be equal.

${F}_{g}=mg\mathrm{sin}\gamma $ ; ${F}_{g\u0442\u0440}=fmg\mathrm{cos}\gamma $ ; ${F}_{\u043a\u0442\u0440}=2mf\omega \stackrel{\dot{}}{r}\mathrm{cos}\alpha $ (4)

We now compose the second-order Lagrange equations

$\frac{\text{d}}{\text{d}t}\left(\frac{\sigma T}{\partial r}\right)-\frac{\sigma T}{\partial r}={F}_{g}-{F}_{g\u0442\u0440}-{F}_{\u043a\u0442\u0440}$ (5)

Substituting the expression T, F_{g}, F_{gtr}, F_{ctr}, from (3) in (2.10) we obtain the equation for determining the distance r:

$\begin{array}{l}\stackrel{\xa8}{r}+2f\omega \mathrm{cos}\alpha \cdot \stackrel{\dot{}}{r}-{\omega}^{2}r\\ ={\omega}^{2}R\mathrm{cos}\alpha +g\left[\mathrm{sin}\left(\alpha +\omega t-{\phi}_{0}\right)-f\mathrm{cos}\left(\alpha +\omega t-{\phi}_{0}\right)\right]\end{array}$ (6)

Equation (5) is integrated under zero initial conditions:

$r=0;\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{\dot{}}{r}=0$ при $t=0$ (7)

The solution of Equation (6) is presented in the form:

$r={A}_{1}{\text{e}}^{{\lambda}_{1}\omega t}+{A}_{2}{\text{e}}^{{\lambda}_{2}\omega t}+{A}_{0}+{B}_{0}\mathrm{sin}\left(\omega t+\alpha -{\phi}_{0}\right)+{C}_{0}\mathrm{cos}\left(\omega t+\alpha -{\phi}_{0}\right)$ (8)

where ${\lambda}_{1}=\sqrt{{f}^{2}{\mathrm{cos}}^{2}\alpha +1}-f\mathrm{cos}\alpha $, ${\lambda}_{2}=-\sqrt{{f}^{2}{\mathrm{cos}}^{2}\alpha +1}-f\mathrm{cos}\alpha $.

We substitute (7) into Equation (6) and equate the coefficients with constant and functions $\mathrm{sin}\left(\alpha +\omega t-{\phi}_{0}\right)$ and $\mathrm{cos}\left(\alpha +\omega t-{\phi}_{0}\right)$ to the right and left of the equality, then we get

${A}_{0}=-R\mathrm{cos}\alpha $, ${B}_{0}=-\frac{g\left(1-{f}^{2}\mathrm{cos}\alpha \right)}{2{\omega}^{2}\left(1+{f}^{2}{\mathrm{cos}}^{2}\alpha \right)}$, ${C}_{0}=-\frac{gf\left(1+\mathrm{cos}\alpha \right)}{2{\omega}^{2}\left(1+{f}^{2}{\mathrm{cos}}^{2}\alpha \right)}$,

To determine the constants A_{1} and A_{2}, we use the initial conditions (7), which give:

$\begin{array}{c}{A}_{1}+{A}_{2}=-{A}_{0}-{B}_{0}\mathrm{sin}\left(\alpha -{\phi}_{0}\right)-{C}_{0}\mathrm{cos}\left(\alpha -{\phi}_{0}\right)+{\lambda}_{1}\omega {A}_{1}+{\lambda}_{2}\omega {A}_{2}\\ =-{B}_{0}\omega \mathrm{cos}\left(\alpha -{\phi}_{0}\right)+{C}_{0}\mathrm{sin}\left(\alpha -{\phi}_{0}\right)\end{array}$

From the last system we find the constants A_{1} and A_{2}:

${A}_{1}=\frac{{\lambda}_{2}{A}_{0}+\left({\lambda}_{2}{B}_{0}+{\u0421}_{0}\right)\mathrm{sin}\left(\alpha -{\phi}_{0}\right)+\left({\lambda}_{2}{\u0421}_{0}-{B}_{0}\right)\mathrm{cos}\left(\alpha -{\phi}_{0}\right)}{{\lambda}_{1}-{\lambda}_{2}}$

${A}_{2}=\frac{{\lambda}_{1}{A}_{0}+\left({\lambda}_{1}{B}_{0}+{\u0421}_{0}\right)\mathrm{sin}\left(\alpha -{\phi}_{0}\right)+\left({\lambda}_{1}{\u0421}_{0}-{B}_{0}\right)\mathrm{cos}\left(\alpha -{\phi}_{0}\right)}{{\lambda}_{2}-{\lambda}_{1}}$

The movement of the fiber along the BM line starts from the moment it is captured by the front face of the tooth, if the condition

$\stackrel{\xa8}{r}>0$ at $t=0$ putting in the Equation (6) z, we find

${\omega}^{2}R\mathrm{cos}\alpha +g\left[\mathrm{sin}\left(\alpha -{\phi}_{0}\right)-f\mathrm{cos}\left(\alpha -{\phi}_{0}\right)\right]\ge 0$

or

${\omega}^{2}\mathrm{cos}\alpha +g\mathrm{cos}\alpha \left[\text{tg}\left(\alpha -{\phi}_{0}\right)-f\right]\ge 0$ (9)

If you assume

$\alpha >{\phi}_{0}$, $\alpha <\frac{\pi}{2}$, $\alpha -{\phi}_{0}<\frac{\pi}{2}$,

then the last inequality is always satisfied, if required $\text{tg}\left(\alpha -{\phi}_{0}\right)>f$ or $\alpha >{\phi}_{0}+\text{arctg}\text{\hspace{0.05em}}f$. If $\text{tg}\left(\alpha -{\phi}_{0}\right)<f$, or ${\phi}_{0}<\alpha <{\phi}_{0}+\text{arctg}\text{\hspace{0.05em}}f$, then the condition of the drum speed satisfies the inequality

${\omega}^{2}>\frac{g\mathrm{cos}\left(\alpha -{\phi}_{0}\right)}{R\mathrm{cos}\alpha}\left(f-\text{tg}\left(\alpha -{\phi}_{0}\right)\right)$ (10)

Let now

$\alpha <{\phi}_{0}$, $\alpha <\frac{\pi}{2}$, ${\phi}_{0}-\alpha <\frac{\pi}{2}$.

Then we rewrite inequality (9):

${\omega}^{2}\mathrm{cos}\alpha -g\mathrm{cos}\left({\phi}_{0}-\alpha \right)\left[f+\text{tg}\left({\phi}_{0}-\alpha \right)\right]\ge 0$

From the last inequality we have

${\omega}^{2}>\frac{g\mathrm{cos}\left({\phi}_{0}-\alpha \right)\left(f+\text{tg}\left({\phi}_{0}\alpha \right)\right)}{R\mathrm{cos}\alpha}$ (11)

If the angular velocity of the drum ω does not satisfy inequalities (10) and (11), then the movement of the fiber along the face of the tooth begins from the time determined from the equation

${\omega}^{2}R\mathrm{cos}\alpha +g\left[\mathrm{sin}\left(\alpha +\omega {t}_{0}-{\phi}_{0}\right)-f\mathrm{cos}\left(\alpha +\omega {t}_{0}-{\phi}_{0}\right)\right]=0$

where, depending on the values ω, α or φ_{0} we find
${t}_{0}=0$ if
${\omega}^{2}\succ \frac{g\sqrt{1+{f}^{2}}}{R\mathrm{cos}\alpha}$,

${t}_{0}=\frac{1}{\omega}\left[{\beta}_{0}+{\phi}_{0}-\alpha -\mathrm{arcsin}\frac{{\omega}^{2}R\mathrm{cos}\alpha}{g\sqrt{1+{f}^{2}}}\right]$

at

${\omega}^{2}<\frac{g\sqrt{1+{f}^{2}}}{R\mathrm{cos}\alpha}$ $\alpha <{\phi}_{0}+\beta -\mathrm{arcsin}\frac{{\omega}^{2}R\mathrm{cos}\alpha}{g\sqrt{1+{f}^{2}}}$

where, ${\beta}_{0}=\mathrm{arcsin}\frac{1}{\sqrt{1+{f}^{2}}}$.

If we denote by t_{k} the time spent by the fiber in the sampling zone, then the fiber moves along the effects of the tooth, if the conditions are met:

$r\left({t}_{k}\right)\le {r}_{0}$ (12)

where ${r}_{0}$ -the length of the front edge of the tooth.

If condition (12) is not satisfied, then the fiber can leave the tooth surface before it enters the transportation zone.

Now we simulate the process of dragging the fiber with the teeth of the headset in the sampling zone. We accept that fibers as a deformable thread of finite length have an arcuate shape, the surface of which interacts with the surrounding pulp. When leaving the feed zone, all fiber points have a speed and a tension of T_{0}.

We accept that the fiber is deformed according to Hooke’s law:

$T=ES\cdot \epsilon $ (13)

where Е-modulus, S-cross-sectional area of the fiber, is the relative deformation. If we denote by V_{0} and V the initial and current volume of the fiber, then the relative change in its volume will be equal

${\epsilon}_{\vartheta}=\frac{V-{V}_{0}}{{V}_{0}}$ (14)

Moreover

${\epsilon}_{\vartheta}={\epsilon}_{1}+{\epsilon}_{2}+{\epsilon}_{3}$

where ${\epsilon}_{1},{\epsilon}_{2},{\epsilon}_{3}$ -changes in the size of the thread in the axial and transverse directions. If we consider a thread of cylindrical shape, то ${\epsilon}_{2}={\epsilon}_{3}=\mu \cdot \epsilon $ $\mu $ - Poisson. In view of the last equality, we have

${\epsilon}_{\vartheta}=\left(1+2\mu \right)\epsilon $

Then from Formula (14) we find

$\stackrel{\xa8}{\epsilon}=\frac{1}{1+2\mu}\cdot \frac{{V}_{0}-V}{{V}_{0}}$.

On the other hand

${V}_{0}=\frac{1}{{\rho}_{0}}$, $V=\frac{1}{\rho}$,

where ${\rho}_{0}$ and $\rho $ -initial (when leaving the feed zone) and the current density of a single fiber, thus we have

$\epsilon =\frac{1}{1+2\mu}\cdot \frac{-\frac{1}{{\rho}_{0}}+\frac{1}{\rho}}{\frac{1}{{\rho}_{0}}}=\frac{{\rho}_{0}-\rho}{\rho}=\left(\frac{{\rho}_{0}}{\rho}-1\right)\frac{1}{1+2\mu}$

From the last equality we express the change in density through deformation $\frac{{\rho}_{0}}{\rho}=\left(1+2\mu \right)\epsilon +1$.

Formula (13) is written as

$T=\frac{E}{1+2\mu}\left(\frac{{\rho}_{0}}{\rho}-1\right)S$ (15)

The process of dragging the fiber is considered stationary; in this case, the law of conservation of mass

${\rho}_{0}{S}_{0}{\vartheta}_{0}=\rho S\vartheta $ (16)

Here S_{0} and S fiber cross section,
$\vartheta $ -section speed. From the Equation (15) taking into account (16), half

$\frac{\vartheta}{{\vartheta}_{0}}=\frac{{\rho}_{0}}{\rho}\cdot \frac{{S}_{0}}{S}=\frac{{S}_{0}}{S}\left[1+\frac{T\left(1+2\mu \right)}{ES}\right]$

In the future we get S_{0} and S we have put

$\vartheta ={\vartheta}_{0}\left[1+\frac{T\left(1+2V\right)}{E{S}_{0}}\right]$ (17)

Accordingly, for density we have:

$\frac{\rho}{{\rho}_{0}}=\frac{{\vartheta}_{0}}{\vartheta}=\frac{1}{1+\frac{T\left(1+2\mu \right)}{E{S}_{0}}}$ (18)

Formulas (17) and (18) establish the relationship between the tension of the thread and the speed as well as the density in an arbitrary section of the fiber.

Let us create the equation of the balance of the thread taking into account the resistance of its environment [10] [11]:

$\frac{\text{d}T}{\text{d}l}-\tau =0$, $\frac{T}{R}-q=0$ (19)

where, l- the length of the arc to be reported from the point of the fiber exit from the feed zone; $\tau $ and q intensity of tangential normal forces on the surface of the thread.: We accept, then the values $\tau $ and q are bound according to the Kulon law by dependence $\tau =f\cdot q$ (f-coefficient of friction). With this in mind...

$\frac{\text{d}T}{\text{d}l}-f\frac{T}{R}=0$, $0\le l\le \omega t\cdot R$ (20)

where $\omega $ and R the angular velocity and radius of the drum. We introduce the moving coordinate $x=l-\omega t\cdot R$ and integrate the Equation (20).

$T=A{\text{e}}^{\frac{fl}{R}}$ (21)

where A is an arbitrary constant determined from the condition of capturing (dragging) the fiber by the teeth of the headset.

$\vartheta =\omega \cdot R$ при $x=0$

Then from (2.23) we find the tension value in the section $S=0$

${T}_{0}=\left(\frac{\omega R}{{\vartheta}_{0}}-1\right)\frac{E{S}_{0}}{1+2\mu}$ (22)

Taking into account (22), the Formula (21) takes the form

$T=\left(\frac{\omega R}{{V}_{0}}-1\right)\frac{E{S}_{0}}{1+2V}{\text{e}}^{f\frac{l-\omega t}{R}}$

The speed and density in any section is calculated by formulas:

$V={V}_{0}\left[1+\left(\frac{\omega R}{{\vartheta}_{0}}-1\right){\text{e}}^{f\frac{l-\omega t}{R}}\right]$

$\rho ={\rho}_{0}\frac{1}{1+\left(\frac{\omega R}{{\vartheta}_{0}}-1\right){\text{e}}^{f\frac{l-\omega t}{R}}}$

3. Analysis of the Results

The obtained equations are processed on a computer using the Maple program with the natural parameters of technological machines and graphical dependencies are obtained fiber movement when moving along the tooth of the headset at different angles from time to time (Figure 2); the distribution of fiber tension along the arc of the sampling zone at different speeds of rotation of the drum (Figure 3), the distribution of fiber speed along the arc of the sampling zone at different speeds of rotation of the drum (Figure 4), according to which it can be noted that all the dependencies have an upward character of change .

The dependence of the movement of the fiber when moving along the tooth of the headset at different angles from time to time (Figure 2) also shows that, in hundredths of a second (0.001 sec), the movement of the fiber does not exceed 3 - 7 mm at angles- $\alpha ={10}^{\circ}$, $\alpha ={30}^{\circ}$, $\alpha ={45}^{\circ}$, $\alpha ={60}^{\circ}$.

A, the dependence of the distribution of fiber tension along the arc of the sampling zone at different speeds of rotation of the drum (Figure 3) shows that at angular speeds 600 - 750 sek^{−1}, a change in distance from 0 to 0.1 m leads to an increase in fiber tension in the range 19 - 24 cN.

Also, the dependence of the distribution of the fiber velocity along the arc of the sampling zone at different speeds of rotation of the drum (Figure 4) shows that at angular speeds of 600 - 750 sec^{−1}, a change in the distance from 0 to 0.1 m contributes to an increase in the fiber speed from 37 - 49 m/s.

Figure 2. The time-dependant $t\left(\text{\u0441\u0435\u043a}\right)$ movement of the fiber $r\left(\text{\u043c\u043c}\right)$ along the headset tooth at different angles $\alpha $ 1- $\alpha ={10}^{\circ}$, 2- $\alpha ={30}^{\circ}$, 3- $\alpha ={45}^{\circ}$, 4- $\alpha ={60}^{\circ}$.

Figure 3. Distribution of fiber tension $\u0422\left(\text{\u0441\u041d}\right)$ along the arc of the sampling zone at different rotational speeds of the drum $\omega \left(1/\mathrm{sec}\right)$ 1- $\omega =600$, 2- $\omega =650$, 3- $\omega =700$, 4- $\omega =750$.

Figure 4. Distribution of fiber velocity $v\left(\text{\u043c}/\text{c}\right)$ along the arc of the sampling zone at different rotational speeds of the drum (1/sec). 1- $\omega =600$, 2- $\omega =650$, 3- $\omega =700$, 4- $\omega =750$.

4. Conclusions

1) The regularity was studied and the equation of fiber motion on the tooth surface of the headset of the sampling drum was obtained.

2) As a result of the analysis of the tension and speed of the fibers, it was found that with an increase in the sampling zone, the tension force and the speed of the fiber increase.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

[1] | Roglin, B., Bouzhek, A., Wilfert, M., et al. (1981) Spinning-Free Spinning. Light and Food Industry, 86-121. |

[2] | Yuldashev, Zh. (2011) Improving the Quality of Pneumomechanical Yarn Using Multi-Thread Winding of the Headset. Dissertation of the Candidate of Technical Sciences, Tashkent Institute of Textile and Light Industry, Tashkent, Uzbekistan, 52-85. |

[3] |
Abdukarimovich, M.O., Ibragimovich, A.K. and Sharipjanovich, S.O. (2018) Designing a New Design of a Loading Cylinder for Pneumomechanical Spinning Machines. Engineering, 10, 345-356. https://doi.org/10.4236/eng.2018.106025 |

[4] | Abbazov, I., Sarimsakov, O., Khodjiev, M. and Mardonov, B. (2018) Waste Produced at Cotton Waste Factories. American Journal ASCIT Communications, 5, 22-28. |

[5] | Kholmirzaev, F., Azimov, S., Abdurahimov, K. and Sarimsakov, O. (2019) Investigation of the Loss of Air Pressure in the Pipeline of the Сotton Pneumatic Conveying. Saudi Journal of Engineering and Technology, 4, 23-27. |

[6] | Muradov, R., Burnashev, R.Z. and Sarimsakov, O.Sh. (2001) The Dynamic Task of the Interaction of Raw Cotton with the Working Bodies of Processing Machines. Journal of Mechanic Problems, No. 34, 55-58. |

[7] |
Obidov, A., Akhmedkhodjaev, Kh., Sarimsakov, O. and Holikov, Q. (2018) Investigation of the Properties of Fibrous Cotton Seeds, for Sorting on a Mesh Surface. Engineering, 10, 572-578. https://doi.org/10.4236/eng.2018.109041 |

[8] | Sarimsakov, O. (2016) The Possibility of Reducing Cotton Consumption in Cotton. American Journal of Science and Technology, 4, 68-72. |

[9] | Sarimsakov, O. (2019) Improving the Process of Filing Cotton in Pneumatic Transport and its Transportation by Air. Monography, Namangan, Uzbekistan, p. 252. |

[10] | Sarimsakov, O. and Gaybnazarov, E. (2017) About Energy Consumption in Pneumatic Conveying of Raw Cotton. American Journal of Energy and Power Engineering, 3, 26-29. |

[11] | Sarimsakov, O., Xusanov, C. and Muradov, R. (2016) The Change in Air Pressure along the Length of the Pipeline Installation for Pneumatic Conveying of Raw Cotton. Journal of Engineering and Technology, 3, 89-92. |

Journals Menu

Copyright © 2021 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.