ansys解函數方程
關注創建者：王靖雯 創建時間：2023-03-07
ansys解函數方程的實例教程
如題，《從形函數與函數的連續可導性到ansys結果中的節點解與單元解的差異》，形函數對結果的影響大部分人都能聯想到二次單元比線性單元求得的結果更精確，但該文要表達的不僅如此，而是從更一般地討論怎么從單元的形函數來理解節點解與單元解之間的差異。
首先討論單元的階次。作為基礎我們應該明白網格與單元的區別，網格是將幾何體離散化后的結構，即組成幾何體的微元，單元是這些微元的幾何、物理或數學屬性（這里我們并不打算詳細討論單元的這些屬性，但是這些知識會方便對本文的理解）。我們經常在使用ansys或其他CAE軟件時經常會遇到單元的選擇以及單元階次的選擇，一般一種單元包括線性單元和二次單元甚至更高級的單元，比如在ansys中經常被使用的shell181（左）和shell281（右），線性單元使用的形函數是一次的多項式，高次單元使用的形函數是高次的多項式，形函數用于描述相鄰節點之間的位移場，所以高次的單元可以更好的描述形狀復雜的幾何體。
不同于常規材料力學中通過平衡方程求解（首先求得的解是力解），有限元方式求解的特點是首先求解出的結果是節點的位移解，即displacement of nodes，所有的節點位移形成了位移場，在空間上位移場一定是連續的，但是不一定是平滑的。哎哎，是不是特別熟悉的感覺，正是和高數中函數的連續性和可導性兩個性質非常相似，不用奇怪，位移場本來就是用函數描述的，所以自然就存在函數的性質，所以用函數的性質來理解就可以方便解釋一些現象了，下圖分別是用兩種形函數描述的位移場，在有限元求解后得到的首先是節點位移解，即圖中5個節點的位移，假如每個節點的位移用坐標x\y\z的函數來表示，然后通過形函數插值得到相鄰節點之間的位移（也是xyz的函數），上圖是用一次形函數插值，下圖是用二次形函數插值。
展開 由此推廣到一般情形，任意兩個函數卷積時，其積分表達式表示其中一個函數在積分點τ的值乘以另一個函數在t-τ的值。任意兩個函數f(t)*g(t)的卷積對應著一個系統，該系統在單位脈沖激勵下的響應為f(t)（或g(t)），而卷積f(t)*g(t)就表示該系統在g(t)（或f(t)）為激勵時下響應。
參考資料：
王新敏《ANSYS結構動力分析與應用》人民交通出版社，2014.
鄭君里《信號與系統上》第三版，高等教育出版社。
已知一階常微分方程，
g＋0.047*du/dt＋ u/6.7＝(13-u)/18
其中，當sin(10/pi*t)>=0時，g=1.18sin(10/pi*t)
當sin(10/pi*t)<0時， g=0
u的初值為u(0)=0，求t>=0時的解。畫出圖像，并求u的極大值。
曾經將g寫成，0.59 sin(10/pi*t)+abs(0.59*sin(10/pi*t)),用dsolve求解，可惜出錯了。
請大家幫忙，謝謝！（用其他函數解也可以。）
程序如下
fun=inline(['((13-u)/18-(sin(10*t/pi)>0)*',...
'1.18*sin(10*t/pi)-u/6.7)/0.047'],'t','u');
[t,u]=ode45(fun,[0,10],[0]);
plot(t,u)
說明g這樣表示的:
gt=(sin(10*t/pi)>0)*1.18*sin(10*t/pi);
感謝蘿卜網友
展開 *dim,a,,10
*dim,b,,10
*do,i,1,10
a(i)=i
b(i)=sin(i/5)
*enddo
/prep7
*do,i,1,10
k,i,a(i),b(i),0
*enddo
*do,i,1,9
l,i,i+1
ansys解函數方程的相關專題、標簽、搜索
ansys解函數方程的最新內容
個人學習總結，懇請指出錯誤。
參考資料見文后，文中引用格式為“作者+頁碼”、“作者名年份+頁碼”等。
-----前言-----
單位脈沖函數（Dirac函數）在一般的數學物理方法書籍中有詳細的介紹。對于該函數的工程應用，在自動控制原理中，可以通過一個系統對單位脈沖激勵的響應（脈沖響應）的表現，來判斷系統的時域穩定性等性質。但是直接求一個系統的脈沖響應不那么容易，往往借助拉普拉斯變換及其逆變換
如題，《從形函數與函數的連續可導性到ansys結果中的節點解與單元解的差異》，形函數對結果的影響大部分人都能聯想到二次單元比線性單元求得的結果更精確，但該文要表達的不僅如此，而是從更一般地討論怎么從單元的形函數來理解節點解與單元解之間的差異。
首先討論單元的階次。作為基礎我們應該明白網格與單元的區別，網格是將幾何體離散化后的結構，即組成幾何體的微元，單元是這些微元的幾何
*dim,a,,10
*dim,b,,10
*do,i,1,10
a(i)=i
b(i)=sin(i/5)
*enddo
/prep7
*do,i,1,10
k,i,a(i),b(i),0
*enddo
*do,i,1,9
l,i,i+1
已知一階常微分方程，
g＋0.047*du/dt＋ u/6.7＝(13-u)/18
其中，當sin(10/pi*t)>=0時，g=1.18sin(10/pi*t)
當sin(10/pi*t)<0時， g=0
u的初值為u(0)=0，求t>=0時的解。畫出圖像，并求u的極大值。
曾經將g寫成，0.59 sin(10/pi*t)+abs(0.59*sin(10/pi*t)),用dsolve
