| 35 | 1/1 | 返回列表 |
| 查看: 5582 | 回復: 34 | |||||||
| 【獎勵】 本帖被評價30次,作者御劍江湖增加金幣 23.6 個 | |||||||
[資源]
站里關(guān)于MS Analysis計算材料的應力應變與Cubic以及Trigonal cell彈性常數(shù)的討論
|
|||||||
|
請教MS Analysis如何得到材料的應力應變曲線? 最近幾天在學習MS,看分析材料的模量,CASTEP和Discover都可以 我要計算的體系比較大,用CASTEP肯定不行,想用Discover 現(xiàn)在有個問題是 Discover里面,Analysis里面可以有Mechanical properties可以給出材料的模量,如果我想得到應力應變曲線,請問應該怎么做呢? 我看MS手冊里倒是說,計算模量的時候,逐步施加小應變,然后得出應力應變曲線,可以擬合曲線得到模量。 不知道具體怎么可以得到這個應力應變曲線 (一) i always calculate a set of point and then fit the Murnaghan equation. For example, the lattice constant is 10 A, then I take the 9.6, 9.7,9.8, 9.9, 10, 10.1,10.2,10,3,10.4, then get the total enery. using the volume-energy, it is quite easy to get the modulous that you want. (二) 一般是作NVT 再作NVE吧 但是也有用NPT的 這個東西可以去查查相關(guān)的文獻 比較多的 (三) 手動構(gòu)建原胞,可以直接加應變 (矩陣)給原胞,結(jié)構(gòu)優(yōu)化,從輸出文件中查最后的 Stress Tensor,找到相應的應力。 用discover如何計算應力應變關(guān)系? 用discover的dy namics可以計算應力應變關(guān)系,需要選擇系綜NPT或NPH。在輸出設(shè)置上選擇每500step輸出1frame,可是怎么運行了5000step還是得不到多張frame?得到文件xtd如何確定有多張frame呢?繼續(xù)用discover analysis卻不能進行應力應變分析,這是為何呢? (一)xtd一般是包含多張frames的 將動畫功能顯示出來,放映一下即可知道是不是多幀了 discover的analys 只能分析每個frame所對應狀態(tài)的力學常數(shù)。 (二)可能后面需要自己再處理 經(jīng)常在別人做分子動力學的文獻中看到有計算彈性常數(shù), 用自己計算的結(jié)果和別人的比較。 驗證自己的勢函數(shù)的好壞。我想問一下能用lammps計算彈性常數(shù)嗎?能給出具體的操作步驟嗎? 我只知道一點皮毛,彈性常數(shù)是通過應力和應變的關(guān)系得來,這個lammps中可以有哪些命令實現(xiàn)? (一)就應力應變曲線求導數(shù)即可,比如用作圖軟件Origin。 (二)在origin里,比如把A(X)列設(shè)為應變,B(Y)列設(shè)為應力,選中B(Y)列,用ananysis-calculus-differentiate,選擇一階,就可以生成一階導數(shù)了,生成的一階導數(shù)在C列~ BTW,diffenrentiate具體的路徑記不清了,好像是上面那個,又好像是ananysis-mathematics-differentiate。 (三)剛查了下,好像是這樣的: 一般說的彈性模量,至少我是這樣叫的,是指楊氏模量(Young's modulus)和剪切模量(Shear modulus) 而彈性常數(shù)包括5個: Lam's constant:拉梅系數(shù),有兩個 Shear modulus:剪切模量 Poisson's ratio:泊松比 Young's modulus:楊氏模量 Bulk modulus:體積模量 在各向同性材料中,上述只有2個是獨立的,所以只需求出兩個即可,我覺得容易求出來的就是楊氏模量和泊松比吧~ 英文定義見: * Young's Modulus, E: a proportionality constant between uniaxial stress and strain. * Poisson's Ratio, n: the ratio of radial strain to axial strain * Bulk modulus, K: a measure of the incompressibility of a body subjected to hydrostatic pressure. * Shear Modulus, G: also called rigidity, a measure of a substance's resistance to shear. * Lame's Constants, l and m: material constants that are derived from Young's Modulus and Poisson's Ratio. 它們之間的關(guān)系可參見:http://geophysics.ou.edu/geomech ... st/elast_const.html 我在2樓說用origin求導數(shù)只是求楊氏模量~ 初學者,上述不妥之處,還望指出~ Relationships Between Elastic Constants (from Fung, A First Course in Continuum Mechanics) Although only 2 independent elastic constants are required to completely the properties of an isotropic, elastic medium, there are many different constants from which to choose. Most of these constants arise as constants of proportionality between stress and strain for various loading conditions. For example, the ratio of hydrostatic stress (pressure) to volumetric strain (dilatation) is , which is given the name bulk modulus and the symbol . The constants referred to below are: • Has no name, as far as I know, so I call it "Lam?s constant" • shear modulus; sometimes G is used instead of • Poisson's ratio • Young's modulus • Bulk modulus To these we may add the following combinations that appear frequently. (四)http://www.mdbbs.org/viewthread. ... table&tid=19827 (五) Cubic以及Trigonal cell彈性常數(shù)的計算 Cubic和Trigonal晶體的彈性常數(shù)個數(shù)比較少,計算擬核也比較方便,目前計算彈性常數(shù)主要有兩個思路,一個是應力-應變曲線關(guān)系,此外就是應變能-應變關(guān)系,在兩種情況下彈性常數(shù)都是曲線的一階導數(shù)。在計算過程采用那種關(guān)系來擬核彈性常數(shù)沒有什么確定標準,目前實際上最大的限制是很多DFT計算軟件實際上不能把應變結(jié)構(gòu)應變能轉(zhuǎn)換成應力張量形式,如DMOL,Crystal,Wine2K等都不具備這個計算功能,CASTEP是少數(shù)直接可以輸出應力的軟件,因此在CASTEP中采用了Stress-Strain關(guān)系來擬和彈性常數(shù): 首先給出MS擬核彈性常數(shù)文件的格式: (六)http://cyh.xjtu.edu.cn/bbs/viewt ... &extra=page%3D1 想問下大家,做高壓下的彈性常數(shù),怎樣才能知道所加壓強超過物質(zhì)的范圍(也就是所加壓強讓物質(zhì)塑變了)?之前有人跟我說,通過應變應力曲線來判斷,這種方法行嗎?我的用的是CASTEP,如果要做這個曲線,在CASTEP中該怎樣設(shè)置呢?知道的能幫下忙嗎? (一)回答者:wanglianli136 回答時間:2009-10-19 3:30:58 2樓 在計算晶格優(yōu)化的task的more后面有一個加靜壓力的,還有你可以選擇的壓力方向,你可以看看 回答者:gavinliu7390 回答時間:2009-10-19 16:08:23 3樓 如果你壓強是一點一點的加,那么物質(zhì)是不會塑變的! 回答者:77585218358 回答時間:2009-10-20 2:04:30 4樓 同問一下,通過應變應力曲線判斷,是不是曲線存在一個屈服點?屈服點的出現(xiàn)代表塑變? 回答者:Huanhuan5360 回答時間:2009-10-20 4:38:19 5樓 你所說的一點 具體應該是多少呢? 回答者:jiangxue 回答時間:2009-10-21 0:07:17 6樓 這個應該是做不出來屈服點的,如果采用分子動力學的話應該可以得到這個屈服點 回答者:wuchenwf 回答時間:2009-10-21 1:16:16 7樓 單獨用MD的話,屈服點可能會給出來但是強度不是特別的準。 用第一性原理的話,屈服點有些問題,但是強度還是可以的。所以兩者各有利弊。 http://openshu.com/wenba/45243354896.html 高階彈性常數(shù)[轉(zhuǎn)] http://abinitio.yo2.cn/ 高階彈性常數(shù) 4分 • 收 藏分 享到頂編輯詞條 新 知社 • 人 人網(wǎng) • 開 心001 • 新 浪微博 • 網(wǎng) 易微博 目錄 • 1 導論(Introduction) • 2 1 廣義Hooker定理以及晶體的Neumann原理 • 3 2 應變自由能與彈性常數(shù)的關(guān)系 • 4 3高階彈性常數(shù)的獲得與計算 • 5 4參考文獻 • 展 開全部 • 1 導論(Introduction) 2 1 廣義Hooker定理以及晶體的Neumann原理 o 2.1 1.1 廣義Hooker定理 o 2.2 1.2 Neumann 原理 • 3 2 應變自由能與彈性常數(shù)的關(guān)系 4 3高階彈性常數(shù)的獲得與計算 o 4.1 3.1 實驗測量 o 4.2 3.2 理論計算 • 5 4參考文獻 • 收 起 摘要請用一段簡單的話描述該詞條,馬上添加摘要。 高階彈性常數(shù) 導論(Introduction) |
仿真建模與計算 | 科研相關(guān)資料 | 材料計算與模擬 | 微觀尺度 |
動力學模擬 |
| 35 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 302求調(diào)劑 +4 | 小賈同學123 2026-03-15 | 8/400 |
|
|
[考研] 289求調(diào)劑 +6 | 步川酷紫123 2026-03-11 | 6/300 |
|
|
[考研] 267一志愿南京工業(yè)大學0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
[文學芳草園] 伙伴們,祝我生日快樂吧 +17 | myrtle 2026-03-10 | 26/1300 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 070300化學學碩求調(diào)劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
[考研] 326求調(diào)劑 +4 | 上岸的小葡 2026-03-15 | 5/250 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 復試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy! 2026-03-13 | 6/300 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 復試調(diào)劑 +9 | Copy267 2026-03-10 | 9/450 |
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
[考研] 工科,求調(diào)劑 +3 | 我887 2026-03-11 | 3/150 |
|
|
[考研] 求調(diào)劑 資源與環(huán)境 285 +3 | 未名考生 2026-03-10 | 3/150 |
|
|
[考研] 321求調(diào)劑(食品/專碩) +3 | xc321 2026-03-12 | 6/300 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|