A Novel Geometric Error Compensation Method for Gantry-Moving CNC Machine Regarding Dominant Errors

Size: px
Start display at page:

Download "A Novel Geometric Error Compensation Method for Gantry-Moving CNC Machine Regarding Dominant Errors"

Transcription

1 processes Article A Novel Geometric Error Compensation Method for Gantry-Moving CNC Machine Regarding Dominant Errors Hong Lu 1, Qian Cheng 1, Xinbao Zhang 2, *, Qi Liu 1, Yu Qiao 1 Yongquan Zhang 1 1 School of Mechanical Electronic Engineering, Wuhan University of Technology, Wuhan , China; lzh@whut.edu.cn (H.L.); nhqian@whut.edu.cn (Q.C.); liunuli@whut.edu.cn (Q.L.); xiaoqiao456@whut.edu.cn (Y.Q.); zhangyongquan@whut.edu.cn (Y.Z.) 2 School of Mechanical Science Engineering, Huazhong University of Science Technology, Wuhan , China * Correspondence: zhangxinbao1@hust.edu.cn Received: 7 July 2020; Accepted: 27 July 2020; Published: 31 July 2020 Abstract: Gantry-type computer numerical control (CNC) machines are widely used in manufacturing industry. A novel structure with moveable gantry is proposed to improve traditional gantry-type machine structure s disadvantage of taking up too much space. Geometric s have direct impacts on actual position of tool, which significantly reduces accuracy of machines. Errors of different components are always coupled have uncertain effects on total geometric. Thus, it is essential to find an effective way to identify dominant s do targeted. First, a novel identification method using value leaded global sensitivity analysis (VLGSA) is proposed to find dominant s. In VLGSA, weighting factors which show influence of range are used to modify multi-body system (MBS) model. Results show that dominant s in three directions respectively contribute 80%, 86% 85% of total in ir directions. n, s identified by VLGSA are modeled by least-square linear fitting B-spline interpolation methods, respectively, according to feature of data. Finally, models are applied in a new real-time system developed on Beckhoff TwinCAT servo system. Experimental results from gantry-type CNC engraving milling machine show proposed method can help figure out most dominant s reduce around 90% of total. Keywords: gantry-type CNC machine; geometric ; multi-body system; VLGSA; method 1. Introduction Gantry-type manufacturing equipment has already played an indispensable role in modern manufacturing industry. In traditional gantry-type machines, X-direction feed movement is realized by motion of workbench, while gantry is fixed on machine bed. A vast room of more than twice feed stroke should be reserved to ensure processing capability, which causes waste of space. Thus, a novel structure with a moveable gantry is proposed. Different from traditional structure, workbench is fixed two ball-screw feed systems are symmetrically arranged to drive gantry. Only half space is needed comparing with traditional one. However, machines with moveable gantries require higher accuracy. Many factors are affecting accuracy of gantry-type machine tools, such as geometric, rmal, servo system so forth. rmal s widely appear in engineering machinery devices many researchers have done great jobs on optimizing heat transformation [1 4]. Geometric is anor common component usually Processes 2020, 8, 906; doi: /pr

2 Processes 2020, 8, of 21 donates most in final inaccuracy [5]. geometric s directly influence final position of tool nose s are generally related to tiny defect of machine itself, like wrong structure of feed system. re are two ways to resolve problem. One is adjusting physical structure, which is complicated expensive may be limited by skill of technical staff. or one is compensating for existing s by numerical control (NC) code or software, which is much more flexible can achieve high quality [6]. However, needs a clear understing of existing s. Error modeling is firstly required to reach goal. technique of geometric modeling has been persistently developed in recent decades. In application of robot industry, Paul, R.P. [7] Veitschegger, W.K. [8] gave methods using homogeneous transformation matrix differential transformation. precise location of robot s terminal point with existence of geometric size s kinematic s was found in ir research. se works established basis of analysis. Chen, J.X. et al. [5] took furr research based on fundamental principle of Paul Veitschegger, by analyzing differential movement differential transformation. y proposed a comprehensive method of both discovering how geometric propagates through every single axis. Zhao, Q.Q. et al. [9] established topology structure of machine tool using multi-body system ory proposed two conceptions, which were relative motion matrices of axis connectors connection matrix of machine tool connectors. Based on homogeneous coordinate transformation ory, Zhou, X.T. [10] built model of classic three-axis machine tool, decomposed final into s existing in components in whole motion chain. se works showed development application of different modeling methods, gave inspirations to furr research. Still, many of m were limited in modeling. Models based on se methods are usually complex mamatical expressions of a vast amount of items, making it a sophisticated work to compensate for all items. However, towards several specific items, which can be defined as dominant s, may provide enough accuracy improvement. Thus, more research could be done to simplify. Sensitivity analysis (SA) is a kind of method to quantify uncertainties of an unknown system using great data, which may help in finding most valuable items to compensate. For a specific model, a vast amount of different input output data can be obtained through all kinds of scientific sampling methods powerful calculation tools. Thus sensitivity analysis can help to identify weight each input has in affecting system outputs. SA consists of two main branches, one is local sensitivity analysis (LSA) or is called global sensitivity analysis (GSA). [11,12] GSA is more commonly used in manufacturing industry with its ability to find out result differences under changes of parameters. It can emphasize interaction power of parameters by using additional methods. Usually, SA is used in economic analysis. However, it is also found suitable in engineering field recently. Fan, J.W. et al. [13] established a prediction model of component tolerance proposed a sensitivity analysis based on component motion. proposed sensitivity analysis was experimentally verified several main parameters that donate to s most were found out. Yan, Y. et al. [14] established model of multi-process assembling built sensitivity analysis model of parameters towards assembling using method of matrix differentiation. Li, J. Xue, F.G. [15] proposed general local sensitivity indices (GLSIs), general global sensitivity indices (GGSIs) general global sensitivity fluctuation indices (GGSFIs) based on traditional sensitivity analysis. A five-axis machine was analyzed using above indices. One common limit of above works is that sensitivity indices calculation is carried out on original model without considering impact of value range towards sensitivity analysis. Besides, sensitivity analysis of machine tools was only done for performance assessment, targeted was missing validity was not verified. Sensitivity analysis tells what to concentrate on; measurement fitting tell how to describe deal with m. With all kinds of advanced measurement equipment, precision data can be obtained. According to data, model can be established for. However, appropriate fitting methods are needed to ensure accuracy. To improve

3 Processes 2020, 8, of 21 manufacturing performance, many researchers have done great works. Fang, K.G. Yang, J.G. [16] considered rmal factors, used orthogonal polynomial method to build total model, compensated 90% of spindle geometric. Wang, W. et al. [17] used Newton interpolation method to construct rmal geometric model of computer numerical control (CNC) milling machine, established a real-time system. experiment based on Fanuc CNC machine was carried out to show system works well. Tuo, Z.Y. et al. [18] using cubic spline interpolation established rmal geometric model improved position accuracy of CNC machine. B-spline curve is commonly used in computer aided design (CAD) areas for its ability to describe continual material. It is also used in trajectory planning for its smooth character. Li, Y.H. et al. [19] used quantic B-spline curve to optimize motion path of pick--place robots. His work showed that B-spline curve is appropriate for smooth machine motion, which is also suitable for relative movement of machine tools. In this paper, multi-body system ory was used to establish geometric model of gantry-moving engraving milling machine. n translation s in all three directions are extracted from matrix to analyze significant items which impact precision most. Value leaded global sensitivity analysis (VLGSA) is proposed used in above analysis process, in which measurement value of each item is considered to modify total geometric model. B-spline curve was developed into B-spline interpolation method used on result of machine measurement. A real-time system based on Beckhoff servo system, TwinCAT control software C# language was built to eliminate s. experiment verified validity of above crucial finding method. 2. Geometric Modeling Based on Multibody System multi-body system (MBS) illustrates a kind of mechanical system in which a series of rigid or flexible bodies interconnect [20]; every single body may undergo translation rotation in different scales. It is a general systematic method for geometric modeling is widely used in analysis of robots, machine tools coordinate measurement. core of MBS is using low-order-body arrays to describe topology structure of system [20], using transformation matrices to calculate motion relationship. Usually, inertial coordinate system is chosen as body B 0, basic component is selected as body B 1, or elements are divided into several branches numbered along direction away from B 1, branch by branch Coordinate System Definition Topology Structure Description of Gantry-Moving CNC Machine In this paper, a gantry-moving engraving milling machine with three-dimension movements is taken to describe ideal of presented geometric modeling method. structure of chosen machine Processes 2020, is shown 8, x FOR PEER in Figure REVIEW 1. 4 of 24 Figure 1. Figure Structure 1. Structure of of gantry-moving engraving milling milling machine. machine. B0-bed, B1-moveable B gantry, 0 -bed, B 1 -moveable B2-Y B3-Z direction slide, B4-spindle with tool, B5-workbench, B6-workpiece. T01 gantry, B 2 -Y direction slide, B 3 -Z direction slide, B 4 -spindle with tool, B 5 -workbench, B 6 -workpiece. T56 represent coordinate transformation between components. T 01 T 56 represent coordinate transformation between components. gantry-moving CNC machine has two servo motors in X direction. Two motors are symmetrically arranged to drive gantry through ball-screw feed systems. X1 X2 represent movement of two feed systems. synchronization between two movements is an important cause of X-direction yaw. Since machine locates on a stable ground coordinate system, machine bed itself can be directly chosen as body B0. As orange lines show in Figure 1, engraving mill machine can be seen as a two-branch topology structure [21]. first branch (tool branch) consists of

4 Processes 2020, 8, of 21 gantry-moving CNC machine has two servo motors in X direction. Two motors are symmetrically arranged to drive gantry through ball-screw feed systems. X1 X2 represent movement of two feed systems. synchronization between two movements is an important cause of X-direction yaw. Since machine locates on a stable ground coordinate system, machine bed itself can be directly chosen as body B 0. As orange lines show in Figure 1, engraving mill machine can be seen as a two-branch topology structure [21]. first branch (tool branch) consists of moveable gantry (B 1 ), Y direction slide (B 2 ), Z direction slide (B 3 ) spindle with tool on it (B 4 ). second branch (workpiece branch) consists of workbench (B 5 ) workpiece (B 6 ). Table 1 shows main sizes of se components. Table 1. Size of main components. Main Parameters of Components Value (mm) Size of workbench (width length) Stroke of gantry (X-direction stroke) 2400 Stroke of Y slide (Y-direction stroke) 1200 Stroke of Z slide (Z-direction stroke) 240 final topology structure coordinate definition with position motion relationships are shown in Figure 2, where vectors between bodies show ir relations. Take B 0 B 1 as an example, T 01p means position transformation between B 0 B 1 ; T 01s means motion transformation between B 0 B 1 ; T 01pe T 01se represent of position motion transformations respectively. T 01p, T 01pe, T 01s T 01se constitute total coordinate transformation betweenprocesses B , B8, x 1 (TFOR 01 ). PEER REVIEW 5 of 24 Figure 2. Figure 2. topology topology structure coordinate definition of of machine. Orange Orange lines show lines show coordinate coordinate transformation transformation between between different bodies, which are are marked with with T01p, T01p, T01pe, T01pe, etc. etc. Table 2 shows low-order-body array of machine, in which L is low-order-body Table 2 shows low-order-body array of machine, in which L is low-order-body operator. operator. For body Bj, its n-th order low-order-body number can be defined as Equation (1): For body B j, its n-th order low-order-body number can be defined as Equation (1): n L(j) = i, (1) where n means that Bi is n-th order low-order-body L n (j) = i of body Bj (1) Similarly, or equation can be inferred as Equations (2) (3): where n means that B i is n-th order low-order-body n n 1 of body B L(j) = L(L (j)) j (2) 0 L(j) = j. (3) Table 2. Low-order-body Array of Chosen Machine. Body L1(j)

5 Processes 2020, 8, of 21 Similarly, or equation can be inferred as Equations (2) (3): L n (j) = L(L n 1 (j)) (2) L 0 (j) = j. (3) Table 2. Low-order-body Array of Chosen Machine. Body L1(j) L2(j) L3(j) 0 1 L4(j) Coordinate Transformation Model Establishment actual position motion between bodies can be expressed by matrix, more specifically by position matrix, motion matrix, position matrix motion matrix. After binging every single body to ir fixed coordinate, position motion transformation can be demonstrated by multiplication of a series of 4 4 homogeneous eigenmatrix [22]. For a three-dimensional Cartesian coordinate system, translation can be expressed as Equation (4), in which x, y z are translation along three axes, respectively. Rotation around three axes can be expressed as Equations (5) (7) x y T = (4) z Rot(x, θ x ) = Rot(y, θ y ) = Rot(z, θ z ) = cosθ x sin θ x 0 0 sin θ x cos θ x 0 cos θ y 0 sin θ y sin θ y 0 cos θ y 0 cos θ z sin θ z 0 0 sin θ z cos θ z However, both structure motion exist in manufacturing process. Take X-direction movement as an example (shown in Figure 3). Errors between two coordinates can be divided into six types, including three linear three angle. Linear s are position straightness along with two directions. Angle s are roll, pitch yaw. Similarly, geometric s in or directions can be defined as Table 3. (5) (6) (7)

