[發明專利]一種快穩簡的方槽體外作用寬度評定方法

在審

• 申請號: 201711489099.6
• 申請日: 2017-12-30
• 公開（公告）號: CN108121876A
• 公開（公告）日: 2018-06-05
• 發明（設計）人: 唐哲敏;黃美發
• 申請（專利權）人: 唐哲敏
• 主分類號: G06F17/50
• 分類號: G06F17/50
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種快穩簡的方槽體外作用寬度評定方法，其特征在于，由以下步驟組成：

步驟1：獲取測點集{p_i}，并根據{p_i}建立特征行向量集{A_i}、邊界元素集{b_i}和狀態元素集{t_i}，其中：

i=1, 2, 3, …, N；i為測點序號，N為測點總數；

p_i={x_i, y_i , z_i}是測點i的平面直角坐標，并且被測方槽的中心平面接近坐標系的XOY平面，測點在坐標系原點附近；

t_i=| z_i |，所有的狀態元素t_i的集合為狀態元素集{t_i}；

A_i=([1, y_i, -x_i])z_i/t_i，是一個特征行向量，所有的特征行向量A_i的集合為特征行向量集{A_i}；

b_i=b，是一個大于0的實數，所有的邊界元素b_i的集合為邊界元素集{b_i}；

步驟1結束后進行步驟2；

步驟2：取t_i最小值t_min對應的測點p_l1為關鍵點，并將其測點序號l₁加入到關鍵點集{l}中；

步驟2結束后進行步驟3；

步驟3：根據關鍵點集{l}建立分析矩陣A和分析列向量b，其中：

A=[…, A_j^T, …, A_k^T, …]^T，是個L行3列的矩陣，L為關鍵點集{l}中的元素個數，j, k為關鍵點集{l}中的元素；

b=[…, b, …]^T，是個L行的列向量；

步驟3結束后進行步驟4；

步驟4：對分析矩陣A及增廣分析矩陣[A, b]進行秩分析；

計算r_A=rank(A)，r_Ab=rank（[A, b]），并比較r_A和r_Ab，只有以下兩種情況：

情況一：如果r_A=r_Ab，那么，應當繼續尋優，跳到步驟5；

情況二：如果r_A< r_Ab，那么，嘗試從分析矩陣A和分析列向量b中刪掉關鍵點集{l}中的某一個元素l對應的行，得到縮小矩陣A_l-和縮小列向量b_l-，求線性方程A_l-v_l-= b_l-的解v_l-=v_l-0，然后計算b_l-=A_lv_l-0；如果關鍵點集{l}中的元素都嘗試過了，并且沒有得到任何一個b_l->b，那么，應當結束尋優，跳到步驟7；如果在嘗試關鍵點集{l}中的元素l時，得到b_l->b，那么，將縮小矩陣A_l-和縮小列向量b_l-分別作為A矩陣及分析列向量b，將元素l移出關鍵點集{l}，并跳到步驟5；其中，v_l-=[v_l-,1, v_l-,2, v_l-,3]^T，v_l-0=[v_l-0,1, v_l-0,2, v_l-0,3]^T；

步驟5：求線性方程Av= b的解v=v₀，其中，v=[v₁, v₂, v₃]^T，v₀=[v_0,1, v_0,2, v_0,3]^T；

步驟5結束后進行步驟6；

步驟6：計算v_i=A_iv₀，然后計算τ_i=(t_i– t_min)÷(b - v_i)；取τ_i中大于零的那部分中的最小值τ_min對應的測點p_l2為新的關鍵點，并將其測點序號l₂加入到關鍵點集{l}中；

將所有t_i更新為t_i+ τ_min? v_i，t_min更新為t_i的最小值；

步驟6結束后完成一次尋優，進行步驟3；

步驟7：計算t=2 t_min就是所求的方槽體外作用寬度。