稳定性分析
Chapter 6 Stability
混沌系統與訊號處理實驗室We deal only with the stability of linear SISO time-invariant systems.
-Absolutely stable stable or not -Relatively stable how stable -Zero-state response -Zero-input response
-Total response= Zero-state response
+ Zero-input response
混沌系統與訊號處理實驗室
()()()()
()()22..2.. 1
:320,01,21.Absolutely stable 2.Relative stability ,
113.:22
341
:2
.22s t
t
t t z s t
t z i more stable
Ex x
x x u x x u t u t e e x s Zero state response Zero input res tea t e e x t pons dy state
e e e
??????++=====?+=???→ ()()()()
22....12homegenerous solutions
115.:22
1
2t t
t t z s z i Tot x t al resp x t x t particular sol K o e K e
n ution se e e ????=+=++=?+→
混沌系統與訊號處理實驗室Bounded-input Bounded-output stability (BIBO stability )
()()With zero initial conditions, the system is said to be BIBO stable or simply stable, if the output is bounded to bounded inp ev ut ery .
y t u t ()()()()()()()()()()()()0
1
0 ,, and is a positive finite number. , y t u t g d y t u t g d u t g d If u t is bounded u t M t M y L
M g d g t τττ
ττττττ
ττ∞
∞
∞
∞
=?∴=
?≤?≤?∈∴≤∫∫∫∫
混沌系統與訊號處理實驗室
()()()()
,
?If g d Q then y t M Q N g d BIBO stable ττττ∞
∞≤<∞≤?<∞∴<∞??∫∫
()()()()()()()()()()000000
0If is a right-half plane or axis pole , then , where and 01
st
st
s t
t
st
t
s s G s L g t g t e dt G s j of G s g t e dt g t e
dt
G s s j e G s g t d j e
e
t
s σσωσωσσω∞
??∞
∞
????∞
??=∞=+≥?≤∴∞==????=
≤?=≤+∫
∫∫∫
混沌系統與訊號處理實驗室()()1
Show that if a system with the impulse response is stable, then .
g t g t L ∈Pf ：()
()()()1
1
~ ~ g L
g L BIBO every bou Want to prove BIBO a boun nde ded input u t unbo d input u t y is boun unde ded y d
is ∈≡????????()()()()000
000. is a right-half s-plane pole not BIBO stable.. is on th sin is of t e he f , orm n .
si If u t t y t t i s ii s j axis s j uns ble t ta ωωωω????=?=
混沌系統與訊號處理實驗室()()()()
()
()()()()
()()()()0
S in c e ,
, 0C h o o se 0 , 0
1 .
.y t g u t d g if g g u if g u b o u n d e d y g u d g d is u n b o u n d e d τττττττττ
τττττ∞∞∞
=????≠?
???
?=??≤=?==∞∫
∫∫ For BIBO stability, the poles of G (s ) cannot be located in the right-half s -plane (or simply called right-half plane RHP ) or on the -axis. They must all lie in the left-half plane LHP .
j ω
混沌系統與訊號處理實驗室Zero-input and Asymptotic stability
()()
()0-Def: zero-input stable
If the zero-input response, , subject to the initial conditions, , as , then the system
is calle finite reaches zero zero-input st d ; otherwise, the system is unsta le bl ab e k y t y
t t →∞()()()
()
1
00
.n k k k y t g t y
t ?==∑()()()()()()()
()()01
1
0 1., 2.lim 0
,0
t n n k k k zero input stable y t M t t y t y t g t y
t g t t →∞
??∴??≤<∞?>=≤∴<∞?≥∑∑∵(or asymptotically stable)
混沌系統與訊號處理實驗室
-Let the n poles be expressed as , .Then if m of n poles are simple and the rest of multiple order,
i i i s j σω=+1,2,,i n
= ()1
1
i i m n m s t
s t
i i i i i y t K e L t e
??===+
∑∑
Ex:
()()()()
()
0121
2
2
1
12312334 211
L K K L G s s s s s s s s n m n m =
=++++++++++==???=For the condition , then,.
a system is BIBO stable zero-input stable or asymptotically
stable.()y t <∞0 i LHP poles σ
Marginally stable: simple roots on j -axis, none in the RHP.
Ex:
Ex:
()
()
()
111
11
=11
t
G s
s s s s
g t e?
==?
++
?→
ω
()
()
222
11
11111
cos sin
cos sin
t
A Bs C
G s
s s s s s
g t Ae B t C t
B t
C t as t
?
=?=++
+++++
=++
→+→∞
≠
won t diverage but
′→
混沌系統與訊號處理實驗室
混沌系統與訊號處理實驗室()()()()
()()
()()()()
()()()()()
()
()2
2
2
432
20
1.
1232012.
1222013. 2410
4.
10410
5. 3010stable s s s s unstable
s s s s Marginally stable s s stable multiple order pole s s Marginally stable s s s s
++++?++?+++++++：：：：：
混沌系統與訊號處理實驗室()1
110The criterion can determine whether the characteristic eq. of a linear time-invariant SISO system
where , has LHP roots.
0,n
n n n i F s a s a s
a s a a ??=+++∈?+=
Ex:
()65432
65432
6
10
6
5
420 53
F s a s a s a s a s a s a s a
s
s
a a a a
a a
=++++++=
1
546352615060
555
35150
4
3
2
1
a
a a a a a a a a a a a a
a
a a a
Aa a B Aa a a
A A
CB AD Ca A
a
C C
E
A B
C
s
s
s
s
s
D Ca
E
Fa
a
F
D
E
F
E
???
=
??
???
=
?
??
=
Routh tabulation
混沌系統與訊號處理實驗室
混沌系統與訊號處理實驗室
1.The coefficient of the first column are all positive LHP poles.
2.# of change of signs in the first col. # of roots in the RHP.
??Ex:
()()()()()
()3
2
32
1
0213460
2,3 1146
4116 2.5 #
42.2unstabl of e poles chan s s s s
s s s s s s ges of sign a s t ?????=?+?=?+++==+?
?()564062.5
????=
混沌系統與訊號處理實驗室
Ex:
()()
432
4
321235100
231015
132511020710
11
75110 6.43 7 s s s s unstable s s s s p ??????=?=????=?
?++++=()0
0.7555 1.4444# 21 0
s j of changes s o of sign les at s =±=
混沌系統與訊號處理實驗室
Special cases
1. Ex:
()2-3
0.4057 1.2928 # 2
unstable poles at s j and of changes of signs εε
4
3212230 1231203?
s s s s s s s s ++++=4
3
21
2-3
0123
23
3
1s s s s s ε
εε>?
nonzero
混沌系統與訊號處理實驗室2. Ex:?
No RHP poles and one pair of poles at
j
±543254
321488740
187 48 46644 0 0
s s s s s s s s s s +++++=543210
18748 4664 48 4
s s s s s s 2:4408Auxiliary equation s dA
s ds
+==zero row
LHP
()()()()()
54322322
2
488741474 1134s s s s s s s s s s s s s +++++=++++=++++No RHP roots
-axis
j ω
混沌系統與訊號處理實驗室
Ex:
()7
7
327
3
2
1073408.31204000 1.5100 112040003408.3 1.510410.36-1.5
0 and 410.36 1.50
K 100
3408.3
1.510s s K s K s s K
s
s K K +++×=∴>×??>×()2
9
0273.57roots are in the LH P
273.573408.3 4.1036100
1097.27
K If K A s s s j ∴<
混沌系統與訊號處理實驗室Ex:
()()3
2
3
2
10
2
3240
12343 2434
0, 36400.528
s Ks K s s s s s K K K K K K K
K K ++++=∴>+?>?>++?
稳定性方法评价
边坡稳定性评价方法概述 （辽宁工程技术大学土木与交通学院辽宁阜新123000 作者：张媛）对边坡稳定性评价方法进行了综述,有：极限平衡法、有限元法、离散单元 法、快速拉格朗日分析法、DDA法、流行元法、块体理论法、可靠度方法、模 糊综合评价法、灰色系统评价法、聚类分析法、神经网络、遗传算法和专家系统。在概要地叙述了各个方法的理论基础上，对各个方法的优缺点进行了叙述，指出了各自的适合条件以及目前的应用状况。其中极限平衡法、块体理论法很多时候 与实际情况不相符合，快速拉格朗日法具有随意性，DDA法在数学收敛上的实 现有一定的难度，有限元法需要定义合适的系数，模糊综合评价法和聚类分析法不能全面、最优，专家系统对于知识的获取具有一定的难度，综合各个方法，其中的离散单元法、流行元法、神经网络、遗传算法的适用性较好。 关键词：边坡稳定性；研究进展；评价方法 Prospect Methods of the Research on Slope Stability Zhang Yuan （ liaoning Technical University Civil Engineering and Transportation Department, Liaoning Fuxin 123000 ） Abstract: The paper reviews the prospect methods of the research on slope stability. There are Limit Equilibrium Method， Finite Element Method, Distinct Element Method, Fast Lagrangion Analysis of Method, Discontinuous Deformation Analysis, Manifold Element Method, Block Theory, Reliability Method, Comprehensive Fuzzy Evaluation, Grey system Evaluation, Clustering Analysis Method, Neural Network, Genetic Algorithm, Expert System. On the base of the theory summary about every method, the paper relate the advantages and disadvantages of these methods,points their suiting conditions and using state. In the outline, Limit Equilibrium Method and Block Theory cannot agree with the fact at the most time. Fast Lagrangion Analysis of Method is at its ease, There is a difficulty of math converge about Discontinuous Deformation Analysis, Finite Element Method needs to definite suitable coefficient, Comprehensive Fuzzy Evaluation and Clustering Analysis Method cannot give a overall result, or often it is not the best, Expert System has a
性能稳定性分析
性能稳定性分析 1功角的具体含义。 电源电势的相角差，发电机q轴电势与无穷大系统电源电势之间的相角差。 电磁功率的大小与δ密切相关，故称δ为“功角”或“功率角”。电磁功率与功角的关系式被称为“功角特性”或“功率特性”。 功角δ除了表征系统的电磁关系之外，还表明了各发电机转子之间的相对空间位置。 2功角稳定及其分类。 电力系统稳态运行时，系统中所有同步发电机均同步运行，即功角δ是稳定值。系统在受到干扰后，如果发电机转子经过一段时间的运动变化后仍能恢复同步运行，即功角δ能达到一个稳定值，则系统就是功角稳定的，否则就是功角不稳定。 根据功角失稳的原因和发展过程，功角稳定可分为如下三类： 静态稳定（小干扰） 暂态稳定（大干扰） 动态稳定（长过程） 3电力系统静态稳定及其特点。 定义：指电力系统在某一正常运行状态下受到小干扰后，不发生自发振荡或非周期性失步，自动恢复到原始运行状态的能力。如果能，则认为系统在该正常运行状态下是静态稳定的。不能，则系统是静态失稳的。 特点：静态稳定研究的是电力系统在某一运行状态下受到微小干扰时的稳定性问题。系统是否能够维持静态稳定主要与系统在扰动发生前的原始运行状态有关，而与小干扰的大小、类型和地点无关。 4电力系统暂态稳定及其特点。 定义：指电力系统在某一正常运行状态下受到大干扰后，各同步发电机保持同步运行并过渡到新的或恢复到原来的稳态运行状态的能力。通常指第一或第二振荡周期不失步。如果能，则认为系统在该正常运行状态下该扰动下是暂态稳定的。不能，则系统是暂态失稳的。 特点：研究的是电力系统在某一运行状态下受到较大干扰时的稳定性问题。系统的暂态稳定性不仅与系统在扰动前的运行状态有关，而且与扰动的类型、地点及持续时间均有关。 作业2 5发电机组惯性时间常数的物理意义及其与系统惯性时间常数的关系。 表示在发电机组转子上加额定转矩后，转子从停顿状态转到额定转速时所经过的时间。TJ=TJG*SGN/SB 6例题6-1 (P152) （补充知识：当发电机出口断路器断开后，转子做匀加速旋转。汽轮发电机极对数p=1。额定频率为50Hz。要求列写每个公式的来源和意义。）题目:已知一汽轮发电机的惯性时间常数Tj=10S,若运行在输出额定功率状态，在t=0时其出口处突然断开。试计算（不计调速器作用） （1）经过多少时间其相对电角度（功角）δ=δ0+PAI.(δ0为断开钱的值）（2）在该时刻转子的转速。 解:（1）Tj=10S,三角M*=1,角加速度d2δ/dt2=三角M*W0/Tj=W0/10=31.4RAD/S2 δ=δ0+0.5dd2δ/dt2 所以PI=0.5*2PI*f/10t方 t=更号10/50=0.447 （2）t=0.447时，
二阶瞬态响应特性与稳定性分析资料报告
广西大学实验报告纸 组长： 组员： 指导老师： 成绩： 学院：电气工程学院 专业：自动化 班级：163 实验容：实验五 二阶瞬态响应特性与稳定性分析 2018年5月11日 【实验时间】 2018年 5月 11日 【实验地点】 综合808 【实验目的】 1、以实际对象为基础，了解和掌握典型二阶系统的传递函数和模拟电路图。 2、观察和分析典型二阶系统在欠阻尼、临界阻尼、过阻尼的响应曲线。 3、学会用MATLAB 分析系统稳定性。 【实验设备与软件】 1、Multisim 10电路设计与仿真软件 2、labACT 试验台与虚拟示波器 3、MATLAB 数值分析软件 【实验原理】 1、被模拟对象模型描述 永磁他励电枢控制式直流电机如图1（a ）所示。根据Kirchhoff 定律和机电转换原理，可得如下方程 u k Ri dt di L e =++ω （1） l t T i k b dt d J -=+ωω （2） ωθ =dt d （3） 式中，各参数如图1（a ）所示：L 、R 为电机和负载折合到电机轴上的转动惯量，Tl 是折合到电机轴上的总的负载转矩，b 是电机与负载折合到电机轴上的粘性摩擦系数；kt 是转矩系数（Nm/A ），k e 是反电动势 系数（Vs/rad ）。令R L /e =τ(电磁时间常数)，b J /m =τ（机械时间常数） ，于是可由这三个方程 画出如图1（b ）的线性模型框图。 将Tl 看成对控制系统的扰动，仅考虑先行模型框图中()()s s U Θ→的传递函数为 ()()()()()s Rb k k s s Rb k s U s s G t e m e t 1 /11/?+++=Θ= ττ （4） 考虑到电枢电感L 较小，在工程应用中常忽略不计，于是上式转化为
边坡稳定性分析方法及其适用条件资料
边坡稳定性分析方法及其适用条件 摘要：边坡是一种自然地质体，在外力的作用下，边坡将沿其裂隙等一些不稳定结构面产生滑移，当土体内部某一面上的滑动力超过土体抗滑动的能力，将导致边坡的失稳。边坡稳定性分析是岩土工程的一个重要研究内容，并已经形成一个应用研究课题，本文对目前边坡稳定性分析中所采用的各种方法进行了归纳，并阐述了其适用条件。 关键词：边坡稳定性分析方法适用条件 正文： 一、工程地质类比法 工程地质类比法，又称工程地质比拟法，属于定性分析，其内容有历史分析法、因素类比法、类型比较法和边坡评比法等。该方法主要通过工程地质勘察，首先对工程地质条件进行分析，如对有关地层岩性、地质构造、地形地貌等因素进行综合调查和分类，对已有的边坡破坏现象进行广泛的调查研究，了解其成因、影响因素和发展规律等；并分析研究工程地质因素的相似性和差异性；然后结合所要研究的边坡进行对比，得出稳定性分析和评价。其优点是综合考虑各种影响边坡稳定的因素，迅速地对边坡稳定性及其发展趋势作出估计和预测；缺点是类比条件因地而异，经验性强，没有数量界限。 适用条件：在地质条件复杂地区，勘测工作初期缺乏资料时，都常使用工程地质类比法，对边坡稳定性进行分区并作出相应的定性评价，因此，需要有丰富实践经验的地质工作者，才能掌握好这种方法。
二、极限分析法 应用理想塑性体或刚塑性体处于极限状态的极小值原理和极大 值原理来求解理想塑性体的极限荷载的一种分析方法。它在土坡稳定分析时，假定土体为刚塑性体，且不必了解变形的全过程，当土体应力小于屈服应力时，它不产生变形，但达到屈服应力，即使应力不变，土体将产生无限制的变形，造成土坡失稳而发生破坏。其最大优点是考虑了材料应力—应变关系，以极限状态时自重和外荷载所做的功等于滑裂面上阻力所消耗的功为条件，结合塑性极限分析的上、下限定理求得边坡极限荷载与安全系数。 三、极限平衡法 该法将滑体作为刚体分析其沿滑动面的平衡状态,计算简单。但由于边坡体的复杂性,计算时模型的建立与参数的选取不可避免地使计算结果与实际结果不吻合。常用的方法有如下几种。 1瑞典条分法。基本假定:A边坡稳定为平面应变问题;B滑动面为圆弧;C计算圆弧面安全系数时,将条块重量向滑面法向分解来求法向力。该方法不考虑条间力的作用,仅能满足滑动体的力矩平衡条件,产生的误差使安全系数偏低。 优缺点：在不能给出应力作用下的结构图像的情况下，仍能对结构的稳定性给出较精确的结论，分析失稳边坡反算的强度参数与室内试验吻合度较好，使分析程序更加可信;但需要先知道滑动面的大致位置和形状，对于均质土坡可以通过搜索迭代确定其危险滑动面，但是对于岩质边坡，由于其结构和构造比较复杂，难以准确确定其滑动
稳定性分析答案
稳定性分析 2009-10-14 14:18 1功角的具体含义。 电源电势的相角差，发电机q轴电势与无穷大系统电源电势之间的相角差。 电磁功率的大小与δ密切相关，故称δ为“功角”或“功率角”。电磁功率与功角的关系式被称为“功角特性”或“功率特性”。 功角δ除了表征系统的电磁关系之外，还表明了各发电机转子之间的相对空间位置。 2功角稳定及其分类。 电力系统稳态运行时，系统中所有同步发电机均同步运行，即功角δ 是稳定值。系统在受到干扰后，如果发电机转子经过一段时间的运动变化后仍能恢复同步运行，即功角δ 能达到一个稳定值，则系统就是功角稳定的，否则就是功角不稳定。 根据功角失稳的原因和发展过程，功角稳定可分为如下三类： 静态稳定（小干扰） 暂态稳定（大干扰） 动态稳定（长过程） 3电力系统静态稳定及其特点。 定义：指电力系统在某一正常运行状态下受到小干扰后，不发生自发振荡或非周期性失步，自动恢复到原始运行状态的能力。如果能，则认为系统在该正常运行状态下是静态稳定的。不能，则系统是静态失稳的。 特点：静态稳定研究的是电力系统在某一运行状态下受到微小干扰时的稳定性问题。系统是否能够维持静态稳定主要与系统在扰动发生前的原始运行状态有关，而与小干扰的大小、类型和地点无关。 4电力系统暂态稳定及其特点。 定义：指电力系统在某一正常运行状态下受到大干扰后，各同步发电机保持同步运行并过渡到新的或恢复到原来的稳态运行状态的能力。通常指第一或第二振荡周期不失步。如果能，则认为系统在该正常运行状态下该扰动下是暂态稳定的。不能，则系统是暂态失稳的。 特点：研究的是电力系统在某一运行状态下受到较大干扰时的稳定性问题。系统的暂态稳定性不仅与系统在扰动前的运行状态有关，而且与扰动的类型、地点及持续时间均有关。 作业2 5发电机组惯性时间常数的物理意义及其与系统惯性时间常数的关系。 表示在发电机组转子上加额定转矩后，转子从停顿状态转到额定转速时所经过的时间。TJ=TJG*SGN/SB 6例题6-1 (P152) （补充知识：当发电机出口断路器断开后，转子做匀加速旋转。汽轮发电机极对数p=1。额定频率为50Hz。要求列写每个公式的来源和意义。）题目:已知一汽轮发电机的惯性时间常数Tj=10S,若运行在输出额定功率状态，在t=0时其出口处突然断开。试计算（不计调速器作用） （1）经过多少时间其相对电角度（功角）δ=δ0+PAI.(δ0为断开钱的值）（2）在该时刻转子的转速。 解:（1）Tj=10S,三角M*=1,角加速度d2δ/dt2=三角M*W0/Tj=W0/10=S2 δ=δ0+δ/dt2 所以PI=*2PI*f/10t方 t=更号10/50=
两种边坡稳定性分析方法比较研究
第10卷 第10期 中 国 水 运 Vol.10 No.10 2010年 10月 China Water Transport October 2010 收稿日期：2010-06-11 作者简介：马玉岩（1987-），男，黑龙江绥化人，武汉大学水资源与水电工程科学国家重点实验室水利水电工程施工与 管理专业硕士研究生，主要研究方向为岩土边坡工程研究以及结构设计。 两种边坡稳定性分析方法比较研究 马玉岩 （武汉大学 水资源与水电工程科学国家重点实验室，湖北 武汉 430072） 摘 要：以某水电工程岩质高边坡做为实例，将强度折减理论与FLAC3D 软件相结合，通过有限差分程序FLAC3D 软件来模拟分析其稳定性。并与极限平衡方法的分析结果对比，探索两种方法的差异性与结果的可靠性，为确定适合工程建设实际的岩质边坡稳定分析方法提出了有益的参考。 关键词：强度折减法；极限平衡法；边坡稳定性 中图分类号：P642.1 文献标识码：A 文章编号：1006-7973（2010）10-0197-03 一、引言 目前，国内在建和待建的大型水电工程大多坐落在西南、西北高山峡谷地区。我国的水电建设面临着一系列高边坡稳定问题。在现代岩土工程和科学技术的新成就的支持下，确定适合工程建设实际的岩质边坡稳定分析方法，是摆在水利水电工程技术人员面前的任务[1]。 目前工程实践中岩质边坡稳定性定量分析主要有三种方法：解析法（最常用的是极限平衡法）、数值方法和概率法。极限平衡法是最常用的解析法，它是在边坡滑动面确定的情况下，根据滑裂面上抗滑力和滑动力比值直接计算安全系数，此外，关键块理论也属于这样的确定性分析方法。数值方法则是借助计算机进行数值分析（例如有限元、快速拉格朗日分析法、离散元、块体元和DDA 等）从而确定边坡的位移场和应力场，再用超载法、强度折减法等使边坡处于极限状态，从而间接得到安全系数。这种方法同时可以考虑位移协调条件和岩体本构关系等。概率法是将概率统计理论被引用到边坡岩体的稳定性分析中来，它通过现场调查，以获得影响边坡稳性影响因素的多个样本，然后进行统计分析，求出它们各自的概率分布及其特征参数，再利用某种可靠性分析方法，来求解边坡岩体的破坏概率即可靠度[2]。 文中选用某水电工程岩质高边坡做为实例，采用强度折减法和极限平衡法对岩质高边坡的稳定性进行对比分析。 二、边坡工程地质条件 模型宽约为700m，高约为700m。 基岩以中粒结构的灰白色、微红色黑云二长花岗岩为主，并有辉绿岩脉（β）、花岗细晶岩脉、闪长岩脉等各类脉岩穿插发育于花岗岩中，尤以辉绿岩脉分布较多。建模过程中考虑了岩体中对边坡稳定影响较大的几个岩脉。 根据岩体风化特点，岸坡岩体由表向内可划分为全风化带、强风化带、弱风化带、微风化—新鲜岩体。岩体风化的水平、垂直分带性明显。 边坡内无地下水分布。 边坡剖面如图1 所示。 图1 边坡剖面 三、强度折减法 强度折减系数法的基本原理是将坡体强度参数凝聚力c 和内摩擦角f 值同时除以一个安全系数K，得到一组新的c k 、f k 值，然后作为新的资料参数输入，再进行试算，当计算不收敛时，对应的K 被称为坡体的最小稳定安全系数，此时坡体达到极限状态，发生剪切破坏，同时可得到坡体的破坏滑动面。 FLAC3D （Three Dimensional Fast Lagrangian Analysis of Continua）是美国Itasca Consulting Goup lnc 开发的三维快速拉格朗日分析程序。该程序能较好地模拟地质材料在达到强度极限或屈服极限时发生的破坏的力学行为，特别适用于分析渐进破坏和失稳。 文中利用FLAC3D，采用“二分法”[3]实现强度折减法，求解安全系数。 所建计算模型节点为29,646个，单元为24,005个。模型的边界条件：模型四周法向约束，底部固定约束，顶部自由，仅受重力作用。 研究表明，随着剪胀角的增大，安全系数也逐渐增大[4]。不过，Vermeer 和de Borst（1984年）研究证明，一般土体、岩石和混凝土的剪胀角要比它们的摩擦角小得多，且通常在0°～20°内变化[5]。因此，剪胀角对强度折减法计算
原料药稳定性试验报告
L- 腈化物稳定性试验报告 一、概述 L-腈化物是L- 肉碱生产过程中的第一步中间体（第二步中间体： L-肉碱粗品；第三步中间体：L-肉碱潮品），由于L- 肉碱生产工艺为 间歇操作，即每生产一步中间体，生产完毕并出具合格检测报告后，存 入中间体仓库，以备下一步生产投料所需。根据本公司L- 肉碱产品的 整个生产周期，L- 腈化物入库后可能存放的最长时间为4 周（约28 天）。以此周期为时间依据制定了L- 腈化物稳定性试验方案，用于验 证L-腈化物在再试验期限内的各项质量指标数据的稳定性，并且能否符 合L- 腈化物的质量标准，此次稳定性试验的整个周期为28 天，具体 的稳定性试验方案以ICH 药物稳定性指导原则为基础制定，以确保L- 腈化化物稳定性试验的可操作性。 二、验证日期 2010 年1 月13 日- 2010 年2 月10 日 三、验证方案 1）样品储存和包装： 考虑到L- 腈化物今后的贮藏、使用过程，本次用于稳定性试验的样品 批次与最终规模生产所用的L- 腈化物的包装和放置条件相同。 2）样品批次选择：此次稳定性试验共抽取三批样品，且抽取样品的批次与 最终规模生产时的合成路线和生产工艺相同
3）抽样频率和日期：从2010.1.13 起，每隔7 天取样一次，共取五次，具体日期为：2010.1.13 、2010.1.20 、2010.1.27 、 2010.2.3 、2010.2.10 ，以确保试验次数足以满足L- 腈化物的稳 定性试验的需要。。 4）检测项目：根据L- 腈化物的质量标准的规定，此次稳定性试验的检测项目共五项，分别为外观、氯含量、熔点、比旋度、干燥失重。这 些指标在L- 腈化物的储存过程中可能会发生变化，且有可能影响 其质量和有效性。 5）试样来源和抽样：L- 腈化物由公司102 车间生产，经检测合格后储存于中间体仓库，本次稳定性试验的L- 腈化物均取自于该中间体仓 库，其抽样方法和抽样量均按照L- 腈化物抽样方案进行抽样。抽 样完毕后直接进行检测分析，并对检测结果进行登记，保存，作为稳 定性数据评估的依据。 四、稳定性试验数据变化趋势分析及评估 通过对三批L- 腈化物的稳定性试验，对其物理、化学方面稳定性资料进行评价，旨在建立未来相似情况下，大规模生产出的L- 腈化物是否适用 现有的再试验期（28天）。批号间的变化程度是否会影响未来生产的
热稳定性分析方法
版 本 号：0.1 页 码：1/3 发布日期：2009-12-09 实验室程序 编 写： 批 准： 签 发： 文件编号：SHLX\LAB\L2-008 题 目：热稳定性测量方法 1.0 目的 提供了产品热稳定性的测量方法。 2.0 概述 （1）原理 Na 2SO 3 方 法 ： 用 1N 的 Na 2SO 3 溶 液 吸 收 样 品 粒 子 中 释 放 的 甲 醛 ， 生 成HOCH 2SO 3Na 和 NaOH 。 CH 2O ＋Na 2SO 3＋H 2O →HOCH 2SO 3Na ＋NaOH （2）本测量方法是利用聚甲醛树脂在高温熔融，产生甲醛气体，随氮气带出，被亚 硫酸钠溶液吸收，由滴定反应生成的氢氧化钠，得出甲醛含量。 3.0 仪器和试剂 【仪器】 （1） 油浴（容量约为 130L ，并配有样品熔融管） （2） 加热器 （3） 过热保护装置 （4） 搅拌器 （5） 自动滴定装置 （6） 数据处理计算机 【试剂】 （1） 0.005mol/l 硫酸 （2） 福尔马林(36.0～38.0%) （3） 亚硫酸钠(Na 2SO 3) （4） 缓冲液(pH 6.86) （5） 缓冲液(pH 9.18) （6） 0.1mol/l NaOH 4.0 定义 甲醛含量通过以下方式表示： （1）K 0 ：表示从 2 分钟到 10 分钟之间，聚合物中溶解的甲醛，不稳定端基和聚合 物主链分解出来的甲醛量。转化为每分钟的甲醛含量。 （2）K 1 ：表示从 10 分钟到 30 分钟之间，聚合物中剩余的溶解甲醛，不稳定端基
文件编号：SHLX\LAB\L2-008 和聚合物主链分解出来的甲醛量。转化为每分钟的甲醛含量。 （3）K2：表示从50 分钟到90 分钟之间，聚合物不稳定端基和聚合物主链分解出来的甲醛量。转化为每分钟的甲醛含量。 5.0安全注意事项 （1）搁置和取出样品过程中，要穿戴安全手套，以防被烫伤。 （2）电极容易损坏，使用时防止碰撞。 （3）作业时，穿戴安全眼镜和防护手套。 （4）实验过程中使用氮气作为载气，所以要控制好氮气流量，并确保良好的通风。6.0步骤 6.1准备 (1) 确认油浴温度223±2℃，硫酸溶液的量。 (2) 打开参比液添加孔，检查电极内饱和KCL 的量，确保液位超过甘汞位置。 (3) 打开自动电位滴定仪、打印机及电脑电源。 (4) 打开电脑桌面上AT-WIN，输入密码并确认与自动电位滴定仪联机。 (5) 调整氮气流量到60 l/h。 (6) 分别用pH 为6.86（25℃）、9.18（25℃）的缓冲液，对电极进行校正（根据 电脑提示进行），若显示“OK”，则校正通过，否则进行检查并重复校正步 骤。 (7) 对自动电位滴定仪进行排气，确保滴定管路中无气泡。 (8) 用250ml 的烧杯，取150ml 吸收液（1mol/L 亚硫酸钠溶液，它的配制方法： 将250g 的Na 2SO3溶于2000ml 的水中，充分搅拌。），放入磁性搅拌子、加 盖、并将电极、Ｎ２管、喷嘴插入溶液中，启动搅拌按钮。 (9) 用硫酸溶液（0.1N）将溶液pH 调节至9.10，待稳定后，用0.1mol/l 甲醛溶 液（配制方法：将81g 的福尔马林放入1L 的容量瓶中，然后加水到刻度线， 配成约0.1mol/l 福尔马林），调节pH 至9.21～9.22，并稳定10 分钟以上。 (10) 电极浸泡液的配制方法：PH=4 的缓冲试剂250ml 一包溶于250ml 水中， 再加入56gKCL，适当加热，搅拌至完全溶解。 6.2步骤 (1) 用铝皿取3.000±0.003g，将其放到小金属底部，然后用钩子，将准备好的 样品放入油浴的熔融管中。 (2) 盖紧硅胶塞，快速按下START，开始试验，试验过程控制pH 值为9.20。 (3) 当实验进行到设定的时间后，自动结束。（按“RESET”键，可手动停止实 验。）测定结束，打印机自动打印结果。 (4) 取出金属筒冷却，取出电极，并将电极放入浸泡液中。
结构动力稳定性的分析方法与进展_何金龙
结构动力稳定性的分析方法与进展 何金龙1,法永生2 (1.卓特建筑设计有限公司,广东佛山528322;2.上海大学土木工程系,上海200074) 【摘 要】 就目前结构动力稳定性问题这一研究领域的若干基本问题,常用的处理方法,判别准则与实验研究方法以及目前取得的主要成果作了简要总结和综述,并且对结构动力稳定性分析与研究今后的发展方向进行了展望。 【关键词】 结构; 动力稳定性; 处理方法; 判别准则; 实验研究 【中图分类号】 T U311.2 【文献标识码】 A 根据结构承受荷载形式的不同,可以将结构稳定问题分为静力稳定和动力稳定两大类。动力载荷作用下结构的稳定性问题是一个动态问题,由于时间参数的引入,使问题变得极为复杂。对于结构动力稳定性的定义一直难以确切给出,这是因为结构自身动力特性具有复杂性使得其在数学意义上的定义很难予以准确表达[1]。长期以来,力学工作者致力于结构稳定性问题的研究,在发展了经典稳定性理论的同时也极大地推动了动力稳定理论研究的前进。如稳定性判定准则的建立、临界载荷的确定、初缺陷的影响或后分叉分析等。理论分析和实验研究逐渐增多,使得这门学科不仅在理论上形成了一个庞大而复杂的体系,而且具有重要的实用价值。可以说,现在的结构动力稳定性研究分析已经是结构动力学、有限元法、数值计算方法及程序设计等诸多学科相互交叉、有机结合的产物,属于现代工程结构研究领域中的一个重要分支。 1 结构动力稳定性的分类及主要的研究问题 结构动力稳定性就其承载的动力形式大致可以分为三类。 (1)结构在周期性荷载作用下的动力稳定性。在简谐荷载等周期性荷载作用下,当结构的自振频率与外载荷的强迫振动频率非常接近时,结构将产生强烈的共振现象;当结构的横向固有振动频率与外荷载的扰动频率之间的比值形成某种特定的关系时,结构将产生强烈的横向振动,即参数振动。对于这类问题,前苏联学者符华·鲍络金(Bolito n)在其著作《弹性体系的动力稳定》中给出了较全面的分析和论述。他们导出的区分稳定区和不稳定区的临界状态方程是一个周期性方程,即M athieu-Hill方程。在周期相同的解之间存在着不稳定区域,便把问题归结为确定微分方程具有周期解的条件,从而解决了稳定的判别问题。但是对于大变形的几何非线形结构,结构的刚度矩阵需要经过迭代,微分方程非常复杂,这些理论将难以成立。 (2)结构在冲击荷载作用下的动力稳定性。在这种情况下,结构的动力稳定性与冲击类型密切相关,而且首要问题在于合理、实用的判别准则,它不仅要在逻辑上站得住脚,又要在实际上可行,遗憾的是这个问题至今未能形成一致的看法。目前对结构承受瞬态冲击作用下的冲击稳定性的试验和理论研究主要集中在理想脉冲以及阶跃荷载下的动力稳定性。在脉冲荷载作用下发生的动力屈曲称为脉冲屈曲,已有的研究表明[2][3][4],脉冲屈曲是一类响应式屈曲或者动力发展型屈曲。阶跃荷载是一类具有恒定幅值和无限长持续时间的载荷形式。在试验或者实际当中,固体与固体之间的冲击引起的屈曲就可看作脉冲冲击。 (3)结构在随动荷载作用下的动力稳定性。所谓随动荷载是指随着时间的变化荷载的幅值保持不变而方向发生变化的作用力,它是非保守力。它的分析将极其复杂,目前还难以见到可借鉴的动力稳定性分析文献。因此,许多学者通常采用结构动力学响应分析常用的手段,将这类荷载作为确定性荷载进行分析。通过对结构的动力平衡路径全过程进行跟踪,根据结构的各参数在动力平衡路径中的变化特性,对结构的动力稳定性进行有效的判定[5]。 综上所述,目前国内外动力稳定性研究的现状大致为:对周期荷载下的参数动力稳定性问题、在冲击荷载作用下的冲击动力稳定性问题和阶跃荷载下的参数阶跃动力稳定性问题研究较多,并取得了满意的效果[6][7][8]。恒幅阶跃载荷及矩形脉冲载荷或其它冲击载荷作用下杆的动力稳定问题也有很多研究,并从不同的角度建立了一些稳定性判定准则。但冲击载荷作用下板的动力稳定问题还没有获得广泛和深入的研究。对于较为复杂的冲击荷载作用下结构的动力稳定性问题,目前的研究主要集中于理想脉冲载荷和阶跃载荷作用下结构的动力稳定问题。在这类问题的分析中,最常采用的屈曲准则有B-R准则、Simitses总势能原理和放大函数法。对非周期激振、参数激振和强迫激振耦合引起的动力稳定问题研究较少;对弹性基本构件和简单模型研究较多(如周期激励下的柱子、梁、拱及壳等已得到了成功的分析),对复杂工程结构研究较少。对于在地震、风荷载等任意动力荷载作用下的具有较强的几何非线性的结构的动力稳定性问题,国内外这方面的文献资料虽然最近几年也有一些,但距离真正地合理解决这类动力稳定性问题还有许多工作要做。 [收稿日期]2006-06-12 [作者简介]何金龙(1962～),男,工学学士,一级注册结构工程师,主要从事工业与民用建筑设计工作。 155 ·工程结构· 四川建筑 第27卷2期 2007.04
控制系统的稳定性分析
精品 实验题目控制系统的稳定性分析 一、实验目的 1．观察系统的不稳定现象。 2．研究系统开环增益和时间常数对稳定性的影响。 二、实验仪器 1．EL-AT-II型自动控制系统实验箱一台 2．计算机一台 三、系统模拟电路图 系统模拟电路图如图3-1 图3-1 系统模拟电路图R3=0～500K； C=1μf或C=0.1μf两种情况。 四、实验报告 1．根据所示模拟电路图，求出系统的传递函数表达式。 G(S)= K=R3/100K,T=CuF/10 2．绘制EWB图和Simulink仿真图。
精品 3．根据表中数据绘制响应曲线。 4．计算系统的临界放大系数，确定此时R3的值，并记录响应曲线。 系统响应曲线 实验曲线Matlab （或EWB）仿真 R3=100K = C=1UF 临界 稳定 （理论值 R3= 200K） C=1UF
精品 临界 稳定 （实测值 R3= 220K） C=1UF R3 =100K C= 0.1UF
精品 临界 稳定 （理论 值R3= 1100 K） C=0.1UF 临界稳定 （实测值 R3= 1110K ） C= 0.1UF
精品 实验和仿真结果 1．根据表格中所给数据分别进行实验箱、EWB或Simulink实验，并进行实验曲线对比，分析实验箱的实验曲线与仿真曲线差异的原因。 对比： 实验曲线中R3取实验值时更接近等幅振荡，而MATLAB仿真时R3取理论值更接近等幅振荡。 原因： MATLAB仿真没有误差，而实验时存在误差。 2．通过实验箱测定系统临界稳定增益，并与理论值及其仿真结果进行比较（1）当C=1uf，R3=200K(理论值)时，临界稳态增益K=2， 当C=1uf，R3=220K(实验值)时，临界稳态增益K=2.2，与理论值相近（2）当C=0.1uf，R3=1100K(理论值)时，临界稳态增益K=11 当C=0.1uf，R3=1110K(实验值)时，临界稳态增益K=11.1，与理论值相近 四、实验总结与思考 1．实验中出现的问题及解决办法 问题：系统传递函数曲线出现截止失真。 解决方法：调节R3。 2．本次实验的不足与改进 遇到问题时，没有冷静分析。考虑问题不够全面，只想到是实验箱线路的问题，而只是分模块连接电路。 改进：在实验老师的指导下，我们发现是R3的取值出现了问题，并及时解决，后续问题能够做到举一反三。 3．本次实验的体会 遇到问题时应该冷静下来，全面地分析问题。遇到无法独立解决的问题，要及时请教老师，
整车操纵稳定性仿真分析报告分析解析
L11整车操纵稳定性仿真分析报告 （HB11A/HB12A） 编制（日期） 校对（日期） 审核（日期） 批准（日期） 简式国际汽车设计（北京）有限公司 L11整车操纵稳定性仿真分析报告（HB11A/HB12A） 1.定半径稳态圆周试验 试验方法 HB11A处于满载状态，沿半径为40m的定半径圆周进行回转运动，开始以最低稳定速度进入圆周，找准方向盘的位置，使汽车可以沿圆周进行回转运动，开始记录，然后缓慢连续而均匀地加速（纵向加速度不超过m/s2），加速的同时调整方向盘转角以维持定半径圆周运动，这个过程中车辆不应超出车道m，直至不能维持稳态定半径圆周运动条件时或受发动机功率限制所能达到的最大侧向加速度为止。记录整个过程，建议使用满足试验条件的最高档位。试验按向左转和向右转两个方向进行，每次试验开始时车身应处于正中位置。 数据处理 “方向盘转角——侧向加速度”拟合曲线线性部分的斜率，取侧向加速度为时的曲线斜率。 图1 方向盘转角—侧向加速度（左转） 从图1 计算得到左转不足转向梯度为137o/g
图2 方向盘转角—侧向加速度（右转） 右转不足转向梯度为g，则HB11A平均不足转向梯度为g。 HB11A的角传动比约为，则不足转向梯度/转向系角传动比为g。 “质心侧偏角——侧向加速度”拟合曲线线性部分的斜率，取侧向加速度为时的曲线斜率。 图3 质心侧偏角——侧向加速度（左转） 左转侧偏角梯度为g。 图4 质心侧偏角——侧向加速度（右转） 右转侧偏角梯度为g，则HB11A平均侧偏角梯度为g。 “车身侧倾角——侧向加速度”拟合曲线线性部分的斜率，取侧向加速度为时的曲线斜率。 图5 车身侧倾角——侧向加速度（左转） 左转侧倾角梯度为g。 图6 车身侧倾角—侧向加速度（右转） 右转侧倾角梯度为g，则HB11A平均侧倾角梯度为g。 2.方向盘转角阶跃输入试验 试验方法 HB11A处于满载状态，以70ｋｍ/ｈ的车速稳定直线行驶，开始记录数据，以尽可能快的速度(阶跃时间为转动方向盘，达到预定的转角，保持方向盘转角不变直至汽车恢复稳定状态，试验过程中油门踏板开度应尽可能保持不变。方向盘转角初始值是10°，每次增加5°，直到车辆达到附着极限，试验分为向左、向右两个方向进行。 数据处理 —方向盘转角滞后时间 横摆角速度达到50%稳态值时相对于方向盘转角达到50%阶跃值时的滞后时间。 图7 时横摆角速度—方向盘转角滞后时间 左转时，横摆角速度——方向盘转角滞后时间为
药物稳定性试验统计分析方法
药物稳定性试验统计分析方法 在确定有效期的统计分析过程中，一般选择可以定量的指标进行处理，通常根据药物含量变化计算，按照长期试验测定数值，以标示量％对时间进行直线回归，获得回归方程，求出各时间点标示量的计算值（y'）,然后计算标示量（y'）95％单侧可信限的置信区间为y'±z ，其中： 2 2 02)()(1X Xi X X N S t z N -∑-+ ??=- （12-21） 式中，t N -2—概率0.05，自由度N-2的t 单侧分布值（见表12-4），N 为数组；X 0—给定自变量；X —自变量X 的平均值； 2 -= N Q S （12-22） 式中，xy yy bL L Q -=；L yy —y 的离差平方和，N y y L yy /)(2 2∑-∑=；L xy —xy 的离差乘 积之和N y x xy L xy /))((∑∑-∑=；b —直线斜率。 将有关点连接可得出分布于回归线两侧的曲线。取质量标准中规定的含量低限（根据各品种实际规定限度确定）与置信区间下界线相交点对应的时间，即为药物的有效期。根据情况也可拟合为二次或三次方程或对数函数方程。 此种方式确定的药物有效期，在药物标签及说明书中均指明什么温度下保存，不得使用“室温”之类的名词。 例：某药物在温度25±2℃，相对温度60±10%的条件下进行长期实验，得各时间的标示量如表12-4。 表12-4 供试品各时间的标示量 时间/月 0 3 6 9 12 18 标示量/% 99.3 97.6 97.3 98.4 96.0 94.0 以时间为自变量（x ），标示量%（y ）为因变量进行回归，得回归方程 y= 99.18－0.26x ，r=0.8970，查T 单侧分布表，当自由度为4，P=0.05得 t N -2=2.132 9279.04 444 .32==-= N Q S 210)(2=-∑X X i
原料药稳定性试验报告
L-腈化物稳定性试验报告 一、概述 L-腈化物是L-肉碱生产过程中的第一步中间体（第二步中间体：L-肉碱粗品；第三步中间体：L-肉碱潮品），由于L-肉碱生产工艺为间歇操作，即每生产一步中间体，生产完毕并出具合格检测报告后，存入中间体仓库，以备下一步生产投料所需。根据本公司L-肉碱产品的整个生产周期，L-腈化物入库后可能存放的最长时间为4周（约28天）。以此周期为时间依据制定了L-腈化物稳定性试验方案，用于验证L-腈化物在再试验期限内的各项质量指标数据的稳定性，并且能否符合L-腈化物的质量标准，此次稳定性试验的整个周期为28天，具体的稳定性试验方案以ICH药物稳定性指导原则为基础制定，以确保L-腈化化物稳定性试验的可操作性。 二、验证日期 2010年1月13日----2010年2月10日 三、验证方案 1）样品储存和包装： 考虑到L-腈化物今后的贮藏、使用过程，本次用于稳定性试验的样品批次与最终规模生产所用的L-腈化物的包装和放置条件相同。 2）样品批次选择：此次稳定性试验共抽取三批样品，且抽取样品的批次与最终规模生产时的合成路线和生产工艺相同
3）抽样频率和日期：从2010.1.13起，每隔7天取样一次，共取五次，具体日期为：2010.1.13、2010.1.20、2010.1.27、2010.2.3、2010.2.10，以确保试验次数足以满足L-腈化物的稳定性试验的需要。。 4）检测项目：根据L-腈化物的质量标准的规定，此次稳定性试验的检测项目共五项，分别为外观、氯含量、熔点、比旋度、干燥失重。这些指标 在L-腈化物的储存过程中可能会发生变化，且有可能影响其质量和有效 性。 5）试样来源和抽样：L-腈化物由公司102车间生产，经检测合格后储存于中间体仓库，本次稳定性试验的L-腈化物均取自于该中间体仓库，其抽 样方法和抽样量均按照L-腈化物抽样方案进行抽样。抽样完毕后直接进 行检测分析，并对检测结果进行登记，保存，作为稳定性数据评估的依 据。 四、稳定性试验数据变化趋势分析及评估 通过对三批L-腈化物的稳定性试验，对其物理、化学方面稳定性资料进行评价，旨在建立未来相似情况下，大规模生产出的L-腈化物是否适用现有的再试验期（28天）。批号间的变化程度是否会影响未来生产的L-腈化物在再试验期内是否仍符合其质量规范。本次试验数据以表格、图解的形式给出，从而对L-腈化物的稳定性数据进行有效的评估。
自动控制实验报告一-控制系统的稳定性分析
实验一控制系统的稳定性分析 一、实验目的 1．观察系统的不稳定现象。 2．研究系统开环增益和时间常数对稳定性的影响。 二、实验仪器 1．自动控制系统实验箱一台 2．计算机一台 三、实验内容 系统模拟电路图如图 系统模拟电路图 其开环传递函数为: G（s）=10K/s(0.1s+1)(Ts+1) 式中K1=R3/R2，R2=100K，R3=0～500K；T=RC，R=100K，C=1f或C=0.1f两种情况。 四、实验步骤
1.连接被测量典型环节的模拟电路。电路的输入U1接A/D、D/A卡的DA1输出，电 路的输出U2接A/D、D/A卡的AD1输入，将纯积分电容两端连在模拟开关上。检查无误后接通电源。 2.启动计算机，在桌面双击图标[自动控制实验系统] 运行软件。 3.在实验项目的下拉列表中选择实验三[控制系统的稳定性分析] 5.取R3的值为50K，100K，200K，此时相应的K=10，K1=5，10，20。观察 不同R3值时显示区内的输出波形(既U2的波形)，找到系统输出产生增幅振荡时相应的R3及K值。再把电阻R3由大至小变化，即R3=200k，100k，50k，观察不同R3值时显示区内的输出波形, 找出系统输出产生等幅振荡变化的R3及K 值，并观察U2的输出波形。 五、实验数据 1模拟电路图
2．画出系统增幅或减幅振荡的波形图。C=1uf时： R3=50K K=5： R3=100K K=10
R3=200K K=20： 等幅振荡：R3=220k:
增幅振荡：R3=220k： R3=260k:
C=0.1uf时：R3=50k:
稳定性评价报告
福鼎市白琳玄武岩矿山北坡地质灾害点治理后斜坡 稳定性评价报告 1、概况 1.1矿区概况 福鼎大嶂山玄武岩矿山位于福鼎城关193°方向，平距20km 处，隶属福鼎市白琳镇山后山村管辖。地理坐标：东经120°09′48.3″--120°10′24.6″，北纬27°9′16.3″--27°9′39″。矿山到白琳镇约5公里。由白琳镇到福鼎八尺门约10公里可与国道主干线沈海高速福鼎至宁德段高速公路相连；温州至福州铁路经过白琳；交通便利（详见交通位置图1）。 福鼎市 27° 省 20km 寿宁 泰顺 柘荣 周宁 往福州 福安市 宁德市 120° 120° 霞浦江 浙 交 通 位 置 图 图1 10 溪潭 南阳 三沙 下白石赛岐 溪南 沙江 长春 下浒 27° 三都澳 福 宁 高 速 路 福安连接线 湾坞 往古田 往屏南 白琳 秦屿 沙埕 苍南 往政和 嵛山 白岩 东海 弃渣场位置 温福 铁路
1.2矿山北坡地质灾害点概况 福鼎白琳玄武岩矿山开发建设始于20世纪80年代初期，由3家公司于不同位置分别对白琳玄武岩体进行掠夺性开采。采区按地理位置分为北坡采场、东坡采场和南坡采场。1997年以前，由于无序开采和监管缺失，北坡采场剥离层剥离后形成的大量废石土就地堆弃于邻近采场的北坡冲沟内。随着时间的推移，无序开采造成白琳玄武岩矿山北坡的废石土超量排放。期间最大排放的废石土总量超过200万m3，大大超出北坡地质环境承载能力。由于北坡废石土的超量排放，致使北坡内及边缘曾多次发生小规模滑坡地质灾害。最为严重是于1998年2月18日受强降雨影响，北坡地质灾害点发生大面积的山体滑坡，滑坡规模在100万m3以上，由于大规模滑坡堵塞沟谷，影响场地内大气降水的自然排泄，并由于进一步引发大规模的泥石流地质灾害，造成18人员死亡、村落毁灭和公路毁坏交通中断的重大事故。泥石流的流通区长度达1km以上，堆积区长度达1km。此后，通过福鼎市政府干预，对矿山无序开采进行整顿，对3个采场进行整合，由福建白琳玄武石材有限公司通过组织白琳玄武岩的开采、经营，并择址建设南坡排土场，集中排放矿山建设、开采所形成的废石土。由于北坡弃碴系历史原因形成，福鼎玄武石材有限公司成立后未对北坡碴进行根本性治理。 2010年12月，受持续强降雨影响，白琳玄武岩矿山北坡临近采场的陡坡坡顶面以及矿山道路路面等出路弃碴的地段出现多道长30～50m，宽度5～15cm，深度0.3～1.5m的裂缝，局部裂缝下错约0.2～0.3m。陡坡坡底的缓坡地段也出现多道长20～30m，宽度5～10cm，深度0.3～1.5m的裂缝，局部裂缝下错约0.1～0.3m。随后裂缝灾害的空间进一步发展，于北坡西侧的冲