ansys正弦函數(shù)的案例
正弦函數(shù)曲線曲面的建立
我其實有很多資料想和大家分享，只是目前還沒有完全分好類別，那么，這一小節(jié)我來和大家分享一下正弦函數(shù)或者其他函數(shù)在曲線曲面中的應(yīng)用方法。
工作情景模式中又很多時候是需要使用到函數(shù)功能的，比如正余弦函數(shù)曲線，當(dāng)然如上一節(jié)我們說的漸開線一樣，還有的時候是需要按照自己畫的線條形成的曲線偏移在某些特殊線面的處理方面能給自己最大的自由度和靈活性。
比如上圖中這個曲面就是被正弦函數(shù)線切割的，那么我們接下來就看一看如何實現(xiàn)上訴需求。
1， 首先在草圖中創(chuàng)建一個園，然后用拉伸面命令拉伸一個圓柱面，當(dāng)然這個可以給高一點
然后我們找到law命令，創(chuàng)建我們需要的正選函數(shù)
鼠標(biāo)點擊fog（law）創(chuàng)建一個新的函數(shù)
然后輸入公式 rad表示弧度一定要加上
另外，如果有些公式不知道的話可以在字典（dictonary）里面找
創(chuàng)建好了參數(shù)之后會在結(jié)構(gòu)樹上關(guān)系（relations）出現(xiàn)：
接下來我們需要做一步split裁切，因為封閉曲線是無法使用平行的law的
按照abcdefg的順序，a,
首先選擇平行線，然后選擇切割的一半，support面選擇拉伸曲面。
關(guān)鍵在于law，點開之后選擇高級（advance）之后選取結(jié)構(gòu)樹上的law，如下
就快看到希望了對吧。別急，關(guān)鍵還要看然后做另外一半，哈哈，同樣做法哦
然后就簡單了，join兩條曲線，之后切割就好了。
希望大家都能靈活應(yīng)用，get到了這個技能你會發(fā)現(xiàn)你會再上一個新臺階。
展開 高熵合金施加循環(huán)應(yīng)力（正弦，三角函數(shù)）的分子動力學(xué)
疲勞失效是工程結(jié)構(gòu)件的主要破壞形式之一，通常由循環(huán)應(yīng)力（如正弦波載荷）作用下的微觀缺陷（如位錯聚集、裂紋萌生與擴展）逐漸累積所致。分子動力學(xué)（MD）模擬能夠在原子尺度揭示高熵合金在循環(huán)載荷下的微觀過程，為理解其抗疲勞機理提供重要依據(jù)。然而，目前針對高熵合金在正弦波循環(huán)應(yīng)力下的MD研究仍較為有限，尤其是不同成分、溫度及加載頻率對疲勞行為的影響仍需深入探索。本研究擬通過分子動力學(xué)模擬，對其開展研究。
1：建立長寬高均為150埃米的正方形盒子，在內(nèi)部填充Ni、Fe、Cr三種原子：
建立的模型如下圖所示：
初始模型在NPT系綜平衡后，在溫度為800K、周期為50ps，拉伸速率以正弦函數(shù)變化，最大拉伸速率為0.05s-1的條件下，使用loop命令循環(huán)10次，使用 fix 3 all deform 100 x erate ${speed} remap x units box命令，在x方向進行拉伸。
在lammps中拉伸的命令設(shè)置如下
模擬結(jié)束之后，在origin中畫出x方向應(yīng)變隨時間的變化情況：從圖中可以看出應(yīng)變符合正弦函數(shù)。
2：在上述條件下，將正弦函數(shù)可調(diào)整為三角形，同樣拉伸10次結(jié)果下圖所示，同樣驗證良好。
本次模擬主要更改了應(yīng)變的函數(shù)形式，溫度，拉伸頻率，周期都是固定的，后續(xù)可通過更改參數(shù)，進行更廣泛的研究，如峰值，循環(huán)次數(shù)，合金成分以及尺寸的影響，也可進一步增加缺陷，探究缺陷對拉伸的影響。
最后，歡迎通過公眾號“320科技工作室”與我們聯(lián)絡(luò)。
展開 由淺入深宏程序，車床旋轉(zhuǎn)正弦函數(shù)宏程序
下圖中正弦曲線如果以其左邊的端點為參考原點，則此條正弦曲線順時針旋轉(zhuǎn)了16度，即b=-16
此正弦曲線周期為24，對應(yīng)直角坐標(biāo)系的360
對應(yīng)關(guān)系 【0,360】 y=sin（x）
【0,24】 y=sin(360*x/24)
可理解為：
360/24是單位數(shù)值對應(yīng)的角度
360*x/24是當(dāng)變量在【0,24】范圍取值為x時對應(yīng)的角度
sin(360*x/24)是當(dāng)角度為360*x/24時的正弦函數(shù)值
旋轉(zhuǎn)正弦函數(shù)曲線粗精加工程序如下：
T0101
M3S800
G0X52Z5
#6=26 工件毛坯假設(shè)為50mm，#6為每層切削時向+X的偏移量。
N5 G0X[#6+18.539]
G1Z0F0.1
#1=48
N10 #2=sin【360*#1/24】
#4=#1*COS[-16]- #2*SIN[-16] 旋轉(zhuǎn)30度之后對應(yīng)的坐標(biāo)值
#5=#1*SIN[-16]+ #2*COS[-16]
#7=#4-【50-3.875】 坐標(biāo)平移后的坐標(biāo)。
#8=45+2*#5+#6
G1X[#8]Z[#7]F0.1 沿小段直線插補加工
#1=#1-0.5 遞減0.5，此值越小，工件表面越光滑。
IF [#1 GE 0] GOTO 10 條件判斷是否到達終點。
Z-50
G1X52 直線插補切到工件外圓之外
G0Z5
#6=#6-2
IF [#6 GE 0] GOTO 5
G0X150Z150
M5
M30
展開 abaqus蠕變基本設(shè)置及雙曲正弦函數(shù)損傷蠕變本構(gòu)CREEP子程序 ￥59.9
該部分為abaqus蠕變計算基本流程
ABAQUS蠕變問題計算流程.pdf
付費部分為使用CREEP子程序建立雙曲正弦函數(shù)蠕變損傷子程序，含到達預(yù)設(shè)損傷值（假設(shè)為1.0）后終止計算，和USDFLD子程序控制材料參數(shù)（該子程序可用于損傷后的材料退化，如蠕變第三階段或者蠕變疲勞分析，若不需要場變量控制可對該部分代碼進行刪除），相關(guān)理論請參考附件sci文獻?？商峁╆P(guān)于CREEP子程序的幫助文件學(xué)習(xí)的相關(guān)指導(dǎo)
Ansys workbench正弦駐頻轉(zhuǎn)隨機窄帶PSD譜的方式 ￥10
問題：
在Ansys Workbench進行隨機振動分析時，有時為了評估結(jié)構(gòu)共振條件下是否可以滿足要求。需要將環(huán)境PSD譜，疊加共振頻率的駐頻進行振動仿真。當(dāng)使用Ncode進行計算時可以實現(xiàn)同時輸入環(huán)境PSD譜和正弦駐頻。但是在Ansys Workbench進行隨機振動分析時，確不能同時輸入PSD譜和正弦駐頻。此時需要將正弦駐頻轉(zhuǎn)為窄帶隨機PSD譜，再將環(huán)境PSD與窄帶PSD的疊加譜輸入到Ansys Workbench進行隨機振動分析。
實現(xiàn)方法：
將正弦駐頻轉(zhuǎn)為窄帶隨機，可以依據(jù)1、能量等效原則。通過正弦信號的均方值等于窄帶隨機信號的均方值來換算。2、也可以通過兩種激勵狀態(tài)下結(jié)構(gòu)的最大加速度響應(yīng)幅值相等來換算。本文參考周炬老師《Ansys workbench有限元分析實例詳解-動力學(xué)》中給出的公式進行轉(zhuǎn)換。具體講解請參考教程。這里僅是將教材的轉(zhuǎn)換方法結(jié)合工作需求轉(zhuǎn)化為可以方便使用的excel工具。
應(yīng)用介紹：
Excel工具表如下。
以下是進行PSD換算所需的輸入信息：
? 首先環(huán)境PSD譜線信息。
? 然后根據(jù)結(jié)構(gòu)的模態(tài)仿真結(jié)果，確定結(jié)構(gòu)固有頻率為駐頻點。
? 正弦激勵幅值：通常依據(jù)頻率值所在范圍有相對應(yīng)的激勵幅值要求。
? 窄帶帶寬：通常由指定寬度、共振頻率的百分比等。
完成以上輸入信息后，點擊左上角“組合”按鈕即可得到，正弦駐頻轉(zhuǎn)窄帶隨機PSD+環(huán)境PSD的疊加結(jié)果。
將疊加后的PSD譜直接復(fù)制到Ansys Workbench中，再進行輸入Improved fit后即可進行正常隨機振動仿真。
示例：
1.模態(tài)疊加法隨機振動分析，計算結(jié)構(gòu)模態(tài)。
展開 從形函數(shù)與函數(shù)的連續(xù)可導(dǎo)性到ansys結(jié)果中的節(jié)點解與單元解的差異
如題，《從形函數(shù)與函數(shù)的連續(xù)可導(dǎo)性到ansys結(jié)果中的節(jié)點解與單元解的差異》，形函數(shù)對結(jié)果的影響大部分人都能聯(lián)想到二次單元比線性單元求得的結(jié)果更精確，但該文要表達的不僅如此，而是從更一般地討論怎么從單元的形函數(shù)來理解節(jié)點解與單元解之間的差異。
首先討論單元的階次。作為基礎(chǔ)我們應(yīng)該明白網(wǎng)格與單元的區(qū)別，網(wǎng)格是將幾何體離散化后的結(jié)構(gòu)，即組成幾何體的微元，單元是這些微元的幾何、物理或數(shù)學(xué)屬性（這里我們并不打算詳細討論單元的這些屬性，但是這些知識會方便對本文的理解）。我們經(jīng)常在使用ansys或其他CAE軟件時經(jīng)常會遇到單元的選擇以及單元階次的選擇，一般一種單元包括線性單元和二次單元甚至更高級的單元，比如在ansys中經(jīng)常被使用的shell181（左）和shell281（右），線性單元使用的形函數(shù)是一次的多項式，高次單元使用的形函數(shù)是高次的多項式，形函數(shù)用于描述相鄰節(jié)點之間的位移場，所以高次的單元可以更好的描述形狀復(fù)雜的幾何體。
不同于常規(guī)材料力學(xué)中通過平衡方程求解（首先求得的解是力解），有限元方式求解的特點是首先求解出的結(jié)果是節(jié)點的位移解，即displacement of nodes，所有的節(jié)點位移形成了位移場，在空間上位移場一定是連續(xù)的，但是不一定是平滑的。哎哎，是不是特別熟悉的感覺，正是和高數(shù)中函數(shù)的連續(xù)性和可導(dǎo)性兩個性質(zhì)非常相似，不用奇怪，位移場本來就是用函數(shù)描述的，所以自然就存在函數(shù)的性質(zhì)，所以用函數(shù)的性質(zhì)來理解就可以方便解釋一些現(xiàn)象了，下圖分別是用兩種形函數(shù)描述的位移場，在有限元求解后得到的首先是節(jié)點位移解，即圖中5個節(jié)點的位移，假如每個節(jié)點的位移用坐標(biāo)x\y\z的函數(shù)來表示，然后通過形函數(shù)插值得到相鄰節(jié)點之間的位移（也是xyz的函數(shù)），上圖是用一次形函數(shù)插值，下圖是用二次形函數(shù)插值。
展開 ansys的取值函數(shù)
有關(guān)實體狀態(tài)的取值函數(shù)
NSEL(N)
ESEL(E)
KSEL(K)
LSEL(L)
ASEL(A)
VSEL(V)
表示某個實體狀態(tài)，其返回值-1，沒有選中，0，沒有定義，1，被選中
有關(guān)下一個被選實體的取值函數(shù)
NDNEXT(N)
ELNEXT(E)
KPNEXT(K)
LSNEXT(L)
ARNEXT(A)
VLNEXT(V)
表示編號大于N,E,K,L,A,V的下一個被選實體
有關(guān)實體位置的取值函數(shù)
CENTRX(E)
CENTRY(E)
CENTRZ(E)
單元E在中心位置的X,Y,Z的坐標(biāo)系（直角坐標(biāo)系），有所選的節(jié)點決定
NX(N)
NY(N)
NZ(N)
KX(K)
KY(K)
KZ(K)
節(jié)點N或關(guān)鍵點K在激活坐標(biāo)系中X,Y,Z的坐標(biāo)值
LX(L,LFRAC)
LY(L,LFRAC)
LZ(L,LFRAC)
線段L在長度比率為LFRAC(0~1)時的X,Y,Z的坐標(biāo)值
有關(guān)最靠近某位置的節(jié)點或關(guān)鍵點編號的取值函數(shù)
NODE(X,Y,Z)
KP(X,Y,Z)
被選擇的節(jié)點嘴靠近X,Y,Z位置的節(jié)點或關(guān)鍵點編號（在激活的坐標(biāo)系下，如果存在多個節(jié)點或關(guān)鍵點，那么取其最小值）
有關(guān)距離的取值函數(shù)
DISTND(N1,N2)
DISTKP(K1,K2)
節(jié)點或關(guān)鍵點兩點之間的距離
DISTEN(E,N)
單元E的中心點與節(jié)點N之間的距離，中心點將由單元上被選擇的節(jié)點確定
有關(guān)角度的取值函數(shù)
ANGLEN(N1,N2,N3)
ANGLEK(K1,K2,K3)
節(jié)點或關(guān)鍵點兩條邊之間的夾角，缺省時單位為弧度，其中所選擇的3個節(jié)點中，N1或K1是頂點
有關(guān)最靠近實體的節(jié)點，關(guān)鍵點和單元的取值函數(shù)
NNEAR(N)
最靠近節(jié)點N的被選節(jié)點
KNEAR(K)
最靠近關(guān)鍵點K的被選關(guān)鍵點
ENEARN(N)
最靠近節(jié)點N的被選單元，單元的位置將由被選節(jié)點確定
有關(guān)面積的取值函數(shù)
展開 ANSYS 內(nèi)部函數(shù)
VLNEXT(N)
Next higher volume number above N in selected set (or zero if none
found).
30.距離函數(shù)
DISTND(N1,N2) Distance between nodes N1 and N2.
DISTKP(K1,K2) Distance between keypoints K1 and K2.
DISTEN(E,N) Distance between the centroid of element E and node N.
Centroid is determined from the selected nodes on the element.
31.角度函數(shù) (缺省單位為弧度，單位變換用 *AFUN 命令)
ANGLEN(N1,N2,N3) Subtended angle between two lines (defined by
three nodes where N1 is the vertex node). Default is in
radians.
ANGLEK(K1,K2,K3) Subtended angle between two lines (defined by
three keypoints where K1 is the vertex keypoint). Default is in
radians.
32.最近實體函數(shù)
NNEAR(N) Selected node nearest node N.
KNEAR(K) Selected keypoint nearest keypoint K.
ENEARN(N) Selected element nearest node N.
展開 Ansys Zemax | 什么是點擴散函數(shù)（ PSF ）
附件下載
聯(lián)系工作人員獲取附件
本文討論了如何在 OpticStudio 中對點擴散函數(shù)進行建模和解釋。使用的分析特征是 Spot Diagram、FFT PSF 和 Huygens PSF。將討論每種工具的優(yōu)點，以及用于最準確分析的有用特征設(shè)置。
介紹
光學(xué)系統(tǒng)的點擴散函數(shù) （PSF） 是單個點光源產(chǎn)生的輻照度分布。（望遠鏡拍攝遙遠恒星的圖像就是一個很好的例子。盡管源可能是一個點，但圖像不是。有兩個主要原因：首先系統(tǒng)中的像差會將圖像傳播到有限的區(qū)域；其次衍射效果也會擴散圖像，即使在沒有像差的系統(tǒng)中也是如此。
OpticStudio 有三種基本類型的 PSF 計算：幾何（無衍射）點列圖、基于衍射的 FFT 和 Huygens PSF。本文將討論基本理論，并就正確使用每種類型的 PSF 提供一些指導(dǎo)。
點列圖
OpticStudio 中最基本的分析功能之一是點列圖。此功能從物空間中的單視場點發(fā)射許多光線，通過光學(xué)系統(tǒng)追跡所有光線，并繪制所有光線相對于某個公共參考的 （x，y） 坐標(biāo)。因此，點列圖本身就可以看作一個幾何 PSF。
這里使用的示例光學(xué)系統(tǒng)是一個焦距為 50 mm 的單拋物面 F/5 反射鏡，物位于無窮遠處。該系統(tǒng)是一個簡化的牛頓望遠鏡，包含的示例文件為 PSF_Newtonian.ZMX。以下是光學(xué)系統(tǒng)的外觀：
兩個視場點（一個在軸上，另一個呈 2 度角）的點列圖如下所示。
請注意，點列圖是光線落點的集合，每個點表示一條光線。光線之間沒有相互作用或干擾。點列圖在顯示望遠鏡的幾何或光線像差的影響方面非常有效。離軸幾何 PSF 清楚地顯示了系統(tǒng)的彗差和像散。然而在軸上，點列圖預(yù)測了完美的成像。但這是否準確代表了光學(xué)系統(tǒng)的性能？為了回答點列圖結(jié)果的這個問題，我們需要將點列分布與衍射極限響應(yīng)進行比較。
展開 ANSYS 支持的函數(shù)列表
ANSYS 支持的函數(shù)列表，備用與共享，以后不要老再去找了
SIN(X) Sine
COS(X) Cosine
TAN(X) Tangent
ASIN(X) Arcsine
ACOS(X) Arccosine
ATAN(X) Arctangent
ATAN2(Y,X) Arctangent (Y/X) with the sign of each component considered
SINH(X) Hyperbolic sine
COSH(X) Hyperbolic cosine
TANH(X) Hyperbolic tangent
SQRT(X) Square root
ABS(X) Absolute value
SIGN(X,Y) Absolute value of X with sign of Y.
展開 ANSYS的“get”函數(shù)列表
自己收藏并與大家分享，來自于ＡＮＳＹＳ的ｈｅｌｐ
“get函數(shù)”可用于某些項，并可用于代替*get命令。函數(shù)返回值并在函數(shù)被輸入的地方使用它，繞過了用參數(shù)名存儲值和在要使用值的地方輸入?yún)?shù)名的需要。
例如，假設(shè)要計算兩個節(jié)點的平均X位置。使用*GET命令，參數(shù)L1可以指定節(jié)點1的X位置，參數(shù)L2可以指定節(jié)點2的X位置。然后，可以從mid=（L1+L2）/2計算mid位置：
*GET,L1,NODE,1,LOC,X
*GET,L2,NODE,2,LOC,X
MID=(L1+L2)/2
但是，使用返回節(jié)點N的X位置的節(jié)點位置“get ”函數(shù)NX（N），可以直接計算MID，而不需要中間參數(shù)L1和L2：
MID=(NX(1)+NX(2))/2
除非另有說明，否則Get函數(shù)返回活動坐標(biāo)系中的值。
Get函數(shù)參數(shù)本身可能是參數(shù)或其他Get函數(shù)。get函數(shù)NELEM（E，NPOS）返回元素編號E的NPOS位置的節(jié)點號。組合函數(shù)NX（NELEM（E，NPOS））返回該節(jié)點的X位置。
下表列出了按功能分組的可用get函數(shù)。*GET命令還列出GET函數(shù)作為*GET items的替代項（如果適用）
Table 1: *GET - Get Function Summary
"Get Function" Summary
Entity Status Get Function Description
NSEL(N) Status of node N: -1=unselected, 0=undefined, 1=selected.
ESEL(E) Status of element E: -1=unselected, 0=undefined, 1=selected.
展開 
Ansys Wrokbench分段復(fù)雜函數(shù)載荷，加載方式記錄 ￥10
問題：
Ansys Workbench的載荷加載形式有三種，constant/table/function。Constant是在載荷步內(nèi)給定恒定值；table形式較為便捷，可以在定義每個子步的載荷大?。?function形式可以輸入以time/X/Y/Z為變量的簡單方程。
但是仍有某些形式的載荷較難輸入，例如分段復(fù)雜函數(shù)載荷等。
解決方法：
需要使用Ansys經(jīng)典界面的function功能編輯分段載荷獲得ADPL載荷命令；再利用Workbench中command的形式施加載荷。
操作方式：
1. Ansys經(jīng)典中function公式編輯器輸入分段函數(shù)。
在function頁卡中選著變量time，在Regime頁卡中逐個定義分段函數(shù)；
定義完成后點擊保存，并輸入函數(shù)名“TEST3.func”
2. 再次點擊標(biāo)題欄的Parameters>Functions>Read From files>找到剛才保存的TEST3.func。并在Table Parameter Name中給編輯導(dǎo)入的分段函數(shù)命名PForce。此后分段函數(shù)即被公式編輯器編譯為表格數(shù)組形式，數(shù)組的名稱為：PForce。
3. 提取分段函數(shù)數(shù)值的ADPL命令形式，用于Workbench使用。
完成分段函數(shù)導(dǎo)入和命名后，在下拉列表中的File>List>Log file中可以查看經(jīng)典界面GUI操作對應(yīng)的ADPL命令。在這里可以將上述function公式編輯器導(dǎo)入的分段函數(shù)數(shù)組對應(yīng)ADPL命令顯示出來。（有時log file顯示不及時，再重復(fù)一次即可）
4. 在Workbench內(nèi)創(chuàng)建加載remote point點，并設(shè)定加載點的ADPL name為“LoadPoint“，用于加載。
展開 ANSYS 中添加窗函數(shù)
比如說我要施加一個周期性的沖擊 希望他作用2個周期后停止
ANSYS知識普及4——如何施加函數(shù)變化的表面載荷 （ANSYS專家編輯，非原創(chuàng)，歡迎轉(zhuǎn)摘）
本人準備出一個ANSYS知識普及系列，將有用的網(wǎng)上資料歸攏，由于知識水平有限，不對之處請諒解。也歡迎各位網(wǎng)友提供好的資料分享，讓我們共同完成這個ANSYS知識普及系列。
編輯人：技術(shù)鄰ANSYS專家
業(yè)務(wù)咨詢網(wǎng)址：http://www.yqgqt.org.cn/content/other/402981
（打個小廣告）
聲 明:1、ANSYS知識普及系列中所有資料均來自網(wǎng)上；
2、如侵犯知識產(chǎn)權(quán)，請聯(lián)系ANSYS專家本人或者技術(shù)鄰，我將第一時間刪除。
小技巧：加本人關(guān)注，可以及時觀看本人發(fā)布的技術(shù)貼
ANSYS具有函數(shù)加載功能，可以很方便地在模型表面施加函數(shù)變化的各種載荷，在ANSYS中，也可以通過變通的方式來實現(xiàn)此功能，其思路是：
首先選定所要施加函數(shù)變化表面載荷的表面上的節(jié)點，利用ANSYS的參數(shù)數(shù)組和嵌入函數(shù)知識寫一簡單的命令流，定義好相應(yīng)節(jié)點位置的面載荷值，然后通過在節(jié)點上施加面載荷來完成。
下面以在一圓柱表面施加函數(shù)變化載荷為例：
/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的函數(shù)加載
最近論壇里很多人問我，如何施加函數(shù)載荷。
我今天給大家奉獻的是任意函數(shù)加載的操作步驟詳解，手把手教大家操作自定義的函數(shù)加載。
如果覺得還不錯，頂下帖子，也算對我的鼓勵了！
大家有什么ANSYS 或 Workbench Mechanical 相關(guān)的問題，可以隨時**我 ansys123@qq.com
手把手教你ansys函數(shù)加載.doc
