Identification of Robot Arm s Joints Time-Varying Stiffness Under Loads
|
|
- Edith McCarthy
- 5 years ago
- Views:
Transcription
1 TELKOMNIKA, Vol.10, No.8, December 2012, pp. 2081~2087 e-issn: X ccredted by DGHE (DIKTI), Decree No: 51/Dkt/Kep/ Identfcton of Robot Arm s Jonts Tme-Vryng Stffness Under Lods Ru Xu 1, Qng Chen* 2, Guol Yng 1 1 School of Mechncl Engneerng, Nnjng Unversty of Scence nd Technology 2 School of Mngement nd Engneerng, Nnjng Unversty No. 22, Hnkou Rod, Nnjng, Chn, , +86(25) *correspondng uthor, e-ml: chenq0772@vp.sn.com Abstrct Ths pper s concerned wth the dentfcton problem of two degree of freedom robot rm s jonts tme-vryng stffness under tme-vryng torque lods. The ndustrl robot rms cn be smplfed nto seres hnged lnkge mechnsm. The unrestrnt dynmc equton of two degrees of freedom robot rm cn be obtned by usng nlytcl mechncs method. The model whch ws estblshed n Pro/E hd been mported nto ADAMS for smulton. And the smulton conssted of free vbrton nd forced vbrton. Then by choosng lmted memory lest squre method nd usng the post-processng mesured vrbles, tme-vryng stffness could hve been dentfed. Fnlly, the clcultve stffness ws compred to the rel stffness whch ws smulted n ADAMS. The whole process shows tht the robot rm s dynmc model nd the method of dentfcton re both effectve. Keywords: tme-vryng stffness, robot rm, dentfcton, lmted memory lest squre method Copyrght 2012 Unversts Ahmd Dhln. All rghts reserved. 1. Introducton As becomng new knd of producton pttern, robot rm hs nfluenced the ndustrl producton lot. Most of the robot rms consst of two prts: motor nd reducer. For ts complex structure, the robot rm s jonts stffness s tme-vryng. Therefore, t s dffcult but necessry to dentfy some prmeters of the robot rm s dynmc model whch s the key problem to estblsh ths model. Becuse of lest squre method s smple prncple nd fst convergence, t s wdely used n dentfyng system prmeters. As one of lest squre method s mproved lgorthms, lmted memory lest squre method cn effectvely overcome dt sturton phenomenon when t s ppled to tme-vryng dentfcton. Xnzh C nd Chong C proved ts effect n ther pper [1]. Xopng L, Dh J nd We Z provde mny dentfyng methods of tme-vryng prmeters n ther ppers [2-4]. But, ther robot rm s dynmc equtons lck physcl sgnfcnce nd t s hrd to exmne the result of dentfcton. In ths pper, tkng two degree of freedom robot rm s the reserch object, we cn dentfy the tme-vryng stffness effectvely when combnng robot rm dynmc model wth lmted memory lest squre method. Becuse the robot rm s stffness hs physcl sgnfcnce, the result of dentfcton cn be compred to the sttc mesurement n order to exmne the method s effect on tme-vryng prmeters. 2. Problem Formulton 2.1. Dynmc equtons Anlytcl mechncs use nlytcl mth method to solve mechncl problems. By nlyzng from energy nd work, the reltonshp between moton nd ntercton cn be found. After mkng use of d lembert s prncple nd vrtul dsplcement prncple, the generl dynmcs equtons re estblshed, nd we substtute Lgrnge multpler n order to get unrestrnt Lgrnge equtons [5-6]. Unrestrnt Lgrnge equtons re shown n Eq. 1: d T T V dt q q q ( 1,2,, ) Q N + = = & (1) Receved September 29, 2012; Revsed November 5, 2012; Accepted November 13, 2012
2 2082 e-issn: X Where T s the knetc energy of system, V s the potentl energy of system nd conssts of grvttonl potentl energy nd elstc potentl energy, Q s the work cused by resultnt of non-conservtve force cused by electrc motors, q s generlzed dsplcement nd N s system s degree of freedom. The structure of robot rm s jonts s smplfed nto hnges wth nonlner torson sprngs. So tht, the potentl energy of system whch s cused by these torson sprngs nd the power of grvty cn be expressed n Eq. 2: ( 1, 2) V = K + P = (2) Where K s the elstc potentl energy of the system, P s the grvttonl potentl energy of system. Then we substtute Eq. 2 nto Eq. 1, so robot rm s dynmc equton cn be obtned n Eq. 3: ( K P) d T T + dt q q q ( 1,2 ) + = Q = & The structure of two degree of freedom robot rm s llustrted n Fgure 1: Where φ 1 s the bg rm s generlzed dsplcement, φ 2 s the smll rm s generlzed dsplcement, l 1 s the bg rm s length, l 2 s the smll rm s length, l 1C s the dstnce between the center of mss nd the bg rm s jont, l 2C s the dstnce between the center of mss nd the smll rm s jont, J 1C s bg rm s moment of nert n center of mss, J 2C s the smll rm s moment of nert n center of mss, m 1 s the bg rm s mss, m 2 s the smll rm s mss. Set the ground s the zero potentl energy locton. After beng substtuted ll these prmeters nto Eq. 3, two degree of freedom robot rm s dynmc equtons cn be expressed n Eq. 4 nd Eq. 5: (3) ( J + m l + m l )&& ϕ + m l l && ϕ cos( ϕ ϕ ) m l l & ϕ sn( ϕ ϕ ) C 1 1C C C k ϕ k ϕ ϕ + m gl cosϕ + m gl cosϕ = Q C ( J + m l )&& ϕ + m l l && ϕ cos( ϕ ϕ ) + m l l & ϕ sn( ϕ ϕ ) 2 2 2C 2 2C C C k ϕ ϕ + m gl cosϕ = Q C 2 2 Where k 1 s the bg rm s jont s stffness, k 2 s the smll rm s jont s stffness. (4) (5) Fgure 1. Structure of two degree of freedom robot rm Fgure 2. Mult-rgd-body system s dynmc model 2.2. Dynmc smulton The robot rms cn be smplfed nto seres hnged lnkge mechnsm. Then the 3D model cn be estblshed n Creo Prmetrc. By usng the functon of mss nlyss n Creo Prmetrc, both the mss, moment of nert nd center of mss cn be obtned n Eq. 6: TELKOMNIKA Vol. 10, No. 8, December 2012 :
3 TELKOMNIKA e-issn: X m1 = kg, J1 C = kg/m, l1 C = m 9 2 m2 = kg, J 2C = kg/m, l2c = 0.15 m (6) Then we mport ths model nto ADAMS for smulton. The mult-rgd-body system s dynmc model s bult nd shown n Fgure 2: In Fgure 2, two degree of freedom robot rm conssts of foundton support, bg rm, bg xs, smll rm nd smll xs. All these subssembles re supposed to be mde from 6061 lumnum. The foundton support s fxed on the ground. The bg rm s hnged on the foundton support by the bg xs, nd the bg torson sprng s loded. The smll rm s hnged on the bg rm by the smll xs, nd the smll torson sprng s loded. The bg motor s fxed on bg rm, nd the smll motor s fxed on smll rm. These motors mss wll be equvlent to the robot rms. In order to dentfy the tme-vryng stffness of torson sprngs, the prmeters of torque sprngs re set up s splnes tht descrbe the reltonshp between torson nd deflecton, nd the prmeters of electrc motors re set up s functons tht descrbe the reltonshp between torson nd tme. These reltonshps cn be expressed n Eq. 7: T torque T = f ( t ) motor = g defo Where T torque s the torque cused by torson sprngs, defo s the jonts deflectons, T motor s the torque cused by electrc motors, t s the tme. The splnes tht descrbe the tme-vryng stffness of torson sprngs re llustrted n Fgure 3 nd Fgure 4. The functons splnes tht descrbe the tme-vryng torque cused by electrc motors re llustrted n Fgure 5 nd Fgure 6. (7) Fgure 3. Splne descrbng the bg torson sprng Fgure 4. Splne descrbng the smll torson sprng Fgure 5. Splne descrbng the bg torque cused by motor Fgure 6. Splne descrbng the smll torque cused by motor Identfcton of Robot Arm s Jonts Tme-Vryng Stffness Under Lods (Ru Xu)
4 2084 e-issn: X The dynmc smulton conssts of two dfferent prts: free vbrton nd forced vbrton. Free vbrton: we ssume tht there s no grvty nd no generlzed force. Intl compresson of two torson sprngs s set up s 0.5 degree (vector) from the blnce poston. So tht the robot rm system vbrtes freely. The smulton tme s set s second nd the whole process conssts of 8000 steps. By ADAMS post-processng progrm, the dt of dsplcement, velocty, ccelerton nd vbrton cn be obtned. Then record ll these smulton dt for dentfcton clculton. Forced vbrton: the drecton vector of grvty s set perpendculr to the ground. Intl compresson of two torson sprngs s set up s zero from the blnce poston. The smulton tme s set s second nd the whole process conssts of 8000 steps, too. Smlrly, the dt of dsplcement, velocty, ccelerton, vbrton nd cn be obtned, nd the tme-vryng torque s well. Then record ll these smulton dt for dentfcton clculton Identfcton method When beng dentfed the system prmeters, f the object cn be descrbed s the controlled uto regressve model, the model cn be expressed n Eq. 8: 1 1 ξ A z y k = B z y k d + k (8) Where ξ(k) s whte nose, nd the structure prmeters m, n, d re known. And: A z = 1+ z + z + + z B z = b + z + z + + z m 1 2 m n n The trget of dentfcton s to estblsh m+n+1prmeters ccordng to the mesurble nput nd output. The model cn be descrbed n lest squre method nd s expressed n Eq. 10: = 1 ( 1) n ( ) n b b + ξ T = ϕ ( k ) θ + ξ ( k ) y k y k y k n b u k d b u k d n k Where φ(k) s the dt vector, θ s the prmeter need to be dentfed. And: ( 1 ),,,,, ϕ k = y k y k n u k d u k d nb R T ( n + nb + 1) 1 θ 1,, n, b0,, b = n R b T ( n + n + 1) 1 So the model of two degree of freedom robot rm s jonts cn be expressed n Eq. 12 nd Eq. 13: b (9) (10) (11) k1 2 2 ϕ1 ϕ2 + ϕ1 = ( J1 C + m1l 1C + m2l1 ) ϕ1 m2l1l 2Cϕ2 cos( ϕ2 ϕ1) k && && 2 [ ] t 2 & 2 1 2Cϕ2 sn( ϕ2 ϕ1) 1 1C cosϕ1 2 1 cos ϕ1 ( motor1 ϕ1) ϕ1 0 + m l l m gl m gl + d T d d k1 2 0 ϕ2 ϕ1 = ( J 2C + m2l2c ) ϕ2 m2l1l 2Cϕ1 cos( ϕ2 ϕ1) k && && 2 [ ] t 2 & 2 1 2Cϕ1 sn( ϕ2 ϕ1) 2 2C cos ϕ2 ( motor 2 ϕ2 ) ϕ2 0 m l l m gl + d T d d (12) (13) In ths pper, the lmted memory lest squre method s used to dentfy the stffness. The bsc de of ths method s usng the ltest perod of dt to dentfy the stffness. Once ddng set of dt, nother set should be declned so s to keep the number of dt constnt. Ths method reflects the tme-vryng chrcterstc of prmeters, nd cn solve dt sturton problem. TELKOMNIKA Vol. 10, No. 8, December 2012 :
5 TELKOMNIKA e-issn: X Numercl clculton nd results The source of dt dpted n the dentfcton s from the robot rm s stffness whch s cqured by numercl methods. A reltve percentge error (RPE) s defned to compre the result of the dentfcton whch s expressed n Eq. 14: Σ ktrue kdentfy RPE = 100% Σ k true (14) 3.1. Free vbrton The results of dentfcton bout rms stffness wthout lods re llustrted n Fgure 7 nd Fgure 8. The results of dentfcton bout reltve error of rms stffness wthout lods re llustrted n Fgure 9 nd Fgure 10. Fgure 7. Identfcton of bg rm s stffness wthout lods Fgure 8. Identfcton of smll rm s stffness wthout lods Fgure 9. Reltve error of bg rm s stffness wthout lods Fgure 10. Reltve error of smll rm s stffness wthout lods So, the RPE of bg rm wthout lods s bout 1.093%, whle the RPE of smll rm wthout lods s bout 1.114% Forced vbrton The results of dentfcton bout rms stffness under lods re llustrted n Fgure 11 nd Fgure 12. The results of dentfcton bout reltve error of rms stffness under lods re llustrted n Fgure 13 nd Fgure 14. Identfcton of Robot Arm s Jonts Tme-Vryng Stffness Under Lods (Ru Xu)
6 2086 e-issn: X Fgure 11. Identfcton of bg rm s stffness under lods Fgure 12. Identfcton of smll rm s stffness under lods Fgure 13. Reltve error of bg rm s stffness under lods Fgure 14. Reltve error of smll rm s stffness under lods So, the RPE of bg rm under lods s bout 1.021%, whle the RPE of smll rm under lods s bout 0.995%. 4. Conclusons In the pper, the unrestrnt Lgrnge equtons nd lmted memory lest squre method s ppled to dentfy the two degree of freedom robot rm s jonts stffness. The robot rm s stffness hs physcl sgnfcnce so tht the result of dentfcton cn be compred to the sttc mesurement whch s obtned by ADAMS. Lmted memory lest squre method cn overcome dt sturton phenomenon effectvely. So t cn be ppled to tme-vryng dentfcton n ths stuton. The result reflects tht lmted memory lest squre method nd the process of estblshng robot rm s dynmc model re effectve n both free vbrton nd forced vbrton. Future work of ths study would be dentfyng rel prmeters nsted of smulton nd exmnng ths reserch method s performnce n other models. Acknowledgements Ths work ws supported n prt by the Mjor Stte Bsc Reserch Development Progrm of Republc of Chn (No ) nd the Fundmentl Reserch Funds for the Centrl Unverstes. TELKOMNIKA Vol. 10, No. 8, December 2012 :
7 TELKOMNIKA e-issn: X 2087 References [1] Xnzh C, Chong C. Prmeter dentfcton of tme-vryng systems bsed on the lest-squres modfed lgorthm. Journl of Fuzhou Unversty (Nturl Scence). 2005; 33(2): [2] Dh J, Xopng L. On-lne Fuzzy Identfyng Physcl Prmeters of Jont Movng Stte Robots. Journl of Bejng Unversty of Posts nd Telecommunctons. 2008; 31(1): [3] Dh J, Xopng L. On-lne Identfcton of Tme-vryng Physcl Prmeters of Robot Jont Bsed on Hrmonc Propgton. Journl of Mechncl Engneerng. 2009; 45(3): [4] We Z, Xopng L. Jont prmeters dentfcton of robot bsed on trvelng wve pproch nd neurl network. Journl of Vbrton Engneerng. 2010; 23(2): [5] Theoretcl Mechncs Techng nd Reserch Secton of Hrbn Insttute of Technology. Theoretcl Mechncs. Sxth Edton. Bejng: Hgher Educton Press [6] Ymng Z. Mechncl Vbrton. Frst Edton. Bejng: Tsnghu Unversty Press Identfcton of Robot Arm s Jonts Tme-Vryng Stffness Under Lods (Ru Xu)
Haddow s Experiment:
schemtc drwng of Hddow's expermentl set-up movng pston non-contctng moton sensor bems of sprng steel poston vres to djust frequences blocks of sold steel shker Hddow s Experment: terr frm Theoretcl nd
More information4. Eccentric axial loading, cross-section core
. Eccentrc xl lodng, cross-secton core Introducton We re strtng to consder more generl cse when the xl force nd bxl bendng ct smultneousl n the cross-secton of the br. B vrtue of Snt-Vennt s prncple we
More informationDefinition of Tracking
Trckng Defnton of Trckng Trckng: Generte some conclusons bout the moton of the scene, objects, or the cmer, gven sequence of mges. Knowng ths moton, predct where thngs re gong to project n the net mge,
More informationMechanical resonance theory and applications
Mechncl resonnce theor nd lctons Introducton In nture, resonnce occurs n vrous stutons In hscs, resonnce s the tendenc of sstem to oscllte wth greter mltude t some frequences thn t others htt://enwkedorg/wk/resonnce
More informationRank One Update And the Google Matrix by Al Bernstein Signal Science, LLC
Introducton Rnk One Updte And the Google Mtrx y Al Bernsten Sgnl Scence, LLC www.sgnlscence.net here re two dfferent wys to perform mtrx multplctons. he frst uses dot product formulton nd the second uses
More information6 Roots of Equations: Open Methods
HK Km Slghtly modfed 3//9, /8/6 Frstly wrtten t Mrch 5 6 Roots of Equtons: Open Methods Smple Fed-Pont Iterton Newton-Rphson Secnt Methods MATLAB Functon: fzero Polynomls Cse Study: Ppe Frcton Brcketng
More informationSolution of Tutorial 5 Drive dynamics & control
ELEC463 Unversty of New South Wles School of Electrcl Engneerng & elecommunctons ELEC463 Electrc Drve Systems Queston Motor Soluton of utorl 5 Drve dynmcs & control 500 rev/mn = 5.3 rd/s 750 rted 4.3 Nm
More informationPrinciple Component Analysis
Prncple Component Anlyss Jng Go SUNY Bufflo Why Dmensonlty Reducton? We hve too mny dmensons o reson bout or obtn nsghts from o vsulze oo much nose n the dt Need to reduce them to smller set of fctors
More informationChapter Newton-Raphson Method of Solving a Nonlinear Equation
Chpter 0.04 Newton-Rphson Method o Solvng Nonlner Equton Ater redng ths chpter, you should be ble to:. derve the Newton-Rphson method ormul,. develop the lgorthm o the Newton-Rphson method,. use the Newton-Rphson
More informationLOCAL FRACTIONAL LAPLACE SERIES EXPANSION METHOD FOR DIFFUSION EQUATION ARISING IN FRACTAL HEAT TRANSFER
Yn, S.-P.: Locl Frctonl Lplce Seres Expnson Method for Dffuson THERMAL SCIENCE, Yer 25, Vol. 9, Suppl., pp. S3-S35 S3 LOCAL FRACTIONAL LAPLACE SERIES EXPANSION METHOD FOR DIFFUSION EQUATION ARISING IN
More informationChapter Newton-Raphson Method of Solving a Nonlinear Equation
Chpter.4 Newton-Rphson Method of Solvng Nonlner Equton After redng ths chpter, you should be ble to:. derve the Newton-Rphson method formul,. develop the lgorthm of the Newton-Rphson method,. use the Newton-Rphson
More informationWork and Energy (Work Done by a Varying Force)
Lecture 1 Chpter 7 Physcs I 3.5.14 ork nd Energy (ork Done y Vryng Force) Course weste: http://fculty.uml.edu/andry_dnylov/techng/physcsi Lecture Cpture: http://echo36.uml.edu/dnylov13/physcs1fll.html
More informationWorkspace Analysis of a Novel Parallel Robot Named 3-R2H2S with Three Freedoms
Reserch Journl of Appled Scences, Engneerng nd Technology 6(0: 3847-3851, 013 ISS: 040-7459; e-iss: 040-7467 Mxwell Scentfc Orgnzton, 013 Submtted: Jnury 17, 013 Accepted: Februry, 013 Publshed: ovember
More informationGAUSS ELIMINATION. Consider the following system of algebraic linear equations
Numercl Anlyss for Engneers Germn Jordnn Unversty GAUSS ELIMINATION Consder the followng system of lgebrc lner equtons To solve the bove system usng clsscl methods, equton () s subtrcted from equton ()
More informationTwo Coefficients of the Dyson Product
Two Coeffcents of the Dyson Product rxv:07.460v mth.co 7 Nov 007 Lun Lv, Guoce Xn, nd Yue Zhou 3,,3 Center for Combntorcs, LPMC TJKLC Nnk Unversty, Tnjn 30007, P.R. Chn lvlun@cfc.nnk.edu.cn gn@nnk.edu.cn
More information829. An adaptive method for inertia force identification in cantilever under moving mass
89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,
More informationQuiz: Experimental Physics Lab-I
Mxmum Mrks: 18 Totl tme llowed: 35 mn Quz: Expermentl Physcs Lb-I Nme: Roll no: Attempt ll questons. 1. In n experment, bll of mss 100 g s dropped from heght of 65 cm nto the snd contner, the mpct s clled
More informationNUMERICAL MODELLING OF A CILIUM USING AN INTEGRAL EQUATION
NUEICAL ODELLING OF A CILIU USING AN INTEGAL EQUATION IHAI EBICAN, DANIEL IOAN Key words: Cl, Numercl nlyss, Electromgnetc feld, gnetton. The pper presents fst nd ccurte method to model the mgnetc behvour
More informationStatistics and Probability Letters
Sttstcs nd Probblty Letters 79 (2009) 105 111 Contents lsts vlble t ScenceDrect Sttstcs nd Probblty Letters journl homepge: www.elsever.com/locte/stpro Lmtng behvour of movng verge processes under ϕ-mxng
More informationINTERPOLATION(1) ELM1222 Numerical Analysis. ELM1222 Numerical Analysis Dr Muharrem Mercimek
ELM Numercl Anlss Dr Muhrrem Mercmek INTEPOLATION ELM Numercl Anlss Some of the contents re dopted from Lurene V. Fusett, Appled Numercl Anlss usng MATLAB. Prentce Hll Inc., 999 ELM Numercl Anlss Dr Muhrrem
More informationInvestigation phase in case of Bragg coupling
Journl of Th-Qr Unversty No.3 Vol.4 December/008 Investgton phse n cse of Brgg couplng Hder K. Mouhmd Deprtment of Physcs, College of Scence, Th-Qr, Unv. Mouhmd H. Abdullh Deprtment of Physcs, College
More informationLesson 2. Thermomechanical Measurements for Energy Systems (MENR) Measurements for Mechanical Systems and Production (MMER)
Lesson 2 Thermomechncl Mesurements for Energy Systems (MEN) Mesurements for Mechncl Systems nd Producton (MME) 1 A.Y. 2015-16 Zccr (no ) Del Prete A U The property A s clled: «mesurnd» the reference property
More informationInternational Journal of Pure and Applied Sciences and Technology
Int. J. Pure Appl. Sc. Technol., () (), pp. 44-49 Interntonl Journl of Pure nd Appled Scences nd Technolog ISSN 9-67 Avlle onlne t www.jopst.n Reserch Pper Numercl Soluton for Non-Lner Fredholm Integrl
More information? plate in A G in
Proble (0 ponts): The plstc block shon s bonded to rgd support nd to vertcl plte to hch 0 kp lod P s ppled. Knong tht for the plstc used G = 50 ks, deterne the deflecton of the plte. Gven: G 50 ks, P 0
More informationSameh Zribi Research Laboratory LARA in Automatic Control, ENIT, University of Tunis El Manar, Tunis, Tunisia
(JACSA) nterntonl Journl of Advnced Computer Scence Applctons, Vol. 8, No. 0, 07 mpct of Externl Dsturbnce Dscontnuous nput on e Redunt Mnpultor Robot Behvour usng e Lner Prmeter Vryng Modellng Approch
More informationPhysics 121 Sample Common Exam 2 Rev2 NOTE: ANSWERS ARE ON PAGE 7. Instructions:
Physcs 121 Smple Common Exm 2 Rev2 NOTE: ANSWERS ARE ON PAGE 7 Nme (Prnt): 4 Dgt ID: Secton: Instructons: Answer ll 27 multple choce questons. You my need to do some clculton. Answer ech queston on the
More informationTorsion, Thermal Effects and Indeterminacy
ENDS Note Set 7 F007bn orson, herml Effects nd Indetermncy Deformton n orsonlly Loded Members Ax-symmetrc cross sectons subjected to xl moment or torque wll remn plne nd undstorted. At secton, nternl torque
More informationChemical Reaction Engineering
Lecture 20 hemcl Recton Engneerng (RE) s the feld tht studes the rtes nd mechnsms of chemcl rectons nd the desgn of the rectors n whch they tke plce. Lst Lecture Energy Blnce Fundmentls F 0 E 0 F E Q W
More informationKINETICS Pipe & Duct Seismic Application Manual
KINETIC pe & Duct esmc Applcton Mnul CODE BAED EIMIC DEIGN FORCE 5.1 Introducton: The code bsed horzontl sesmc force requrements for ppe nd duct re ether clculted by the sesmc restrnt mnufcturer s prt
More informationRemember: Project Proposals are due April 11.
Bonformtcs ecture Notes Announcements Remember: Project Proposls re due Aprl. Clss 22 Aprl 4, 2002 A. Hdden Mrov Models. Defntons Emple - Consder the emple we tled bout n clss lst tme wth the cons. However,
More informationANALYTICAL AND EXPERIMENTAL STUDY OF CONICAL TELESCOPING SPRINGS WITH NON-CONSTANT PITCH
ANALYTICAL AND EXPEIMENTAL STUDY OF CONICAL TELESCOPING SPINGS WITH NON-CONSTANT PITCH Mnuel Predes Unversté de Toulouse; INSA, UPS, Mnes Alb, ISAE; ICA (Insttut Clément Ader) 5, venue de nguel, F-077
More informationStrong Gravity and the BKL Conjecture
Introducton Strong Grvty nd the BKL Conecture Dvd Slon Penn Stte October 16, 2007 Dvd Slon Strong Grvty nd the BKL Conecture Introducton Outlne The BKL Conecture Ashtekr Vrbles Ksner Sngulrty 1 Introducton
More informationApplied Statistics Qualifier Examination
Appled Sttstcs Qulfer Exmnton Qul_june_8 Fll 8 Instructons: () The exmnton contns 4 Questons. You re to nswer 3 out of 4 of them. () You my use ny books nd clss notes tht you mght fnd helpful n solvng
More informationCENTROID (AĞIRLIK MERKEZİ )
CENTOD (ĞLK MEKEZİ ) centrod s geometrcl concept rsng from prllel forces. Tus, onl prllel forces possess centrod. Centrod s tougt of s te pont were te wole wegt of pscl od or sstem of prtcles s lumped.
More informationAttribute reduction theory and approach to concept lattice
Scence n Chn Ser F Informton Scences 2005 Vol48 No6 713 726 713 Attrbute reducton theory nd pproch to concept lttce ZHANG Wenxu 1, WEI Lng 1,2 & QI Jnun 3 1 Insttute for Informton nd System Scences, Fculty
More informationSEISMIC SOIL-STRUCTURE INTERACTION ANALYSIS IN TIME-DOMAIN BASED ON COMBINATION COMPUTATION METHOD OF INTERACTION FORCE
Octoer -7, 8, Beng, Chn SEISMIC SOIL-STRUCTURE INTERACTION ANALYSIS IN TIME-DOMAIN BASED ON COMBINATION COMPUTATION METHOD OF INTERACTION FORCE Jn-feng Zho nd Xu-l Du Lecture, Dept.of Cvl Engneerng, Qngdo
More informationChemical Reaction Engineering
Lecture 20 hemcl Recton Engneerng (RE) s the feld tht studes the rtes nd mechnsms of chemcl rectons nd the desgn of the rectors n whch they tke plce. Lst Lecture Energy Blnce Fundmentls F E F E + Q! 0
More informationThe Schur-Cohn Algorithm
Modelng, Estmton nd Otml Flterng n Sgnl Processng Mohmed Njm Coyrght 8, ISTE Ltd. Aendx F The Schur-Cohn Algorthm In ths endx, our m s to resent the Schur-Cohn lgorthm [] whch s often used s crteron for
More informationDCDM BUSINESS SCHOOL NUMERICAL METHODS (COS 233-8) Solutions to Assignment 3. x f(x)
DCDM BUSINESS SCHOOL NUMEICAL METHODS (COS -8) Solutons to Assgnment Queston Consder the followng dt: 5 f() 8 7 5 () Set up dfference tble through fourth dfferences. (b) Wht s the mnmum degree tht n nterpoltng
More informationJens Siebel (University of Applied Sciences Kaiserslautern) An Interactive Introduction to Complex Numbers
Jens Sebel (Unversty of Appled Scences Kserslutern) An Interctve Introducton to Complex Numbers 1. Introducton We know tht some polynoml equtons do not hve ny solutons on R/. Exmple 1.1: Solve x + 1= for
More informationPerformance analysis of a novel planetary speed increaser used in single-rotor wind turbines with counter-rotating electric generator
OP Conference Seres: Mterls Scence nd Engneerng PAPE OPEN ACCESS Performnce nlyss of novel plnetry speed ncreser used n sngle-rotor wnd turbnes wth counter-rottng electrc genertor elted content - Fult
More informationCALIBRATION OF SMALL AREA ESTIMATES IN BUSINESS SURVEYS
CALIBRATION OF SMALL AREA ESTIMATES IN BUSINESS SURVES Rodolphe Prm, Ntle Shlomo Southmpton Sttstcl Scences Reserch Insttute Unverst of Southmpton Unted Kngdom SAE, August 20 The BLUE-ETS Project s fnnced
More informationResearch Article On the Upper Bounds of Eigenvalues for a Class of Systems of Ordinary Differential Equations with Higher Order
Hndw Publshng Corporton Interntonl Journl of Dfferentl Equtons Volume 0, Artcle ID 7703, pges do:055/0/7703 Reserch Artcle On the Upper Bounds of Egenvlues for Clss of Systems of Ordnry Dfferentl Equtons
More informationIdentification of modal parameters from nonstationary ambient vibration data using the channel-expansion technique
Journl of Mechncl Scence nd Technology 25 (5) (2011) 1307~1315 www.srngerln.com/content/1738-494x DOI 10.1007/s12206-011-0229-y Identfcton of modl rmeters from nonsttonry mbent vbrton dt usng the chnnel-exnson
More informationTHE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD
Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS
More informationLinear and Nonlinear Optimization
Lner nd Nonlner Optmzton Ynyu Ye Deprtment of Mngement Scence nd Engneerng Stnford Unversty Stnford, CA 9430, U.S.A. http://www.stnford.edu/~yyye http://www.stnford.edu/clss/msnde/ Ynyu Ye, Stnford, MS&E
More informationLecture 4: Piecewise Cubic Interpolation
Lecture notes on Vrtonl nd Approxmte Methods n Appled Mthemtcs - A Perce UBC Lecture 4: Pecewse Cubc Interpolton Compled 6 August 7 In ths lecture we consder pecewse cubc nterpolton n whch cubc polynoml
More informationPlease initial the statement below to show that you have read it
EN0: Structural nalyss Exam I Wednesday, March 2, 2005 Dvson of Engneerng rown Unversty NME: General Instructons No collaboraton of any nd s permtted on ths examnaton. You may consult your own wrtten lecture
More informationSmart Motorways HADECS 3 and what it means for your drivers
Vehcle Rentl Smrt Motorwys HADECS 3 nd wht t mens for your drvers Vehcle Rentl Smrt Motorwys HADECS 3 nd wht t mens for your drvers You my hve seen some news rtcles bout the ntroducton of Hghwys Englnd
More informationPhysics for Scientists and Engineers I
Phscs for Scentsts nd Engneers I PHY 48, Secton 4 Dr. Betr Roldán Cuen Unverst of Centrl Flord, Phscs Deprtment, Orlndo, FL Chpter - Introducton I. Generl II. Interntonl Sstem of Unts III. Converson of
More informationModel Fitting and Robust Regression Methods
Dertment o Comuter Engneerng Unverst o Clorn t Snt Cruz Model Fttng nd Robust Regresson Methods CMPE 64: Imge Anlss nd Comuter Vson H o Fttng lnes nd ellses to mge dt Dertment o Comuter Engneerng Unverst
More informationAnalytical and Simulation Method Validation for Flutter Aeroservoelasticity Studies
Vldton for Flutter Aeroservoelstcty Studes Ruxndr Mhel Botez, Professor, PhD, RTO Member Teodor Lucn Grgore, Postdoctorl fellow Adrn Hlut, PhD ETS, Unversty of Québec Lbortory of Actve Controls, Avoncs
More informationMany-Body Calculations of the Isotope Shift
Mny-Body Clcultons of the Isotope Shft W. R. Johnson Mrch 11, 1 1 Introducton Atomc energy levels re commonly evluted ssumng tht the nucler mss s nfnte. In ths report, we consder correctons to tomc levels
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationLumped Parameter Models Representing Impedance. Functions at the Interface of a Rod on a Viscoelastic Medium
Lumped Prmeter Models Representng Impednce Functons t the Interfce of Rod on Vscoelstc Medum Msto Stoh Grdute School of Scence nd Engneerng, Stm Unversty 55 Shmo-Oubo, Sur-Ku, Stm, Jpn Tel.: +8-48-858-356
More informationANALOG CIRCUIT SIMULATION BY STATE VARIABLE METHOD
U.P.B. Sc. Bull., Seres C, Vol. 77, Iss., 25 ISSN 226-5 ANAOG CIRCUIT SIMUATION BY STATE VARIABE METHOD Rodc VOICUESCU, Mh IORDACHE 22 An nlog crcut smulton method, bsed on the stte euton pproch, s presented.
More informationState Estimation in TPN and PPN Guidance Laws by Using Unscented and Extended Kalman Filters
Stte Estmton n PN nd PPN Gudnce Lws by Usng Unscented nd Extended Klmn Flters S.H. oospour*, S. oospour**, mostf.sdollh*** Fculty of Electrcl nd Computer Engneerng, Unversty of brz, brz, Irn, *s.h.moospour@gml.com
More informationRobust algorithm for calibration of robotic manipulator model
Robust lgorthm for clbrton of robotc mnpultor model Alexndr Klmchk, Yer Wu, Gbrel Abb, Benoît Furet, Antol Pshkevch o cte ths verson: Alexndr Klmchk, Yer Wu, Gbrel Abb, Benoît Furet, Antol Pshkevch. Robust
More informationSolubilities and Thermodynamic Properties of SO 2 in Ionic
Solubltes nd Therodync Propertes of SO n Ionc Lquds Men Jn, Yucu Hou, b Weze Wu, *, Shuhng Ren nd Shdong Tn, L Xo, nd Zhgng Le Stte Key Lbortory of Checl Resource Engneerng, Beng Unversty of Checl Technology,
More informationDecomposition of Boolean Function Sets for Boolean Neural Networks
Decomposton of Boolen Functon Sets for Boolen Neurl Netorks Romn Kohut,, Bernd Stenbch Freberg Unverst of Mnng nd Technolog Insttute of Computer Scence Freberg (Schs), Germn Outlne Introducton Boolen Neuron
More informationStudy and modeling on saponification dynamics of the mixture of insect wax and oil-tea camellia seed oil
Avlble onlne www.jocpr.com Journl of Chemcl nd Phrmceutcl Reserch, 04, 6(4):568-574 Reserch Artcle ISSN : 0975-7384 CODEN(USA) : JCPRC5 Study nd modelng on sponfcton dynmcs of the mxture of nsect wx nd
More informationReproducing Kernel Hilbert Space for. Penalized Regression Multi-Predictors: Case in Longitudinal Data
Interntonl Journl of Mthemtcl Anlyss Vol. 8, 04, no. 40, 95-96 HIKARI Ltd, www.m-hkr.com http://dx.do.org/0.988/jm.04.47 Reproducng Kernel Hlbert Spce for Penlzed Regresson Mult-Predctors: Cse n Longudnl
More informationResearch on the Fuzzy Control for Vehicle Semi-active Suspension. Xiaoming Hu 1, a, Wanli Li 1,b
Advanced Materals Research Onlne: 0-0- ISSN: -9, Vol., pp -9 do:0.0/www.scentfc.net/amr.. 0 Trans Tech Publcatons, Swterland Research on the Fuy Control for Vehcle Sem-actve Suspenson Xaomng Hu, a, Wanl
More informationMEASUREMENT OF MOMENT OF INERTIA
1. measurement MESUREMENT OF MOMENT OF INERTI The am of ths measurement s to determne the moment of nerta of the rotor of an electrc motor. 1. General relatons Rotatng moton and moment of nerta Let us
More informationTHE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions
THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George
More informationEffects of polarization on the reflected wave
Lecture Notes. L Ros PPLIED OPTICS Effects of polrzton on the reflected wve Ref: The Feynmn Lectures on Physcs, Vol-I, Secton 33-6 Plne of ncdence Z Plne of nterfce Fg. 1 Y Y r 1 Glss r 1 Glss Fg. Reflecton
More informationragsdale (zdr82) HW6 ditmire (58335) 1 the direction of the current in the figure. Using the lower circuit in the figure, we get
rgsdle (zdr8) HW6 dtmre (58335) Ths prnt-out should hve 5 questons Multple-choce questons my contnue on the next column or pge fnd ll choces efore nswerng 00 (prt of ) 00 ponts The currents re flowng n
More informationA Family of Multivariate Abel Series Distributions. of Order k
Appled Mthemtcl Scences, Vol. 2, 2008, no. 45, 2239-2246 A Fmly of Multvrte Abel Seres Dstrbutons of Order k Rupk Gupt & Kshore K. Ds 2 Fculty of Scence & Technology, The Icf Unversty, Agrtl, Trpur, Ind
More informationDESIGN OF MULTILOOP CONTROLLER FOR THREE TANK PROCESS USING CDM TECHNIQUES
DESIGN OF MULTILOOP CONTROLLER FOR THREE TANK PROCESS USING CDM TECHNIQUES N. Kngsb 1 nd N. Jy 2 1,2 Deprtment of Instrumentton Engneerng,Annml Unversty, Annmlngr, 608002, Ind ABSTRACT In ths study the
More informationSlovak University of Technology in Bratislava Institute of Information Engineering, Automation, and Mathematics PROCEEDINGS
lovk Unverst of Technolog n Brtslv Insttute of Inmton Engneerng, Automton, nd themtcs PROCEEDING of the 8 th Interntonl Conference on Process Control Hotel Ttrs, Ttrnská Lomnc, lovk, June 4 7, IBN 978-8-7-357-9
More informationDennis Bricker, 2001 Dept of Industrial Engineering The University of Iowa. MDP: Taxi page 1
Denns Brcker, 2001 Dept of Industrl Engneerng The Unversty of Iow MDP: Tx pge 1 A tx serves three djcent towns: A, B, nd C. Ech tme the tx dschrges pssenger, the drver must choose from three possble ctons:
More informationTwo Activation Function Wavelet Network for the Identification of Functions with High Nonlinearity
Interntonl Journl of Engneerng & Computer Scence IJECS-IJENS Vol:1 No:04 81 Two Actvton Functon Wvelet Network for the Identfcton of Functons wth Hgh Nonlnerty Wsm Khld Abdulkder Abstrct-- The ntegrton
More informationDemand. Demand and Comparative Statics. Graphically. Marshallian Demand. ECON 370: Microeconomic Theory Summer 2004 Rice University Stanley Gilbert
Demnd Demnd nd Comrtve Sttcs ECON 370: Mcroeconomc Theory Summer 004 Rce Unversty Stnley Glbert Usng the tools we hve develoed u to ths ont, we cn now determne demnd for n ndvdul consumer We seek demnd
More informationLAPLACE TRANSFORM SOLUTION OF THE PROBLEM OF TIME-FRACTIONAL HEAT CONDUCTION IN A TWO-LAYERED SLAB
Journl of Appled Mthemtcs nd Computtonl Mechncs 5, 4(4), 5-3 www.mcm.pcz.pl p-issn 99-9965 DOI:.75/jmcm.5.4. e-issn 353-588 LAPLACE TRANSFORM SOLUTION OF THE PROBLEM OF TIME-FRACTIONAL HEAT CONDUCTION
More informationCHOVER-TYPE LAWS OF THE ITERATED LOGARITHM FOR WEIGHTED SUMS OF ρ -MIXING SEQUENCES
CHOVER-TYPE LAWS OF THE ITERATED LOGARITHM FOR WEIGHTED SUMS OF ρ -MIXING SEQUENCES GUANG-HUI CAI Receved 24 September 2004; Revsed 3 My 2005; Accepted 3 My 2005 To derve Bum-Ktz-type result, we estblsh
More informationAn Improved multiple fractal algorithm
Advanced Scence and Technology Letters Vol.31 (MulGraB 213), pp.184-188 http://dx.do.org/1.1427/astl.213.31.41 An Improved multple fractal algorthm Yun Ln, Xaochu Xu, Jnfeng Pang College of Informaton
More informationWeek 9 Chapter 10 Section 1-5
Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,
More informationStudy on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI
2017 2nd Internatonal Conference on Electrcal and Electroncs: echnques and Applcatons (EEA 2017) ISBN: 978-1-60595-416-5 Study on Actve Mcro-vbraton Isolaton System wth Lnear Motor Actuator Gong-yu PAN,
More informationPhysics 207: Lecture 20. Today s Agenda Homework for Monday
Physcs 207: Lecture 20 Today s Agenda Homework for Monday Recap: Systems of Partcles Center of mass Velocty and acceleraton of the center of mass Dynamcs of the center of mass Lnear Momentum Example problems
More informationParameter Estimation for Dynamic System using Unscented Kalman filter
Parameter Estmaton for Dynamc System usng Unscented Kalman flter Jhoon Seung 1,a, Amr Atya F. 2,b, Alexander G.Parlos 3,c, and Klto Chong 1,4,d* 1 Dvson of Electroncs Engneerng, Chonbuk Natonal Unversty,
More informationSpring 2002 Lecture #13
44-50 Sprng 00 ecture # Dr. Jaehoon Yu. Rotatonal Energy. Computaton of oments of nerta. Parallel-as Theorem 4. Torque & Angular Acceleraton 5. Work, Power, & Energy of Rotatonal otons Remember the md-term
More informationNovember 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project
SE 8 Fnal Project Story Shear Frame u m Gven: u m L L m L L EI ω ω Solve for m Story Bendng Beam u u m L m L Gven: m L L EI ω ω Solve for m 3 3 Story Shear Frame u 3 m 3 Gven: L 3 m m L L L 3 EI ω ω ω
More informationName: SID: Discussion Session:
Nme: SID: Dscusson Sesson: hemcl Engneerng hermodynmcs -- Fll 008 uesdy, Octoer, 008 Merm I - 70 mnutes 00 onts otl losed Book nd Notes (5 ponts). onsder n del gs wth constnt het cpctes. Indcte whether
More informationModeling of Dynamic Systems
Modelng of Dynamc Systems Ref: Control System Engneerng Norman Nse : Chapters & 3 Chapter objectves : Revew the Laplace transform Learn how to fnd a mathematcal model, called a transfer functon Learn how
More informationA concept for Vertical Folding Convertible Design
The th IFToMM Worl Congress, Tpe, Twn, Octoer -0, 0 DOI Numer: 0./IFToMM.TH.WC.PS.00 A concept or Vertcl Folng Convertle Desgn Hsng-Hu Hung Chng-Yu Weng Sheng-u T Deprtment o Vehcle Engneerng, Ntonl Pngtung
More informationSequences of Intuitionistic Fuzzy Soft G-Modules
Interntonl Mthemtcl Forum, Vol 13, 2018, no 12, 537-546 HIKARI Ltd, wwwm-hkrcom https://doorg/1012988/mf201881058 Sequences of Intutonstc Fuzzy Soft G-Modules Velyev Kemle nd Huseynov Afq Bku Stte Unversty,
More informationAccurate Instantaneous Frequency Estimation with Iterated Hilbert Transform and Its Application
Proceedngs of the 7th WSEAS Interntonl Conference on SIGAL PROCESSIG, ROBOTICS nd AUTOMATIO (ISPRA '8) Unversty of Cmbrdge, UK, Februry -, 8 Accurte Instntneous Frequency Estmton wth Iterted Hlbert Trnsform
More informationPosition and Speed Control. Industrial Electrical Engineering and Automation Lund University, Sweden
Poton nd Speed Control Lund Unverty, Seden Generc Structure R poer Reference Sh tte Voltge Current Control ytem M Speed Poton Ccde Control * θ Poton * Speed * control control - - he ytem contn to ntegrton.
More informationIntroduction to Numerical Integration Part II
Introducton to umercl Integrton Prt II CS 75/Mth 75 Brn T. Smth, UM, CS Dept. Sprng, 998 4/9/998 qud_ Intro to Gussn Qudrture s eore, the generl tretment chnges the ntegrton prolem to ndng the ntegrl w
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationGravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11)
Gravtatonal Acceleraton: A case of constant acceleraton (approx. hr.) (6/7/11) Introducton The gravtatonal force s one of the fundamental forces of nature. Under the nfluence of ths force all objects havng
More informationStudy on Non-Linear Dynamic Characteristic of Vehicle. Suspension Rubber Component
Study on Non-Lnear Dynamc Characterstc of Vehcle Suspenson Rubber Component Zhan Wenzhang Ln Y Sh GuobaoJln Unversty of TechnologyChangchun, Chna Wang Lgong (MDI, Chna [Abstract] The dynamc characterstc
More informationChapter 11 Angular Momentum
Chapter 11 Angular Momentum Analyss Model: Nonsolated System (Angular Momentum) Angular Momentum of a Rotatng Rgd Object Analyss Model: Isolated System (Angular Momentum) Angular Momentum of a Partcle
More informationThree hours UNIVERSITY OF MANCHESTER SCHOOL OF COMPUTER SCIENCE. Date: Wednesday 16th January 2013 Time: 09:45-12:45
Three hours UNIVERSITY OF MANHESTER SHOOL OF OMPUTER SIENE Prllel Progrms nd ther Performnce Dte: Wednesdy 16th Jnury 2013 Tme: 09:45-12:45 Plese nswer ny TWO Questons from the FOUR Questons provded. Ths
More informationOnline Appendix to. Mandating Behavioral Conformity in Social Groups with Conformist Members
Onlne Appendx to Mndtng Behvorl Conformty n Socl Groups wth Conformst Members Peter Grzl Andrze Bnk (Correspondng uthor) Deprtment of Economcs, The Wllms School, Wshngton nd Lee Unversty, Lexngton, 4450
More informationKatholieke Universiteit Leuven Department of Computer Science
Updte Rules for Weghted Non-negtve FH*G Fctorzton Peter Peers Phlp Dutré Report CW 440, Aprl 006 Ktholeke Unverstet Leuven Deprtment of Computer Scence Celestjnenln 00A B-3001 Heverlee (Belgum) Updte Rules
More informationOperating conditions of a mine fan under conditions of variable resistance
Paper No. 11 ISMS 216 Operatng condtons of a mne fan under condtons of varable resstance Zhang Ynghua a, Chen L a, b, Huang Zhan a, *, Gao Yukun a a State Key Laboratory of Hgh-Effcent Mnng and Safety
More informationElectrical double layer: revisit based on boundary conditions
Electrcal double layer: revst based on boundary condtons Jong U. Km Department of Electrcal and Computer Engneerng, Texas A&M Unversty College Staton, TX 77843-318, USA Abstract The electrcal double layer
More information13 Design of Revetments, Seawalls and Bulkheads Forces & Earth Pressures
13 Desgn of Revetments, Sewlls nd Bulkheds Forces & Erth ressures Ref: Shore rotecton Mnul, USACE, 1984 EM 1110--1614, Desgn of Revetments, Sewlls nd Bulkheds, USACE, 1995 Brekwters, Jettes, Bulkheds nd
More informationThe Study of Teaching-learning-based Optimization Algorithm
Advanced Scence and Technology Letters Vol. (AST 06), pp.05- http://dx.do.org/0.57/astl.06. The Study of Teachng-learnng-based Optmzaton Algorthm u Sun, Yan fu, Lele Kong, Haolang Q,, Helongang Insttute
More information