6 sinθ cosθ 0 0 z z Rot(z, θ ) = z However, both structure motion exist in manufacturing process. Take X- direction movement as an example (shown in Figure 3). Errors between two coordinates can be divided into six types, including three linear three angle. Linear s are position Processes 2020, 8, of 21 straightness along with two directions. Angle s are roll, pitch yaw. Similarly, geometric s in or directions can be defined as Table 3. (7) (a) (b) (c) (d) Figure 3. Error measurement. (a) Principle of position measurement; (b) Principle of straightness measurement; (c) Position measurement experiment; (d) Straightness measurement experiment. Figure 3. Error measurement. (a) Principle of position measurement; (b) Principle of straightness measurement; (c) Position measurement experiment; (d) Straightness measurement experiment. Table 3. Geometric s in chosen caving milling machine. Error Type Error Item Error Symbol Position x x Y-direction straightness y x X-direction feed s Z-direction straightness z x Roll α x Pitch β x Yaw γ x Position y y X-direction straightness x y Y-direction feed s Z-direction straightness z y Roll β y Pitch α y Yaw γ y Position z z X-direction straightness x z Z-direction feed s Y-direction straightness y z Roll γ z Pitch β z Yaw α z X-Y perpendicularity γ xy Structure s X-Z perpendicularity β xz Y-Z perpendicularity α yz According to principle of homogeneous matrix transformation, general matrix between B 0 B 1 which containing all six s can be calculated by multiplying six homogeneous eigenmatrices, as Equation (8) shows (where c for cos s for sin ). T 01e = T x x T y x T z x T θx x T θy x T θz x c θy x c θz x c θy x s θz x s θy x x x c θx x s θz x + s θx x s θy = x c θz x c θx x c θz x s θx x s θy x s θz x c θy x s θx x y x s θx x s θz x c θx x c θz x s θy x c θz x s θx x + c θx x s θy x s θz x c θx x c θy x z x (8)

7 Processes 2020, 8, of 21 In translation movement, three angle s are usually tiny, so approximation that sin θ θ cos θ 1 can be introduced. Putting approximation into Equation (8) ignoring second-order high-order elements, equation can be simplified as Equation (9). T 01e = 1 θz x θy x x θz x 1 θx x y θy x θx x 1 z According to method above, transformation matrices of chosen engraving mill machine are listed in Table 4. Table 4. transformation matrix of chosen caring mill machine. Neighbor Bodies Ideal Position Motion Matrix Position Motion Error Matrix 1 γ xy β xz 0 γ T 01p = I 4 4 T pe = xy 1 α yz 0 β xz α yz x 1 γ x β x x x γ T 01s = T se = x 1 α x y x β x α x 1 z x (9) T 12s = T 23s = T 12p = I 4 4 T 12pe = I y T se = 1 γ y β y x y γ y 1 α y y y β y α y 1 z y T 23p = I 4 4 T 23p = I y T se = 1 γ z β z x z γ z 1 α z y z β z α z 1 z z T 34p = I 4 4 T 34pe = I 4 4 T 34s = I 4 4 T 34se = I 4 4 T 05p = I 4 4 T 05pe = I 4 4 T 05s = I 4 4 T 05se = I 4 4 T 56p = I 4 4 T 56pe = I 4 4 T 56s = I 4 4 T 56se = I 4 4 coordinate system eigenmatrix of two neighboring bodies can be written as Equation (10). Eigenmatrix with s can be written as Equation (11). T 01 = T 01p T 01pe T 01s T 01se (10) T 01e = T 01 T 01i. (11) Coordinate system eigenmatrix of any two bodies can be written as Equation (12). T SG = K=1 K=m,L m (S)=Q T K L (S)LK 1 (S) 1 J=1 J=n,L n (G)=Q T J L (G)LJ 1 (G) (12) As for topology structure of chosen engraving milling machine mentioned above, coordinate transformation eigenmatrix in both ideal situation actual situation (with s) should be calculated as Equation (13) T 64 = ( k=1 1 j=1 T L k (6)L (6)) k 1 T L j (4)L j 1 (4) k=2 j=4 (13) = (T 05 T 56 ) 1 T 01 T 12 T 23 T 34

8 Processes 2020, 8, of 21 geometric of machine can n be calculated by subtracting ideal transformation eigenmatrix actual transformation eigenmatrix between tool workpiece, like Equation (14). result T 64e is final matrix containing angle s linear s between tool coordinate system workpiece coordinate system. In planar milling, angle s represent small posture fault linear s along three axes have most significant impact on misalignment between tool nose target point on workpiece. 3. Dominant Errors Identification with VLGSA T 64e = T 64 T 64i (14) From perspective of black box, system models can always be demonstrated as a function Y = f(x), in which x is vector of several uncertain inputs of system y is one of system outputs which is chosen to be analyzed. In this paper, model targets total geometric of engraving milling machine. Y represents final X equals vector of each item like X-direction position. By doing VLGSA to model function, effects of inputs on a particular output can be quantified into indices. In or words, impact of each item on total geometric can be described with numbers dominant s, which cause most of total, can be found Decomposition of Variance Indices Calculating In system, model function Y = f(x) is under assumption that inputs distribute independently within a unit uniform space. refore, actual input scales should be transformed into unit space [23]. function can be decomposed as Equation (15). Y = f 0 + d f i (X i ) + i=1 1 i<j d f ij (X i, X j ) f 12...d (X 1, X 2,..., X d ), (15) where f 0 is a component of constant; f i is function of X i, which represents effect X i donates alone ( geometric caused by one specific component); f ij is function with two inputs X i X j, representing effect to output given by X i X j toger. Similarly, or higher parts can be defined as above. This composition must be in condition that sub-components are independent of each or. In or words, all components in equation are orthogonal. It can be shown in following Equation (16). 1 f i1 i 2...i s (X i1, X i2,..., X is )dx ik = 0, (k = 1, 2,..., s). (16) 0 It makes result that components of decomposition function can be defined in terms of expected values, as Equations (17) (19) show. f 0 = E(Y) (17) f i (X i ) = E(Y Xi ) f 0 (18) f ij (X i, X j ) = E(Y Xi, X j ) f 0 f i f j (19) Or components containing more inputs can be defined as equations above. f i which shows effect of change input X i can be known as primary interaction or first-order interaction. f ij shows effect of changing both X i X j at same time. It can be known as second-order interaction. higher-order interactions can be defined similarly.

9 Processes 2020, 8, of 21 Square integrate decomposition Equation (15), new equation is got can be organized as Equations (20) (21). 1 0 f2 (X)dX = f d 1 0 i=1 1 0 f i 2 dx i + 1 i<j d f 2 (X)dX f 0 2 = f ij 2 dx i dx j f 12...d 2 dx 1 dx 2... dx d (20) d j=1 1 i<j d f 2 i...j dx i... dx j (21) 0 It is found that left part of equation is variance of Y(X) right part is variance terms. This leads to decomposition with variance expression, shown as Equations (22) (24). Var(Y) = d V i + i=1 1 i<j d V ij V 12...d (22) V i = Var Xi (E X i (Y Xi )) (23) V ij = Var Xij (E X ij (Y X i, X j )) V i V j, (24) where X i represents whole series of variables except X i. n first-order sensitivity index, which is direct sensitivity factor based on variance, can be defined as Equation (25). S i = V i Var(Y). (25) After estimating, component Var Xi (E X i (Y Xi )) can be written as Equation (26). Var Xi (E X i (Y X i )) 1 N N f(b) j (f(ab i ) j f(a) j ) (26) j=1 second-order higher-order sensitivity index can be defined in same way, which represents effect of output variance causing by change of several inputs, simultaneously. sensitivity indices satisfy following principle. d S i + i=1 1 i<j d S ij S 12...d = 1 (27) For some functions which are simple in structure, first-order sensitivity indices can be obtained by calculating integrals of decomposition function. However, more common way for most cases is estimating, which provides an easy solution to all kinds of functions it is often finished by Monte Carlo method. Monte Carlo method needs to generate a set of rom points in unit hypercube space, n use points to calculate use statistic method to estimate character of model. In actual work, rom points are usually pseudorom. most used points sequence is Latin hypercube sequence Sobol sequence. Several steps should be taken as follow to calculate sensitivity index using estimation method. 1. Generate a matrix of N 2d, in which N is number of sample points d is dimension of model (number of inputs). Points in se matrices should be romly distributed. 2. Generate matrix A using first d columns of above matrix matrix B using rest d columns. se two matrices distribute first two samples, which have n points in d-dimension hypercube space.

10 Processes 2020, 8, of Generate d more matrices AB i by replacing i-th column of matrix A with same column in matrix B. 4. Use above d + 2 matrices as inputs of model get d + 2 result matrices. ses matrices contain f(a) j, f(b) j f(ab i ) j mentioned in Equation (26). 5. Use matrices above to calculate sensitivity index with Equation (27) Machine Geometric Error Measurement for Value Leaded Factors Construction sample points of global sensitivity analysis distribute in unit hypercube space. However, items are different in scale usual sample method may cause inaccuracy in estimating actual impact of each item. Thus, a value-leaded factor based on measurement is proposed to modify geometric model, which normalizes s effects to same level. Keysight 5530 dual-frequency laser interferometer measurement system is used in measurement, cooperating with several kinds of mirror groups, 21 geometric s of engraving milling machine can be detected. Measurements are done with stard of ISO230-1 processes are shown in Figure 3. This work includes many time-consuming steps. For each component, measurement is done in different feed speeds, with range from 500 mm/min to 6000 mm/min interval of 500 mm/min. measurement stroke is x = 2000 mm, y = 1000 mm, z = 200 mm, which covers main processing area. start points are set in x = 200 mm, y = 100 mm z = 20 mm in machine coordinate system to leave enough space for installing measurement instruments. All measurements are done for three times. scale of each ir ratios are obtained from measurements. y are used for normalization of sample of global sensitivity analysis. measured s are shown in Table 5. Table 5. Max value in whole stroke. No. Item Max Value No. Item Max Value 1 x x mm 12 γ y mm/1000 mm 2 y x mm 13 z z mm 3 z x mm 14 x z mm 4 α x mm/1000 mm 15 y z mm 5 β x mm/1000 mm 16 γ z mm/1000 mm 6 γ x mm/1000 mm 17 β z mm/1000 mm 7 y y mm 18 α z mm/1000 mm 8 x y mm 19 γ xy mm/1000 mm 9 z y mm 20 β xz mm/1000 mm 10 β y mm/1000 mm 21 α yz mm/1000 mm 11 α y mm/1000 mm mm/1000 mm means deviation caused by α x in stroke of 1000 mm is mm. Or angle s are described in same way Indices Calculation Global Sensitivity Analysis analysis program is made in Matlab according to steps of proposed VLGAS method. geometric model of engraving milling machine built with multi-body system ory is adjusted for normalization based on measurement. Putting model into analysis program setting analysis parameters, first-order sensitivity indices in three directions are calculated listed in Table 6. angle s are at a very low level, high-order elements linking to angle s can be eliminated from MBS model. Thus, some elements sensitivity indices are zero in a particular direction. For X-direction movement, as shown in Figure 4, with geometric s changing in a certain range, sensitivity indices of x x, x y x z are largest contribute 80% of total value toger. For Y direction, most dominant element is y y, which has highest sensitivity index value of 0.424, or two dominant elements are y x y z, se three elements contribute about

11 Processes 2020, 8, of 21 86% of total value. Similarly, three linear s z x, z y z z have largest impact on Z direction, with proportion of around 85%. Number Table 6. First-order sensitivity index of each item. Item First-Order Sensitivity X-Direction Y-Direction Z-Direction 1 x x y x z x α x β x γ x y y x y z y β y α y γ y z z x z y z γ z β z α z γ xy β xz α yz data is drawn shown in Figure 4 to demonstrate clearly. In conclusion, for chosen engraving milling machine, dominant s are mainly Processes 2020, 8, x FOR PEER REVIEW 14 of 24 linear items, specifically speaking position two straightness s from or two feed directions. feed directions. In orin words, or words, aiming at at se se three three s s can provide can provide remarkable remarkable improvement improvement to chosen to chosen machine machine with with minimum cost. (a) (b) (c) Figure 4. Figure Image4. Image of first-order of sensitivity index. index. (a) First-order (a) First-order sensitivity of sensitivity X-direction ; of (b) X-direction Firstorder sensitivity of Y-direction ; (c) First-order sensitivity of Z-direction. ; (b) First-order sensitivity of Y-direction ; (c) First-order sensitivity of Z-direction. 4. Error Compensation In this section, data in whole stroke are listed modeled based on ir character. Different models are applied to a novel real-time system compared. linear in X direction is chosen as an example. Errors in or directions can be dealt same way. As shown in last section, main components, which have largest impact on X-direction linear, are X-direction position, X-direction straightness along with Y- direction movement X-direction straightness along with Z-direction movement. y all have a close relation to feed movement, which means y can be compensated by

