參數(shù)取值的案例
abaqus混凝土損傷模型參數(shù)取值
誰能把a(bǔ)baqus混凝土損傷模型參數(shù)取值這個問題說清楚呢？
LS-DYNA 土壤MAT147 土壤本構(gòu)參數(shù)取值范圍
https://www.fhwa.dot.gov/publications/research/safety/04094/04.cfm
單位制 kg-mm-ms-Gpa
基于ansys的梁單元、實體單元徐變精細(xì)化分析（含各參數(shù)解釋） ￥25
徐變應(yīng)變可表達(dá)為：
其中， ?(t,τ)為徐變系數(shù)，需通過規(guī)范公式或?qū)嶒灁?shù)據(jù)擬合確定
Ansys程序中內(nèi)置金屬蠕變規(guī)律如下：
命令中詳細(xì)解釋了改公式的具體用法，以及參數(shù)意義。
二者除個別參數(shù)外形式具有異曲同工之妙，因此本案例給出用ansys精確分析混凝土徐變的方法，案例背景模擬了一個混凝土PK梁特定工況下的徐變發(fā)生過程。
案例文件中包含：
1. 00-ConcreteCreep-benchmark.mac【徐變標(biāo)定文件，開箱即用，可以用來和手算對比是否正確】
2. 01-ConcreteCreep-solid.mac【分輸入模塊的參數(shù)化徐變計算文件【詳細(xì)解釋了各參數(shù)取值】。只需要改文件和計算邊界荷載即可計算實體徐變。】
3. ansa文件，用來生成網(wǎng)格
4. .cdb文件，網(wǎng)格文件
5. excel轉(zhuǎn)apdl命令流文件，用來輸入徐變系數(shù)。
進(jìn)一步白話闡述一下：
1、什么是徐變？別看公式一大堆，理論一大推，簡單講就是：受力的結(jié)構(gòu)，啥邊界條件、荷載不變的情況下，結(jié)構(gòu)還是慢慢變形了。將這種慢慢變形的變形結(jié)果以及應(yīng)力重分配準(zhǔn)確分析出來就是徐變分析。機(jī)理一大堆，教科書上都比較詳盡，在此不做贅述，只講應(yīng)用，而且是拿到案例開箱即用。
白話闡述要點(diǎn)：
1、案例是ansys apdl（命令流）分析的，給出了全套參數(shù)化命令流，材料模型定義、材料參數(shù)定義、求解，拿過來可以直接運(yùn)行。
2、機(jī)理是用了ansys中關(guān)于金屬蠕變的材料模型。（細(xì)想蠕變和徐變的現(xiàn)象，表征都是一樣的。至于機(jī)理，各有各的理論，但不影響材料模型使用。）
具體使用：
1、，先跑一遍，看看到底徐變是怎么個事兒。
展開 OpenSEES減隔震常用單元本構(gòu)命令介紹
課程介紹了減隔震分析中常用四種本構(gòu)單元的參數(shù)取值和使用注意事項：
多線性彈性本構(gòu)ElasticMultiLinear
滯回材料本構(gòu)Hysteretic
摩擦擺支座單元singleFPBearing
粘性阻尼材料本構(gòu)Viscous和ViscousDamper
附件里有4種材料案例的完整命令流和課程說明文檔。
https://www.yqgqt.org.cn/college/video/c194583?SUID=4aa4997d-e3cf-4959-8863-1162548b72f2&n=OpenSEES%E5%87%8F&c=%E6%9C%AC%E8%AF%BE%E7%A8%8B%E4%BB%8B%E7%BB%8D%E4%BA%86%E5%87%8F%E9%9A%94%E9%9C%87
展開 
新課程：精細(xì)化軌道-橋梁耦合振動模型建模與分析
主要知識點(diǎn)：
橫向節(jié)點(diǎn)數(shù)目不匹配的兩種處理方式
矩陣奇異原因：約束不足
MinMax材料本構(gòu)
單向受壓材料ENT
TwoNodeLink單元
Concrete02材料本構(gòu)參數(shù)取值
Steel02材料本構(gòu)
理想彈塑性本構(gòu)ElasticPP
PS：由于本課程介紹的模型為本人碩士畢業(yè)論文中所用案例，后續(xù)可能用于發(fā)表文章，故不提供完整命令流和Word文檔，僅提供涉及知識點(diǎn)的代碼，介意勿拍。
基于磁路法與等效熱網(wǎng)絡(luò)法的航天永磁同步電機(jī)設(shè)計與仿真
給出了定子內(nèi)徑、定子外徑、鐵芯長度、匝數(shù)等關(guān)鍵參數(shù)的取值準(zhǔn)則，建立了包含36個節(jié)點(diǎn)集總參數(shù)熱網(wǎng)絡(luò)模型，并以端部繞組為例給出了熱平衡方程的詳細(xì)推導(dǎo)過程。通過與成熟商業(yè)軟件對比，其電磁計算最大誤差出現(xiàn)在電流有效值上，偏差值為6.07%；與樣機(jī)實測值對比，繞組溫升最大誤差為7.3%，滿足方案設(shè)計階段預(yù)示精度要求，為方案設(shè)計階段航天永磁同步電機(jī)快速性能預(yù)估提供有力支撐。
引言
機(jī)電伺服系統(tǒng)作為伺服機(jī)構(gòu)的重要一員，越來越多地應(yīng)用到航天領(lǐng)域中。航天機(jī)電伺服系統(tǒng)具有短時高功率、長時低功率、制動負(fù)功率的特性，由于永磁同步電機(jī)(Permanent Magnet Synchronous Motor, PMSM)高效率、高功率因數(shù)和高功率密度的特點(diǎn)受到航空航天領(lǐng)域研究者的廣泛關(guān)注。電機(jī)作為航天機(jī)電伺服系統(tǒng)的能量轉(zhuǎn)換裝置，其電磁、磁力轉(zhuǎn)化效率及熱損耗估算的快速性和準(zhǔn)確性是產(chǎn)品研制的關(guān)鍵。現(xiàn)有電機(jī)電磁設(shè)計方法包括磁路法、解析法和有限元法等；熱設(shè)計方法包括熱路法、等效熱網(wǎng)絡(luò)法、有限元法和流體力學(xué)方法等。解析法在工程上常無法獲得精確解；有限元法依賴于詳細(xì)的幾何參數(shù)；在信息較少的方案設(shè)計階段均無法使用；磁路法和等效熱網(wǎng)絡(luò)法原理清晰、便于理解，常常用于性能的初步預(yù)估，但在航天領(lǐng)域，尚無文獻(xiàn)給出磁路法與等效熱網(wǎng)絡(luò)法用于電磁熱分析的詳細(xì)流程及各關(guān)鍵參數(shù)的取值準(zhǔn)則。
基于此,本文提出了基于磁路法和等效熱網(wǎng)絡(luò)法的電機(jī)磁熱快速設(shè)計仿真方法。以表貼式永磁同步電機(jī)為例，給出了基于磁路法進(jìn)行PMSM電磁設(shè)計關(guān)鍵參數(shù)的取值標(biāo)準(zhǔn)；同時建立了包含36個節(jié)點(diǎn)的PMSM集總參數(shù)熱網(wǎng)絡(luò)分析模型，并選取其中關(guān)鍵節(jié)點(diǎn)推導(dǎo)了熱平衡方程，最后應(yīng)用該方法對某航天電機(jī)產(chǎn)品進(jìn)行了磁熱仿真分析，并分別與商軟計算數(shù)據(jù)和實物試驗數(shù)據(jù)對比。
展開 分類預(yù)測 | MATLAB實現(xiàn)WOA-CNN-LSTM-Attention數(shù)據(jù)分類預(yù)測
SearchAgents_no = 8; % 數(shù)量Max_iteration = 5; % 最大迭代次數(shù)dim = 3; % 優(yōu)化參數(shù)個數(shù)lb = [1e-3,10 1e-4]; % 參數(shù)取值下界(學(xué)習(xí)率，隱藏層節(jié)點(diǎn)，正則化系數(shù))ub = [1e-2, 30,1e-1]; % 參數(shù)取值上界(學(xué)習(xí)率，隱藏層節(jié)點(diǎn)，正則化系數(shù))
fitness = @(x)fical(x,num_dim,num_class,p_train,t_train,T_train);
[Best_score,Best_pos,curve]=WOA(SearchAgents_no,Max_iteration,lb ,ub,dim,fitness)Best_pos(1, 2) = round(Best_pos(1, 2)); best_hd = Best_pos(1, 2); % 最佳隱藏層節(jié)點(diǎn)數(shù)best_lr= Best_pos(1, 1);% 最佳初始學(xué)習(xí)率best_l2 = Best_pos(1, 3);% 最佳L2正則化系數(shù) %% 建立模型lgraph = layerGraph(); % 建立空白網(wǎng)絡(luò)結(jié)構(gòu)tempLayers = [ sequenceInputLayer([num_dim, 1, 1], "Name", "sequence") % 建立輸入層，輸入數(shù)據(jù)結(jié)構(gòu)為[num_dim, 1, 1] sequenceFoldingLayer("Name", "seqfold
展開 深基坑設(shè)計計算時涉及8個問題的探討
將目錄列出，可根據(jù)需要閱讀自己感興趣的部分~
一、土層強(qiáng)度參數(shù)取值與水土分算與合算；二、土壓力分布模式；三、超載計算；四、被動區(qū)留土與被動區(qū)加固；五、樁、墻嵌固深度；六、鉆工墻底部抗隆起穩(wěn)定性驗算；七、樁、墻剛度界定；八、分階支護(hù)問題
（點(diǎn)擊圖片可查看高清大圖哦~）
【來源：筑龍巖土】
基于灰狼算法優(yōu)化支持向量機(jī)的matlab算法
c和g，number of your variables
lb=[0.01,0.01]; % 參數(shù)取值下界
ub=[100,100]; % 參數(shù)取值上界
% v = 5; % SVM Cross Validation參數(shù),默認(rèn)為5
Convergence_curve=zeros(1,Max_iteration);
l=0; % Loop counter循環(huán)計數(shù)器
% Main loop主循環(huán)
while l<Max_iteration % 對迭代次數(shù)循環(huán)
for i=1:size(Positions,1) % 遍歷每個狼
if fitness>Alpha_score && fitness<Beta_score % 如果目標(biāo)函數(shù)值介于于Alpha狼和Beta狼的目標(biāo)函數(shù)值之間
Beta_score=fitness; % 則將Beta狼的目標(biāo)函數(shù)值更新為最優(yōu)目標(biāo)函數(shù)值，Update beta
Beta_pos=Positions(i,:); % 同時更新Beta狼的位置
end
if fitness>Alpha_score && fitness>Beta_score && fitness<Delta_score % 如果目標(biāo)函數(shù)值介于于Beta狼和Delta狼的目標(biāo)函數(shù)值之間
Delta_score=fitness; % 則將Delta狼的目標(biāo)函數(shù)值更新為最優(yōu)目標(biāo)函數(shù)值，Update delta
Delta_pos=Positions(i,:); % 同時更新Delta狼的位置
end
end
展開 Amesim電池教程 電池?zé)崾Э啬Ｐ?在線閱讀
定義其反應(yīng)速率為RSEI，則有：
由該反應(yīng)產(chǎn)生的放熱量為：
上述反應(yīng)公式中的參數(shù)取值參考值如下：
3.1.2. 負(fù)極發(fā)生與電解液溶劑反應(yīng)
電解液溶劑的消耗速率與負(fù)極上穩(wěn)態(tài)SEI生成的速率一致，此消彼長。其中，用ZSEI表示穩(wěn)態(tài)SEI所占比重，xneg表示負(fù)極電解液溶劑消耗量所占比重。以-Rneg表示反應(yīng)速率。
這里通過項來表示反應(yīng)速率受到擴(kuò)散濃度的制約，zSEI越大，則反應(yīng)越慢。鋰離子的傳輸阻力與穩(wěn)態(tài)SEI的厚度成指數(shù)相關(guān)。
因此，該部分的反應(yīng)放熱率為：
上述反應(yīng)公式中的參數(shù)取值參考值如下：
3.1.3. 正極發(fā)生與電解液溶劑反應(yīng)
把正極與電解液溶劑反應(yīng)的轉(zhuǎn)化因子定義為xpos，則它的導(dǎo)數(shù)即為反應(yīng)速率。
因此，該部分的反應(yīng)放熱率為：
由于正極反應(yīng)導(dǎo)致的電池容量損失：
用比例因子FQres來表示剩余容量與初始容量的關(guān)系，該因子的計算方式如下：
參考取值如下：
xpos在上表中的值為電池充滿狀態(tài)下的取值，它與放熱反應(yīng)開始時的SOC相關(guān)，用戶可以調(diào)整。其中為wvpos為密度，因此代入到放熱率計算公式還需乘以單體體積。
3.1.4. 電解液溶質(zhì)分解反應(yīng)
高溫狀態(tài)下，電解液不再穩(wěn)定，電解液溶質(zhì)（例如：LiPF6）將與空氣發(fā)生氧化反應(yīng)。以無量綱數(shù) θe表示參加反應(yīng)的電解質(zhì)。則反應(yīng)速率為（此處xe應(yīng)為θe）：
因此，該部分的反應(yīng)放熱率為：
電解質(zhì)LiPF6對應(yīng)的參考取值如下，其他類型的數(shù)值需用戶自行修改：
3.1.5. 自放電（隔膜熔化）
高溫下，隔膜熔融后電池內(nèi)部將出現(xiàn)大量短路，這將直接導(dǎo)致鋰電池端電壓輸出下降到0V。
展開 修正劍橋模型的應(yīng)用
參數(shù)取值
劍橋模型涉及的幾個參數(shù)，比較難取的無非就是κ、λ、a0。具體到上海軟土，可以參考上海華東院的參考取值（來自于網(wǎng)絡(luò)，見附件）。在沒有實驗條件的情況下，根據(jù)地勘報告，可分別將Cc、Cs除以ln10來求λ、κ，具體原因就不用說了吧。對于M，可以根據(jù)有效內(nèi)摩擦角按公式求出。對于a0，在abaqus里可以通過給定e1來設(shè)定，如果已給定e1，a0一欄默認(rèn)即可。關(guān)于e1的取值，我是通過地勘報告里的各地層高壓固結(jié)曲線來求的，具體做法是根據(jù)最后一個（e，lnp）坐標(biāo)值以及正常固結(jié)曲線斜率λ反算出截距便是e1。
2. 初始條件
包括初始地應(yīng)力、孔隙比、孔壓。
初始地應(yīng)力：
*initial conditions, type=stress, geostatic
ele-soil-1, -89689.6, 124.8, 0, 130, 0.47, 0.47
各數(shù)據(jù)含義就不用廢話了。注意前面要用單元集（element set）注：我這里沒有考慮固結(jié)，沒有采用有效應(yīng)力
abaqus里默認(rèn)的地應(yīng)力方向是Z軸負(fù)方向，如果是在CAE中建模且開始將地應(yīng)力方向設(shè)為了Y軸負(fù)方向，可以在Assembly中將整個模型旋轉(zhuǎn)90°，如果彈出“...converted to absolute positions....”，yes就是了。
注意，劍橋模型中不支持在無應(yīng)力條件下施加重力進(jìn)行地應(yīng)力平衡的方法。
孔隙比：
*initial conditions, type=ratio
soil-1, 1.112
這個是假設(shè)孔隙比與高程無關(guān)。
展開 
纏繞條帶加固柱體
具體為Part模塊中通過Create Part選擇三維空間，基于shell通過旋轉(zhuǎn)建模（圖2），然后在圖3所示的繪圖空間中繪制條帶的起始邊，確認(rèn)后調(diào)出土4所示的對話框，參數(shù)含義如圖所示，可隨意調(diào)整查看效果。
圖1
圖2
圖3
圖4
條帶單元采用膜單元，所以在創(chuàng)建Section時，要選擇shell下的membrane（圖5）。
圖5
柱體可再創(chuàng)建一個實體part，條帶與柱之間通過設(shè)置接觸，本例中采用Cohesive Behavior模擬所可能的粘結(jié)（膠水）作用，也可以設(shè)置一個初始拉應(yīng)力模擬實際情況。
a 有條帶
b無條帶
圖6對比了有無條帶，柱頂軸向加壓后的豎向變形，在條帶約束下，側(cè)向變形受限，軸向變形相應(yīng)減小。需要注意，本例中的參數(shù)取值較為隨意，不反映實際情況，只是提供一種可能的分析思路。
作者： 費(fèi)康
來源ABAQUS在巖土工程中的應(yīng)用
展開 Abaqus子程序HETVAL模擬混凝土水化熱溫度場
Q0為混凝土水化放熱的最大值，取為364000J/kg（一般試驗或文獻(xiàn)里給的都是kJ/kg這種單位，但Abaqus里的標(biāo)準(zhǔn)單位為J，這里要統(tǒng)一單位），而b、c這兩個參數(shù)取值可在文獻(xiàn)中查到，不同水泥類型的參數(shù)取值不一樣，對于52.5普通硅酸鹽水泥，這里b取為0.36，c為0.74。注意到式1中的t的單位為天。
式1：
然后對式1求導(dǎo)，可得到單位立方米混凝土單位時間的水化放熱率，qv (J/m^3/hr)，見式2。這里將上式中的時間單位改為了hr，因為一般混凝土溫度場監(jiān)測以小時為單位進(jìn)行記錄。這里的qv就是HETVAL中的FLUX(1)。具體HETVAL子程序如下：
式2：
在該子程序中最重要為熱源放熱率FLUX(1)，并且及時用statev結(jié)果狀態(tài)變量來存儲，這樣方便在后處理中輸出查看，我稱之為實時調(diào)試，也就是可以實時監(jiān)控自己定義或想要輸出變量在有限元計算中隨時間的變化規(guī)律。這里statev(1)記錄的是熱源放熱率flux(1)（時間單位為小時），statev(2)記錄的是時變的水化放熱總量（注意到時間單位也換為了小時，因為模型中的混凝土熱工性能參數(shù)，比如熱傳導(dǎo)系數(shù)用的時間單位也是小時，這里統(tǒng)一對應(yīng)起來），而statev(3)、statev(4)和statev(5)分別為Q0，b，c，這樣可以驗證下子程序的正確性，看下輸出這些確定值是否在計算過程中發(fā)生了變化。
另外，特別需要注意的是，下列程序中的Cemass這個變量一定要出現(xiàn)在FLUX(1)的計算表達(dá)式中，否則計算出來的溫度變化極其小。
展開 降溫法+追蹤法模擬巖土開挖與支護(hù)中的shell面殼單元與實體solid element合理連接
如果采用接觸對定義接觸，個人感覺很難被工程接受，原因是接觸參數(shù)取值的合理性和可信性，接觸狀態(tài)的合理性等。
提出以上開挖支護(hù)中遇到的一個問題，希望大家發(fā)表自己的心得。
首先感謝Simwe論壇，也Robert_Su的熱心關(guān)注與指導(dǎo)。
通過幾天的摸索嘗試，應(yīng)力集中問題已經(jīng)基本搞定：采用surface to surface 綁定約束，面殼網(wǎng)格與其支護(hù)的實體網(wǎng)格密度完全相同，非共節(jié)點(diǎn)綁定約束的響應(yīng)和共節(jié)點(diǎn)模型的響應(yīng)幾乎沒有差別。
現(xiàn)將采用降溫法+追蹤法模擬洞室開挖和支護(hù)的共節(jié)點(diǎn)與非共節(jié)點(diǎn)小模型的inp上傳，希望對巖土開挖與支護(hù)的朋友有用。簡單介紹下開挖實現(xiàn)的過程：前提---無論是共節(jié)點(diǎn)還是非共節(jié)點(diǎn)模型，面殼shell單元都要建立兩組相同節(jié)點(diǎn)不同單元號的襯砌shell單元組，其中一組給定很小的材料參數(shù)，用來追蹤降溫過中襯砌-圍巖交界面處的幾何位置（參考幫助或本論壇相關(guān)帖子）。（1）地應(yīng)力平衡，這應(yīng)該沒有什么說的，并殺死將來用來支護(hù)的shell單元，保留用來追蹤幾何位置的shell單元。（2）利用溫度或產(chǎn)變量對開挖巖體進(jìn)行降溫，以達(dá)到應(yīng)力釋放的目的。（3）殺死開挖的巖體和追蹤的shell單元組，并激活支護(hù)的shell單元組，計算平衡后，開挖過程完畢。
**如果開挖洞室比較長時，如200m，模型會出現(xiàn)收斂問題，其原因是shell面殼單元rotation自由度過大，調(diào)節(jié)收斂參數(shù)也不好用，尤其是tie綁定非共節(jié)點(diǎn)情況，不知道大家是否遇到，有什么好的解決方法。
excave.rar
展開 考慮企業(yè)異質(zhì)性的外部知識對流程創(chuàng)新作用機(jī)制的博弈分析
并用mathematica中用提取系數(shù)的函數(shù)將求偏導(dǎo)結(jié)果整理，得到：
其中，參數(shù)為：
解得最優(yōu)引入外部知識量分別為：
代入回式(5)(6)，得
結(jié)果分析
最后進(jìn)行數(shù)值分析，驗證本篇論文復(fù)現(xiàn)的正確性。Mathematica即可進(jìn)行符號計算公式推導(dǎo)，也可進(jìn)行數(shù)值分析。前文在mathematica進(jìn)行公式推導(dǎo)的結(jié)果可直接應(yīng)用于后續(xù)的數(shù)值分析，避免了切換軟件造成的麻煩。參數(shù)取值參考[1]。
企業(yè)異質(zhì)性體現(xiàn)在初始邊際生產(chǎn)成本差異系數(shù)ε，初始邊際生產(chǎn)成本差異系數(shù)對均衡外部知識引入量和均衡產(chǎn)量影響的數(shù)值模擬結(jié)果如圖1(上圖為復(fù)現(xiàn)結(jié)果，右圖為論文原文結(jié)果)。
圖1 ε 與均衡引入外部知識量的關(guān)系
圖2 k1與均衡利潤的關(guān)系
參考文獻(xiàn)
[1] 考慮企業(yè)異質(zhì)性的外部知識對流程創(chuàng)新作用機(jī)制的博弈分析[J]. 中國管理科學(xué), 2023, 31(12): 175-184.
最后，有相關(guān)需求歡迎通過公眾號“320科技工作室”與我們聯(lián)絡(luò)
展開 