位移時(shí)程曲線的案例
基于ABAQUS的建筑結(jié)構(gòu)時(shí)程分析
本文選擇帝谷波ELCENTRO作為地震動(dòng)輸入，其圖形如下圖所示：
分別輸入不同方向下的地震波數(shù)據(jù)后得到各個(gè)方向的計(jì)算指標(biāo)：
X向底部剪力時(shí)程曲線如下圖所示：
X向頂點(diǎn)（下圖）的位移時(shí)程曲線如下圖所示：
位移角數(shù)值如下圖所示
以上分析可以看出，位移角為1/2225，滿足1/1000的要求；
最大水平位移為0.08m，遠(yuǎn)小于0.3的要求。
Y向的計(jì)算結(jié)果與X向一樣的處理方法：輸入地震波數(shù)據(jù)，提取基底剪力與位移時(shí)程曲線，后續(xù)工作主要為與振型分解法的結(jié)果進(jìn)行比較，當(dāng)大于反應(yīng)譜法結(jié)果時(shí)需要進(jìn)行內(nèi)力及配筋的調(diào)整。
四、 小結(jié)
本文從ABAQUS的隱式分析算法原理，軟件設(shè)置中的關(guān)鍵參數(shù)理解，國家規(guī)范及實(shí)際案例這四大方面介紹了H.H.T算法在高層建筑結(jié)構(gòu)分析中的應(yīng)用，對實(shí)際復(fù)雜建筑結(jié)構(gòu)分析中具有很好借鑒作用，同時(shí)可以作為超限結(jié)構(gòu)分析的補(bǔ)充方法。
CPU:I7-10750H
內(nèi)存：16384MB
計(jì)算模型的處理技術(shù)：ABAQUS隱式分析
計(jì)算機(jī)耗時(shí)：30min
基于ABAQUS的建筑結(jié)構(gòu)時(shí)程分析20210824 V2.0.pdf
展開 沖擊碰撞時(shí)程曲線
我的沖擊碰撞時(shí)程曲線數(shù)據(jù)不是那種連續(xù)的平滑曲線，除了峰值以外，其他數(shù)據(jù)都是零，附圖一張，附表一份
沖擊力時(shí)程數(shù)據(jù).xlsx
關(guān)于ABAQUS動(dòng)力時(shí)程分析相對位移處理的問題
大家好，我是剛剛學(xué)ABAQUS的新手，前面做了一些關(guān)于靜態(tài)分析的問題，感覺ABAQUS用起來比較好，現(xiàn)在在學(xué)習(xí)用他作動(dòng)力時(shí)程分析，但是在后處理階段發(fā)現(xiàn)一個(gè)問題：做出的動(dòng)畫和想像中的不一樣，經(jīng)查找證明是結(jié)果輸出的絕對位移值，現(xiàn)在我想讓它輸出相對位移值的動(dòng)畫，怎么辦？
請前輩們指點(diǎn)，謝謝了！
時(shí)程分析初位移的施加,振動(dòng)衰減和固有頻率
刪除初位移約束
D,1,all !施加固定端約束
KBC,0 !斜坡
TIME,20*1/fre1 !計(jì)算時(shí)間20個(gè)周期
TIME,4 !計(jì)算時(shí)間4s
!ALPHAD,afa1
SOLVE !求解
FINISH
/POST26 !時(shí)間歷程后處理
FILE,'file','rst','.'
NSOL,2,2,U,Y,
XVAR,1
/YRANGE,-uy1,uy1,1
PLVAR,2, !撓度曲線
轉(zhuǎn)載自好學(xué)ANSYSA。詳細(xì)模型及文件請移步公眾號~鏈接：https://mp.weixin.qq.com/s/EO_StXY70wSSMDNraOc3Gg
展開 
無網(wǎng)格劃分新技術(shù)midas MeshFree - 瞬態(tài)響應(yīng)分析案例
⑤施加邊界條件和載荷
定義剛性連接
定義約束
定義瞬態(tài)速度載荷
⑥分析設(shè)置
瞬態(tài)響應(yīng)分析控制
定義時(shí)間步
⑦求解后查看結(jié)果
MeshFree位移動(dòng)畫
MeshFree應(yīng)力動(dòng)畫
MeshFree瞬態(tài)響應(yīng)分析—總體隨時(shí)間的響應(yīng)
同樣，利用ANSYS軟件對該模型進(jìn)行了瞬態(tài)動(dòng)力學(xué)分析，也給出了整體結(jié)構(gòu)的總體響應(yīng)云圖。
ANSYS動(dòng)畫-位移
ANSYS動(dòng)畫-應(yīng)力
隨后，提取模型中桿件中部一點(diǎn)的變形和應(yīng)力的時(shí)程數(shù)據(jù)，進(jìn)行對比
提取點(diǎn)所在位置
MeshFree中該點(diǎn)位移時(shí)程曲線
MeshFree中該點(diǎn)應(yīng)力時(shí)程曲線
ANSYS中該點(diǎn)位移時(shí)程曲線
ANSYS中該點(diǎn)應(yīng)力時(shí)程曲線
經(jīng)過上述云圖及曲線的對比，兩款軟件計(jì)算時(shí)程響應(yīng)是一致的，也說明MeshFree的瞬態(tài)響應(yīng)分析模塊結(jié)果是準(zhǔn)確的。
可靠
ANSYS分析過程
①選擇分析類型
進(jìn)入ANSYS Workbench，將瞬態(tài)響應(yīng)分析模塊（Transient Structural）拖入工程面板。
②導(dǎo)入幾何
③設(shè)置材料模型
為保證兩款軟件數(shù)據(jù)對比的合理性，材料的選擇與MeshFree一致。
④基礎(chǔ)設(shè)置
接觸采用程序自動(dòng)創(chuàng)建的接觸，接觸類型為綁定接觸。
⑤網(wǎng)格劃分
由于物體較規(guī)則，只需設(shè)置網(wǎng)格的總體尺寸即可。
特備注意的是，下方長塊部件只做顯示導(dǎo)軌的作用，不參與計(jì)算，因此設(shè)置為剛體。
展開 Abaqus采用Python腳本保存定制化格式時(shí)程曲線圖 ￥19.89
在之前的一篇帖子里采用Python腳本輸出了某個(gè)值相對于時(shí)間的變化曲線的截圖
https://www.yqgqt.org.cn/content/post/1272569
但是輸出的曲線圖是Abaqus本身默認(rèn)的格式，包括坐標(biāo)軸字體、顏色、曲線的線型等，美觀性上比較差勁。當(dāng)然在用腳本做這個(gè)事情的時(shí)候是可以定制化圖表中的所有格式的。比如要想實(shí)現(xiàn)如下圖片中的效果：修改了legend位置和字體、title字體、x/y坐標(biāo)軸字體和坐標(biāo)軸title的字體、曲線的線型和Marker，附件中是相應(yīng)的腳本代碼。
展開 一致輸入和多點(diǎn)輸入下超長鋼框架結(jié)構(gòu)動(dòng)力彈塑性時(shí)程分析
考慮行波效應(yīng)影響，分別輸入一致位移激勵(lì)和多點(diǎn)位移激勵(lì)。一致位移激勵(lì)輸入時(shí)，各柱底位移時(shí)程輸入相同，如圖2（c）所示。擬設(shè)定以X正方向?yàn)榈卣鸩▊鞑シ较颍暡ㄋ偃?00m/s，共劃分6組位移激勵(lì)，每組間隔為20m，即時(shí)間差為0.2s。多點(diǎn)位移激勵(lì)輸入如圖3所示，C1~C6分別表示6組位移時(shí)程曲線輸入，D1、D2表示角柱柱底位置，D3、D4表示結(jié)構(gòu)中部的邊柱柱底位置。
（a）加速度時(shí)程曲線
（b）速度時(shí)程曲線
（c）位移時(shí)程曲線
圖2 EL-centro波（南北）時(shí)程曲線圖
3.算例結(jié)果
3.1 基底剪力
選取Ｘ向最大基底剪力分析，兩種位移輸入模式下的各基底剪力最大值對比如表2所示。在一致激勵(lì)中，C1~C5位置基底最大剪力相差較小，平均約為650kN；最大值出現(xiàn)在中部C3位置，為652.921kN；最小值出現(xiàn)在端部第C6組，其的最大基底剪力最小，平均和前五組相差約100kN；整體由C1~C6先增大（C1~C3）后減小（C3~C6）。由于各組基底位移時(shí)程同步，各組會(huì)同時(shí)出現(xiàn)最大基底剪力，因此，其最大的總基底剪力（3814.861kN）與6個(gè)組的總的最大基底剪力之和相同。
圖3 位移激勵(lì)輸入示意圖
在多點(diǎn)激勵(lì)中，最大值出現(xiàn)在地震波傳播的起始輸入端第C1組位置，為2335.78kN；最小值出現(xiàn)在結(jié)構(gòu)中部第C4組位置，為1764.1kN；整體由C1~C6先減小（C1~C4）后增大（C4~C6）。由于各組基底位移時(shí)程曲線有相位差，各組的最大基底剪力不會(huì)同時(shí)出現(xiàn)，且各組基底剪力相互影響，因此，總和的基底剪力反而較小，在同一時(shí)刻最大的總基底剪力僅為730.01kN，遠(yuǎn)小于一致激勵(lì)輸入。
展開 請教各位仁兄,在求解空間框架結(jié)構(gòu)抗震時(shí)程分析中的位移時(shí),不知道用什么軟件可以求出?我先前用MATLAB編程,
無
ABAQUS UEL二次開發(fā)（動(dòng)力隱式C3D8自定義單元）
荷載幅值為：1e6，即：
*Cload, amplitude=Amp-1
Part-1-1.forced, 2, -1e+06
懸臂端角點(diǎn)加載向（U2）位移時(shí)程曲線如下圖：
懸臂端角點(diǎn)軸向（U3）位移時(shí)程曲線如下圖：
第3s加載向加速度云圖如下圖：
第7s加載向加速度云圖如下圖：
第10s加載向加速度云圖如下圖：
（五）結(jié)語
該UEL的BBAR修正部分，參考了帖子：https://www.yqgqt.org.cn/post/05a24d5b-ca1d-4206-861d-633f702a80d7。
得到了發(fā)帖者的幫助。向發(fā)帖者表示感謝！
展開 地震時(shí)程分析、滯回曲線分析、（鋼管）混凝土本構(gòu)成和數(shù)據(jù)提取等軟件大集合 ￥99
<p>地震時(shí)程分析軟件合集：①人工波生成軟件simqke（施加地震波）；②中國規(guī)范反應(yīng)譜生成軟件GB_Spectrum；③太平洋地震工程研究中心PEER網(wǎng)站（<a href="https://ngawest2.berkeley.edu/）" rel="noopener noreferrer" target="_blank">https://ngawest2.berkeley.edu/）</a> 地震動(dòng)數(shù)據(jù)格式轉(zhuǎn)換器SeismoSingal 。</p><p>滯回分析軟件合集：①求解屈服點(diǎn)，極限荷載/位移，延性系數(shù)的軟件→Yielding Point and Ductility/YDP；②求解滯回環(huán)能量、退化剛度、耗能系數(shù)軟件（Hysteristic Loop Analysis/HLA）；③骨架曲線提取軟件（Skeleton）；④往復(fù)加載制度生成軟件(Cyclic-Loading-Process-Generator)；⑤提取論文圖片中滯回曲線數(shù)據(jù)到excel軟件</p><p>pushover分析軟件（可求解需求譜曲線、能力譜曲線和性能點(diǎn)）、中國規(guī)范混凝土骨架曲線生成軟件、鋼管混凝土塑性<a href="https://www.yqgqt.org.cn/qa/4700" rel="noopener noreferrer" target="_blank">應(yīng)力</a>應(yīng)變計(jì)算軟件（Python腳本插件）、ABAQUS纖維離散生成器、Rayleigh阻尼與頻率的計(jì)算、混凝土CDP模型與鋼筋本構(gòu)模型計(jì)算表格、約束混凝土Mander本構(gòu)計(jì)算表格等</p><p>PS：如果有小伙伴想單獨(dú)買某個(gè)軟件也可以QQ聯(lián)系我</p><p><br></p>
展開 基于ANSYS APDL的車橋耦合振動(dòng)分析程序 ￥299
主要內(nèi)容包括：
（1）120m連續(xù)鋼混組合梁橋模型（實(shí)體單元+殼單元+梁單元+栓釘建模細(xì)節(jié)、支座建模細(xì)節(jié)、橋墩建模細(xì)節(jié)）；
（2）空間整車模型，可考慮車體豎向，俯仰和側(cè)傾振動(dòng)加速度；
（3）車橋耦合振動(dòng)分析程序（可以修改車速，車重和路面不平整度）；
（4）結(jié)果提取程序，可以提取橋梁任意節(jié)點(diǎn)位移時(shí)程曲線，加速度時(shí)程曲線，車輛多個(gè)方向動(dòng)力響應(yīng)。
（使用該程序已發(fā)表sci論文3篇，1篇檢索，1篇已錄用，1篇返修中，可提供檢索論文）
buildings-13-01109-v2.pdf
Driving adaptability of highway steel-concrete composite beam bridge with multiple damages theory technology and practice.pdf
4304704.pdf
展開 
COMSOL黏彈性動(dòng)力邊界及地震動(dòng)輸入
持續(xù)時(shí)間為0.2s，計(jì)算時(shí)長為1s，計(jì)算時(shí)間步為0.001s，瞬態(tài)隱式求解，時(shí)間進(jìn)步方法為向后差分。
圖1 二維土體計(jì)算模型
圖2 入射波位移時(shí)程曲線圖
圖3 入射波速度時(shí)程曲線圖
計(jì)算結(jié)果如圖4、5所示，入射波在經(jīng)過0.25s之后到達(dá)自由表面與反射波疊加，變成入射波位移的2倍，0.4s之后自由地表停止振動(dòng)（圖中藍(lán)色部分為數(shù)值震蕩），說明入射波在底部黏彈性邊界處被吸收，沒有二次反射。與文獻(xiàn)上的結(jié)果也是一樣的。
圖4 頂部位移檢測點(diǎn)時(shí)程曲線
圖5 底部監(jiān)測點(diǎn)時(shí)程曲線
文獻(xiàn)結(jié)果
從上述結(jié)果可看出，本貼方法對比ABAQUS模擬結(jié)果效果相當(dāng)。當(dāng)然該方法在數(shù)值結(jié)果中精確度也有待提高，歡迎給位聯(lián)系我互相探討！想要學(xué)習(xí)本算例的歡迎各位前來咨詢。
展開 四層RC框架開洞，八度多遇地震下底部剪力頂層位移?模態(tài)分析，push over，時(shí)程分析,反應(yīng)譜計(jì)算 ￥100
計(jì)算八度多遇地震下底部剪力及頂層位移?(方法：push over，時(shí)程分析，反應(yīng)譜計(jì)算作對比)
提共模型(DAT,CAE,INP,ODB文件)，計(jì)算小插件，計(jì)算結(jié)果，以及計(jì)算截圖!
橋梁箱梁渦激振動(dòng)（渦振）仿真 ￥245
udf導(dǎo)入：2dof，龍格庫塔法
監(jiān)測：x、y向位移，三分力系數(shù)等
結(jié)果：速度云圖
結(jié)果：位移時(shí)程曲線
ABAQUS 二次開發(fā)（UEL+UMAT C3D8 單元?jiǎng)恿﹄[式計(jì)算）
以方塊的受動(dòng)力簡諧荷載為例，采用上述程序，應(yīng)用動(dòng)力隱式計(jì)算分析步，最終計(jì)算的位移、應(yīng)變等時(shí)程曲線均與ABAQUS保持一致。
（）模型信息
模型尺寸為10x10x10，彈性模量1e10，密度2000，泊松比0.25，荷載和邊界條件示意圖為：
一面的所有節(jié)點(diǎn)均固定。另一面所有的節(jié)點(diǎn)施加簡諧荷載，簡諧荷載曲線為（詳細(xì)的參數(shù)見附件）：
計(jì)算的網(wǎng)格圖為：
網(wǎng)格尺寸為1，共計(jì)10x10x=1000個(gè)單元。總的計(jì)算時(shí)長為1s，這只固定增量步長為0.01s，所以總增量步數(shù)為100。
（）計(jì)算結(jié)果
以上面網(wǎng)格圖中中間角點(diǎn)為例，提取加載向位移時(shí)程曲線如下圖：
0.89s時(shí)刻x向應(yīng)力云圖為：
注：左一為ABAQUS計(jì)算結(jié)果，中間為umat實(shí)現(xiàn)線彈性各向同性本構(gòu)計(jì)算結(jié)果，右一為UEL+UMAT計(jì)算結(jié)果，后面的云圖也是這樣排布。
0.89s時(shí)刻y向應(yīng)力云圖為：
寫文字好麻煩，不寫了！放個(gè)動(dòng)圖算了，有興趣的移步附件。
展開 