12 Processes 2020, 8, of Error Compensation In this section, data in whole stroke are listed modeled based on ir character. Different models are applied to a novel real-time system compared. linear in X direction is chosen as an example. Errors in or directions can be dealt same way. As shown in last section, main components, which have largest impact on X-direction linear, are X-direction position, X-direction straightness along with Y-direction movement X-direction straightness along with Z-direction movement. y all have a close relation to feed movement, which means y can be compensated by adjusting feed trajectory Position Error Model Based on Least-Square Method Position s are measured three times on each point average values are calculated to reduce measurement. values are listed in Tables 7 9. Position (mm) Table 7. Position in X direction. Error (µm) First Second Third Average Origin of position located on x = 200 mm in machine coordinate system. Position (mm) Table 8. Position in Y direction. Error (µm) First Second Third Average Origin of position located on y = 100 mm in machine coordinate system. According to tables, position s show an obverse trend of linear distribution in whole stroke. value accumulates through feed movement maximum value appears at end of stroke. maximum values in three directions are , , respectively. Thus, least square method is chosen to establish mamatical model of position s to describe linear trend.

13 Processes 2020, 8, of 21 Position (mm) Table 9. Position in Z direction. Error (µm) First Second Third Average Processes 2020, 8, x FOR PEER REVIEW 16 of 24 Origin of position located on z = 20 mm in machine coordinate system. end of stroke. maximum values in three directions are , , respectively. Thus, least square method is chosen to establish mamatical model of least-square method is a commonly used fitting method, in which fitting curve not position s to describe linear trend. necessarily cross least-square sample points method but is a makes commonly se used points fitting as method, closein aswhich possible fitting to curve curve. not For sample points necessarily (x i, y i ), cross a series sample of equations points but can makes bese built points as used close toas create possible to function curve. For of points, sample points ( x residual sum of squares n i, y i ), a series δ 2 of equations can be built used to create function of points, is chosen to decide best coefficients of unknown function. i residual sum of squares i=1 n 2 δ is chosen to decide best coefficients of unknown function. i i= 1 fitting functions are shown as Equations (28) (30). fitting curves of position s are fitting functions are shown as Equations (28) (30). fitting curves of position s are shown in Figure 5. shown in Figure 5. y = x (28) y = x (28) y = 0.071x (29) y = 0.071x (29) y = 0.108x y = (30) (30) (a) (b) (c) Figure 5. Image Figure of 5. fitting Image of curves. fitting curves. (a) Fitting (a) Fitting curve of X-direction position position ; ; (b) Fitting (b) curve Fitting of Y-curve of direction position ; (c) Fitting curve of Z-direction position. Y-direction position ; (c) Fitting curve of Z-direction position Straightness Error Model Based on B-Spline Interpolation. straightness s are also measured three times to calculate average value. curve of s is drawn in Figure 6 to show intuitively.

14 Processes 2020, 8, of Straightness Error Model Based on B-Spline Interpolation. straightness s are also measured three times to calculate average value. curve of Processes 2020, 8, x FOR PEER REVIEW 17 of 24 s is drawn in Figure 6 to show intuitively. (a) (b) Figure 6. Figure Image 6. of Image straightness of straightness. (a). Straightness (a) Straightness in X-direction in X-direction movement; movement; (b) Straightness (b) Straightness in Y-direction movement; (c) Straightness in Z-direction movement. in Y-direction movement; (c) Straightness in Z-direction movement. According to Figure 6, max value of straightness s is lower than position s. According to Figure 6, max value of straightness s is lower than position s. curves fluctuate with total trend of increasing first n decreasing, except Y-direction curves fluctuate along with with X-direction totalmovement. trend of increasing maximum first values respectively n decreasing, appear in except point Y-direction of along x1 = 1600 withmm, x2 X-direction = 1000 mm, movement. y1 = 400 mm, y2 = 500 maximum mm, z1 = values 100 mm respectively z2 = 100 mm appear for insix point of x 1 = straightness 1600 mm, xs. 2 = 1000 mm, distortion y 1 = 400 of mm, guides y 2 = during 500 mm, manufacturing z 1 = 100 mm is leading z 2 = 100 cause mmof for six straightness s. distortion of guides during manufacturing is leading cause of B-spline curve is a kind of unique curve, which is developed from Bezier curve is straightness s. constructed by linear combination of B-spline base curves. B-spline curves have many excellent properties, B-spline such curve as geometric is a kind invariance, of unique convex curve, hull, which convexity is developed retention, from variation reduction Bezier curve is constructed local support by linear so forth. combination B-spline surface of B-spline reconstruction base curves. is one of B-spline research curves hotspots have many critical excellent properties, technologies such asof geometric reverse engineering. invariance, refore, convexin hull, this paper, convexity B-spline retention, interpolation variationis reduction used to local support establish so mamatical forth. B-spline surface model of reconstruction straightness is, onewhich of can research describe hotspots continual critical technologies material of property reverse of engineering. guide. B-spline refore, base in function this paper, curve construction B-spline interpolation function are listed is used to in Appendix A. establish mamatical model of straightness, which can describe continual Exping measurement point number doing simulation with material interpolation property of result, guide. it is shown B-spline that values base function of residual curve keep construction at a low level function after are listed in Appendix interpolation A.. Interpolation results are shown in Figure 7. Exping measurement point number doing simulation with interpolation result, it is shown that values of residual keep at a low level after interpolation. Interpolation results are shown in Figure 7. (c)

15 Processes 2020, 8, of 21 Processes 2020, 8, x FOR PEER REVIEW 18 of 24 (a) (b) (c) (d) (e) (f) Figure Figure 7. Compensation 7. Compensation simulation simulation with with B-spline B-spline interpolation model. model. (a) (a) Y-direction straightness in X-direction in X-direction movement; movement; (b) Z-direction (b) Z-direction straightness straightness in X-direction in X-direction movement; movement; (c) X-direction (c) direction straightness in Y-direction movement; (d) Z-direction straightness in Y-direction straightness in Y-direction movement; (d) Z-direction straightness in Y-direction movement; movement; (e) X-direction straightness in Z-direction movement; (f) Y-direction straightness (e) X-direction straightness in Z-direction movement; (f) Y-direction straightness in Z-direction movement. in Z-direction movement Compensation 4.3. Compensation System System Experiment Experiment Aiming to, many researchers in previous works try to change G code Aiming to, many researchers in previous works try to change G code according to modeling result, which performs well in reducing value [24]. according However, to this kind of modeling method need result, pre-work which for performs each program well in cannot reducing be real-time. With value [24]. However, advantages this kind of Beckhoff s of method open need CNC pre-work system, for eachmodel program can help to cannot improve bemachine real-time. accuracy With advantages in a more of Beckhoff s effective open real-time CNCway. system, can help to improve machine accuracy in a more effective real-time way. strategy combines virtual axis electronic cam technology based on TwinCAT control system of Beckhoff. core of strategy is uploading model points into cam table to help change actual tool path offsetting trajectory deviation. Firstly,

16 Processes 2020, 8, x FOR PEER REVIEW 19 of 24 strategy combines virtual axis electronic cam technology based on Processes 2020, 8, of 21 TwinCAT control system of Beckhoff. core of strategy is uploading model points into cam table to help change actual tool path offsetting trajectory deviation. Firstly, a virtual axis is established to accept motion comm. virtual axis is a kind of electronic symbol in control software. It It does not connect to any physical component only has digital motion simulation. refore, it itrepresents represents ideal ideal motion motion state. state. Secondly, Secondly, set set virtual virtual axis axis as as master master create create anor anor axis axis symbol symbol connecting connecting to to physical physical axis axis motor motor as as slave. slave. Thus, Thus, slave slave axis axis can can move move with with master master axis, axis, which whichis is virtual virtualaxis axis has hasno physical output, according to to specific specific rules. rules. Thirdly, Thirdly, set set rule rule between between two axes two as axes electronic as electronic cam, n cam, n physical axis physical can do axis superimposed can do superimposed movements movements of both ideal of both movements ideal movements cam-curve cam-curve movements movements based on based on models. models. cam cam table table can can be created be created in in two two ways: ways: cam cam design design tool tool assembled in TwinCAT software software user user interface interface linking linking to softto programmable soft programmable logic controller logic controller (PLC) program. (PLC) Soft program. PLC method Soft PLC has method better real-time has better property real-time property table value can table bevalue changed can in be changed feed process. in feed process. system can be described system in can Figure be described 8. hardware in Figure of 8. hardware system of is shown in Figure system 9. is shown in Figure 9. Figure 8. structure of system. Figure 9. hardware of system. Suppose that x is ideal position when feed system moves along X direction xδx is is position generated in in feed feed stroke, stroke, n n cam cam table table value value in this in this position position should should be set be as set x. as Δx Give. Give motion motion comm comm to to virtual virtual axis, axis, electronic electronic cam cam offset offset will will compensate compensate in position x. With help of Beckhoff ErCAT communication system, massive data can be transmitted in real-time. refore, coordinate information in all directions can be read cam

17 Processes 2020, 8, x FOR PEER REVIEW 20 of 24 Processes 2020, 8, x FOR PEER REVIEW 20 of 24 Processes 2020, data of 21 in 8,position x. With help of Beckhoff ErCAT communication system, massive can transmitted refore, coordinate information in all directions be read be in positioninx.real-time. With help of Beckhoff ErCAT communication system,can massive data cam table can be modified with calculated value based on models, can can be transmitted real-time. refore, coordinate information in all can be read table be modifiedinwith calculated value based on directions models, without delay. without process expressed in Figure 10.value based on models, camdelay. table can be modified with is calculated process is expressed in Figure 10. without delay. process is expressed in Figure 10. Figure 10. process of proposed method. Figure processof of proposed proposed method. Figure process method. system is composed of both hardware software. hardware platform composedof ofboth bothhardware hardware software. software. hardware platform system system is is composed hardware platform of system is provided by Beckhoff industrial personal computer (PC) ofof by Beckhoff Beckhoffindustrial industrialpersonal personal computer (PC) system system is provided provided by computer (PC) software is written in C# language. With Beckhoff ADS communication protocol, C# language. language.with With Beckhoff Beckhoff ADS communication protocol, software software is is written written in C# ADS communication protocol, data in PLC Servo controller areare easy toto bebe dealt with algorithm inin software data PLC Servo controller easy dealt with algorithm software data inin PLC Servo controller are easy to be dealt with algorithm in software n transmitted back. n ntransmitted transmittedback. back. Straightness is measured verify effect. Straightness after after is effect. Straightness after ismeasured measuredtoto toverify verify effect. models built by three or usual fitting methods (Polynomial Fitting, Newton Interpolation Cubic models built by three or usual fitting methods (Polynomial Fitting, Newton Interpolation models built by three or usual fitting methods (Polynomial Fitting, Newton Interpolation Spline are usedare to compare with with chosen method ofmethod B-spline interpolation. As for CubicInterpolation) SplineInterpolation) Interpolation) are used to chosen of of B-spline interpolation. Cubic Spline used to compare compare with chosenmethod B-spline interpolation. total effect, a space movement of x = 2000 mm, y = 1000 mm z = 200 mm carried for total total effect, effect, aa space mm, y =y 1000 mm z = z200 mmmm AsAsfor spacemovement movementofofx x= = mm, = 1000 mm =is200 out X-direction total is measured. values are measured by absolute grating is carried out X-direction total is measured. values are measured by absolute is carried out X-direction total is measured. values are measured by absolute grating rulers. total isatcurbed at a low level, both modeling method rulers. Ifrulers. total is curbed a low level, both method sensitivity grating IfIf total is curbed at a low level, bothmodeling modeling method sensitivity analysis method can be verified to be effective, for only dominant components analysis method can be verified to be effective, for only dominant components analyzed by sensitivity analysis method can be verified to be effective, for only dominant components analyzed by method are compensated. Compensation resultsin are shown11. in Figure 11. method are compensated. Compensation results are shown Figure analyzed by method are compensated. Compensation results are shown in Figure 11. (a) (a) (b) Figure 11. Cont. (b)

