实干、实践、积累、思考、创新。
- 程序图标 ( Program Icon )
- 程序介绍 ( Program Introduction)
GMP ( A tool for Calculating Ground Motion Parameters for Seismic Analysis of Structures) 是一款结构抗震分析地震动强度指标/地震动参数计算工具,程序将地震波积分+反应谱分析+地震动参数分析等功能集合于一身,使用方便。GMP的主要功能包括:
(1) 丰富的地震波数据格式导入支持
(2) 地震动加速度、速度、位移时程
(3) 地震加速度谱计算功能,包括伪加速度谱、速度谱、位移谱
(4) 提供多达 59个地震动参数 及 Husid plot、Energy Flux plot
(5) 地震动参数批量分析功能
程序经历了多个版本的更新和维护,一直在持续改进:
GMP v2024 (当前版): [软件][地震动][更新] GMP v2024: A tool for Calculating Ground Motion Parameters for Seismic Analysis of Structures [结构抗震分析地震动强度指标/地震动参数计算工具]
- 软件界面 ( Screent Short)
- 地震动加速度格式支持 ( Supported Ground Motion Format)
目前支持的地震波导入格式:
a. PEER数据格式(NGA WEST, NGA WEST2, NGA EAST)
b. GML or GMS (from www.jdcui.com)
c. SPECTR (from www.jdcui.com)
d. YJK (盈建科)
e. SeismoSignal
d. 通用的单列格式
e. 双列地震波格式
- 地震动参数列表 ( Ground Motion Paramters )
目前GMP总共提供 59 个地震动参数的计算,其中主要的一些参数介绍如下:
(1)PGA,Peak ground values of acceleration , 峰值地震动加速度。
$$PGA = \max \left| {{{\ddot u}_g}\left( t \right)} \right|$$
(2)PGV,Peak ground values of velocity , 峰值地震动速度。
$$PGV{\rm{ = }}\max \left| {\mathop {\dot u}\nolimits_g (t)} \right|$$
(3)PGD,Peak ground values of displacement, 峰值地震动位移。
$$PGD = \max \left| {\mathop u\nolimits_g (t)} \right|$$
(4)FR1 = PGV /PGA , 地面速度峰值与加速度峰值比。
(5)FR2 = PGD /PGV , 地面位移峰值与速度峰值比。
(6)SI , Housner Intensity (HI), Housner谱强度。
$$SI\left( \zeta \right) = {I_H} = \int_{0.1}^{2.5} {{S_v}} \left( {T,\zeta } \right)dT = {\pi \over {2g}}\int_{0.1}^{2.5} {{S_a}} \left( {T,\zeta } \right)dT$$
(7)aRMS, Root-mean-square (RMS) of acceleration, 均方根加速度。
$${a_{RMS}} = \sqrt {{{\int_0^{{t_d}} {\ddot u_g^2dt} } \over {{t_d}}}} $$
(8)vRMS, Root-mean-square (RMS) of velocity, 均方根速度。
$${v_{RMS}} = \sqrt {{{\int_0^{{t_d}} {\dot u_g^2dt} } \over {{t_d}}}} $$
(9)dRMS, Root-mean-square (RMS) of displacement, 均方根加位移。
$${d_{RMS}} = \sqrt {{{\int_0^{{t_d}} {u_g^2dt} } \over {{t_d}}}} $$
(10)Td, Significant duration (TIA95% – TIA5%), The interval of time over which a proportion (percentage) of the total Arias Intensity is accumulated (default is the interval between the 5% and 95% thresholds). Arias强度分别占95%和5%的对应时刻之间的时间,目前使用最广泛的地震强震持时。
$${T_d} = t\left( {0.95{I_A}} \right) – t\left( {0.05{I_A}} \right)$$
(11)Td, Significant duration (TIA75% – TIA5%), The interval of time over which a proportion (percentage) of the total Arias Intensity is accumulated (default is the interval between the 5% and 75% thresholds). Arias强度分别占75%和5%的对应时刻之间的时间。
$${T_d} = t\left( {0.75{I_A}} \right) – t\left( {0.05{I_A}} \right)$$
(12)Td, Significant duration (TIA15% – TIA5%), The interval of time over which a proportion (percentage) of the total Arias Intensity is accumulated (default is the interval between the 5% and 15% thresholds). Arias强度分别占15%和5%的对应时刻之间的时间。
$${T_d} = t\left( {0.15{I_A}} \right) – t\left( {0.05{I_A}} \right)$$
(13)IA, Arias Intensity, Arias强度。
$$\mathop I\nolimits_A = {\pi \over {2g}}\int_0^{\mathop T\nolimits_d } {\mathop {\ddot u}\nolimits_g^2 } dt$$
(14)Characteristic Intensity (Ic),特征强度。
$${I_c} = A_{_{rms}}^{1.5}T_d^{0.5} = {\left( {\sqrt {{{\int_{\mathop t\nolimits_5 }^{\mathop t\nolimits_{95} } {\mathop {\ddot u}\nolimits_g^2 dt} } \over {\left( {\mathop t\nolimits_{95} – \mathop t\nolimits_5 } \right)}}} } \right)^{1.5}}{\left( {\mathop t\nolimits_{95} – \mathop t\nolimits_5 } \right)^{0.5}}$$
(15)A95 Parameter (g)。
For a given acceleration record and for a given acceleration level A, we can determine the area Ex bounded bythe curve a2(t) and the horizontal line at A2 level. Such areas are computed for various levels of acceleration. When the acceleration level A is the zero acceleration, the area Es, gives the Arias Intensity. The acceleration level which gives Ex/Es, equal to 0.05 is defined as the new parameter A95.
( Sarma S K, Yang K S. An evaluation of strong motion records and a new parameter A95[J]. Earthquake engineering & structural dynamics, 1987, 15(1): 119-132. )
(16)Specific Energy Density (SED),能量密度。
$$SED = \int_0^{{t_d}} {\dot u_g^2dt} $$
(17)Cumulative Absolute Velocity (CAV),累积绝对速度。
$$CAV{\rm{ = }}\int_0^{{t_d}} {\left| {{{\ddot u}_g}\left( t \right)} \right|} dt$$
(18)Cumulative Absolute Displacement (CAD),累积绝对位移。
$$CAD = \int_0^{{t_d}} {\left| {{{\dot u}_g}\left( t \right)} \right|} dt$$
(19)Cumulative Absolute Energy (CAI),累积绝对能量。
$$CAI = \int_0^{{t_d}} {\left| {{u_g}\left( t \right)} \right|} dt$$
(20)Ea, NAU和HALL指标。
$${E_a} = \int_0^{\mathop t\nolimits_d } {\ddot u_{\rm{g}}^2\left( t \right)} dt$$
(21)Ars, Sqrt(Ea) , NAU和HALL指标。
$$Ars = \sqrt {{E_a}} $$
(22)Ev, NAU和HALL指标。
$${E_v} = \int_0^{\mathop t\nolimits_d } {\dot u_{\rm{g}}^2\left( t \right)} dt$$
(23)Vrs, NAU和HALL指标。
$$Vrs = \sqrt {{E_v}} $$
(24)Ed, NAU和HALL指标。
$${E_d} = \int_0^{\mathop t\nolimits_d } {u_{\rm{g}}^2\left( t \right)} dt$$
(25)Drs, NAU和HALL指标。
$$Drs = \sqrt {{E_d}} $$
(26)Pa ,Housner 均方根加速度
$$\mathop P\nolimits_a = {{\int_{\mathop t\nolimits_5 }^{\mathop t\nolimits_{95} } {\mathop {\ddot u}\nolimits_g^2 dt} } \over {\mathop T\nolimits_d }} = {{\int_{\mathop t\nolimits_5 }^{\mathop t\nolimits_{95} } {\mathop {\ddot u}\nolimits_g^2 dt} } \over {\mathop t\nolimits_{95} – \mathop t\nolimits_5 }}$$
(27)Arms = sqrt(Pa ),Housner 均方根指标。
$$\mathop A\nolimits_{rms} = \sqrt {\mathop P\nolimits_a } $$
(28)Pv,Housner 均方根速度
$$\mathop P\nolimits_v = {{\int_{\mathop t\nolimits_5 }^{\mathop t\nolimits_{95} } {\mathop {\dot u}\nolimits_g^2 } dt} \over {\mathop T\nolimits_d }}$$
(29)Vrms = sqrt(Pv ),Housner 均方根指标。
$$\mathop V\nolimits_{rms} = \sqrt {\mathop P\nolimits_v } $$
(30)Pd,Housner 均方根位移。
$$\mathop P\nolimits_d = {{\int_{\mathop t\nolimits_5 }^{\mathop t\nolimits_{95} } {\mathop u\nolimits_g^2 dt} } \over {\mathop T\nolimits_d }}$$
(31)Drms = sqrt(Pd ),Housner 均方根指标。
$$\mathop D\nolimits_{rms} = \sqrt {\mathop P\nolimits_d } $$
(32)PSA,Peak acceleration,谱加速度峰值 。
$$PSA = \max \left( {{S_a}\left( {{T_{\rm{i}}}} \right)} \right)$$
(33)PSV,Peak velocity, 谱速度峰值 。
$$PSV = \max \left( {{S_v}\left( {{T_{\rm{i}}}} \right)} \right)$$
(34)PSD,Peak displacement,谱位移峰值
$$PSD = \max \left( {{S_d}\left( {{T_{\rm{i}}}} \right)} \right)$$
(35)SAavg,加速度谱均值。
$$S{A_{avg}} = \sum\limits_{{T_{\rm{i}}}}^{{T_{\rm{n}}}} {{{{S_a}\left( {{T_{\rm{i}}}} \right)} \over n}} $$
(36)SVavg,速度谱均值。
$$S{V_{avg}} = \sum\limits_{{T_{\rm{i}}}}^{{T_{\rm{n}}}} {{{{S_v}\left( {{T_{\rm{i}}}} \right)} \over n}} $$
(37)SDavg,位移谱均值。
$$S{D_{avg}} = \sum\limits_{{T_{\rm{i}}}}^{{T_{\rm{n}}}} {{{{S_{\rm{d}}}\left( {{T_{\rm{i}}}} \right)} \over n}} $$
(38)ASI,Acceleration Spectrum Intensity 加速度谱烈度。
$$ASI = \int_{0.1}^{0.5} {{S_{\rm{a}}}\left( {T,\zeta } \right)} dT$$
(39)VSI,Velocity Spectrum Intensity , 速度谱烈度。
$$VSI = \int_{0.1}^{2.0} {{S_v}\left( {T,\zeta } \right)} dT$$
(40)DSI,Displacement Spectrum Intensity , 位移谱烈度。
$$DSI = \int_{2.5}^{4.0} {{S_d}\left( {T,\zeta } \right)} dT$$
(41)EPA,Effective Peak Acceleration,有效峰值加速度
$$EPA = {{\int_{0.1}^{0.5} {{S_{\rm{a}}}\left( {T,\zeta } \right)} dT} \over {2.5}}$$
(42)EPV,Effective Peak Velocity ,有效峰值速度
$$EPV = {{\int_{0.8}^{2.0} {{S_v}\left( {T,\zeta } \right)} dT} \over {2.5}}$$
(43)EPD,Effective Peak Displacement ,有效峰值位移。
$$EPD = {{\int_{2.5}^{4.0} {{S_d}\left( {T,\zeta } \right)} dT} \over {2.5}}$$
(44)v0,The number of times the seismic acceleration curve passes through zero point per unit time. 加速度单位时间内通过零点次数
(45)IAM, Modified Arias Intensity, 修正Arias强度
$${I_{AM}} = {{{I_A}} \over {v_0^2}}$$
(46)Ia, Riddell强度指标
$${I_a} = {a_{\max }}t_d^{1/3}$$
(47)Id, Riddell强度指标
$${I_d} = {d_{\max }}t_d^{1/3}$$
(48)Iv, Riddell强度指标
$${I_v} = v_{\max }^{2/3}t_d^{1/3}$$
(49)If, Fajfar强度指标
$${I_f} = {v_{\max }}t_d^{1/4}$$
(50)Beta1 ( Ti>=2s and Ti<=10s)
$${\beta _1} = {{\sum {T_i^2\left( {{{{S_a}\left( {{T_i}} \right)} \over {PGA}}} \right)} } \over {\sum {T_i^2} }}$$
(51)PPV,PPV = max(vg) – min(vg) ,最大正负速度峰值差
(52)SMA, The third-largest absolute peak in acceleration
(53)SMV, The third-largest absolute peak in velocity
(54)SMD, The third-largest absolute peak in displacement
(55)Sa (Ti), 指定周期Ti下的 伪加速度谱
(56)Sv (Ti), 指定周期Ti下的 伪速度谱
(57)Sd (Ti), 指定周期Ti下的 位移谱
(58)MIV,Maximum Incremental Velocity,最大增量速度
(59)MID,Maximum Incremental Displacement,最大增量位移
- 批量分析功能介绍 ( Batch Analysis )
- 勾选需要分析的地震波
- 设置地震波是否需要极限修正,设置地震波的反应谱分析参数
- 点击批量分析,软件即可自动对勾选的地震波进行批量分析。
- 分析完毕,软件输出各组地震波的地震动参数计算结果及分析结果汇总。
- 下载 ( Download )
如果您需要这个软件,请在这给网站捐助:http://www.jdcui.com/?page_id=4813,捐助金额不低于300RMB,捐助后在评论区留下您的评论,或者联系管理员QQ:3014479529,给您配置软件,对软件有任何疑问也欢迎联系,敬请尊重劳动成果 !!!PS. 如果原软件功能无法满足您的需求,您可以和我们联系,可以针对您的需求,增加参数、做修改和调整。
- 相关博文( Related Topics)
[01]. [软件][工具] SPECTR: A program for Response Spectra Analysis [SPECTR地震波反应谱计算程序 v1.0]
[02]. [地震波][软件]GMS: Ground Motion Selection System [强震记录选取系统]
[03]. GML: Ground Motion Library Management System [强震记录管理系统]
[04]. [工具][软件][地震动] AEEG: A Program for Artificial Earthquake Accelerograms Generation [人工地震波合成软件]
[05]. ATC-63 Ground Motion Record Sets [ATC-63推荐的地震动记录集]
[06]. [软件][Tool] NGA Converter: PEER NGA Ground Motion Convertor [PEER地震波数据格式转换器]
[07]. Frequency Analysis of Ground Motion by MATLAB [MATLAB地震动频域分析]
[09]. [Tool][编程][软件] 水电工程水工建筑物标准设计反应谱计算程序(NB 35047-2015)
[10]. GB18306-2015地震动参数确定步骤与实例 (Ground motion parameter determination steps and examples of GB18306-2015)
[11]. 设计基本地震加速度与水平地震影响系数最大值
[12]. Tripartite Plot of Response Spectra [三联反应谱的基本概念与绘制]
[13]. Earthquake Response Spectra Analysis by MATLAB [MATLAB地震反应谱分析]
[14]. YJK地震波反应谱分析与地震波选取注意事项
[15]. 地震反应谱的特性(The Characteristics of Earthquake Response Spectrum)
[16]. Basics of Earthquake Response Spectra[地震反应谱的基本概念]
[17]. [下载][软件][规范]GB-SPECT: Chinese Code’s Design Response Spectrum[中国规范反应谱生成程序]
[18]. ETABS 2015人工波功能初步测试(一)
[19]. ETABS 2015人工波功能初步测试(二)
[21]. Dynamic Analysis of SDOF System by MATLAB [MATLAB单自由度系统动力响应分析]
[22]. 地震波基线是否修正对结构分析结果的影响(一)
[24]. [Software][科研小工具] FDEnergy: 力-位移累积耗能计算 [Force-Displacement Cumulative Dissipated Energy Calculator ]
[25]. [Tool][更新][Updates] SPECTR (v1.1) – A program for Response Spectra Analysis [SPECTR反应谱计算程序更新]
[27]. [编程][软件][PERFORM-3D] PERFORM-3D地震波时程逆向转换工具
[29]. [软件][地震动][Tool] 地震动参数计算软件GMP:补充介绍及说明
[30]. [地震工程][科研]地震动能量反应谱 (Earthquake Energy Demand Spectra)
[32]. [减隔震][地震工程][软件][Tool] 软钢·摩擦阻尼器减震结构的性能曲线计算
[33]. [地震工程][科研][软件][Tool] 等延性反应谱计算工具 (Construction of Constant-Ductility Response Spectrum)
[34]. [选波][地震波][Tool] GMS: Ground Motion Selection Program Updates (GMS选波软件更新)
[35]. [选波][地震波] 记录一个实际工程有趣的选波问题
[36]. [PEER][地震波] 关于PEER强震数据库地震波的单位及数据格式
[37]. [科研][更新][Update] 地震动参数计算软件(GMP)更新 (支持更多数据导入格式)
[38]. YJK1.7人工波功能测试
[40]. [数学][地震动][软件] FOUR_TRAN: Fourier Analysis Tool [傅里叶分析工具]
[44]. [软件][工具][地震波] KNETConvert: Kik-net Ground Motion Converter [K-Net, Kik-net 地震波数据转换工具]
[45]. [地震工程][科研][软件] IRSA 2020: Inelastic Response Spectra Analysis Program (弹塑性反应谱及单自由度非线性地震分析工具)
[49]. [地震][动力学] 对称结构的地震剪力规律 (Seismic Shear Law of Symmetrical Structures)
[50]. [软件][地震波][抗震] GMS_DESIGN: Ground Motion Selection Program for Practicing Engineers [地震波选波工具 工程师版]
[51]. [地震][结构] 双向地震作用效应,【先振型组合,再方向组合】及【先方向组合再方向组合】的差异?(实际案例测算)
[52]. [地震计算][反应谱][动力学][CQC] 振型叠加法随着组合振型数量的增加各种响应量是怎么变化的?
[53]. [抗震][动力学] 对于整体结构,X向地震作用下有Y向剪力吗?有!!
[54]. [地震波][选波][教程] 时程分析地震波选波介绍 ( Introduction of Seismic Wave Selection in Time History Analysis)
[55]. [抗震][减震][笔记] 黏弹性阻尼器的减震性能曲线 Response Reduction Curve of Viscoelastic (VE) System
[56]. [地震][动力学][Dynamics] 将阻尼矩阵的非对角线元素取为0计算结果会怎么样?
[57]. [工具][软件][规范] 广东省标准《高层建筑混凝土结构设计规范》反应谱计算工具
[58]. [资料] 1940 El Centro (Imperial Valley) Earthquake 地震波
[59] [动力学][Dynamics][SAP2000] SAP2000中振型向量的标准化方法
[61] [Dynamics][动力学][SAP2000] 梁的振动形态及振型质量 (Vibration Modes and Modal Mass of Beams)
[63] [地震动][动力学] 地震波的基线修正是什么意思?什么时候要基线修正?
[64] [SPECTR][SeismoSignal] SPECTR与SeismoSignal反应谱计算有差异?
[65] [动力学][舒适度] 增加结构的刚度可以降低的结构的加速度吗?
[67] [软件][动力学][Dynamics] NSDOF算例1——单自由度体系弹性动力时程分析
[68] [软件][动力学][Dynamics] NSDOF算例2——单自由度体系非线性动力时程分析
[69] [软件][动力学][Dynamics] NSDOF算例3——非线性粘滞阻尼单自由度体系动力时程分析
[70] [软件][动力学][Dynamics] NSDOF算例4——摩擦阻尼单自由度体系动力时程分析
[71] [软件][动力学][Dynamics] NSDOF算例5——非线性粘滞阻尼器+材料非线性 单自由度体系动力时程分析
[72] [软件][动力学][Dynamics] NSDOF算例6——设置黏弹性阻尼器单自由度体系动力时程分析
[73] [下载][软件]GB-SPECT V2021: 中国规范反应谱生成程序 [Chinese Code’s Design Response Spectrum]
[74] [结构][软件] TLDPC: 调谐液体阻尼器参数计算器 [TLDPC: Tuned Liquid Damper (TLD) Parameter Calculator]
[75] [下载][软件] 黏滞阻尼系数单位转换工具 [Viscous Damping Coefficient Unit Conversion Tool]
[76] [笔记][算例][减振] 调谐液体阻尼器减震算例 [Example of Vibration Reduction of Tuned Liquid Damper (TLD) ]
[77] [动力学][地震] 振型分解反应谱法构件地震力的计算过程?
[78] [下载][软件][地震工程] Spectr_Evolution: Seismic response spectrum evolution [地震波反应谱演化][反应谱随积分时间长短的变化]
[79] [软件][数学][地震动] FPSA: Fourier and Power Spectra Analysis [地震波频谱分析工具]
[80] [动力学][地震工程] 一个有趣的问题: SPECTR中的Newmark-Beta法计算反应谱发散?
[81] [结构设计][动力学] YJK中CQC振型组合地震力的复核
[82] [工具][软件][地震动] AEEG: A Program for Artificial Earthquake Accelerograms Generation [人工地震波合成软件]
[84] [动力学][地震动] SPECTR与SeismoSignal反应谱计算有差异?
[92] [软件][地震动] AEEG人工地震波合成 —— 案例1(拟合自定义反应谱)
[93] [软件][地震动] AEEG人工地震波合成 —— 案例2(如何提高人工波拟合精度)
[95]. [软件][地震工程][科研][更新] IRSA 2022: Inelastic Response Spectra Analysis Program (弹塑性反应谱及单自由度非线性地震分析工具)
[96]. [软件][地震工程] IRSA 案例1 —— 等延性反应谱分析算例 ( Constant Ductility Response Spectra Analysis Examples of IRSA)
[97]. [软件][地震工程] IRSA 案例2——单自由度非线性地震分析算例 ( Inelastic SDOF Earthquake Analysis Examples of IRSA)
[100]. [软件][地震工程] IRSA 案例5——地震波能量谱分析 ( Seismic Wave Energy Spectrum Analysis Examples of IRSA)
[101]. [软件][地震工程] IRSA 案例6——地震波基线修正 ( Seismic Wave Baseline Correction Examples of IRSA)
[102]. [地震工程][动力学] 地震作用下结构的能量分析 [Energy analysis of structures under earthquake]
[105]. [视频][Video] FPSA: Fourier and Power Spectra Analysis Batch Analysis Video [FPSA频谱分析工具批量分析录屏]
[106]. [程序][Web开发] 国标反应谱计算工具 网络版 [Chinese Design Code Response Spectrum Curve Generator —— Online version]
[107]. [地震动][研究][软件] RSFS: Response Spectra and Fourier Spectra [反应谱及傅里叶谱对比分析工具]
[108]. [软件][地震动] GMRS_Scale: Ground Motion Response Spectrum Scaling Tool [地震波&反应谱缩放工具]
[109]. [动力学][振动控制][编程] SDOF_FRE 案例1 —— 动力时程响应分析 [SDOF_FRE Example 1: Dynamic Force Time History Analysis]
[110]. [动力学][振动控制][软件] SDOF_FRE 案例 2 —— 地震时程响应分析 [SDOF_FRE Example 2: Earthquake Time History Analysis]
[112]. [软件][地震动][编程] GM_Truncation: A Program for truncating ground motion records [地震波截断工具]
[113]. [软件][编程][地震动] NFAGM: Near‐Field Artificial Ground Motions Generation Program [近场脉冲型地震动人工合成程序]
[114]. [软件][编程][地震动] NFAGM近场脉冲型地震动人工合成程序 案例1—— 单向脉冲波生成 (永久地面位移现象)
[115]. [软件][编程][地震动] NFAGM近场脉冲型地震动人工合成程序 案例2—— 正弦脉冲波生成与拟合
[116]. [编程][软件] PDDVA: Parameter design of dynamic vibration absorber [PDDVA: 动力吸振器参数设计软件]
[121]. [软件][地震动] AEEG人工地震波合成 —— 案例3(拟合自定义反应谱)
[122]. [软件][编程][研究] CASD: Calculate Average and Standard Deviation Curves [平均与标准偏差谱曲线计算工具]
[123]. [获奖][研究] “细腰型平面结构设计方法与关键技术研究” 获广东省工程勘察设计行业协会科学技术奖
[124]. [软件][编程][地震动] NFAGM近场脉冲型地震动人工合成程序 案例3—— 多峰值速度脉冲波
[125]. [动力学][结构] 大震弹塑性顶点位移时程为何”不收敛”?
[126]. [振动控制] 质量阻尼器分类 [Passive, semi-active, active and hybrid mass dampers]
[127]. [地震工程] 功率谱系列2(一): 设计反应谱及其功率谱 (Design response spectrum and its power spectrum)
[134]. [软件][地震工程] GMS_DESIGN选波软件——多频段选波案例 (Selection of earthquake ground motion using GMS_DESIGN)
结构|超限|软件开发|编程|参数化|优化|算法|振动控制|减隔震|有限元|技术培训
追求卓越 脚踏实地 致力于探索和拓展行业设计前沿
https://www.jdcui.com
合作及技术咨询
COOPERATION & CONTACT
E-mail:jidong_cui@163.com
WeChat & Tel: 13450468449