CN108229045A - A kind of high speed pantograph key parameter discrimination method based on sensitivity analysis - Google Patents
A kind of high speed pantograph key parameter discrimination method based on sensitivity analysis Download PDFInfo
- Publication number
- CN108229045A CN108229045A CN201810039385.0A CN201810039385A CN108229045A CN 108229045 A CN108229045 A CN 108229045A CN 201810039385 A CN201810039385 A CN 201810039385A CN 108229045 A CN108229045 A CN 108229045A
- Authority
- CN
- China
- Prior art keywords
- high speed
- speed pantograph
- pantograph
- sensitivity
- bow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Abstract
The invention discloses a kind of high speed pantograph key parameter discrimination methods based on sensitivity analysis, using Sobol ' Global sensitivity analysis methods, it derives and calculates one order value, Second Order Sensitivity value and total sensitivity value of the high speed pantograph structural parameters to bow head movement locus, and analyze different high speed pantograph structural parameters and the influence of displacement is lifted to high speed pantograph dynamic, so that it is determined that influencing the key structural parameters of high speed pantograph-catenary current collection quality；And then determine which kind of design parameter of adjustment is maximally efficient, reasonably select design parameter.The present invention can greatly simplify process of optimization compared with traditional high speed pantograph optimum design method, to improve the sliding contact characteristic of high speed bow net to greatest extent.
Description
Technical field
The present invention relates to high speed pantograph parameter identification technique field, it is specially a kind of based on the high speed of sensitivity analysis by
Pantograph key parameter discrimination method.
Background technology
High-speed railway pantograph contact line relation be influence railway security reliability service an important factor for one of, and contact net and at a high speed by
Pantograph Parameters Optimal Design is always the research emphasis of high speed pantograph contact line relation.To improve high speed pantograph-catenary current collection quality, part document
The prioritization scheme of centenary design parameter is proposed, such as Song Yang proposes a kind of contact net three-dimension modeling method, is considering
Non-linear solution under the conditions of windage yaw.However, these schemes need mostly on the implementation rebuild contact net, need consumption it is higher into
This.Compared with catenary's parameters prioritization scheme, high speed pantograph parameter design optimization scheme is more convenient for implementing, therefore with higher
Feasibility and realistic meaning.Such as Zhouning County uses multi-mass model, has studied single Parameters variation, analyzes high speed pantograph
Influence of the structural parameters to bow net contact power；Horse fruit build etc. using dynamics software Simpack establish at a high speed by
Pantograph-contact net Coupling Simulation Model has studied it to bow net dynamic characteristic using high speed pantograph mass block model parameter
It influences.But in existing high speed pantograph parameter optimization scheme, do not consider that high speed pantograph different designs are joined in analysis
Several influence degrees to bow net sliding contact characteristic.During Parameters Optimal Design, high speed pantograph univers parameter is carried out
Optimization design no doubt can be improved high speed Pantograph-OCS system current carrying quality, but such scheme optimization target is excessively general,
Easily increase cost in implementation process, and high speed pantograph performance degree of optimization can not be made to reach best.
Invention content
In view of the above-mentioned problems, the purpose of the present invention is to provide a kind of kinematic accuracies that can improve high speed pantograph；Together
When consider parameters influence and parameter between influence each other, make high speed pantograph Optimal Structure Designing accuracy and can
By the Research on Identification method of high speed pantograph key parameter of the property more preferably based on sensitivity analysis.Technical solution is as follows：
A kind of high speed pantograph key parameter discrimination method based on sensitivity analysis, which is characterized in that including following step
Suddenly：
Step 1：The random error of high speed pantograph pair clearance is included in, the bow head movement for establishing high speed pantograph is non-
Linear model derives bow head Movement Locus Equation according to the frame model of high speed pantograph；
Step 2：Using Latin hypercube samplings, the design variable random sample of high speed pantograph is simulated；
Step 3：Using Sobol ' Global sensitivity analysis methods, derive and calculate high speed pantograph design variable to bow head
One order value, Second Order Sensitivity value and the total sensitivity value of movement locus；
Step 4：High speed pantograph normalized wave function model is derived, calculates frame mass and frame under the influence of different designs variable
Frame damps；
Step 5：The different high speed pantograph frame mass of gained will be derived and frame damping substitutes into simulation model, analysis is different
High speed bow net contact pressure change under the influence of design variable verifies the sensitivity analysis knot of high speed pantograph key Design variable
Fruit.
Further, bow head balancing pole and horizontal direction angle in the bow head kinematic nonlinearities model of the high speed pantograph
β is：
Wherein, (Ex, Ey)、(Hx, Hy) it is respectively bow head balancing pole x, y-coordinate component；
The bow head movement locus is：
Wherein, Ex(i) it is the x coordinate component of bow head E locus of points curve discrete point, Ey(i) for bow head E locus of points curve from
Scatterplot y-coordinate component, β (i) are bow head counter-jib in i-th of position and the angle of horizontal direction, and n is is taken discrete point sum；
x1,x2,...,x11Represent the design variable of high speed pantograph, x1Length l for lower arm rod ACAC；x2For upper armed lever CD sections of length
lCD；x3For BG sections of length l of push rodBG；x4For BD sections of length l of push rodBD；x5For upper armed lever DE sections of length lDE；x6For upper ledge
The length l of frame control-rod GHGH；x7Length l for bow head balancing pole EHEH；x8For two two dot center's distances of fixed-hinged support A and B
lAB；x9For upper frame CD bars and the angle of DE bars；x10For BG bars on push rod and the angle of BD bars, x11For A the and B lines of centres and x
The angle of axis.
Further, the method that the design variable random sample to high speed pantograph is simulated is：Modulus is intended
Number is N, and the uniformly distributed function of all design variables is divided into the subinterval of N number of non-overlapping copies, in each subinterval respectively
Independent sampling with equal probability is carried out, each subinterval only generates a random number, then using inverse transformation method, by N number of sub-district
Between the random number that generates obtain N number of sample of random variable value, the serial number in the affiliated section of the sample value of each stochastic variable is carried out
Random alignment.
Further, the high speed pantograph design variable that calculates is to the one order value of bow head movement locus and total
The method of Sensitirity va1ue is：
Bow head E point path curves are decomposed into：
Ωn=(xi|ximin<xi<ximax, i=1,2 ..., n, n=11)
Wherein, Ey0It is a constant, the integration of any variable that other addend items include it is zero, i.e.,
It is orthogonal between knowing each addend item by formula (3) and formula (4), that is, if
(i1,i2,…is)≠(j1,j2,…jl)
So
Therefore
Wherein, dx represents dx1；…；dxn；
Decomposition in formula (3) is unique, and addend item is acquired by multiple integral：
Wherein, x~i, x~ijIt is represented respectively except xiAnd except xiAnd xjExcept variable；
Ey(x) population variance D is：
Partial variance is acquired by each addend item of formula (3)：
Wherein, 1≤i1<…is≤ n and s=1,2 ..., n；To formula (3) in entire ΩnDomain square, integration, obtain：
Then sensitivityIt is expressed as：
According to definition, obtained according to equation (11)：
Wherein, SiThe x of expressioniOne order；SijRepresent xiAnd xjSecond Order Sensitivity, i ≠ j；
The sensitivity of total system is the sum of each rank sensitivity coefficient of a variable, is expressed as
Wherein,
Wherein, Ey0、D、DiAnd D~iIt is acquired by monte carlo integration method：
According to equation (12), (16), (17) and (18), xiFirst-order sensitivity coefficient be：
According to equation (14), (16), (17) and (19), xiTotal sensitivity coefficient be：
Wherein, hits of the k for Monte Carlo Method, xmIn ΩnThe sampled point in domain；In formula (19) and formula (20) subscript (1) and
(2) two k × m dimension sampling arrays of X are represented；
Finally E is derived in aforementioned mannersx(i)=Ex(x1,x2,...,x11) sensitivity.
Further, the detailed process of the step 4 is：
According to lagrange equation of motion and the pull-type establishing equation frame part equation of motion：
Wherein, α represents rising bow angle,yeFor bow head E point ordinates, obtained by geometrical relationship：
ye=csin γ+x1sinα (24)
Wherein, c is CE pole lengths,
s2=x1 2+x8 2+2x1x8cos(α+x11)；
k10For according to intermediate variable obtained by virtual displacement variation principle, and k10=k4ccosγ+x1Cos α,
E, f is respectively B point ordinates and abscissa；
Taylor series expansion is carried out in bow head working depth, obtains equivalent linear of the pantograph frame in operating position attachment
Differential equation of motion：
Wherein, MFFor high speed pantograph frame mass, KH、KFRespectively head of high-speed pantograph and frame stiffness, Bh、BFPoint
Not Wei head of high-speed pantograph and frame frictional force, UH、UFThe respectively damping of head of high-speed pantograph and frame, ye、yhPoint
It Wei not bow head E points and H point ordinates；
Wherein, yhIt is obtained by geometrical relationship：
yh=e+x3sinξ+x6sinλ (26)
Wherein, ξ is GB bars and negative direction of the x-axis angle, and λ is HG bars and x-axis angle；
Pantograph inearized model reduction parameter is obtained as a result,：
Wherein, M be equivalent framework quality, MxFor rising bow torque.
The beneficial effects of the invention are as follows：The present invention is when carrying out high speed pantograph R ＆ D design, using the global spirits of Sobol '
Basis of sensitivity analysis method, to high speed pantograph Different structural parameters to one order value of the bow head track in x and y directions, second order
Sensitirity va1ue and total sensitivity value are calculated.Compared with traditional high speed pantograph optimum design method, this method is using at a high speed
The actual structure parameters model of pantograph can improve the kinematic accuracy of high speed pantograph；The influence of parameters is considered simultaneously
And influencing each other between parameter, make high speed pantograph Optimal Structure Designing accuracy and reliability more preferably, effectively reduce excellent
Change design cost.
Description of the drawings
Fig. 1 is high speed pantograph key structural parameters discrimination method flow chart.
Fig. 2 is high speed pantograph geometrical relationship schematic diagram.
Fig. 3 is the total sensitivity value (x directions) of high speed pantograph structural parameters.
Fig. 4 is the total sensitivity value (y directions) of high speed pantograph structural parameters.
Fig. 5 is high speed pantograph-contact net coupling model.
Fig. 6 is 350km/h bow net contact pressure change results under Different structural parameters.
Specific embodiment
The present invention is described in further details in the following with reference to the drawings and specific embodiments.The present embodiment is complete using Sobol '
Office's Sensitivity Analysis Method, to one order value of the high speed pantograph structural parameters to bow head movement locus, Second Order Sensitivity
Value and total sensitivity value carry out derivation calculating, and analyze different high speed pantograph structural parameters to high speed pantograph dynamic lifting position
The influence of shifting, so that it is determined that influencing the key structural parameters of high speed pantograph-catenary current collection quality.Method flow is as shown in Figure 1, in specific
Hold and design method step is as follows：
Step 1：Traditional quality block models are replaced, and introduce movement wherein using high speed pantograph nonlinear motion model
Auxiliary air gap error, high speed pantograph geometrical relationship schematic diagram are as shown in Figure 2.On this basis, according to each structure of high speed pantograph
Geometrical relationship between parameter carries out the derivation of head of high-speed pantograph Movement Locus Equation.Bow head balancing pole is pressed from both sides with horizontal direction
Shown in angular dependence such as formula (1), shown in bow head track such as formula (2).
Wherein, β is bow head balancing pole and horizontal direction angle；(Ex,Ey)、(Hx,Hy) it is respectively bow head balancing pole x, y is sat
Mark component.In formula (2), Ex(i) it is the x coordinate component of bow head E locus of points curve discrete point；Ey(i) it is bow head E locus of points curves
Discrete point y-coordinate component；β (i) is bow head counter-jib in i-th of position and the angle of horizontal direction；N is total to be taken discrete point
Number.x1,x2,...,x11Represent the design variable of high speed pantograph, wherein, x1Length l for lower arm rod ACAC；x2For upper armed lever CD
The length l of sectionCD；x3For BG sections of length l of push rodBG；x4For BD sections of length l of push rodBD；x5For upper armed lever DE sections of length lDE；
x6Length l for upper frame control-rod GHGH；x7Length l for bow head balancing pole EHEH；x8For two 2 points of fixed-hinged support A and B
Centre distance lAB；x9For upper frame CD bars and the angle of DE bars；x10For BG bars on push rod and the angle of BD bars, x11For in A and B
The angle of heart line and x-axis.
Step 2：Using the Latin hypercube method of samplings, for high speed pantograph x1,x2... ..., x1111 designs altogether
Variable carries out the sampling of random sample：
(1) it is 20 that modulus, which intends times N, and the normal distyribution function of all design variables is divided into the sub-district of 20 non-overlapping copies
Between, carry out independent sampling with equal probability respectively in each subinterval；
(2) in order to ensure that the random number extracted belongs to each subinterval, the random number V in i-th of sectioniIt should meet：
(3) each subinterval only generates a random number, then using inverse transformation method, by N number of subinterval generate with
Machine number obtains N number of sample of random variable value, and random alignment is carried out to the serial number in the affiliated section of the sample value of each stochastic variable.
Step 3：Using Sobol ' Global sensitivity analysis methods, derive and calculate high speed pantograph structural parameters to bow head
One order value, Second Order Sensitivity value and the total sensitivity value of movement locus：
It is by bow head E locus of points curve separatings：
Ωn=(xi|ximin<xi<ximax, i=1,2 ..., n, n=11)
Wherein, Ey0It is a constant, the integration of any variable that other addend items include it is zero；
It is orthogonal between each addend item can be obtained by formula (5) and formula (6), that is to say, that if
(i1,i2,…is)≠(j1,j2,…jl)
So
Therefore
Wherein, dx represents dx1；…；dxn。
Decomposition in formula (5) is unique, and addend item can be acquired by multiple integral
Wherein, x~i, x~ijIt is represented respectively except xiAnd except xiAnd xjExcept variable, higher order term can similarly be obtained.
Ey(x) population variance D is：
Partial variance can be acquired by each addend item of formula (5)：
Here 1≤i1<…is≤ n and s=1,2 ..., n.To formula (5) in entire ΩnDomain square, integration, obtain：
In this way, sensitivityIt is expressed as：
It can be obtained according to definition and equation (13)：
Wherein, SiThe x of expressioniOne order, Sij(i ≠ j) represents xiAnd xjSecond Order Sensitivity, and so on.Always
The sensitivity of system is the sum of each rank sensitivity coefficient of a variable, is expressed as：
Wherein,
Ey0、D、DiAnd D~iIt can be acquired by monte carlo integration method, it can be deduced that：
According to equation (14), (18), (19) and (20), xiFirst-order sensitivity coefficient be：
According to equation (16), (18), (19) and (21), xiTotal sensitivity coefficient be：
Wherein, hits of the k for Monte Carlo Method, xmIn ΩnThe sampled point in space.Subscript (1) in formula (21) and formula (22)
(2) two k × m dimension sampling arrays of X are represented.
E is derived in aforementioned mannersx(i)=Ex(x1,x2,...,x11) sensitivity.
Step 4：According to the sensitivity analysis of high speed pantograph structural parameters as a result, carrying out high speed pantograph normalized wave function model
Derivation, to the frame mass under the influence of Different structural parameters and frame damping calculate：
According to lagrange equation of motion and the pull-type establishing equation frame part equation of motion：
Wherein, α represents rising bow angle,yeFor bow head E point ordinates, can be obtained by geometrical relationship：
ye=csin γ+x1sinα (26)
Wherein, c is CE pole lengths,s2=x1 2+x8 2+
2x1x8cos(α+x11)。
k10For the intermediate variable according to obtained by virtual displacement variation principle：
k10=k4ccosγ+x1cosα (27)
Wherein,
E, f is respectively B point ordinates and abscissa.
Taylor series expansion is carried out in bow head working depth, obtains equivalent linear of the pantograph frame in operating position attachment
Differential equation of motion：
Wherein, MFFor high speed pantograph frame mass, KH、KFRespectively head of high-speed pantograph and frame stiffness, Bh、BFPoint
Not Wei head of high-speed pantograph and frame frictional force, UH、UFThe respectively damping of head of high-speed pantograph and frame, ye、yhPoint
It Wei not bow head E points and H point ordinates；
Wherein, yhIt can be obtained by geometrical relationship：
yh=e+x3sinξ+x6sinλ (29)
Wherein, ξ is GB bars and negative direction of the x-axis angle, and λ is HG bars and x-axis angle.
Pantograph inearized model reduction parameter is obtained as a result,：
Wherein, M be equivalent framework quality, MxFor rising bow torque.
Step 5：The different high speed pantograph frame mass of gained will be derived and frame damping substitutes into high speed pantograph-contact net
Coupling model, calculate Different structural parameters under the influence of bow net dynamic touch pressure, and to gained contact simulation result into
Row analysis, analysis result is compared with high speed pantograph structural parameters Calculation of Sensitivity result, to verify sensitivity analysis
As a result.
It in this example, can be obtained through Calculation of Sensitivity, high speed pantograph structural parameters are distinguished in x, the distribution situation in y directions
As shown in table 1, table 2, table 3, table 4, Fig. 3 and Fig. 4, strong affecting parameters (l can be divided intoAC lCD lBD), medium influence parameter (lDE
lEH) and weak affecting parameters (lAB lBG lGH)。
The one order assay value (x directions) of 1 pantograph parameters of table
The Second Order Sensitivity value (x directions) of 2 pantograph parameters of table
The one order assay value (y directions) of 3 pantograph parameters of table
The Second Order Sensitivity value (y directions) of 4 pantograph parameters of table
The analysis result of contact during 5 350km/h of table
It is damped later according to the high speed pantograph frame mass under the influence of step 4 and 5 calculating Different structural parameters and frame,
And it is substituted into high speed pantograph shown in Fig. 5-contact net coupling model and carries out subsequent simulation, simulation result such as table 5 and Fig. 6 institutes
Show.Simulation result, which is analyzed, to be obtained, and in this example, it is upper armed lever l to influence the bigger structural parameters of contactCDWith
lDE, contact variation is respectively 6.81%, 6.69%, is secondly balancing pole lEH(4.87%) and push rod lBD(4.63%), with
Calculation of Sensitivity result is basically identical.This shows using high speed pantograph structure obtained by Sobol ' Global sensitivity analysis methods
The sensitivity analysis result of parameter substantially conforms to practical dynamic analysis as a result, with tradition high speed pantograph parameter optimization scheme phase
Than that can realize the specific aim of high speed pantograph structure parameter optimizing using such scheme, high speed pantograph optimization can be effectively improved
Design efficiency.
Claims (5)
1. a kind of high speed pantograph key parameter discrimination method based on sensitivity analysis, which is characterized in that include the following steps：
Step 1：The random error of high speed pantograph pair clearance is included in, establishes the bow head kinematic nonlinearities of high speed pantograph
Model derives bow head Movement Locus Equation according to the frame model of high speed pantograph；
Step 2：Using Latin hypercube samplings, the design variable random sample of high speed pantograph is simulated；
Step 3：Using Sobol ' Global sensitivity analysis methods, derive and calculate high speed pantograph design variable to bow head movement
One order value, Second Order Sensitivity value and the total sensitivity value of track；
Step 4：Derive high speed pantograph normalized wave function model, calculate different designs variable under the influence of frame mass and frame resistance
Buddhist nun；
Step 5：The different high speed pantograph frame mass of gained will be derived and frame damping substitutes into simulation model, analyze different designs
High speed bow net contact pressure change under the influence of variable verifies the sensitivity analysis result of high speed pantograph key Design variable.
2. the high speed pantograph key parameter discrimination method according to claim 1 based on sensitivity analysis, feature exist
In bow head balancing pole is with horizontal direction angle β in the bow head kinematic nonlinearities model of the high speed pantograph：
Wherein, (Ex, Ey)、(Hx, Hy) it is respectively bow head balancing pole x, y-coordinate component；
The bow head movement locus is：
Wherein, Ex(i) it is the x coordinate component of bow head E locus of points curve discrete point, Ey(i) it is bow head E locus of points curve discrete point y
Coordinate components, β (i) are bow head counter-jib in i-th of position and the angle of horizontal direction, and n is is taken discrete point sum；x1,
x2,...,x11Represent the design variable of high speed pantograph, x1Length l for lower arm rod ACAC；x2For upper armed lever CD sections of length
lCD；x3For BG sections of length l of push rodBG；x4For BD sections of length l of push rodBD；
x5For upper armed lever DE sections of length lDE；x6Length l for upper frame control-rod GHGH；x7Length for bow head balancing pole EH
lEH；x8For two two dot center distance l of fixed-hinged support A and BAB；x9For upper frame CD bars and the angle of DE bars；x10For BG on push rod
The angle of bar and BD bars, x11For A the and B lines of centres and the angle of x-axis.
3. the high speed pantograph key parameter discrimination method according to claim 1 based on sensitivity analysis, feature exist
In the method that the design variable random sample to high speed pantograph is simulated is：Number realization is taken as N, all designs
The uniformly distributed function of variable is divided into the subinterval of N number of non-overlapping copies, carries out independent equiprobability respectively in each subinterval
Sampling, each subinterval only generate a random number, and then using inverse transformation method, the random number generated by N number of subinterval obtains
To N number of sample of random variable value, random alignment is carried out to the serial number in the affiliated section of the sample value of each stochastic variable.
4. the high speed pantograph key parameter discrimination method according to claim 2 based on sensitivity analysis, feature exist
In the calculating high speed pantograph design variable is to the one order value of bow head movement locus, Second Order Sensitivity value and total spirit
The method of sensitivity value is：
Bow head E point path curves are decomposed into：
Ωn=(xi|ximin<xi<ximax, i=1,2 ..., n, n=11)
Wherein, Ey0It is a constant, the integration of any variable that other addend items include it is zero, i.e.,
It is orthogonal between knowing each addend item by formula (3) and formula (4), that is, if
(i1,i2,…is)≠(j1,j2,…jl)
So
Therefore
Wherein, dx represents dx1；…；dxn；
Decomposition in formula (3) is unique, and addend item is acquired by multiple integral：
Wherein, x~i, x~ijIt is represented respectively except xiAnd except xiAnd xjExcept variable；
Ey(x) population variance D is：
Partial variance is acquired by each addend item of formula (3)：
Wherein, 1≤i1<…is≤ n and s=1,2 ..., n；To formula (3) in entire ΩnDomain square, integration, obtain：
Then sensitivityIt is expressed as：
According to definition, obtained according to equation (11)：
Wherein, SiThe x of expressioniOne order；SijRepresent xiAnd xjSecond Order Sensitivity, i ≠ j；
The sensitivity of total system is the sum of each rank sensitivity coefficient of a variable, is expressed as
Wherein,
Wherein, Ey0、D、DiAnd D~iIt is acquired by monte carlo integration method：
According to equation (12), (16), (17) and (18), xiFirst-order sensitivity coefficient be：
According to equation (14), (16), (17) and (19), xiTotal sensitivity coefficient be：
Wherein, hits of the k for Monte Carlo Method, xmIn ΩnThe sampled point in domain；Subscript (1) and (2) in formula (19) and formula (20)
Represent two k × m dimension sampling arrays of X；
Finally E is derived in aforementioned mannersx(i)=Ex(x1,x2,...,x11) sensitivity.
5. the high speed pantograph key parameter discrimination method according to claim 4 based on sensitivity analysis, feature exist
In the detailed process of the step 4 is：
According to lagrange equation of motion and the pull-type establishing equation frame part equation of motion：
Wherein, α represents rising bow angle,yeFor bow head E point ordinates, obtained by geometrical relationship：
ye=csin γ+x1Wherein, c is CE pole lengths to sin α (24),
s2=x1 2+x8 2+2x1x8cos(α+x11)；
k10For according to intermediate variable obtained by virtual displacement variation principle, and k10=k4ccosγ+x1Cos α,
E, f distinguishes
For B point ordinates and abscissa；
Taylor series expansion is carried out in bow head working depth, the equivalent linear for obtaining pantograph frame in operating position attachment moves
The differential equation：
Wherein, MFFor high speed pantograph frame mass, KH、KFRespectively head of high-speed pantograph and frame stiffness, Bh、
BFThe respectively frictional force of head of high-speed pantograph and frame, UH、UFThe respectively resistance of head of high-speed pantograph and frame
Buddhist nun, ye、yhRespectively bow head E points and H point ordinates；
Wherein, yhIt is obtained by geometrical relationship：
yh=e+x3sinξ+x6sinλ (26)
Wherein, ξ is GB bars and negative direction of the x-axis angle, and λ is HG bars and x-axis angle；
Pantograph inearized model reduction parameter is obtained as a result,：
Wherein, M be equivalent framework quality, MxFor rising bow torque.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810039385.0A CN108229045B (en) | 2018-01-16 | 2018-01-16 | High-speed pantograph key parameter identification method based on sensitivity analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810039385.0A CN108229045B (en) | 2018-01-16 | 2018-01-16 | High-speed pantograph key parameter identification method based on sensitivity analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108229045A true CN108229045A (en) | 2018-06-29 |
| CN108229045B CN108229045B (en) | 2021-05-04 |
Family
ID=62641508
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810039385.0A Expired - Fee Related CN108229045B (en) | 2018-01-16 | 2018-01-16 | High-speed pantograph key parameter identification method based on sensitivity analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108229045B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109214139A (en) * | 2018-11-16 | 2019-01-15 | 中车株洲电力机车有限公司 | A kind of pantograph Strength Simulation calculating fast modeling method |
| CN109766627A (en) * | 2019-01-08 | 2019-05-17 | 西南交通大学 | A kind of analysis method of the pantograph unsteady characteristic based on slide plate spacing |
| CN110059415A (en) * | 2019-04-22 | 2019-07-26 | 西南交通大学 | A kind of high speed pantograph multi-subject design method based on Cooperative Optimization Algorithm |
| CN110287562A (en) * | 2019-06-14 | 2019-09-27 | 西南交通大学 | A kind of high speed pantograph structure parameter optimizing method |
| CN110362786A (en) * | 2019-06-14 | 2019-10-22 | 西南交通大学 | A kind of decoupling modeling method of high-speed railway bow net coupled system dynamics response analysis |
| CN113656900A (en) * | 2021-07-06 | 2021-11-16 | 扬州大学 | Sensitivity correction-based method for analyzing influence factors of structural strength of harvester |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060064288A1 (en) * | 2004-09-22 | 2006-03-23 | Liem Ferryanto | System and method of interactive design of a product |
| US20130194268A1 (en) * | 2012-01-27 | 2013-08-01 | Nvidia Corporation | System and process for improved sampling for parallel light transport simulation |
| CN106202734A (en) * | 2016-07-14 | 2016-12-07 | 中国人民解放军国防科学技术大学 | global sensitivity analysis method based on Gaussian radial basis function |
| CN106503362A (en) * | 2016-10-26 | 2017-03-15 | 西南交通大学 | A kind of high speed pantograph Multi-objective Robust H ∞ controller design methods |
| CN106599340A (en) * | 2016-10-19 | 2017-04-26 | 北京交通大学 | Method for optimizing parameters of contact line and pantograph based on sensitivity analysis |
| CN106709264A (en) * | 2017-01-11 | 2017-05-24 | 合肥工业大学 | Double-crane system variable amplitude angle response modeling algorithm and random response domain prediction method |
| CN107315859A (en) * | 2017-06-07 | 2017-11-03 | 西南交通大学 | A kind of high ferro car pessimistic concurrency control method for building up of the secondary arcing influence of meter and bow net |
-
2018
- 2018-01-16 CN CN201810039385.0A patent/CN108229045B/en not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060064288A1 (en) * | 2004-09-22 | 2006-03-23 | Liem Ferryanto | System and method of interactive design of a product |
| US20130194268A1 (en) * | 2012-01-27 | 2013-08-01 | Nvidia Corporation | System and process for improved sampling for parallel light transport simulation |
| CN106202734A (en) * | 2016-07-14 | 2016-12-07 | 中国人民解放军国防科学技术大学 | global sensitivity analysis method based on Gaussian radial basis function |
| CN106599340A (en) * | 2016-10-19 | 2017-04-26 | 北京交通大学 | Method for optimizing parameters of contact line and pantograph based on sensitivity analysis |
| CN106503362A (en) * | 2016-10-26 | 2017-03-15 | 西南交通大学 | A kind of high speed pantograph Multi-objective Robust H ∞ controller design methods |
| CN106709264A (en) * | 2017-01-11 | 2017-05-24 | 合肥工业大学 | Double-crane system variable amplitude angle response modeling algorithm and random response domain prediction method |
| CN107315859A (en) * | 2017-06-07 | 2017-11-03 | 西南交通大学 | A kind of high ferro car pessimistic concurrency control method for building up of the secondary arcing influence of meter and bow net |
Non-Patent Citations (4)
| Title |
|---|
| O. VO VAN .ETAL: "Sensitivity analysis of catenary geometry on current collection quality", 《THIRD INTERNATIONAL CONFERENCE ON RAILWAY TECHNOLOGY: RESEARCH, DEVELOPMENT AND MAINTENANCE》 * |
| 张剑: "基于代理模型技术的高速列车性能参数设计及优化", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
| 郭汝坤，冯春，周东，李世海: "基于破裂度可靠性的边坡灾变状态评价方法研究", 《岩石力学与工程学报》 * |
| 齐伟 等: "Sobol"方法分析TOPMODEL水文模型参数敏感性", 《水文》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109214139A (en) * | 2018-11-16 | 2019-01-15 | 中车株洲电力机车有限公司 | A kind of pantograph Strength Simulation calculating fast modeling method |
| CN109766627A (en) * | 2019-01-08 | 2019-05-17 | 西南交通大学 | A kind of analysis method of the pantograph unsteady characteristic based on slide plate spacing |
| CN110059415A (en) * | 2019-04-22 | 2019-07-26 | 西南交通大学 | A kind of high speed pantograph multi-subject design method based on Cooperative Optimization Algorithm |
| CN110059415B (en) * | 2019-04-22 | 2022-09-02 | 西南交通大学 | High-speed pantograph multidisciplinary design method based on collaborative optimization algorithm |
| CN110287562A (en) * | 2019-06-14 | 2019-09-27 | 西南交通大学 | A kind of high speed pantograph structure parameter optimizing method |
| CN110362786A (en) * | 2019-06-14 | 2019-10-22 | 西南交通大学 | A kind of decoupling modeling method of high-speed railway bow net coupled system dynamics response analysis |
| CN110287562B (en) * | 2019-06-14 | 2022-07-01 | 西南交通大学 | High-speed pantograph structure parameter optimization method |
| CN110362786B (en) * | 2019-06-14 | 2023-03-14 | 西南交通大学 | Decoupling modeling method for dynamic response analysis of high-speed railway bow-net coupling system |
| CN113656900A (en) * | 2021-07-06 | 2021-11-16 | 扬州大学 | Sensitivity correction-based method for analyzing influence factors of structural strength of harvester |
| CN113656900B (en) * | 2021-07-06 | 2023-06-20 | 扬州大学 | Harvester structural strength influence factor analysis method based on sensitivity correction |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108229045B (en) | 2021-05-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108229045A (en) | A kind of high speed pantograph key parameter discrimination method based on sensitivity analysis | |
| Cao et al. | Research on dynamic nonlinear input prediction of fault diagnosis based on fractional differential operator equation in high-speed train control system | |
| CN111475596B (en) | Sub-segment similarity matching method based on multi-level track coding tree | |
| CN104807039B (en) | Variable dimensionality reduction modeling method for boiler combustion optimization | |
| CN103995937B (en) | Precision machine tool mass matching design method based on response surface and genetic algorithm | |
| CN103714166B (en) | laser radar power line point cloud data fitting method | |
| CN106779384B (en) | Iron and steel industry blast furnace gas long-term interval prediction method based on optimal information granularity distribution | |
| CN104679956A (en) | High-speed press base reliability robust design method considering dynamic characteristics | |
| CN105844015B (en) | A kind of automobile glass lifter preferred arrangement design method | |
| Cui et al. | Shape optimization of high-speed train with the speed of 500kph | |
| CN106951657A (en) | One kind abrasion steel rail grinding target profile Fast design method | |
| CN111980688B (en) | Integrated learning algorithm-based inclination angle prediction method | |
| Yu | Understanding regional development mechanisms in Greater Beijing Area, China, 1995–2001, from a spatial–temporal perspective | |
| CN109918793A (en) | A kind of bogie suspension parameter fast Optimization | |
| CN114509991B (en) | Numerical control machine tool cutting stability prediction and optimization method considering parameter uncertainty | |
| Etesami et al. | Optimal transmission angle and dynamic balancing of slider-crank mechanism with joint clearance using Pareto Bi-objective Genetic Algorithm | |
| CN109255141A (en) | A kind of body of a motor car forward direction conceptual design cross sectional shape optimization method | |
| CN105512427B (en) | Consider that bow net contacts the high-speed railway bow net dynamic-simulation method of momentum impact again | |
| CN111325441B (en) | Quantitative evaluation method for shale gas target preservation conditions | |
| CN103810324B (en) | A kind of build-up tolerance Optimization Design optimized based on cost objective | |
| Shi et al. | Efficient method for calculating the composite stiffness of parabolic leaf springs with variable stiffness for vehicle rear suspension | |
| CN106201990A (en) | A kind of Mathematical Modeling Methods utilizing computer hyperspace | |
| CN115374536A (en) | Method for constructing four-degree-of-freedom numerical model of cab suspension system | |
| Cheng et al. | HIL simulation of high-speed train active stability employing active inertial actuators | |
| CN108919008B (en) | Online power quality disturbance identification method and system based on time-frequency database |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210504 Termination date: 20220116 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |