Abaqus材料參數(shù)設(shè)置之密度動(dòng)態(tài)定義

2017-05-24 by:CAE仿真在線來(lái)源:互聯(lián)網(wǎng)

問(wèn)題:有時(shí)我們需要在材料的不同區(qū)域定義不同的密度,比如當(dāng)巖體采用了Gravity(重力)作為荷載施加類(lèi)型時(shí),不同的深度的巖體密度會(huì)有所不同。這時(shí)就不能簡(jiǎn)單將密度定義為一常數(shù)。

一、CAE當(dāng)中先前設(shè)置

路徑:在材料設(shè)置窗口選擇【General/Density】

1、Number of field variables:表示和密度相關(guān)的變量個(gè)數(shù)。如果是1個(gè),比如我們提出的問(wèn)題只與豎向位置有關(guān),則這個(gè)值就設(shè)置為1;如果有多個(gè)相關(guān),那有幾個(gè)值相關(guān)就設(shè)置幾個(gè)值。

2、Mass Density:密度值

3、Field1:對(duì)應(yīng)密度的第1個(gè)變量的值

二、后期設(shè)置

這里主要分兩種情況將密度具體設(shè)置到指定節(jié)點(diǎn)。

1、inp文件設(shè)置

如果數(shù)據(jù)量不是很大,且操作簡(jiǎn)單,可以使用這種方式,如下進(jìn)行定義

格式:*Field, variable=n

Node, N1, N2,……Nn

n:表示變量個(gè)數(shù)

Node:節(jié)點(diǎn)號(hào)或集合

N1:對(duì)應(yīng)的第1個(gè)變量值

N2:對(duì)應(yīng)的第2個(gè)變量值

Nn:對(duì)應(yīng)的第n個(gè)變量值

2、用戶子程序自定義

在使用這種方式時(shí),依然需要先修改inp文件

格式:*Field, variable=n, User

Node

n:表示變量個(gè)數(shù)

User:表示提示程序?qū)⑹褂糜脩糇远x程序

Node:節(jié)點(diǎn)號(hào)或集合

這里用到的用戶子程序?yàn)閁Field,其格式如下,其中函數(shù)里傳遞變量的具體含義可以查閱Abaqus幫助文檔.

SubRoutineUField(Field,KField,Nsecpt,KStep,Kinc,Time,Node,Coords,Temp,DTemp,NField)

Include'ABA_PARAM.inc'

DimensionField(Nsecpt,NField),Time(2),Coords(3),Temp(Nsecpt),DTemp(Nsecpt)

Field(1,1)=Coords(2)!如果這里表示設(shè)置第1個(gè)變量的值為此節(jié)點(diǎn)2方向(深度方向)的值(當(dāng)然如果第1個(gè)變量的值是節(jié)點(diǎn)1方向與2方向共同作用,類(lèi)似這種定義 Field(1,1)=(Coords(1)-Coords(2))*3.1.如果有多個(gè)變量,可分別在這里進(jìn)行賦值

Return

End SubRoutine

注意:

1、Abaqus的密度是指的干密度,即不包含水分時(shí)的密度

2、Field variable必須在*Step語(yǔ)句之后

參考:

[1] 費(fèi)康,張建偉. Abaqus 在巖土工程中的應(yīng)用[M]. 中國(guó)水利水電出版社, 2013:186-187

開(kāi)放分享：優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才