18 Processes Processes 2020, 2020, 8, 8, 906 x FOR PEER REVIEW 21 of 1824 of 21 (c) Figure Compensation experiment result. (a) Residual Y-direction straightness in in X-direction movement before after ; (b) Residual Z-direction straightness in in X-direction movement before after ; (c) (c) Residual Residual X-direction X-direction total total before before after after. In In Figure 11a,b, straightness s along X-direction movement are are shown. five five data data curves curves are are respectively residual residual of four of chosen four chosen models, models, including including proposed B-spline proposed model B-spline three model commonly three used commonly fitting used models fitting models without anywithout. any Three. comparisonthree models comparison are Newton models are interpolation Newton model, interpolation cubicmodel, spline interpolation cubic spline interpolation polynomial model. polynomial From image, model. From it is obvious image, that it is obvious that ability of ability B-spline interpolation of B-spline model interpolation is most significant. model is Newton most significant. interpolation Newton model interpolation has good performance model has good in performance in middle part of stroke but it also shows Runge s phenomenon, which leads middle part of stroke but it also shows Runge s phenomenon, which leads to huge residual to huge residual in beginning end of movement. polynomial model has a in beginning end of movement. polynomial model has a good ability to good ability to follow trend of curve had has a smooth residual curve. With follow trend of curve had has a smooth residual curve. With polynomial model, polynomial model, process avoids steep fluctuation of velocity acceleration. process avoids steep fluctuation of velocity acceleration. cubic spline cubic spline interpolation model has same capability of eliminating with B-spline interpolation model has same capability of eliminating with B-spline interpolation model interpolation model in most positions but still has defects in beginning end of in prediction. most positions B-spline but still interpolation has defects in method beginning is a kind of sectional end interpolation, of prediction. which means B-spline interpolation shape of method interpolation is a kindcurve of sectional only depends interpolation, on nearby which points. means Thus, shape curve of will interpolation not be curve influenced only depends by Runge s on nearby phenomenon. points. Thus, massive curve jump willof not both be velocity influenced byacceleration Runge s phenomenon. does not exist massive in jump interpolation of both velocity curve because acceleration of local doesadjustable not exist inproperty, interpolation B-spline curve interpolation because ofcan local fit irregular adjustable property, data without B-spline influence interpolation fitting canof fitor irregular flat areas. data In conclusion, without influence fitting of orbased flat areas. on In conclusion, B-spline interpolation model based improves on B-spline machine s interpolation straightness model improves accuracy to a machine s remarkable straightness level. Y-direction accuracy to a remarkable Z-direction level. straightness Y-direction s in X-direction Z-direction straightness movement s reduce into 2.92 X-direction μm 5.27 movement μm from reduce to μm 2.92 µm80.62 μm µm total fromstraightness µm accuracy µm. has total been straightness improved by accuracy around 95%. has been improved by around 95%. Figure Figure 11c 11c shows total total of of X-direction position in space movement with a stroke of x = of 2000 x = 2000 mm, mm, y = 1000 y = 1000 mm mm z = 200 z = 200 mm. mm. Two Two curves curves respectively respectively show show before before after. after. max value max value before before appears appears at at end end of of stroke, stroke, with with around around 110 µm. 110 After μm., After, total total shows shows a trend a of trend fluctuating of fluctuating enlargement, enlargement, with with max of 9.88 μm. In conclusion, after aiming at dominant max of 9.88 µm. In conclusion, after aiming at dominant components components calculated by VLGSA method, total of X direction in space movement calculated by VLGSA method, total of X direction in space movement reduces around reduces around 90%. remarkable accuracy improvement of machine shows validity of 90%. remarkable accuracy improvement of machine shows validity of using VLGSA using VLGSA in finding key s. Meanwhile, considering VLGSA method is applied to in finding key s. Meanwhile, considering VLGSA method is applied to MBS model, MBS model, model accuracy is also verified. model accuracy is also verified. 5. Conclusions

