ansys中函數如何編寫的案例
提高代碼質量：如何編寫函數
函數是實現程序功能的最基本單位，每一個程序都是由一個個最基本的函數構成的。寫好一個函數是提高程序代碼質量最關鍵的一步。本文就函數的編寫，從函數命名，代碼分布，技巧等方面入手，談談如何寫好一個可讀性高、易維護，易測試的函數。
命名
首先從命名說起，命名是提高可讀性的第一步。如何為變量和函數命名一直是開發者心中的痛點之一，對于母語非英語的我們來說，更是難上加難。下面我來說說如何為函數命名的一些想法和感受：
采用統一的命名規則
在談及如何為函數取一個準確而優雅的名字之前，首先最重要的是要有統一的命名規則。這是提高代碼可讀性的最基礎的準則。
帕斯卡命名法和駝峰命名法是目前比較流行的兩種規則，不同語言采用的規則可能不一樣，但是要記住一點：保持團隊和個人風格一致。
1、帕斯卡命名法
帕斯卡命名法簡單地說就是：多個單詞組成一個名稱時，每個單詞的首字母大寫。比如：
在C#中，這種命名法常用于類、屬性，函數等等，在JS中，構造函數也推薦采用這種方式命名。
2、駝峰命名法
駝峰命名法和帕斯卡命名法很類似，多個單詞組成一個名稱時，第一個單詞全部小寫，后面單詞首字母大寫。比如：
駝峰命名法一般用于字段、局部變量、函數參數等等。在JS中，函數也常用此方法命名。
采用哪種命名規則并不絕對，最重要的是要遵守團隊約定，語言規范。
盡可能完整地描述函數所做的所有事情
有的開發者可能覺得相較于長函數名來說，短函數名看起來可能更簡潔，看起來也更舒服。但是通常來說，函數名稱越短其描述的意思越抽象。函數使用者對函數的第一印象就是函數名稱，進而了解函數的功能，我們應該盡可能地描述到函數所做的所有事情，防止使用者不知道或誤解造成潛在的錯誤。
舉個例子，假設我們做一個添加評論的功能，添加完畢后并返回評論總數量，如何命名比較合適呢？
展開 如何在ANSYS WORKBENCH中施加分段函數激勵
本篇回答一位朋友提出來的問題，說明如何在ANSYS WOKRBENCH中施加分段函數激勵。
假設分段的分布載荷如下
該載荷施加在一長方體的頂面上，作為分布力系施加。
下面說明操作方法。
1. 創建一個瞬態動力學分析系統
2.創建一長方體，尺寸任意。
3.劃分網格
4.分析設置
設置兩個時間步，
第一步終止時間為1秒，打開自動時間步長，通過載荷步來定義載荷子步，初始子步10步，最小5步，最多20步。
再定義第二步如下
其含義是
第2步終止時間為2秒，打開自動時間步長，通過載荷步來定義載荷子步，初始子步10步，最小5步，最多20步.
5.固定左端
6.在上面施加分布載荷1
首先定義第一個載荷步內的函數載荷
接著休眠期第二段（1-2秒內的部分）
得到結果如下
7.在上面施加分布載荷2
接著休眠期第1段（0-1秒內的部分）
得到結果如下
這就可以了。
至于后面的求解就不再贅述了。
來源：宋博士的博客，版權歸作者所有。
展開 『分享』在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
從形函數與函數的連續可導性到ansys結果中的節點解與單元解的差異
如題，《從形函數與函數的連續可導性到ansys結果中的節點解與單元解的差異》，形函數對結果的影響大部分人都能聯想到二次單元比線性單元求得的結果更精確，但該文要表達的不僅如此，而是從更一般地討論怎么從單元的形函數來理解節點解與單元解之間的差異。
首先討論單元的階次。作為基礎我們應該明白網格與單元的區別，網格是將幾何體離散化后的結構，即組成幾何體的微元，單元是這些微元的幾何、物理或數學屬性（這里我們并不打算詳細討論單元的這些屬性，但是這些知識會方便對本文的理解）。我們經常在使用ansys或其他CAE軟件時經常會遇到單元的選擇以及單元階次的選擇，一般一種單元包括線性單元和二次單元甚至更高級的單元，比如在ansys中經常被使用的shell181（左）和shell281（右），線性單元使用的形函數是一次的多項式，高次單元使用的形函數是高次的多項式，形函數用于描述相鄰節點之間的位移場，所以高次的單元可以更好的描述形狀復雜的幾何體。
不同于常規材料力學中通過平衡方程求解（首先求得的解是力解），有限元方式求解的特點是首先求解出的結果是節點的位移解，即displacement of nodes，所有的節點位移形成了位移場，在空間上位移場一定是連續的，但是不一定是平滑的。哎哎，是不是特別熟悉的感覺，正是和高數中函數的連續性和可導性兩個性質非常相似，不用奇怪，位移場本來就是用函數描述的，所以自然就存在函數的性質，所以用函數的性質來理解就可以方便解釋一些現象了，下圖分別是用兩種形函數描述的位移場，在有限元求解后得到的首先是節點位移解，即圖中5個節點的位移，假如每個節點的位移用坐標x\y\z的函數來表示，然后通過形函數插值得到相鄰節點之間的位移（也是xyz的函數），上圖是用一次形函數插值，下圖是用二次形函數插值。
展開 
ANSYS 中添加窗函數
比如說我要施加一個周期性的沖擊 希望他作用2個周期后停止
ANSYS中可以使用的數學函數
在ANSYS幫助系統中關于*SET命令的注釋下列出了ANSYS中可以使用的數學函數。所有這些數學函數均可以在ANSYS環境中使用，這些數學函數包括：
ABS(X) 求絕對值
ACOS(X) 反余弦
ASIN(X) 反正弦
ATAN(X) 反正切
ATAN2(X,Y) 反正切, ArcTangent of (Y/X) , 可以考慮變量X,Y 的符號
COS(X) 求余弦
COSH(X) 雙曲余弦
EXP(X) 指數函數
GDIS(X,Y) 求以X為均值,Y為標準差的高斯分布,在使用蒙地卡羅法研究隨機荷載和隨機材料參數時,可以用該函數處理計算結果
LOG(X) 自然對數
LOG10(X) 常用對數(以10為基)
MOD(X,Y) 求 X/Y的余數. 如果 Y=0, 函數值為 0
NINT(X) 求最近的整數
RAND(X,Y) 取隨機數,其中X 是下限, Y是上限
SIGN(X,Y) 取 X的絕對值并賦予Y的符號. Y>=0, 函數值為|X|, Y<0, 函數值為-|X|,.
SIN(X) 正弦
SINH(X) 雙曲正弦
SQRT(X) 平方根
TAN(X) 正切
TANH(X) 雙曲正切
展開 如何利用ANSYS的隨機分布函數功能
作者：水哥ANSYS
來源：本文源于ANSYS結構院，上海安世亞太授權轉載
隨機分布在材料微觀力學分析中扮演著重要角色，例如混凝土骨料力學、新型材料纖維力學分析等內容，提及隨機分布，更多的同學可能會聯想到采用第三方軟件如Matlab來生成，并導入ANSYS計算，其實ANSYS本身自帶隨機分布功能，只是功能略有限制。
ANSYS中產生隨機分布的一個重要函數是 *VFILL，該函數主要的作用是對數組進行填充賦值，而在賦值的過程中，用戶既可以選擇自定義數據內容，也可以選擇利用隨機函數產生數值，ANSYS Help中*VFILL說明如下：
該函數主要輸入參數為數組名稱以及輸入數據的函數，當選擇為data時，表示用戶自定義數據進行填充，當選擇其他選項時，則根據函數類型進行填充。
*VFILL隨機數生成支持均勻分布(Rand)、高斯分布(GDIS)、三角分布(TRIA)、貝塔分布（BETA）、伽馬分布（GRMM），*VFILL用于批量生成，如果需要單獨生成數據，則可以分別使用函數：
1) num=Rand(con1,con2)
2) num=Gdis(con1,con2)
3) num=Tria(con1,con2,con3)
4) num=Beta(con1,con2,con3,con4)5) num=Gram(con1,con2,con3)
上述con1~con4分別表示函數參數，例如針對均布分布，con1和con2分別表示分布的下限和上限。
下面分別以均布分布、高斯分布、伽馬分布為例進行演示。
1、均布分布
APDL代碼：
finish
/clear
/prep7
numA=1000
!
展開 Ansys Zemax光學設計軟件技術教程：如何編寫ZPL宏：計算環帶垂軸色差
本文介紹如何使用Zemax編程語言 (ZPL) 創建宏來計算和繪制用戶自定義的性能指標。 在本示例中，系統的環帶垂軸色差 (Zonal Transverse Chromatic Aberration, ZTCA) 將被分析。在解釋了需要使用的基本算法之后，本文將描述如何將宏通用化以允許用戶進行更加靈活的使用，包括用戶輸入和錯誤捕獲。作者 Alessandra Croce附件下載文章附件簡介假設我們要計算系統的環帶垂軸色差 (ZTCA)。ZTCA作為光瞳區域的函數，定義為兩個波長在像平面上的光線高度差ZTCA = YF - YC基本算法如下：
定義步進式光瞳區域
對應追跡每個波長的邊緣光線
減去y坐標
輸出結果
我們建議在編寫ZPL宏時將采用以下編程規范，使宏更易于閱讀：
將所有變量寫成小寫形式
將所有ZPL函數和關鍵字寫成大寫形式
對IF和FOR循環采用文本縮進
賦予變量有意義的名稱 （例如 ray_angle 而不是 a）
使用描述性注釋（例如使用 ! 或 #）
為了使宏更易于閱讀，OpticStudio具有內置的文本編輯器，將語法突出顯示：初版宏的編寫對于此示例，我們將使用示例文件<Documents>\Zemax\Samples\Sequential\Objectives\ Doublet.zmx。
展開 ANSYS知識普及4——如何施加函數變化的表面載荷 （ANSYS專家編輯，非原創，歡迎轉摘）
本人準備出一個ANSYS知識普及系列，將有用的網上資料歸攏，由于知識水平有限，不對之處請諒解。也歡迎各位網友提供好的資料分享，讓我們共同完成這個ANSYS知識普及系列。
編輯人：技術鄰ANSYS專家
業務咨詢網址：http://www.yqgqt.org.cn/content/other/402981
（打個小廣告）
聲 明:1、ANSYS知識普及系列中所有資料均來自網上；
2、如侵犯知識產權，請聯系ANSYS專家本人或者技術鄰，我將第一時間刪除。
小技巧：加本人關注，可以及時觀看本人發布的技術貼
ANSYS具有函數加載功能，可以很方便地在模型表面施加函數變化的各種載荷，在ANSYS中，也可以通過變通的方式來實現此功能，其思路是：
首先選定所要施加函數變化表面載荷的表面上的節點，利用ANSYS的參數數組和嵌入函數知識寫一簡單的命令流，定義好相應節點位置的面載荷值，然后通過在節點上施加面載荷來完成。
下面以在一圓柱表面施加函數變化載荷為例：
/prep7
et,1,45
cyl4,,,0.5,,,,3
vsweep,all
asel,s,loc,y,0.01,1
nsla
!
*get,nmax,node,,num,max,
*get,nmin,node,,num,min,
*afun,deg
*dim,t1,array,nmax,1,1,
csys,1
*do,k,nmin,nmax
*if,nsel(k),eq,1,then
t1(k)=1000*sin(ny(k))
*else
t1(k)=0
*endif
*enddo
!
sffun,pres,t1(1)
sf,all,pres,0
展開 ANSYS中不同形狀的波函數書寫方法
在ANSYS仿真中經常會遇到一些動態的加載方法，加載的載荷（位移、力、電流、溫度等）隨著時間而變化，表示不同的狀態。而相應的在workbench中可以方便的采用表格方法設置不同時間狀態下加載的位移或受力等載荷。但是又實用需要APDL命令的方式書寫不同時刻的載荷，但是函數庫當中又沒有相應的函數，那么如何書寫呢，下面我們選取幾個有代表性的書寫方法
（1）三角波的使用
一個物體在一個平面上移動，從左到右勻速運動，然后再從右到左的勻速運動。如果次數多，則采用表格方式實現比較麻煩，而采用do循環命令的方式就可以方便的加載。
該方法可以采用三角波的形式，給物體一個位移，使它不斷的左右移動，獲取摩擦熱或者應力應變等過程方法如下：
*do,x,1,25,1
cc=ABS(ASIN(ABS(SIN((x+PI()/2)/2)))-PI()/4)!
展開 如何在ADAMS中基于機械手終端位移確定電機驅動函數？
在機械手運動學問題中，經常給定了機械手終端的位移，要求據此確定轉動副處電機的驅動函數，以完成指定的機械手軌跡。那么如何根據給定的終端位移而計算轉動副處的驅動函數呢？
這屬于一個機器人運動學的逆解問題。本文使用ADAMS，對一個簡單的機械手，首先給定終端位移，然后據此仿真，得到各轉動副的運動方程。然后使用上述運動方程，驅動轉動副，以實現機械手的運動。此例最終說明，使用過ADAMS進行機器人運動學的反解是一種合適的選擇。
問題描述：
如圖所示的簡單機械手，桿件1固定在地面，而桿件2與桿件1，桿件3與桿件2之間通過轉動副連接。現在給定了桿件3終端的位移，要使用ADAMS計算出轉動副1和轉動副2的驅動位移函數。
求解步驟：
1)
創建機構。如下圖所示，創建桿1、2、3，桿1與地面之間創建固定副1，桿1、2和桿2、3之間分別創建轉動副1、2。
2)
創建一般點驅動.如下圖所示，選項為兩個物體一個位置，垂直于網格，其中第一個物體選擇桿3，第二個物體選擇地面，位置選擇桿3末端點，在彈出的對話框中選擇disp(time)位移-時間函數，如圖5所示。
3)
編輯點驅動函數。
展開 
如何在ADAMS中基于機械手終端位移確定電機驅動函數
在機械手運動學問題中，經常給定了機械手終端的位移，要求據此確定轉動副處電機的驅動函數，以完成指定的機械手軌跡。那么如何根據給定的終端位移而計算轉動副處的驅動函數呢？
這屬于一個機器人運動學的逆解問題。本文使用ADAMS，對一個簡單的機械手，首先給定終端位移，然后據此仿真，得到各轉動副的運動方程。然后使用上述運動方程，驅動轉動副，以實現機械手的運動。此例最終說明，使用過ADAMS進行機器人運動學的反解是一種合適的選擇。
問題描述：
如圖所示的簡單機械手，桿件1固定在地面，而桿件2與桿件1，桿件3與桿件2之間通過轉動副連接。現在給定了桿件3終端的位移，要使用ADAMS計算出轉動副1和轉動副2的驅動位移函數。
求解步驟：
1)
創建機構。如下圖所示，創建桿1、2、3，桿1與地面之間創建固定副1，桿1、2和桿2、3之間分別創建轉動副1、2。
2)
創建一般點驅動.如下圖所示，選項為兩個物體一個位置，垂直于網格，其中第一個物體選擇桿3，第二個物體選擇地面，位置選擇桿3末端點，在彈出的對話框中選擇disp(time)位移-時間函數，如圖5所示。
3)
編輯點驅動函數。
展開 Comsol小技巧| 8-在Comsol中如何設置電流隨時間變化的分段函數？
1
在 Comsol中，如何設置電流隨時間變化的分段函數？
可以采用邏輯表達式的方法，將電流寫成類似 I=I1*(t>=0 & t<=600)+I2*(t>600 & t<1200)+I3*(t>=1200 & t<=1800)的形式，I1、I2 和 I3分別表示 3 個階段下輸入的電流值。
2
在Comsol中如何自定義函數？
在設置函數（functions）時，要指定自變量和因變量。function name（函數名）就是因變量名。在函數列表中設定離散數據時，x 表示的是自變量數據，f(x)是對應的因變量數據。其中 x 不指坐標分量，而是用戶要設置的函數的自變量。例如，如果要設置 E_rod 是 H 的函數，就把 function name 設置為 E_rod，在函數列表的 x 列中輸入 H 的數據，在 f(x)列中輸入 E_rod 的數據。
3
Comsol中的變量 s 有何含義？
變量 s 是一個表示弧長的參數化幾何變量，該值是一個相對值，即考察的弧長與總弧長之間的比值。s 的定義與時間無關，僅僅與空間有關，即一個曲線(或直線)從起點開始為 0，到終點為 1，s 就表示測定點距起點的距離與整個弧長之間的相對比值，因此其范圍是[0,1]。詳細說明可參考用戶手冊中幾何變量這章的參數化變量部分。
歡迎關注 ~
編輯 | 電子F430
文案 | 小蘇
審核 | 趙佳樂
展開 如何從Ansys APDL中提取剛度矩陣與質量矩陣？ ￥69
首先，我們打開導出的矩陣：
從這個文件中我們可以看到第二行有五個數字，其對應的意義分別是：總行數、列指針個數、矩陣行索引總行數、矩陣元素數總行數、結點力向量總行數。
這些數據具體的意義與使用方法，都可以在下面的矩陣轉換文件”Transformer.m”中找出。
接下來我們直接打開matlab，將工作路徑設置為”Transformer.m”文件所在的路徑，
然后我們便可以開始使用矩陣提取m文件：Transformer.m，復制Stiffness_mat.dat文件（或者通過GUI導出的txt格式矩陣）的路徑，只需在matlab中輸入如下語句，便能直接得到我們想要的矩陣K1：
點擊K1
便可以看到最終的總剛度、質量矩陣了，其為大型稀疏矩陣，且為對稱矩陣。
到此，我們便成功通過Transformer.m函數完整提取出了我們想要的總剛度以及質量矩陣，整體過程十分簡單，易上手，不管你是要在工作還是學習中應用都有著不錯的可操作性，使用matlab中的m文件還會加深你對有限元程序設計的理解。
4.實戰應用與范例講解
接上一個矩陣的例子，其實際為Ansys中的一個應力集中問題模型所導出的剛度矩陣，那么我們如何來驗證其結果的準確性呢，這時我們就要用到結點力矩陣來進行驗證了，只要所解出來的位移與Ansys中可展示的結點位移相同，那么就證明我們的結果是準確無誤的。
展開 如何在ANSYS中模擬非線性三維隔震支座 ￥299
最近有很多同學聯系我，問到如何數值模擬三維隔震支座。假期加個班，做個算例分析。
1. 包含的內容
（1）算例模型命令流
（2）三維隔震支座命令流
（3）計算過程excel文件
（4）建筑隔震橡膠支座規范
（5）常用隔震支座的設計參數
2. 進階內容（需另付費，有需要可聯系）
（1）隔震支座在ANSYS中的批量建模方法，預計時間2024年02月
（2）如何在ABAQUS中模擬非線性單位隔震支座（連接器單元），預計時間2024年03月
3. 解決的問題
（1）如何在ANSYS中模擬橡膠隔震支座？
（2）如何確定隔震模型的力學參數與隔震支座設計參數的定量對應關系？
（3）如何模擬隔震支座的非線性特性？
（4）如何驗證隔震支座模擬的正確性？
4. 隔震模型的力學參數與隔震支座設計參數的定量對應關系
我們知道，實際應用中，我們可以采用廠家提供的標準型號的隔震支座，也可以訂制特殊類型的隔震支座，不管采用那種形式，在仿真模擬時，我們都要將設計參數與隔震模型的力學參數對應起來，從而進行力學分析。
ANSYS中并沒有特定的隔震單元，但提供了一系列的彈簧-阻尼器單元，可以通過組合單元模擬隔震支座的力學特性。采用COMBIN14單元模擬隔震支座的豎向剛度，COMBIN14又稱彈簧-阻尼器單元，具有1D、2D和3D的軸向或扭轉能力。軸向彈簧-阻尼器為單軸拉壓行為，每個單元有2個節點，每個節點有3個自由度，即沿著X、Y和Z方向的三個平動或轉動位移。水平方向上，采用COMBIN40單元模擬隔震支座的水平剛度和阻尼，COMBIN40單元將彈簧、滑塊和阻尼器并聯，再用串聯的方式與間隙耦合形成組合體，適用于多種情況的分析。該單元可以引入雙線性強化模型，并考慮粘滯阻尼的影響。詳細參考《ANSYS結構分析單元與應用》。
展開 