CN108763674A - A kind of method for solving of elastic boundary condition downhaul beam frequency - Google Patents
A kind of method for solving of elastic boundary condition downhaul beam frequency Download PDFInfo
- Publication number
- CN108763674A CN108763674A CN201810465820.6A CN201810465820A CN108763674A CN 108763674 A CN108763674 A CN 108763674A CN 201810465820 A CN201810465820 A CN 201810465820A CN 108763674 A CN108763674 A CN 108763674A
- Authority
- CN
- China
- Prior art keywords
- indicate
- drag
- line
- formula
- value
- 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/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Vibration Prevention Devices (AREA)
Abstract
The invention discloses a kind of method for solving of elastic boundary condition downhaul beam frequency, the performance parameter of drag-line to be detected is obtained first and inhaul cable vibration partial differential equation are built based on the performance parameter of the drag-line got, then the separation of variable is used to carry out abbreviation to inhaul cable vibration partial differential equation, obtain the inhaul cable vibration differential equation, establish the elastic boundary condition of inhaul cable vibration, the model function of vibration of drag-line is unfolded in the form of Chebyshev series, then the model function of vibration based on drag-line and elastic boundary condition obtain the characteristic equation of inhaul cable vibration, the characteristic equation for being finally based on inhaul cable vibration obtains the frequency equation of drag-line, it solves the frequency equation and obtains the vibration frequency of drag-line;Advantage is that accuracy is higher.
Description
Technical field
The present invention relates to a kind of method for solving of drag-line beam frequency, more particularly, under a kind of elastic boundary condition
The method for solving of drag-line beam frequency.
Background technology
Currently, in order to ensure the safe handling of bridge, it usually needs be monitored to the indices of bridge.And drag-line
Important monitoring index one of of the bending vibration as science of bridge building usually assumes that the boundary condition at drag-line both ends is ideal clamped
Or freely-supported condition is calculated the vibration frequency of drag-line.However, in fact, the boundary condition at drag-line both ends and nonideal letter
Branch and clamped, but between freely-supported and it is clamped between, such as the end of cable-stayed bridge cable is by anchorage and girder and Sarasota phase
Connection, anchorage mode is complex in engineering, and the boundary condition at drag-line both ends is assumed to clamped or simple boundary, this
It is assumed that being constrained there are larger error with actual boundary, there are larger mistakes for the drag-line beam frequency thus caused
Difference, accuracy be not high.
Invention content
Technical problem to be solved by the invention is to provide a kind of higher elastic boundary condition downhaul bendings of accuracy
The method for solving of vibration frequency.
Technical solution is used by the present invention solves above-mentioned technical problem:A kind of bending of elastic boundary condition downhaul is shaken
The method for solving of dynamic frequency, includes the following steps:
(1) performance parameter of drag-line to be detected and inclined based on the performance parameter of the drag-line got structure inhaul cable vibration is obtained
The differential equation:
In formula (1), x indicates that the vibration displacement of drag-line to be detected, x ∈ [0, L], L indicate the length of drag-line to be detected, T tables
Show that the design value of the Suo Li of drag-line to be detected, E indicate that the elasticity modulus of drag-line to be detected, I indicate that the section of drag-line to be detected is used
Property square, ρ indicate that the density of drag-line to be detected, A indicate that the sectional area of drag-line to be detected, t indicate the time of vibration of drag-line to be detected,
U (x, t) indicates each point displacement function that t changes at any time of drag-line to be detected,Indicate that u (x, t) asks two to time t
Rank partial derivative;Indicate that u (x, t) seeks second-order partial differential coefficient to x,Indicate that u (x, t) seeks quadravalence partial derivative to x;
(2) separation of variable is used to carry out abbreviation, specific abbreviation mistake to the inhaul cable vibration partial differential equation that step (1) is built
Cheng Wei:
A. u (x, t) is decomposed into the multiplication relationship between the function of time and spatial function, obtained:
Y (x) is the model function of vibration of drag-line in formula (2), and sin (*) is SIN function, and ω is the vibration frequency of drag-line,To draw
The phase angle of Suo Zhendong;
B. formula (2) is substituted into formula (1) and obtains the inhaul cable vibration differential equation:
EIY(4)(x)-TY(2)(x)-ρAω2Y (x)=0 (3)
Wherein, Y(4)(x) indicate that Y (x) asks Fourth-Derivative, Y to x(2)(x) indicate that Y (x) seeks second dervative to x;
(3) elastic boundary condition of inhaul cable vibration is established:
C. it is K to assume that the left end of drag-line is connected separately with elastic stiffness1Extension spring and elastic stiffness be K2Torsion bullet
Spring, it is K that the right end of drag-line, which is connected separately with elastic stiffness,3Extension spring and elastic stiffness be K4Torsionspring, thus
It arrives:
In formula (4), u " ' |X=0Indicate that u (x, t) seeks three rank partial derivatives to x at x=0, u |X=0Indicate at x=0 u (x,
T) functional value, u " |X=0Indicate that u (x, t) asks second-order partial differential coefficient, u'| to x at x=0X=0Indicate that u (x, t) is to x at x=0
First-order partial derivative is sought, u " ' |X=LIndicate that u (x, t) seeks three rank partial derivatives to x at x=L, u |X=LIndicate u (x, t) at x=L
Functional value, u " |X=LIndicate that u (x, t) asks second-order partial differential coefficient, u'| to x at x=LX=LIndicate that u (x, t) seeks x at x=L
First-order partial derivative, k1It is K for elastic stiffness1Extension spring rigidity nondimensionalization coefficient, k2It is K for elastic stiffness2Torsion bullet
Spring rigidity nondimensionalization coefficient, k3It is K for elastic stiffness3Extension spring rigidity nondimensionalization coefficient, k4It is K for elastic stiffness4
Torsionspring rigidity nondimensionalization coefficient, k1, k2, k3And k4Value be respectively:
In formula (5), symbol "/" is division operation symbol;k1, k2, k3And k4Value range be [0,104], and it is full simultaneously
It is enough lower three conditions:One, k1, k2, k3And k4Value cannot simultaneously be 0;Two, k1, k2, k3And k4Value cannot be simultaneously
104;Three, work as k1And k3Value simultaneously be 104, and k2Value be 0 when, k4Value cannot be 104, work as k1And k3Value
It is 10 simultaneously4, k4Value be 0 when, k2Value cannot be 104;
D. it is by the elastic boundary condition for obtaining drag-line in (2) formula substitution (4) formula:
In formula (6):Y"'|X=0Indicate that Y (x) seeks three order derivatives to x at x=0, Y |X=0Indicate the letter of the Y (x) at x=0
Numerical value, Y " |X=0Indicate that Y (x) asks second dervative, Y'| to x at x=0X=0Indicate that Y (x) seeks first derivative to x at x=0,
Y”'|X=LIndicate that Y (x) seeks three order derivatives to x at x=L, Y |X=LIndicate the functional value of the Y (x) at x=L, Y " |X=LIt indicates
Y (x) asks second dervative, Y'| to x at x=LX=LIndicate that Y (x) seeks first derivative to x at x=L;
(4) the model function of vibration Y (x) of drag-line is unfolded in the form of Chebyshev series, obtains the expansion of Y (x)
Formula:
In formula (7):∑ (*) expression sums to *, and n is the most high-order of Chebyshev series, and n is whole more than or equal to 1
Number, i=1,2 ..., n;Φi(x) it is to move letter using the i-th component level for meeting boundary condition formula (6) that Chebyshev series construct
Number, aiFunction phi is moved for the i-th component leveli(x) coefficient is undetermined coefficient, wherein
Φi(x)=h (x) Ti(x)+g(x) (8)
Wherein, Ti(x) it is the i-th rank Chebyshev series, h (x) and g (x) are the auxiliary function of construction, Ti(x) expression formula
For:
In formula (9):Cos (*) indicates that cosine function, arccos (*) indicate inverse cosine function;
The expression formula of h (x) and g (x) indicates as follows respectively:
(5) according to the principle of virtual displacement enable drag-line on virtual displacement δ Y (x) made by work(be zero, wherein δ indicate variation, δ Y
(x) expression formula is:
Drag-line on virtual displacement δ Y (x) made by work(:
According to formula (12) and formula (13), can obtain:
Wherein Φi (4)(x) Φ is indicatedi(x) Fourth-Derivative, Φ are asked to xi (2)(x) Φ is indicatedi(x) second dervative, Φ are asked to xj
(x) jth rank displacement function, j=1,2 ..., n are indicated;Arrangement formula can obtain equation after (14):
In formula (15):
Write formula (15) as matrix form, the as characteristic equation of inhaul cable vibration:
(D-ω2M) (16) H=0
In formula (16), D indicates that stiffness matrix, M indicate that mass matrix, H indicate coefficient matrix, dijThe i-th row of representing matrix D
Jth column element, mijThe i-th row jth column element of representing matrix M, i.e.,:
(6) characteristic equation (16) of inhaul cable vibration is solved, obtaining frequency equation is:
|D-ω2M |=0 (17)
Wherein | * | the determinant of representing matrix * is solved to obtain under different boundary constraints to frequency equation (17)
The vibration frequency ω of drag-line.
Compared with the prior art, the advantages of the present invention are as follows the performance parameters and base first by obtaining drag-line to be detected
Inhaul cable vibration partial differential equation are built in the performance parameter of the drag-line got, then use the separation of variable inclined to inhaul cable vibration
The differential equation carries out abbreviation, obtains the inhaul cable vibration differential equation, establishes the elastic boundary condition of inhaul cable vibration, the elastic boundary item
The boundary condition that drag-line left and right ends are eliminated in part is clamped (k1, k2, k3And k4Value simultaneously be 104) the case where, drawing
Rope left and right ends are in free state (k1, k2, k3And k4Value simultaneously for 0) the case where, drag-line left and right ends boundary condition
It is that freely-supported (works as k1And k3Value simultaneously be 104When, k2And k4Value simultaneously for 0) the case where, drag-line left and right ends side
Boundary's condition is that clamped-freely-supported (works as k1And k3Value simultaneously be 104, k2Value be 104, k4Value be 0) the case where and
The boundary condition of drag-line left and right ends, which is that freely-supported-is clamped, (works as k1And k3Value simultaneously be 104, k2Value be 0, k4Value is
104) the case where, make the left and right ends boundary condition of drag-line be:Freely-supported-elastic boundary, elasticity-simple boundary, clamped-elastic edge
Then the model function of vibration of drag-line is used Chebyshev series by one kind in boundary, elasticity-built-in boundary and inelastic-elastic boundary
Form be unfolded, then the model function of vibration based on drag-line and elastic boundary condition obtain the characteristic equation of inhaul cable vibration, most
The characteristic equation based on inhaul cable vibration obtains the frequency equation of drag-line afterwards, solves the frequency equation and obtains the vibration frequency of drag-line,
In the method for the present invention, the frequency equation of drag-line is the model of elastic boundary constraints downhaul bending vibration, is considered simultaneously
Bending stiffness and different boundary constrain the influence to drag-line beam frequency, and use Chebyshev polynomial sides
Method (i.e. Chebyshev series) approaches drag-line shaping function, further increases the accuracy and convergence of measurement result,
Get the higher inhaul cable vibration frequency of accuracy.
Description of the drawings
Fig. 1 is the realization structure for the elastic boundary condition that the present invention is built.
Specific implementation mode
Below in conjunction with attached drawing embodiment, present invention is further described in detail.
Embodiment:A kind of method for solving of elastic boundary condition downhaul beam frequency, includes the following steps:
(1) performance parameter of drag-line to be detected and inclined based on the performance parameter of the drag-line got structure inhaul cable vibration is obtained
The differential equation:
In formula (1), x indicates that the vibration displacement of drag-line to be detected, x ∈ [0, L], L indicate the length of drag-line to be detected, T tables
Show that the design value of the Suo Li of drag-line to be detected, E indicate that the elasticity modulus of drag-line to be detected, I indicate that the section of drag-line to be detected is used
Property square, ρ indicate that the density of drag-line to be detected, A indicate that the sectional area of drag-line to be detected, t indicate the time of vibration of drag-line to be detected,
U (x, t) indicates each point displacement function that t changes at any time of drag-line to be detected,Indicate that u (x, t) asks two to time t
Rank partial derivative;Indicate that u (x, t) seeks second-order partial differential coefficient to x,Indicate that u (x, t) seeks quadravalence partial derivative to x;
(2) separation of variable is used to carry out abbreviation, specific abbreviation mistake to the inhaul cable vibration partial differential equation that step (1) is built
Cheng Wei:
A. u (x, t) is decomposed into the multiplication relationship between the function of time and spatial function, obtained:
Y (x) is the model function of vibration of drag-line in formula (2), and sin (*) is SIN function, and ω is the vibration frequency of drag-line,To draw
The phase angle of Suo Zhendong;
B. formula (2) is substituted into formula (1) and obtains the inhaul cable vibration differential equation:
EIY(4)(x)-TY(2)(x)-ρAω2Y (x)=0 (3)
Wherein, Y(4)(x) indicate that Y (x) asks Fourth-Derivative, Y to x(2)(x) indicate that Y (x) seeks second dervative to x;
(3) elastic boundary condition of inhaul cable vibration is established, as shown in Figure 1:
C. it is K to assume that the left end of drag-line is connected separately with elastic stiffness1Extension spring and elastic stiffness be K2Torsion bullet
Spring, it is K that the right end of drag-line, which is connected separately with elastic stiffness,3Extension spring and elastic stiffness be K4Torsionspring, thus
It arrives:
In formula (4), u " ' |X=0Indicate that u (x, t) seeks three rank partial derivatives to x at x=0, u |X=0Indicate at x=0 u (x,
T) functional value, u " |X=0Indicate that u (x, t) asks second-order partial differential coefficient, u'| to x at x=0X=0Indicate that u (x, t) is to x at x=0
First-order partial derivative is sought, u " ' |X=LIndicate that u (x, t) seeks three rank partial derivatives to x at x=L, u |X=LIndicate u (x, t) at x=L
Functional value, u " |X=LIndicate that u (x, t) asks second-order partial differential coefficient, u'| to x at x=LX=LIndicate that u (x, t) seeks x at x=L
First-order partial derivative, k1It is K for elastic stiffness1Extension spring rigidity nondimensionalization coefficient, k2It is K for elastic stiffness2Torsion bullet
Spring rigidity nondimensionalization coefficient, k3It is K for elastic stiffness3Extension spring rigidity nondimensionalization coefficient, k4It is K for elastic stiffness4
Torsionspring rigidity nondimensionalization coefficient, k1, k2, k3And k4Value be respectively:
In formula (5), symbol "/" is division operation symbol;k1, k2, k3And k4Value range be [0,104], and it is full simultaneously
It is enough lower three conditions:One, k1, k2, k3And k4Value cannot simultaneously be 0;Two, k1, k2, k3And k4Value cannot be simultaneously
104;Three, work as k1And k3Value simultaneously be 104, and k2Value be 0 when, k4Value cannot be 104, work as k1And k3Value
It is 10 simultaneously4, k4Value be 0 when, k2Value cannot be 104;D. (2) formula is substituted into (4) formula and obtains the elastic edge of drag-line
Boundary's condition is:
In formula (6):Y"'|X=0Indicate that Y (x) seeks three order derivatives to x at x=0, Y |X=0Indicate the letter of the Y (x) at x=0
Numerical value, Y " |X=0Indicate that Y (x) asks second dervative, Y'| to x at x=0X=0Indicate that Y (x) seeks first derivative to x at x=0,
Y”'|X=LIndicate that Y (x) seeks three order derivatives to x at x=L, Y |X=LIndicate the functional value of the Y (x) at x=L, Y " |X=LIt indicates
Y (x) asks second dervative, Y'| to x at x=LX=LIndicate that Y (x) seeks first derivative to x at x=L;
(4) the model function of vibration Y (x) of drag-line is unfolded in the form of Chebyshev series, obtains the expansion of Y (x)
Formula:
In formula (7):∑ (*) expression sums to *, and n is the most high-order of Chebyshev series, and n is whole more than or equal to 1
Number, i=1,2 ..., n;Φi(x) it is to move letter using the i-th component level for meeting boundary condition formula (6) that Chebyshev series construct
Number, aiFunction phi is moved for the i-th component leveli(x) coefficient is undetermined coefficient, wherein
Φi(x)=h (x) Ti(x)+g(x) (8)
Wherein, Ti(x) it is the i-th rank Chebyshev series, h (x) and g (x) are the auxiliary function of construction, Ti(x) expression formula
For:
In formula (9):Cos (*) indicates that cosine function, arccos (*) indicate inverse cosine function;
The expression formula of h (x) and g (x) indicates as follows respectively:
(5) according to the principle of virtual displacement enable drag-line on virtual displacement δ Y (x) made by work(be zero, wherein δ indicate variation, δ Y
(x) expression formula is:
Drag-line on virtual displacement δ Y (x) made by work(:
According to formula (12) and formula (13), can obtain:
Wherein Φi (4)(x) Φ is indicatedi(x) Fourth-Derivative, Φ are asked to xi (2)(x) Φ is indicatedi(x) second dervative, Φ are asked to xj
(x) jth rank displacement function, j=1,2 ..., n are indicated;Arrangement formula can obtain equation after (14):
In formula (15):
Write formula (15) as matrix form, the as characteristic equation of inhaul cable vibration:
(D-ω2M) (16) H=0
In formula (16), D indicates that stiffness matrix, M indicate that mass matrix, H indicate coefficient matrix, dijThe i-th row of representing matrix D
Jth column element, mijThe i-th row jth column element of representing matrix M, i.e.,:
(6) characteristic equation (16) of inhaul cable vibration is solved, obtaining frequency equation is:
|D-ω2M |=0 (17)
Wherein | * | the determinant of representing matrix * is solved to obtain under different boundary constraints to frequency equation (17)
The vibration frequency ω of drag-line.
Claims (1)
1. a kind of method for solving of elastic boundary condition downhaul beam frequency, it is characterised in that include the following steps:
(1) it obtains the performance parameter of drag-line to be detected and inhaul cable vibration partial differential is built based on the performance parameter of the drag-line got
Equation:
In formula (1), x indicates that the vibration displacement of drag-line to be detected, x ∈ [0, L], L indicate that the length of drag-line to be detected, T expressions wait for
The design value of the Suo Li of drag-line is detected, E indicates that the elasticity modulus of drag-line to be detected, I indicate the cross sectional moment of inertia of drag-line to be detected,
ρ indicates the density of drag-line to be detected, and A indicates the sectional area of drag-line to be detected, and t indicates the time of vibration of drag-line to be detected, u (x,
T) each point displacement function that t changes at any time of drag-line to be detected is indicated,Indicate that u (x, t) seeks Second Order Partial to time t
Derivative;Indicate that u (x, t) seeks second-order partial differential coefficient to x,Indicate that u (x, t) seeks quadravalence partial derivative to x;
(2) separation of variable is used to carry out abbreviation, specific abbreviation process to the inhaul cable vibration partial differential equation that step (1) is built
For:
A. u (x, t) is decomposed into the multiplication relationship between the function of time and spatial function, obtained:
Y (x) is the model function of vibration of drag-line in formula (2), and sin (*) is SIN function, and ω is the vibration frequency of drag-line,It shakes for drag-line
Dynamic phase angle;
B. formula (2) is substituted into formula (1) and obtains the inhaul cable vibration differential equation:
EIY(4)(x)-TY(2)(x)-ρAω2Y (x)=0 (3)
Wherein, Y(4)(x) indicate that Y (x) asks Fourth-Derivative, Y to x(2)(x) indicate that Y (x) seeks second dervative to x;
(3) elastic boundary condition of inhaul cable vibration is established:
C. it is K to assume that the left end of drag-line is connected separately with elastic stiffness1Extension spring and elastic stiffness be K2Torsionspring,
It is K that the right end of drag-line, which is connected separately with elastic stiffness,3Extension spring and elastic stiffness be K4Torsionspring, thus obtain:
In formula (4), u " ' |X=0Indicate that u (x, t) seeks three rank partial derivatives to x at x=0, u |X=0Indicate the u (x, t) at x=0
Functional value, u " |X=0Indicate that u (x, t) asks second-order partial differential coefficient, u'| to x at x=0X=0Indicate that u (x, t) asks one to x at x=0
Rank partial derivative, u " ' |X=LIndicate that u (x, t) seeks three rank partial derivatives to x at x=L, u |X=LIndicate the letter of the u (x, t) at x=L
Numerical value, u " |X=LIndicate that u (x, t) asks second-order partial differential coefficient, u'| to x at x=LX=LIndicate that u (x, t) seeks single order to x at x=L
Partial derivative, k1It is K for elastic stiffness1Extension spring rigidity nondimensionalization coefficient, k2It is K for elastic stiffness2Torsionspring it is rigid
Spend nondimensionalization coefficient, k3It is K for elastic stiffness3Extension spring rigidity nondimensionalization coefficient, k4It is K for elastic stiffness4Torsion
Turn spring rate nondimensionalization coefficient, k1, k2, k3And k4Value be respectively:
In formula (5), symbol "/" is division operation symbol;k1, k2, k3And k4Value range be [0,104], and meet simultaneously with
Lower three conditions:One, k1, k2, k3And k4Value cannot simultaneously be 0;Two, k1, k2, k3And k4Value cannot simultaneously be 104;
Three, work as k1And k3Value simultaneously be 104, and k2Value be 0 when, k4Value cannot be 104, work as k1And k3Value simultaneously
It is 104, k4Value be 0 when, k2Value cannot be 104;D. (2) formula is substituted into (4) formula and obtains the elastic boundary item of drag-line
Part is:
In formula (6):Y"'|X=0Indicate that Y (x) seeks three order derivatives to x at x=0, Y |X=0Indicate the functional value of the Y (x) at x=0,
Y”|X=0Indicate that Y (x) asks second dervative, Y'| to x at x=0X=0Indicate that Y (x) seeks first derivative to x at x=0, Y " ' |X=L
Indicate that Y (x) seeks three order derivatives to x at x=L, Y |X=LIndicate the functional value of the Y (x) at x=L, Y " |X=LIt indicates at x=L
Y (x) asks second dervative, Y'| to xX=LIndicate that Y (x) seeks first derivative to x at x=L;
(4) the model function of vibration Y (x) of drag-line is unfolded in the form of Chebyshev series, obtains the expansion of Y (x):
In formula (7):∑ (*) expression sums to *, and n is the most high-order of Chebyshev series, and n is the integer more than or equal to 1, i
=1,2 ..., n;Φi(x) it is the i-th rank displacement function for meeting boundary condition formula (6) constructed using Chebyshev series, aiFor
I-th component level moves function phii(x) coefficient is undetermined coefficient, wherein
Φi(x)=h (x) Ti(x)+g(x) (8)
Wherein, Ti(x) it is the i-th rank Chebyshev series, h (x) and g (x) are the auxiliary function of construction, Ti(x) expression formula is:
In formula (9):Cos (*) indicates that cosine function, arccos (*) indicate inverse cosine function;
The expression formula of h (x) and g (x) indicates as follows respectively:
(5) according to the principle of virtual displacement enable drag-line on virtual displacement δ Y (x) made by work(be zero, wherein δ indicate variation, δ Y's (x)
Expression formula is:
Drag-line on virtual displacement δ Y (x) made by work(:
According to formula (12) and formula (13), can obtain:
Wherein Φi (4)(x) Φ is indicatedi(x) Fourth-Derivative, Φ are asked to xi (2)(x) Φ is indicatedi(x) second dervative, Φ are asked to xj(x)
Indicate jth rank displacement function, j=1,2 ..., n;Arrangement formula can obtain equation after (14):
In formula (15):
Write formula (15) as matrix form, the as characteristic equation of inhaul cable vibration:
(D-ω2M) (16) H=0
In formula (16), D indicates that stiffness matrix, M indicate that mass matrix, H indicate coefficient matrix, dijThe i-th row jth of representing matrix D
Column element, mijThe i-th row jth column element of representing matrix M, i.e.,:
(6) characteristic equation (16) of inhaul cable vibration is solved, obtaining frequency equation is:
|D-ω2M |=0 (17)
Wherein | * | the determinant of representing matrix * solves frequency equation (17) to obtain different boundary constraints downhaul
Vibration frequency ω.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810465820.6A CN108763674B (en) | 2018-05-16 | 2018-05-16 | Method for solving bending vibration frequency of inhaul cable under elastic boundary condition |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810465820.6A CN108763674B (en) | 2018-05-16 | 2018-05-16 | Method for solving bending vibration frequency of inhaul cable under elastic boundary condition |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108763674A true CN108763674A (en) | 2018-11-06 |
CN108763674B CN108763674B (en) | 2021-12-17 |
Family
ID=64008003
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810465820.6A Active CN108763674B (en) | 2018-05-16 | 2018-05-16 | Method for solving bending vibration frequency of inhaul cable under elastic boundary condition |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108763674B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109829211A (en) * | 2019-01-21 | 2019-05-31 | 东南大学 | A kind of thermal environment lower plate structure high frequency partial method of response calculation |
CN109885898A (en) * | 2019-01-28 | 2019-06-14 | 华北水利水电大学 | The measuring method of the eigentone of non-linear rectangular section convexity spring |
CN111460662A (en) * | 2020-03-31 | 2020-07-28 | 浙江大学 | Method for establishing aeroelastic model of flexible photovoltaic support structure |
CN115017681A (en) * | 2022-05-10 | 2022-09-06 | 西北工业大学 | Fine analysis method for dynamic characteristics of small-sag stay cable |
CN115046620A (en) * | 2021-12-22 | 2022-09-13 | 中南大学 | Structural vibration sensing method based on FBG strain sensor array |
CN115452226A (en) * | 2022-10-05 | 2022-12-09 | 大连理工大学 | Inhaul cable force identification algorithm considering semi-rigid constraints at two ends |
CN117271949A (en) * | 2023-11-23 | 2023-12-22 | 华东交通大学 | Suspension cable self-vibration analysis method and system considering elastic boundary and large sag influence |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5289120A (en) * | 1990-11-28 | 1994-02-22 | Skf Industrie S.P.A. | Railway axle bearing and sealing assembly with integrated rotational speed detector |
AU2004320069A1 (en) * | 2004-05-25 | 2005-12-08 | Al-S Technology Bv | Method and system for real-time non-destructive testing of projection welds and system for implementing such a method comprising force sensor means and displacement sensor means |
JP2006180519A (en) * | 2004-12-23 | 2006-07-06 | Gain Net Systems Kk | Hybrid connector and relay system therewith |
CN101201282A (en) * | 2007-12-20 | 2008-06-18 | 宁波大学 | Fundamental frequency identification method for detecting cord force of cable-stayed bridge |
US20100148589A1 (en) * | 2008-10-01 | 2010-06-17 | Hamam Rafif E | Efficient near-field wireless energy transfer using adiabatic system variations |
CN103226076A (en) * | 2013-04-01 | 2013-07-31 | 河海大学 | Method and apparatus for measuring early elastic modulus of cement-based materials |
CN103902504A (en) * | 2014-04-11 | 2014-07-02 | 北京工业大学 | Method for calculating inherent frequency of Euler-Bernoulli beam through improved differential transformation method |
US20150304144A1 (en) * | 2012-11-19 | 2015-10-22 | Orange | Method for determining at least one filter of a filter bank of a transmission or coding system, corresponding device and computer program |
CN106096178A (en) * | 2016-06-24 | 2016-11-09 | 哈尔滨大金工程试验检测有限公司 | A kind of bridge cable flexural rigidity identification method |
CN107808038A (en) * | 2017-10-12 | 2018-03-16 | 宁波大学 | A kind of method for solving of Arbitrary Boundaries constraints drag-line oscillation crosswise frequency |
-
2018
- 2018-05-16 CN CN201810465820.6A patent/CN108763674B/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5289120A (en) * | 1990-11-28 | 1994-02-22 | Skf Industrie S.P.A. | Railway axle bearing and sealing assembly with integrated rotational speed detector |
AU2004320069A1 (en) * | 2004-05-25 | 2005-12-08 | Al-S Technology Bv | Method and system for real-time non-destructive testing of projection welds and system for implementing such a method comprising force sensor means and displacement sensor means |
JP2006180519A (en) * | 2004-12-23 | 2006-07-06 | Gain Net Systems Kk | Hybrid connector and relay system therewith |
CN101201282A (en) * | 2007-12-20 | 2008-06-18 | 宁波大学 | Fundamental frequency identification method for detecting cord force of cable-stayed bridge |
US20100148589A1 (en) * | 2008-10-01 | 2010-06-17 | Hamam Rafif E | Efficient near-field wireless energy transfer using adiabatic system variations |
US20150304144A1 (en) * | 2012-11-19 | 2015-10-22 | Orange | Method for determining at least one filter of a filter bank of a transmission or coding system, corresponding device and computer program |
CN103226076A (en) * | 2013-04-01 | 2013-07-31 | 河海大学 | Method and apparatus for measuring early elastic modulus of cement-based materials |
CN103902504A (en) * | 2014-04-11 | 2014-07-02 | 北京工业大学 | Method for calculating inherent frequency of Euler-Bernoulli beam through improved differential transformation method |
CN106096178A (en) * | 2016-06-24 | 2016-11-09 | 哈尔滨大金工程试验检测有限公司 | A kind of bridge cable flexural rigidity identification method |
CN107808038A (en) * | 2017-10-12 | 2018-03-16 | 宁波大学 | A kind of method for solving of Arbitrary Boundaries constraints drag-line oscillation crosswise frequency |
Non-Patent Citations (2)
Title |
---|
PROF. SIVACOUMAR R等: "Analysis of Radio over fiber System’s Performance after optimizing Delay Interferometer’s parameters and Incorporating Low pass Chebyshev filters at the Receiver", 《2015 INTERNATIONAL CONFERENCE ON APPLIED AND THEORETICAL COMPUTING AND COMMUNICATION TECHNOLOGY (ICATCCT)》 * |
刘占生等: "拉伸弹簧与扭转弹簧模拟动态网格方法分析与改进", 《哈尔滨工业大学学报》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109829211A (en) * | 2019-01-21 | 2019-05-31 | 东南大学 | A kind of thermal environment lower plate structure high frequency partial method of response calculation |
CN109885898A (en) * | 2019-01-28 | 2019-06-14 | 华北水利水电大学 | The measuring method of the eigentone of non-linear rectangular section convexity spring |
CN111460662A (en) * | 2020-03-31 | 2020-07-28 | 浙江大学 | Method for establishing aeroelastic model of flexible photovoltaic support structure |
CN111460662B (en) * | 2020-03-31 | 2022-05-03 | 浙江大学 | Method for establishing aeroelastic model of flexible photovoltaic support structure |
CN115046620A (en) * | 2021-12-22 | 2022-09-13 | 中南大学 | Structural vibration sensing method based on FBG strain sensor array |
CN115046620B (en) * | 2021-12-22 | 2023-08-18 | 中南大学 | Structural vibration sensing method based on FBG strain sensor array |
CN115017681A (en) * | 2022-05-10 | 2022-09-06 | 西北工业大学 | Fine analysis method for dynamic characteristics of small-sag stay cable |
CN115017681B (en) * | 2022-05-10 | 2023-06-23 | 西北工业大学 | Fine analysis method for dynamic characteristics of small-sag inhaul cable |
CN115452226A (en) * | 2022-10-05 | 2022-12-09 | 大连理工大学 | Inhaul cable force identification algorithm considering semi-rigid constraints at two ends |
CN115452226B (en) * | 2022-10-05 | 2023-12-19 | 大连理工大学 | Cable force recognition algorithm considering semi-rigid constraint at two ends |
CN117271949A (en) * | 2023-11-23 | 2023-12-22 | 华东交通大学 | Suspension cable self-vibration analysis method and system considering elastic boundary and large sag influence |
CN117271949B (en) * | 2023-11-23 | 2024-03-12 | 华东交通大学 | Suspension cable self-vibration analysis method and system considering elastic boundary and large sag influence |
Also Published As
Publication number | Publication date |
---|---|
CN108763674B (en) | 2021-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108763674A (en) | A kind of method for solving of elastic boundary condition downhaul beam frequency | |
Tornabene et al. | Higher-order structural theories for the static analysis of doubly-curved laminated composite panels reinforced by curvilinear fibers | |
CN111783201B (en) | Rapid analysis method for dynamic characteristics of three-span self-anchored suspension bridge | |
Eisenberger | Exact static and dynamic stiffness matrices for general variable cross section members | |
Khare et al. | Free vibration of thick laminated circular and annular plates using three-dimensional finite element analysis | |
Banerjee | Explicit frequency equation and mode shapes of a cantilever beam coupled in bending and torsion | |
CN111783198A (en) | Suspension bridge refined power analysis method based on double-beam model | |
CN110287505B (en) | Aircraft stability analysis method | |
CN113076677B (en) | Flexible body structure high-order nonlinear finite element numerical simulation method based on quintic hermite function | |
Wang et al. | Assumed stress quasi‐conforming technique for static and free vibration analysis of Reissner–Mindlin plates | |
Zhang et al. | A Spatial Euler‐Bernoulli Beam Element for Rigid‐Flexible Coupling Dynamic Analysis of Flexible Structures | |
Barbone et al. | Optimal modal reduction of vibrating substructures | |
Ton That et al. | Nonlinear bending analysis of functionally graded plates using SQ4T elements based on twice interpolation strategy | |
CN109948180A (en) | A kind of orthotropy opposite side freely-supported rectangular thin plate vibration analysis method | |
Ritter et al. | Comparison of nonlinear aeroelastic methods for maneuver simulation of very flexible aircraft | |
Wang et al. | Treatment of elastically restrained ends for beam buckling in finite difference, microstructured and nonlocal beam models | |
CN117271949B (en) | Suspension cable self-vibration analysis method and system considering elastic boundary and large sag influence | |
CN112115616A (en) | Dynamic characteristic analysis method and device of power transmission tower | |
Rezaiee-Pajand et al. | Hybrid trefftz formulation for thin plate analysis | |
CN107220678B (en) | Multi-degree-of-freedom beam type structure nonlinear type determination method | |
Lee et al. | Frame elements with mixed formulation for singular section response | |
CN115452226A (en) | Inhaul cable force identification algorithm considering semi-rigid constraints at two ends | |
Carrera et al. | Comparisons between 1d (beam) and 2d (plate/shell) finite elements to analyze thin walled structures | |
CN109726454B (en) | Fluid-solid coupling modeling method and device for pipeline system | |
Morandini et al. | Implementation and validation of a 4-node shell finite element |
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 |