19 Processes 2020, 8, of Conclusions Structural motional geometric s widely exist in CNC machines. y are most common factors reducing precision machine accuracy, affecting quality of products. refore, characterizing trends, finding most dominant items finally decreasing total geometric play an essential role in development of modern industry. However, considering complexity of precision CNC machine tools, process of targeted can be sophisticated. following works have been done in this paper to achieve goal of improving machine accuracy. 1. A geometric model was established with MBS method. By building topological structure relationship, sophisticated transmitting was decoupled into calculation of several matrices, which can be measured separately. 2. Based on MBS model, items dominant most to total geometric were identified by value leaded global sensitivity analysis (VLGSA), in which actual range of s was considered in process of sensitivity analysis. interaction of different parameters was analyzed in proposed analysis method. 3. Appropriate fitting methods were applied to measured data. data was measured by a high-precision dual-frequency laser interferometer. Considering its excellent ability to describe continual material significant usage in CAD, B-spline interpolation method developed on B-spline curve was used to establish straightness model. 4. A real-time online system was built based on Beckhoff TwinCAT servo system. A experiment was carried out using designed system. straightness in X-direction movement was compensated with four models to verify ability of B-spline interpolation model. X-direction in a space movement of x = 2000 mm, y = 1000 mm z = 200 mm was compensated as a sample. remarkable reduction of geometric in X-direction showed above research was effective. However, this research is based on CNC machine operating at a single speed in cold start temperature. working condition is more sophisticated in actual manufacturing process. Thus, furr works need to concentrate more on relationship of geometric s working condition, such as acceleration changing workforce. Author Contributions: Conceptualization, H.L. Q.C.; methodology, X.Z. Q.C.; software, Q.C. Y.Q.; validation data curation, Q.C. Q.L.; writing-original draft preparation, Q.C.; writing-review editing, Q.C. Y.Z. All authors have read agreed to published version of manuscript. Funding: This research was funded by National Natural Science Foundation of China (Grant no ), Hubei Province Intellectual Property Guidance Development Project (Grant no ) Excellent Dissertation Cultivation Funds of Wuhan University of Technology (Grant no YS-035). Conflicts of Interest: authors declare no conflict of interest. Appendix A For n+1 control points P i (i = 0, 1,..., n) a knot vector U = {u 0, u 1,..., u m }, a feature polygon can be formed by connecting control points. A k + 1 order B-spline curve can be expressed as follow (2 k n + 1): n C k (u) = N i,k (u)p i (A1) i=0 where N i,k (u) is k order B-spline base function, whose recursive function is shown in Equation (A1). N i,0 (u) = N i,k (u) = (u u i)n i,k 1 (u) u i+k u i (A2) { 1, ui u u i+1 0, else + (u i+k+1 u)n i+1,k 1 (u) u i+k+1 u, u i+1 k < u < u k+1

20 Processes 2020, 8, of 21 References 1. Sarafraz, M.M.; Peyghambarzadeh, S.M. Experimental study on subcooled flow boiling heat transfer to water diethylene glycol mixtures as a coolant inside a vertical annulus. Exp. rm. Fluid Science 2013, 50, [CrossRef] 2. Goodarzi, M.; D Orazio, A.; Keshavarzi, A. Develop nano scale method of lattice Boltzmann to predict fluid flow heat transfer of air in inclined lid driven cavity with a large heat source inside, Two case studies: Pure natural convection & mixed convection. Phys. Stat. Mech. Appl. 2018, 509, Sarafraz, M.M. Experimental investigation on pool boiling heat transfer to formic acid, propanol 2-butanol pure liquids under atmospheric pressure. J. Appl. Fluid Mech. 2013, 6, Goodarzi, M.; Amiri, A.; Goodarzi, M.S. Investigation of heat transfer pressure drop of a counter flow corrugated plate heat exchanger using MWCNT based nanofluids. Int. Commun. Heat Mass Transf. 2015, 66, [CrossRef] 5. Chen, J.X.; Lin, S.W.; Zhou, X.L. A comprehensive analysis method for geometric of multi-axis machine tool. Int. J. Mach. Tools Manuf. 2016, 106, [CrossRef] 6. Bai, H.; Shen, J.X. Synsis modeling of CNC machine tool based on multi-body system ory. Aeronaut. Manuf. Technol. 2017, 1, Paul, R.P.; Zhang, H. Computationally efficient kinematics for manipulators with spherical wrists based on homogeneous transformation representation. Int. J. Robot. Res. 1986, 5, [CrossRef] 8. Veitschegger, W.K.; Wu, C.H. Robot accuracy analysis based on kinematics. IEEE J. Robot. Autom. 1986, 2, [CrossRef] 9. Zhao, Q.Q.; Jun, H.; Liu, Z.G. Modeling method on motive axes transfer chain for machine tool of arbitrary topological structure. J. Mech. Eng. 2016, 52, [CrossRef] 10. Zhou, X.T. Research on Machining Accuracy Prediction Based on Machine Tool Error Modeling. Master s sis, Shong University, Shong, China, 20 March Chhatre, S.; Francis, R.; Zhou, Y.H.; Titchener, H.N.; King, J.; Keshavarz, M.E. Global sensitivity analysis for determination of parameter importance in bio-manufacturing processes. Biotechnol. Appl. Biochem. 2008, 51, [CrossRef] 12. Cossarini, G.; Solidoro, C. Global sensitivity analysis of a trophodynamic model of Gulf of Trieste. Ecol. Model. 2008, 212, [CrossRef] 13. Fan, J.W.; Tao, H.H.; Wu, C.J.; Tang, Y.H. Method of sensitivity analysis of machine tools. J. Beijing Univ. Technol. 2019, 45, Yan, Y.; Wang, G.; Zhang, F.P. Precision assembly geometric sensitivity analysis based on transformation model for precision assembly. Trans. Beijing Inst. Technol. 2017, 37, Li, J.; Xue, F.G.; Liu, X.J. Geometric modeling sensitivity analysis of a five-axis machine tool. Int. J. Adv. Manuf. Technol. 2015, 82, [CrossRef] 16. Fan, K.G.; Yang, J.G.; Yang, L.Y. Orthogonal polynomials-based rmally induced spindle geometric modeling. Int. J. Adv. Manuf. Technol. 2013, 65, [CrossRef] 17. Wang, W.; Zhang, Y.; Yang, J.G. Geometric rmal for CNC milling machines based on Newton interpolation method. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2012, 227, [CrossRef] 18. Tuo, Z.Y.; Huang, Y.Q.; Shen, M.W. Modeling of geometric rmal complex position of CNC machine tools based on cubic spline interpolation. J. Shanghai Jiao Tong Univ. 2016, 50, Li, Y.H.; Huang, T.; Derek, G.C. An approach for smooth trajectory planning of high-speed pick--place parallel robots using quantic B-splines. Mech. Mach. ory 2018, 126, [CrossRef] 20. Guo, C.; Yang, L.; Li, Q.Y. Modeling of volumetric s for NC machine tools on multi-body system ories. Mach. Des. Manuf. 2005, 3, Liu, Z.F.; Li, D.D.; Liu, Z.Z. Gantry machining tool assembly method predictive optimization based on multi-body system. Comput. Integr. Manuf. Syst. 2014, 20, Fan, J.W.; Tang, Y.H.; Wang, H.L. Error modeling identification technology for a CNC camshaft grinder. J. Vib. Shock 2017, 36,

21 Processes 2020, 8, of Yu, D.J.; Li, R. Application of Sobol method to sensitivity analysis of a non-linear passive vibration isolators. J. Vib. Eng. 2004, 17, Gao, X.; Tong, H.; Zhou, L. Geometric of CNC machine tools based on G-code modification. Manuf. Technol. Mach. Tool 2015, 1, by authors. Licensee MDPI, Basel, Switzerl. This article is an open access article distributed under terms conditions of Creative Commons Attribution (CC BY) license (

2/80 2

2/80 2 2/80 2 3/80 3 DSP2400 is a high performance Digital Signal Processor (DSP) designed and developed by author s laboratory. It is designed for multimedia and wireless application. To develop application

More information

University of Science and Technology of China A dissertation for master s degree Research of e-learning style for public servants under the context of

University of Science and Technology of China A dissertation for master s degree Research of e-learning style for public servants under the context of 中 国 科 学 技 术 大 学 硕 士 学 位 论 文 新 媒 体 环 境 下 公 务 员 在 线 培 训 模 式 研 究 作 者 姓 名 : 学 科 专 业 : 导 师 姓 名 : 完 成 时 间 : 潘 琳 数 字 媒 体 周 荣 庭 教 授 二 一 二 年 五 月 University of Science and Technology of China A dissertation for

More information

WTO

WTO 10384 200015128 UDC Exploration on Design of CIB s Human Resources System in the New Stage (MBA) 2004 2004 2 3 2004 3 2 0 0 4 2 WTO Abstract Abstract With the rapid development of the high and new technique

More information

THE APPLICATION OF ISOTOPE RATIO ANALYSIS BY INDUCTIVELY COUPLED PLASMA MASS SPECTROMETER A Dissertation Presented By Chaoyong YANG Supervisor: Prof.D

THE APPLICATION OF ISOTOPE RATIO ANALYSIS BY INDUCTIVELY COUPLED PLASMA MASS SPECTROMETER A Dissertation Presented By Chaoyong YANG Supervisor: Prof.D 10384 070302 9825042 UDC 2001.6. 2001.7. 20016 THE APPLICATION OF ISOTOPE RATIO ANALYSIS BY INDUCTIVELY COUPLED PLASMA MASS SPECTROMETER A Dissertation Presented By Chaoyong YANG Supervisor: Prof.Dr. Xiaoru

More information

Microsoft Word - TIP006SCH Uni-edit Writing Tip - Presentperfecttenseandpasttenseinyourintroduction readytopublish

Microsoft Word - TIP006SCH Uni-edit Writing Tip - Presentperfecttenseandpasttenseinyourintroduction readytopublish 我 难 度 : 高 级 对 们 现 不 在 知 仍 道 有 听 影 过 响 多 少 那 次 么 : 研 英 究 过 文 论 去 写 文 时 作 的 表 技 引 示 巧 言 事 : 部 情 引 分 发 言 该 生 使 在 中 用 过 去, 而 现 在 完 成 时 仅 表 示 事 情 发 生 在 过 去, 并 的 哪 现 种 在 时 完 态 成 呢 时? 和 难 过 道 去 不 时 相 关? 是 所 有

More information

Thesis for the Master degree in Engineering Research on Negative Pressure Wave Simulation and Signal Processing of Fluid-Conveying Pipeline Leak Candi

Thesis for the Master degree in Engineering Research on Negative Pressure Wave Simulation and Signal Processing of Fluid-Conveying Pipeline Leak Candi U17 10220 UDC624 Thesis for the Master degree in Engineering Research on Negative Pressure Wave Simulation and Signal Processing of Fluid-Conveying Pipeline Leak Candidate:Chen Hao Tutor: Xue Jinghong

More information

國家圖書館典藏電子全文

國家圖書館典藏電子全文 i ii Abstract The most important task in human resource management is to encourage and help employees to develop their potential so that they can fully contribute to the organization s goals. The main

More information

:1949, 1936, 1713 %, 63 % (, 1957, 5 ), :?,,,,,, (,1999, 329 ),,,,,,,,,, ( ) ; ( ), 1945,,,,,,,,, 100, 1952,,,,,, ,, :,,, 1928,,,,, (,1984, 109

:1949, 1936, 1713 %, 63 % (, 1957, 5 ), :?,,,,,, (,1999, 329 ),,,,,,,,,, ( ) ; ( ), 1945,,,,,,,,, 100, 1952,,,,,, ,, :,,, 1928,,,,, (,1984, 109 2006 9 1949 3 : 1949 2005, : 1949 1978, ; 1979 1997, ; 1998 2005,,, :,,, 1949, :, ;,,,, 50, 1952 1957 ; ; 60 ; 1978 ; 2003,,,,,,, 1953 1978 1953 1978,,,, 100,,,,, 3,, :100836, :wulijjs @263. net ;,, :

More information

untitled

untitled LBS Research and Application of Location Information Management Technology in LBS TP319 10290 UDC LBS Research and Application of Location Information Management Technology in LBS , LBS PDA LBS

More information

A VALIDATION STUDY OF THE ACHIEVEMENT TEST OF TEACHING CHINESE AS THE SECOND LANGUAGE by Chen Wei A Thesis Submitted to the Graduate School and Colleg

A VALIDATION STUDY OF THE ACHIEVEMENT TEST OF TEACHING CHINESE AS THE SECOND LANGUAGE by Chen Wei A Thesis Submitted to the Graduate School and Colleg 上 海 外 国 语 大 学 SHANGHAI INTERNATIONAL STUDIES UNIVERSITY 硕 士 学 位 论 文 MASTER DISSERTATION 学 院 国 际 文 化 交 流 学 院 专 业 汉 语 国 际 教 育 硕 士 题 目 届 别 2010 届 学 生 陈 炜 导 师 张 艳 莉 副 教 授 日 期 2010 年 4 月 A VALIDATION STUDY

More information

CL-68x00,00,00,00,00, CL-78x00,00,00,00,6000 Spindle 181mm mm Spindle bore 181mm is standard. 255,5 or 5mm is option. Chuck is optional. You ca

CL-68x00,00,00,00,00, CL-78x00,00,00,00,6000 Spindle 181mm mm Spindle bore 181mm is standard. 255,5 or 5mm is option. Chuck is optional. You ca CNC KINWA FLAT BED CNC LATHE CL-68x00 CL-68/CL-78 Series Swing 11, mm Between centers 00~7000mm Spindle bore 181, 255, 5, 5mm Spindle motor 60HP CL-68x00,00,00,00,00, CL-78x00,00,00,00,6000 Spindle 181mm

More information

/MPa / kg m - 3 /MPa /MPa 2. 1E ~ 56 ANSYS 6 Hz (a) 一阶垂向弯曲 (b) 一阶侧向弯曲 (c) 一阶扭转 (d) 二阶侧向弯曲 (e) 二阶垂向弯曲 (f) 弯扭组合 2 6 Hz

/MPa / kg m - 3 /MPa /MPa 2. 1E ~ 56 ANSYS 6 Hz (a) 一阶垂向弯曲 (b) 一阶侧向弯曲 (c) 一阶扭转 (d) 二阶侧向弯曲 (e) 二阶垂向弯曲 (f) 弯扭组合 2 6 Hz 31 3 Vol. 31 No. 3 218 9 Journal of Shijiazhuang Tiedao University Natural Science Edition Sep. 218 1 1 2 1 2 1 1. 543 2. 543 U462. 3 217-2 - 16 A 295-373 218 3-63 - 6 1-4 5-7 8-11 1 11 11 398 mm 86 mm

More information

<4D6963726F736F667420576F7264202D20B8DFB9B0B0D3B0D3F5E0D3A6C1A6CAB5B2E2D3EBBCC6CBE3BDE1B9FBB2EED2ECD4ADD2F2B7D6CEF62DD5C5B9FAD0C22E646F6378>

<4D6963726F736F667420576F7264202D20B8DFB9B0B0D3B0D3F5E0D3A6C1A6CAB5B2E2D3EBBCC6CBE3BDE1B9FBB2EED2ECD4ADD2F2B7D6CEF62DD5C5B9FAD0C22E646F6378> 高 拱 坝 坝 踵 应 力 实 测 与 计 算 结 果 差 异 原 因 分 析 张 国 新 周 秋 景 ( 中 国 水 利 水 电 科 学 研 究 院, 北 京 100038) 摘 要 : 坝 踵 应 力 是 关 系 到 混 凝 土 坝 是 否 开 裂 和 安 全 的 一 个 重 要 指 标, 不 管 是 用 结 构 力 学 法 还 是 有 限 元 法 都 能 计 算 出 坝 踵 有 一 定 的 拉

More information

BC04 Module_antenna__ doc

BC04 Module_antenna__ doc http://www.infobluetooth.com TEL:+86-23-68798999 Fax: +86-23-68889515 Page 1 of 10 http://www.infobluetooth.com TEL:+86-23-68798999 Fax: +86-23-68889515 Page 2 of 10 http://www.infobluetooth.com TEL:+86-23-68798999

More information

Microsoft PowerPoint - NCBA_Cattlemens_College_Darrh_B

Microsoft PowerPoint - NCBA_Cattlemens_College_Darrh_B Introduction to Genetics Darrh Bullock University of Kentucky The Model Trait = Genetics + Environment Genetics Additive Predictable effects that get passed from generation to generation Non-Additive Primarily

More information

%

% 38 1 2014 1 Vol. 38No. 1 January 2014 51 Population Research 2010 2010 2010 65 100028 Changing Lineal Families with Three Generations An Analysis of the 2010 Census Data Wang Yuesheng Abstract In contemporary

More information

ENGG1410-F Tutorial 6

ENGG1410-F Tutorial 6 Jianwen Zhao Department of Computer Science and Engineering The Chinese University of Hong Kong 1/16 Problem 1. Matrix Diagonalization Diagonalize the following matrix: A = [ ] 1 2 4 3 2/16 Solution The

More information

[1-3] (Smile) [4] 808 nm (CW) W 1 50% 1 W 1 W Fig.1 Thermal design of semiconductor laser vertical stack ; Ansys 20 bar ; bar 2 25 Fig

[1-3] (Smile) [4] 808 nm (CW) W 1 50% 1 W 1 W Fig.1 Thermal design of semiconductor laser vertical stack ; Ansys 20 bar ; bar 2 25 Fig 40 6 2011 6 Vol.40 No.6 Infrared and Laser Engineering Jun. 2011 808 nm 2000 W 1 1 1 1 2 2 2 2 2 12 (1. 710119 2. 710119) : bar 808 nm bar 100 W 808 nm 20 bar 2 000 W bar LIV bar 808 nm : : TN248.4 TN365

More information

Microsoft Word - 24.doc

Microsoft Word - 24.doc 水 陸 畢 陳 晚 明 飲 食 風 尚 初 探 蕭 慧 媛 桃 園 創 新 技 術 學 院 觀 光 與 休 閒 事 業 管 理 系 摘 要 飲 食 是 人 類 維 持 與 發 展 生 命 的 基 礎 之 一, 飲 食 風 尚 會 隨 著 社 會 地 位 物 質 條 件 以 及 人 為 因 素 轉 移, 不 同 階 層 的 飲 食 方 式, 往 往 標 誌 著 他 們 的 社 會 身 分, 甚 至 反

More information

Microsoft Word doc

Microsoft Word doc 中 考 英 语 科 考 试 标 准 及 试 卷 结 构 技 术 指 标 构 想 1 王 后 雄 童 祥 林 ( 华 中 师 范 大 学 考 试 研 究 院, 武 汉,430079, 湖 北 ) 提 要 : 本 文 从 结 构 模 式 内 容 要 素 能 力 要 素 题 型 要 素 难 度 要 素 分 数 要 素 时 限 要 素 等 方 面 细 致 分 析 了 中 考 英 语 科 试 卷 结 构 的

More information

% GIS / / Fig. 1 Characteristics of flood disaster variation in suburbs of Shang

% GIS / / Fig. 1 Characteristics of flood disaster variation in suburbs of Shang 20 6 2011 12 JOURNAL OF NATURAL DISASTERS Vol. 20 No. 6 Dec. 2011 1004-4574 2011 06-0094 - 05 200062 1949-1990 1949 1977 0. 8 0. 03345 0. 01243 30 100 P426. 616 A Risk analysis of flood disaster in Shanghai

More information

致 谢 本 人 自 2008 年 6 月 从 上 海 外 国 语 大 学 毕 业 之 后, 于 2010 年 3 月 再 次 进 入 上 外, 非 常 有 幸 成 为 汉 语 国 际 教 育 专 业 的 研 究 生 回 顾 三 年 以 来 的 学 习 和 生 活, 顿 时 感 觉 这 段 时 间 也

致 谢 本 人 自 2008 年 6 月 从 上 海 外 国 语 大 学 毕 业 之 后, 于 2010 年 3 月 再 次 进 入 上 外, 非 常 有 幸 成 为 汉 语 国 际 教 育 专 业 的 研 究 生 回 顾 三 年 以 来 的 学 习 和 生 活, 顿 时 感 觉 这 段 时 间 也 精 英 汉 语 和 新 实 用 汉 语 课 本 的 对 比 研 究 The Comparative Study of Jing Ying Chinese and The New Practical Chinese Textbook 专 业 : 届 别 : 姓 名 : 导 师 : 汉 语 国 际 教 育 2013 届 王 泉 玲 杨 金 华 1 致 谢 本 人 自 2008 年 6 月 从 上 海 外

More information

10384 199928010 UDC 2002 4 2002 6 2002 2002 4 DICOM DICOM 1. 2. 3. Canny 4. 5. DICOM DICOM DICOM DICOM I Abstract Eyes are very important to our lives. Biologic parameters of anterior segment are criterions

More information

by industrial structure evolution from 1952 to 2007 and its influence effect was first acceleration and then deceleration second the effects of indust

by industrial structure evolution from 1952 to 2007 and its influence effect was first acceleration and then deceleration second the effects of indust 2011 2 1 1 2 3 4 1. 100101 2. 100124 3. 100039 4. 650092 - - - 3 GDP U 20-30 60% 10% TK01 A 1002-9753 2011 02-0042 - 10 Analysis on Character and Potential of Energy Saving and Carbon Reducing by Structure

More information

Microsoft Word - A200810-897.doc

Microsoft Word - A200810-897.doc 基 于 胜 任 特 征 模 型 的 结 构 化 面 试 信 度 和 效 度 验 证 张 玮 北 京 邮 电 大 学 经 济 管 理 学 院, 北 京 (100876) E-mail: weeo1984@sina.com 摘 要 : 提 高 结 构 化 面 试 信 度 和 效 度 是 面 试 技 术 研 究 的 核 心 内 容 近 年 来 国 内 有 少 数 学 者 探 讨 过 基 于 胜 任 特 征

More information

Microsoft Word - 24217010311110028谢雯雯.doc

Microsoft Word - 24217010311110028谢雯雯.doc HUAZHONG AGRICULTURAL UNIVERSITY 硕 士 学 位 论 文 MASTER S DEGREE DISSERTATION 80 后 女 硕 士 生 择 偶 现 状 以 武 汉 市 七 所 高 校 为 例 POST-80S FEMALE POSTGRADUATE MATE SELECTION STATUS STUDY TAKE WUHAN SEVEN UNIVERSITIES

More information

Microsoft PowerPoint _代工實例-1

Microsoft PowerPoint _代工實例-1 4302 動態光散射儀 (Dynamic Light Scattering) 代工實例與結果解析 生醫暨非破壞性分析團隊 2016.10 updated Which Size to Measure? Diameter Many techniques make the useful and convenient assumption that every particle is a sphere. The

More information

國立中山大學學位論文典藏.PDF

國立中山大學學位論文典藏.PDF 國 立 中 山 大 學 企 業 管 理 學 系 碩 士 論 文 以 系 統 動 力 學 建 構 美 食 餐 廳 異 國 麵 坊 之 管 理 飛 行 模 擬 器 研 究 生 : 簡 蓮 因 撰 指 導 教 授 : 楊 碩 英 博 士 中 華 民 國 九 十 七 年 七 月 致 謝 詞 寫 作 論 文 的 過 程 是 一 段 充 滿 艱 辛 與 淚 水 感 動 與 窩 心 的 歷 程, 感 謝 這 一

More information

Settlement Equation " H = CrH 1+ e o log p' o + ( p' p' c o! p' o ) CcH + 1+ e o log p' c + p' f! ( p' p' c c! p' o ) where ΔH = consolidation settlem

Settlement Equation  H = CrH 1+ e o log p' o + ( p' p' c o! p' o ) CcH + 1+ e o log p' c + p' f! ( p' p' c c! p' o ) where ΔH = consolidation settlem Prediction of Compression and Recompression Indices of Texas Overconsolidated Clays Presented By: Sayeed Javed, Ph.D., P.E. Settlement Equation " H = CrH 1+ e o log p' o + ( p' p' c o! p' o ) CcH + 1+

More information

10384 19020101152519 UDC Rayleigh Quasi-Rayleigh Method for computing eigenvalues of symmetric tensors 2 0 1 3 2 0 1 3 2 0 1 3 2013 , 1. 2. [4], [27].,. [6] E- ; [7], Z-. [15]. Ramara G. kolda [1, 2],

More information

A dissertation for Master s degree Metro Indoor Coverage Systems Analysis And Design Author s Name: Sheng Hailiang speciality: Supervisor:Prof.Li Hui,

A dissertation for Master s degree Metro Indoor Coverage Systems Analysis And Design Author s Name: Sheng Hailiang speciality: Supervisor:Prof.Li Hui, 中 国 科 学 技 术 大 学 工 程 硕 士 学 位 论 文 地 铁 内 移 动 通 信 室 内 覆 盖 分 析 及 应 用 作 者 姓 名 : 学 科 专 业 : 盛 海 亮 电 子 与 通 信 导 师 姓 名 : 李 辉 副 教 授 赵 红 媛 高 工 完 成 时 间 : 二 八 年 三 月 十 日 University of Science and Technology of Ch A dissertation

More information

A Study on the Relationships of the Co-construction Contract A Study on the Relationships of the Co-Construction Contract ( ) ABSTRACT Co-constructio in the real estate development, holds the quite

More information

一次辽宁暴雨过程的诊断及风场反演分析

一次辽宁暴雨过程的诊断及风场反演分析 Climate Change Research Letters 气 候 变 化 研 究 快 报, 2013, 2, 139-146 http://dx.doi.org/10.12677/ccrl.2013.24024 Published Online October 2013 (http://www.hanspub.org/journal/ccrl.html) Analysis of the Diagnosis

More information

地質調査研究報告/Bulletin of the Geological Survey of Japan

地質調査研究報告/Bulletin of the Geological Survey of Japan Shigeru Suto, Takayuki Inomata, Hisashi Sasaki and Sakae Mukoyama (2007) Data base of the volcanic ash fall distribution map of Japan. Bull. Geol. Surv. Japan, vol. 58(9/10), p.261-321, 8 figs, 2 tables,

More information

Microsoft PowerPoint - STU_EC_Ch08.ppt

Microsoft PowerPoint - STU_EC_Ch08.ppt 樹德科技大學資訊工程系 Chapter 8: Counters Shi-Huang Chen Fall 2010 1 Outline Asynchronous Counter Operation Synchronous Counter Operation Up/Down Synchronous Counters Design of Synchronous Counters Cascaded Counters

More information

Microsoft Word - 01李惠玲ok.doc

Microsoft Word - 01李惠玲ok.doc 康 寧 學 報 11:1-20(2009) 1 數 位 學 習 於 護 理 技 術 課 程 之 運 用 與 評 值 * 李 惠 玲 ** 高 清 華 *** 呂 莉 婷 摘 要 背 景 : 網 路 科 技 在 教 育 的 使 用 已 成 為 一 種 有 利 的 教 學 輔 助 工 具 網 路 教 學 的 特 性, 在 使 學 習 可 不 分 時 間 與 空 間 不 同 進 度 把 握 即 時 性 資

More information

168 健 等 木醋对几种小浆果扦插繁殖的影响 第1期 the view of the comprehensive rooting quality, spraying wood vinegar can change rooting situation, and the optimal concent

168 健 等 木醋对几种小浆果扦插繁殖的影响 第1期 the view of the comprehensive rooting quality, spraying wood vinegar can change rooting situation, and the optimal concent 第 31 卷 第 1 期 2013 年 3 月 经 济 林 研 究 Nonwood Forest Research Vol. 31 No.1 Mar. 2013 木醋对几种小浆果扦插繁殖的影响 健 1,2 杨国亭 1 刘德江 2 (1. 东北林业大学 生态研究中心 黑龙江 哈尔滨 150040 2. 佳木斯大学 生命科学学院 黑龙江 佳木斯 154007) 摘 要 为了解决小浆果扦插繁殖中生根率及成活率低等问题

More information

1 引言

1 引言 P P 第 40 卷 Vol.40 第 7 期 No.7 计 算 机 工 程 Computer Engineering 014 年 7 月 July 014 开 发 研 究 与 工 程 应 用 文 章 编 号 :1000-348(014)07-081-05 文 献 标 识 码 :A 中 图 分 类 号 :TP391.41 摘 基 于 图 像 识 别 的 震 象 云 地 震 预 测 方 法 谢 庭,

More information

g 100mv /g 0. 5 ~ 5kHz 1 YSV8116 DASP 1 N 2. 2 [ M] { x } + [ C] { x } + [ K]{ x } = { f t } 1 M C K 3 M C K f t x t 1 [ H( ω )] = - ω 2

g 100mv /g 0. 5 ~ 5kHz 1 YSV8116 DASP 1 N 2. 2 [ M] { x } + [ C] { x } + [ K]{ x } = { f t } 1 M C K 3 M C K f t x t 1 [ H( ω )] = - ω 2 10 2016 10 No. 10 Modular Machine Tool & Automatic Manufacturing Technique Oct. 2016 1001-2265 2016 10-0012 - 05 DOI 10. 13462 /j. cnki. mmtamt. 2016. 10. 004 * 116024 MIMO TH166 TG502 A Dynamic Performance

More information

Microsoft Word - 刘藤升答辩修改论文.doc

Microsoft Word - 刘藤升答辩修改论文.doc 武 汉 体 育 学 院 硕 士 学 位 论 文 ( 全 日 制 硕 士 ) 社 会 需 求 视 角 下 武 汉 体 院 乒 乓 球 硕 士 研 究 生 就 业 状 况 研 究 研 究 生 : 刘 藤 升 导 师 : 滕 守 刚 副 教 授 专 业 : 体 育 教 育 训 练 学 研 究 方 向 : 乒 乓 球 教 学 训 练 理 论 与 方 法 2015 年 5 月 分 类 号 : G8 学 校 代

More information

Microsoft PowerPoint - ATF2015.ppt [相容模式]

Microsoft PowerPoint - ATF2015.ppt [相容模式] Improving the Video Totalized Method of Stopwatch Calibration Samuel C.K. Ko, Aaron Y.K. Yan and Henry C.K. Ma The Government of Hong Kong Special Administrative Region (SCL) 31 Oct 2015 1 Contents Introduction

More information

PowerPoint Presentation

PowerPoint Presentation Decision analysis 量化決策分析方法專論 2011/5/26 1 Problem formulation- states of nature In the decision analysis, decision alternatives are referred to as chance events. The possible outcomes for a chance event

More information

PCA+LDA 14 1 PEN mL mL mL 16 DJX-AB DJ X AB DJ2 -YS % PEN

PCA+LDA 14 1 PEN mL mL mL 16 DJX-AB DJ X AB DJ2 -YS % PEN 21 11 2011 11 COMPUTER TECHNOLOGY AND DEVELOPMENT Vol. 21 No. 11 Nov. 2011 510006 PEN3 5 PCA + PCA+LDA 5 5 100% TP301 A 1673-629X 2011 11-0177-05 Application of Electronic Nose in Discrimination of Different

More information

IPCC CO (IPCC2006) 1 : = ( 1) 1 (kj/kg) (kgc/gj) (tc/t)

IPCC CO (IPCC2006) 1 : = ( 1) 1 (kj/kg) (kgc/gj) (tc/t) 2011 5 5 (278 ) China Industrial Economics May 2011 No.5 1 12 (1. 100005; 2. 066004) [ ] : ; ; : ; ; [ ] ; ; ; [ ]F290 [ ]A [ ]1006-480X(2011)05-0047-11 2008 CO 2 ( ) (2009) (GDP) (Binhocker et al. 2008)

More information

10 ( ) ( ) [5] 1978 : [1] (P13) [6] [1] (P217) [7] [1] (P19) : : [1] [4] (P1347) (P18) 1985 : [1] (P343) 1300 : [1] (P12) 1984 :

10 ( ) ( ) [5] 1978 : [1] (P13) [6] [1] (P217) [7] [1] (P19) : : [1] [4] (P1347) (P18) 1985 : [1] (P343) 1300 : [1] (P12) 1984 : 27 3 ( ) Vol.27 No.3 2010 5 Journal of Shenzhen University (Humanities & Social Sciences) May 2010 ( 518060) : ; : ; : ; ; ; ; : F 127.9 :A :1000-260X(2010)03-0009-13 30 [2] : [2] (P381) 1979 : : [3] :1978

More information

Construction of Chinese pediatric standard database A Dissertation Submitted for the Master s Degree Candidate:linan Adviser:Prof. Han Xinmin Nanjing

Construction of Chinese pediatric standard database A Dissertation Submitted for the Master s Degree Candidate:linan Adviser:Prof. Han Xinmin Nanjing 密 级 : 公 开 学 号 :20081209 硕 士 学 位 论 文 中 医 儿 科 标 准 数 据 库 建 设 研 究 研 究 生 李 楠 指 导 教 师 学 科 专 业 所 在 学 院 毕 业 时 间 韩 新 民 教 授 中 医 儿 科 学 第 一 临 床 医 学 院 2011 年 06 月 Construction of Chinese pediatric standard database

More information

畜牧 动物医学 蚕 蜂

畜牧 动物医学 蚕 蜂 中国农学通报 2015,31(20):181-185 Chinese Agricultural Science Bulletin 老河口市砂梨花期的物候特点分析 邓 环 1 马德栗 1 徐洪海 2 洪国平 1 2 武汉区域气候中心 武汉 430074 老河口气象局 湖北老河口 441800 1 摘 要 为掌握湖北省老河口市砂梨花期的物候规律 利用老河口市 20 年的砂梨裂蕾期 初花期 和盛花期 4

More information

m m m ~ mm

m m m ~ mm 2011 10 10 157 JOURNAL OF RAILWAY ENGINEERING SOCIETY Oct 2011 NO. 10 Ser. 157 1006-2106 2011 10-0007 - 0124-05 710043 6 TBM TBM U455. 43 A Structural Calculation and Analysis of Transfer Node of Three

More information

附件1:

附件1: 附 件 1: 全 国 优 秀 教 育 硕 士 专 业 学 位 论 文 推 荐 表 单 位 名 称 : 西 南 大 学 论 文 题 目 填 表 日 期 :2014 年 4 月 30 日 数 学 小 组 合 作 学 习 的 课 堂 管 理 攻 硕 期 间 及 获 得 硕 士 学 位 后 一 年 内 获 得 与 硕 士 学 位 论 文 有 关 的 成 果 作 者 姓 名 论 文 答 辩 日 期 学 科 专

More information

Microsoft PowerPoint - Aqua-Sim.pptx

Microsoft PowerPoint - Aqua-Sim.pptx Peng Xie, Zhong Zhou, Zheng Peng, Hai Yan, Tiansi Hu, Jun-Hong Cui, Zhijie Shi, Yunsi Fei, Shengli Zhou Underwater Sensor Network Lab 1 Outline Motivations System Overview Aqua-Sim Components Experimental

More information

mm 5 1 Tab 1 Chemical composition of PSB830 finishing rolled rebars % C Si Mn P S V 0 38 ~ 1 50 ~ 0 80 ~ ~

mm 5 1 Tab 1 Chemical composition of PSB830 finishing rolled rebars % C Si Mn P S V 0 38 ~ 1 50 ~ 0 80 ~ ~ PSB830 365000 32 mm PSB830 PSB830 TG 335 64 A Productive Practition of PSB830 Finishing Rolled Rebars PAN Jianzhou Bar Steel Rolling Minguang Co Ltd of Fujian Sansteel Sanming 365000 China Abstract High

More information

http / /yxxy. cbpt. cnki. net / % % %

http / /yxxy. cbpt. cnki. net / % % % 2017 3 Mar. 2017 5 2 Chongqing Higher Education Research Vol. 5 No. 2 DOI 10. 15998 /j. cnki. issn1673-8012. 2017. 02. 006 230039 2011 2015 2016 G649. 21 A 1673-8012 2017 02-0037-11 2017-01-03 2015zdjy024

More information

國 史 館 館 刊 第 23 期 Chiang Ching-kuo s Educational Innovation in Southern Jiangxi and Its Effects (1941-1943) Abstract Wen-yuan Chu * Chiang Ching-kuo wa

國 史 館 館 刊 第 23 期 Chiang Ching-kuo s Educational Innovation in Southern Jiangxi and Its Effects (1941-1943) Abstract Wen-yuan Chu * Chiang Ching-kuo wa 國 史 館 館 刊 第 二 十 三 期 (2010 年 3 月 ) 119-164 國 史 館 1941-1943 朱 文 原 摘 要 1 關 鍵 詞 : 蔣 經 國 贛 南 學 校 教 育 社 會 教 育 掃 盲 運 動 -119- 國 史 館 館 刊 第 23 期 Chiang Ching-kuo s Educational Innovation in Southern Jiangxi and

More information

國立中山大學學位論文典藏.PDF

國立中山大學學位論文典藏.PDF I II III The Study of Factors to the Failure or Success of Applying to Holding International Sport Games Abstract For years, holding international sport games has been Taiwan s goal and we are on the way

More information

Microsoft Word - 刘 慧 板.doc

Microsoft Word - 刘  慧 板.doc 中 国 环 境 科 学 2012,32(5):933~941 China Environmental Science 系 统 动 力 学 在 空 港 区 域 规 划 环 境 影 响 评 价 中 的 应 用 刘 慧 1,2, 郭 怀 成 1*, 盛 虎 1, 都 小 尚 1,3, 李 娜 1 1, 杨 永 辉 (1. 北 京 大 学 环 境 科 学 与 工 程 学 院, 北 京 100871; 2.

More information

% % % % % % ~

% % % % % % ~ 1001-5558 2015 03-0021-16 2010 C91 A 2014 5 2010 N. W. Journal of Ethnology 2015 3 86 2015.No.3 Total No.86 2010 2010 2181.58 882.99 40.47% 1298.59 59.53% 2013 2232.78 847.29 37.95% 1385.49 62.05% 1990

More information

JOURNAL OF EARTHQUAKE ENGINEERING AND ENGINEERING VIBRATION Vol. 31 No. 5 Oct /35 TU3521 P315.

JOURNAL OF EARTHQUAKE ENGINEERING AND ENGINEERING VIBRATION Vol. 31 No. 5 Oct /35 TU3521 P315. 31 5 2011 10 JOURNAL OF EARTHQUAKE ENGINEERING AND ENGINEERING VIBRATION Vol. 31 No. 5 Oct. 2011 1000-1301 2011 05-0075 - 09 510405 1 /35 TU3521 P315. 8 A Earthquake simulation shaking table test and analysis

More information

KUKA W. Polini L. Sorrentino Aized Shirinzadeh 6 7 MF Tech Pitbull Fox Taniq Scorpo Scorpo Compositum Windows KUKA 1 P 1 P 2 KU

KUKA W. Polini L. Sorrentino Aized Shirinzadeh 6 7 MF Tech Pitbull Fox Taniq Scorpo Scorpo Compositum Windows KUKA 1 P 1 P 2 KU 23 5 2018 10 Vol. 23 No. 5 JOURNAL OF HARBIN UNIVERSITY OF SCIENCE AND TECHNOLOGY Oct. 2018 1 1 2 2 1. 150080 2. 150080 Matlab AD- AMS CAD DOI 10. 15938 /j. jhust. 2018. 05. 016 TP391. 9 A 1007-2683 2018

More information

10384 200115009 UDC Management Buy-outs MBO MBO MBO 2002 MBO MBO MBO MBO 000527 MBO MBO MBO MBO MBO MBO MBO MBO MBO MBO MBO Q MBO MBO MBO Abstract Its related empirical study demonstrates a remarkable

More information

(baking powder) 1 ( ) ( ) 1 10g g (two level design, D-optimal) 32 1/2 fraction Two Level Fractional Factorial Design D-Optimal D

(baking powder) 1 ( ) ( ) 1 10g g (two level design, D-optimal) 32 1/2 fraction Two Level Fractional Factorial Design D-Optimal D ( ) 4 1 1 1 145 1 110 1 (baking powder) 1 ( ) ( ) 1 10g 1 1 2.5g 1 1 1 1 60 10 (two level design, D-optimal) 32 1/2 fraction Two Level Fractional Factorial Design D-Optimal Design 1. 60 120 2. 3. 40 10

More information

85% NCEP CFS 10 CFS CFS BP BP BP ~ 15 d CFS BP r - 1 r CFS 2. 1 CFS 10% 50% 3 d CFS Cli

85% NCEP CFS 10 CFS CFS BP BP BP ~ 15 d CFS BP r - 1 r CFS 2. 1 CFS 10% 50% 3 d CFS Cli 1 2 3 1. 310030 2. 100054 3. 116000 CFS BP doi 10. 13928 /j. cnki. wrahe. 2016. 04. 020 TV697. 1 A 1000-0860 2016 04-0088-05 Abandoned water risk ratio control-based reservoir pre-discharge control method

More information

9330.doc

9330.doc The research of the ecotourism operated by the cooperative operating system in northern Tapajen Mountain The research of the ecotourism operated by the cooperative operating system in northern Tapajen

More information

從詩歌的鑒賞談生命價值的建構

從詩歌的鑒賞談生命價值的建構 Viktor E. Frankl (logotherapy) (will-to-meaning) (creative values) Ture (Good) (Beauty) (experiential values) (attitudinal values) 1 2 (logotherapy) (biological) (2) (psychological) (3) (noölogical) (4)

More information

Public Projects A Thesis Submitted to Department of Construction Engineering National Kaohsiung First University of Science and Technology In Partial

Public Projects A Thesis Submitted to Department of Construction Engineering National Kaohsiung First University of Science and Technology In Partial Public Projects A Thesis Submitted to Department of Construction Engineering National Kaohsiung First University of Science and Technology In Partial Fulfillment of the Requirements For the Degree of Master

More information

~ ~

~ ~ * 40 4 2016 7 Vol. 40 No. 4 July 2016 35 Population Research 2014 1 2016 2016 9101. 0 40 49. 6% 2017 ~ 2021 1719. 5 160 ~ 470 100872 Accumulated Couples and Extra Births under the Universal Tw o-child

More information

1

1 Activity- based Cost Management: A New Mode of Medical cost Management () 1 Activity - based Cost Management A New Mode of Medical cost Management Abstract With the development of medical market, the defects

More information

三維空間之機械手臂虛擬實境模擬

三維空間之機械手臂虛擬實境模擬 VRML Model of 3-D Robot Arm VRML Model of 3-D Robot Arm MATLAB VRML MATLAB Simulink i MATLAB Simulink V-Realm Build Joystick ii Abstract The major purpose of this thesis presents the procedure of VRML

More information

1 ABSTRACT 3 5 13 14 14 15 15 16 17 18 18 18 18-21 28 33 34 46 47 53 54 2

1 ABSTRACT 3 5 13 14 14 15 15 16 17 18 18 18 18-21 28 33 34 46 47 53 54 2 : - S99017212 O O 1 1 ABSTRACT 3 5 13 14 14 15 15 16 17 18 18 18 18-21 28 33 34 46 47 53 54 2 - A L Aeq CNE 463 3-5 L Aeq L Aeq CNE WHO/ISH 1999 [49] 18.6kPa / 12.0kPa EpiInfo6.04 SPSS10.0 Logistic 10

More information

M M. 20

M M. 20 37 1 Vol. 37 No.1 2 0 1 6 1 TSINGHUA JOURNAL OF EDUCATION Jan. 2 0 1 6 4. 0 100872 1. 0 2. 0 3. 0 4. 0 4. 0 4. 0 G640 A 1001-4519 2016 01-0006 - 10 DOI 10. 14138 /j. 1001-4519. 2016. 01. 000610 11-12 18

More information

IP TCP/IP PC OS µclinux MPEG4 Blackfin DSP MPEG4 IP UDP Winsock I/O DirectShow Filter DirectShow MPEG4 µclinux TCP/IP IP COM, DirectShow I

IP TCP/IP PC OS µclinux MPEG4 Blackfin DSP MPEG4 IP UDP Winsock I/O DirectShow Filter DirectShow MPEG4 µclinux TCP/IP IP COM, DirectShow I 2004 5 IP TCP/IP PC OS µclinux MPEG4 Blackfin DSP MPEG4 IP UDP Winsock I/O DirectShow Filter DirectShow MPEG4 µclinux TCP/IP IP COM, DirectShow I Abstract The techniques of digital video processing, transferring

More information

θ 1 = φ n -n 2 2 n AR n φ i = 0 1 = a t - θ θ m a t-m 3 3 m MA m 1. 2 ρ k = R k /R 0 5 Akaike ρ k 1 AIC = n ln δ 2

θ 1 = φ n -n 2 2 n AR n φ i = 0 1 = a t - θ θ m a t-m 3 3 m MA m 1. 2 ρ k = R k /R 0 5 Akaike ρ k 1 AIC = n ln δ 2 35 2 2012 2 GEOMATICS & SPATIAL INFORMATION TECHNOLOGY Vol. 35 No. 2 Feb. 2012 1 2 3 4 1. 450008 2. 450005 3. 450008 4. 572000 20 J 101 20 ARMA TU196 B 1672-5867 2012 02-0213 - 04 Application of Time Series

More information

UDC Empirical Researches on Pricing of Corporate Bonds with Macro Factors 厦门大学博硕士论文摘要库

UDC Empirical Researches on Pricing of Corporate Bonds with Macro Factors 厦门大学博硕士论文摘要库 10384 15620071151397 UDC Empirical Researches on Pricing of Corporate Bonds with Macro Factors 2010 4 Duffee 1999 AAA Vasicek RMSE RMSE Abstract In order to investigate whether adding macro factors

More information

國立中山大學學位論文典藏.PDF

國立中山大學學位論文典藏.PDF 中 國 文 學 系 國 立 中 山 大 學, 碩 士 論 文 國 立 中 山 大 學 中 國 文 學 系 碩 士 論 文 Department of Chinese Literature 肉 蒲 團 研 究 National Sun Yat-sen University Master Thesis 肉 蒲 團 研 究 The Research of Rou Pu Tuan 研 究 生 : 林 欣 穎

More information

Revit Revit Revit BIM BIM 7-9 3D 1 BIM BIM 6 Revit 0 4D 1 2 Revit Revit 2. 1 Revit Revit Revit Revit 2 2 Autodesk Revit Aut

Revit Revit Revit BIM BIM 7-9 3D 1 BIM BIM 6 Revit 0 4D 1 2 Revit Revit 2. 1 Revit Revit Revit Revit 2 2 Autodesk Revit Aut 60 2 2016 2 RAILWAY STANDARD DESIGN Vol. 60 No. 2 Feb. 2016 1004-2954201602-0071-06 BIM 1 1 2 2 1 1. 7140992. 710054 BIM BIM 3D 4D nd BIM 1 3D 4D Revit BIM BIM U442. 5TP391. 72 A DOI10. 13238 /j. issn.

More information

Microsoft Word - 专论综述1.doc

Microsoft Word - 专论综述1.doc 1 基 于 协 同 过 滤 的 高 考 志 愿 推 荐 系 统 徐 兰 静, 李 珊, 严 钊 ( 南 京 航 空 航 天 大 学 经 济 与 管 理 学 院, 南 京 211100) 摘 要 : 近 年 来 信 息 过 载 问 题 的 出 现 使 得 个 性 化 推 荐 技 术 应 运 而 生, 其 中 协 同 过 滤 推 荐 技 术 通 过 在 用 户 和 信 息 之 间 建 立 联 系, 被

More information

HC50246_2009

HC50246_2009 Page: 1 of 7 Date: June 2, 2009 WINMATE COMMUNICATION INC. 9 F, NO. 111-6, SHING-DE RD., SAN-CHUNG CITY, TAIPEI, TAIWAN, R.O.C. The following merchandise was submitted and identified by the vendor as:

More information

HC20131_2010

HC20131_2010 Page: 1 of 8 Date: April 14, 2010 WINMATE COMMUNICATION INC. 9 F, NO. 111-6, SHING-DE RD., SAN-CHUNG CITY, TAIPEI, TAIWAN, R.O.C. The following merchandise was submitted and identified by the vendor as:

More information

/3 CAD JPG GIS CAD GIS GIS 1 a CAD CAD CAD GIS GIS ArcGIS 9. x 10 1 b 1112 CAD GIS 1 c R2VArcscan CAD MapGIS CAD 1 d CAD U

/3 CAD JPG GIS CAD GIS GIS 1 a CAD CAD CAD GIS GIS ArcGIS 9. x 10 1 b 1112 CAD GIS 1 c R2VArcscan CAD MapGIS CAD 1 d CAD U 1006-3862 2010 05-0059 - 07 361005 1 GIS 2 3 What if 2. 0 1 2 3 4 GIS TU984. 11 A 1 Planning Support System MIS PSS 1989 1 90 23 4-7 GIS Planning Support GIS System SDSS PSS GIS GIS CAD GIS SDSS CAD CAD

More information

untitled

untitled Co-integration and VECM Yi-Nung Yang CYCU, Taiwan May, 2012 不 列 1 Learning objectives Integrated variables Co-integration Vector Error correction model (VECM) Engle-Granger 2-step co-integration test Johansen

More information

致 谢 开 始 这 篇 致 谢 的 时 候, 以 为 这 是 最 轻 松 最 愉 快 的 部 分, 而 此 时 心 头 却 充 满 了 沉 甸 甸 的 回 忆 和 感 恩, 一 时 间 竟 无 从 下 笔 虽 然 这 远 不 是 一 篇 完 美 的 论 文, 但 完 成 这 篇 论 文 要 感 谢

致 谢 开 始 这 篇 致 谢 的 时 候, 以 为 这 是 最 轻 松 最 愉 快 的 部 分, 而 此 时 心 头 却 充 满 了 沉 甸 甸 的 回 忆 和 感 恩, 一 时 间 竟 无 从 下 笔 虽 然 这 远 不 是 一 篇 完 美 的 论 文, 但 完 成 这 篇 论 文 要 感 谢 中 国 科 学 技 术 大 学 博 士 学 位 论 文 论 文 课 题 : 一 个 新 型 简 易 电 子 直 线 加 速 器 的 关 键 技 术 研 究 学 生 姓 名 : 导 师 姓 名 : 单 位 名 称 : 专 业 名 称 : 研 究 方 向 : 完 成 时 间 : 谢 家 麟 院 士 王 相 綦 教 授 国 家 同 步 辐 射 实 验 室 核 技 术 及 应 用 加 速 器 物 理 2006

More information

2012 2 157 No. 2 2012 Jinan Journal Philosophy and Social Sciences Sum No. 157 100871 I127. 41 A 1000-5072 2012 02-0002 - 10 4 1 1 2 1540 1503 2 2011-06 - 27 1969 1 1985 19 2 3 1997 1998 3 2 34 2 3 1540

More information

Untitled-3

Untitled-3 SEC.. Separable Equations In each of problems 1 through 8 solve the given differential equation : ü 1. y ' x y x y, y 0 fl y - x 0 fl y - x 0 fl y - x3 3 c, y 0 ü. y ' x ^ y 1 + x 3 x y 1 + x 3, y 0 fl

More information

课题调查对象:

课题调查对象: 1 大 陆 地 方 政 府 大 文 化 管 理 职 能 与 机 构 整 合 模 式 比 较 研 究 武 汉 大 学 陈 世 香 [ 内 容 摘 要 ] 迄 今 为 止, 大 陆 地 方 政 府 文 化 管 理 体 制 改 革 已 经 由 试 点 改 革 进 入 到 全 面 推 行 阶 段 本 文 主 要 通 过 结 合 典 型 调 查 法 与 比 较 研 究 方 法, 对 已 经 进 行 了 政 府

More information

mode of puzzle-solving

mode of puzzle-solving 91 12 145 174 * * 146 1 1 mode of puzzle-solving 91 12 147 83-105 148 2 3 2 3 151 91 12 149 150 4 4 101-104 91 12 151 identity 5 6 7 5 6 7 100 140 152 8 9 10 8 31-32 9 27-29 10 sense of political efficacy

More information

chr82-06DuanW.hwp

chr82-06DuanW.hwp 上 海 小 三 線 建 設 在 縣 域 分 佈 特 點 的 歷 史 地 理 考 察 - 以 安 徽 省 寧 國 縣 为 例 - 段 偉 ( 中 國 復 旦 大 學 ) Ⅰ. 序 言 Ⅱ. 上 海 小 三 線 建 設 的 選 址 與 寧 國 縣 的 區 位 Ⅲ. 小 三 線 企 事 業 單 位 在 寧 國 縣 的 分 佈 Ⅳ. 改 造 後 的 寧 國 縣 小 三 線 建 設 與 城 鎮 發 展 Ⅴ.

More information

Microsoft Word - KSAE06-S0262.doc

Microsoft Word - KSAE06-S0262.doc Stereo Vision based Forward Collision Warning and Avoidance System Yunhee LeeByungjoo KimHogi JungPaljoo Yoon Central R&D Center, MANDO Corporation, 413-5, Gomae-Ri, Gibeung-Eub, Youngin-Si, Kyonggi-Do,

More information

我国原奶及乳制品安全生产和质量安全管理研究

我国原奶及乳制品安全生产和质量安全管理研究 密 级 论 文 编 号 中 国 农 业 科 学 院 硕 士 学 位 论 文 我 国 原 奶 及 乳 制 品 质 量 安 全 管 理 研 究 Study on Quality and Safety Management of Raw Milk and Dairy Products in China 申 请 人 : 段 成 立 指 导 教 师 : 叶 志 华 研 究 员 张 蕙 杰 研 究 员 申 请

More information

热设计网

热设计网 例 例 Agenda Popular Simulation software in PC industry * CFD software -- Flotherm * Advantage of Flotherm Flotherm apply to Cooler design * How to build up the model * Optimal parameter in cooler design

More information

1990-1997 1980-1997 Abstract The relationship between human resource development and economy increase has been an important research item with the coming of knowledge economy. In this paper, ahead item

More information

Microsoft PowerPoint - talk8.ppt

Microsoft PowerPoint - talk8.ppt Adaptive Playout Scheduling Using Time-scale Modification Yi Liang, Nikolaus Färber Bernd Girod, Balaji Prabhakar Outline QoS concerns and tradeoffs Jitter adaptation as a playout scheduling scheme Packet

More information

硕 士 学 位 论 文 论 文 题 目 : 北 岛 诗 歌 创 作 的 双 重 困 境 专 业 名 称 : 中 国 现 当 代 文 学 研 究 方 向 : 中 国 新 诗 研 究 论 文 作 者 : 奚 荣 荣 指 导 老 师 : 姜 玉 琴 2014 年 12 月

硕 士 学 位 论 文 论 文 题 目 : 北 岛 诗 歌 创 作 的 双 重 困 境 专 业 名 称 : 中 国 现 当 代 文 学 研 究 方 向 : 中 国 新 诗 研 究 论 文 作 者 : 奚 荣 荣 指 导 老 师 : 姜 玉 琴 2014 年 12 月 硕 士 学 位 论 文 论 文 题 目 : 北 岛 诗 歌 创 作 的 双 重 困 境 专 业 名 称 : 中 国 现 当 代 文 学 研 究 方 向 : 中 国 新 诗 研 究 论 文 作 者 : 奚 荣 荣 指 导 老 师 : 姜 玉 琴 2014 年 12 月 致 谢 文 学 是 我 们 人 类 宝 贵 的 精 神 财 富 两 年 半 的 硕 士 学 习 让 我 进 一 步 接 近 文 学,

More information

Improved Preimage Attacks on AES-like Hash Functions: Applications to Whirlpool and Grøstl

Improved Preimage Attacks on AES-like Hash Functions: Applications to Whirlpool and Grøstl SKLOIS (Pseudo) Preimage Attack on Reduced-Round Grøstl Hash Function and Others Shuang Wu, Dengguo Feng, Wenling Wu, Jian Guo, Le Dong, Jian Zou March 20, 2012 Institute. of Software, Chinese Academy

More information

01-0982.doc

01-0982.doc 第 32 卷 第 7 期 岩 土 工 程 学 报 Vol.32 No.7 2010 年 7 月 Chinese Journal of Geotechnical Engineering July 2010 沿 海 碎 石 回 填 地 基 上 高 能 级 强 夯 系 列 试 验 对 比 研 究 年 廷 凯 1,2, 水 伟 厚 3, 李 鸿 江 4, 杨 庆 1,2, 王 玉 立 (1. 大 连 理 工

More information

2015年4月11日雅思阅读预测机经(新东方版)

2015年4月11日雅思阅读预测机经(新东方版) 剑 桥 雅 思 10 第 一 时 间 解 析 阅 读 部 分 1 剑 桥 雅 思 10 整 体 内 容 统 计 2 剑 桥 雅 思 10 话 题 类 型 从 以 上 统 计 可 以 看 出, 雅 思 阅 读 的 考 试 话 题 一 直 广 泛 多 样 而 题 型 则 稳 中 有 变 以 剑 桥 10 的 test 4 为 例 出 现 的 三 篇 文 章 分 别 是 自 然 类, 心 理 研 究 类,

More information

Shanghai International Studies University A STUDY ON SYNERGY BUYING PRACTICE IN ABC COMPANY A Thesis Submitted to the Graduate School and MBA Center I

Shanghai International Studies University A STUDY ON SYNERGY BUYING PRACTICE IN ABC COMPANY A Thesis Submitted to the Graduate School and MBA Center I 上 海 外 国 语 大 学 工 商 管 理 硕 士 学 位 论 文 ABC 中 国 食 品 公 司 的 整 合 采 购 研 究 学 科 专 业 : 工 商 管 理 硕 士 (MBA) 作 者 姓 名 :0113700719 指 导 教 师 : 答 辩 日 期 : 2013 年 12 月 上 海 外 国 语 大 学 二 一 四 年 一 月 Shanghai International Studies

More information

CNC KINWA FLAT BED CNC LATHE CL-58x3000 CL-38/CL-58 Series Swing 660, 900mm Between centers 1000~6000mm Spindle bore 120,186,258,375mm Spindle motor 3

CNC KINWA FLAT BED CNC LATHE CL-58x3000 CL-38/CL-58 Series Swing 660, 900mm Between centers 1000~6000mm Spindle bore 120,186,258,375mm Spindle motor 3 CNC KINWA FLAT BED CNC LATHE CL-8x000 CL-8/CL-8 Series Swing 660, 900mm Between centers 00~6000mm Spindle bore 1,186,8,7mm Spindle motor 0,40HP 14 CL-8x00,00,00,000,4000 CL-8x00,00,000,4000,000, C L8/C

More information

% % 34

% % 34 * 2000 2005 1% 1% 1% 1% * VZDA2010-15 33 2011. 3 2009 2009 2004 2008 1982 1990 2000 2005 1% 1 1 2005 1% 34 2000 2005 1% 35 2011. 3 2000 0. 95 20-30 209592 70982 33. 9% 2005 1% 258 20-30 372301 115483 31.

More information

天 主 教 輔 仁 大 學 社 會 學 系 學 士 論 文 小 別 勝 新 婚? 久 別 要 離 婚? 影 響 遠 距 家 庭 婚 姻 感 情 因 素 之 探 討 Separate marital relations are getting better or getting worse? -Exp

天 主 教 輔 仁 大 學 社 會 學 系 學 士 論 文 小 別 勝 新 婚? 久 別 要 離 婚? 影 響 遠 距 家 庭 婚 姻 感 情 因 素 之 探 討 Separate marital relations are getting better or getting worse? -Exp 天 主 教 輔 仁 大 學 社 會 學 系 學 士 論 文 小 別 勝 新 婚? 久 別 要 離 婚? 影 響 遠 距 家 庭 婚 姻 感 情 因 素 之 探 討 Separate marital relations are getting better or getting worse? -Explore the impact of emotional factors couples do not

More information

LH_Series_Rev2014.pdf

LH_Series_Rev2014.pdf REMINDERS Product information in this catalog is as of October 2013. All of the contents specified herein are subject to change without notice due to technical improvements, etc. Therefore, please check

More information

#4 ~ #5 12 m m m 1. 5 m # m mm m Z4 Z5

#4 ~ #5 12 m m m 1. 5 m # m mm m Z4 Z5 2011 6 6 153 JOURNAL OF RAILWAY ENGINEERING SOCIETY Jun 2011 NO. 6 Ser. 153 1006-2106 2011 06-0014 - 07 300142 ABAQUS 4. 287 mm 6. 651 mm U455. 43 A Analysis of Impact of Shield Tunneling on Displacement

More information