./三向輸入簡化后的單向輸入

首先,將三個方向的地震加速度放到一個文本文件里,如accexyz.txt,在這個數(shù)據(jù)文件里共放三列數(shù)據(jù),每列為一個方向的地震加速度值,這里僅給出數(shù)據(jù)文件中前幾行的數(shù)據(jù)：-0.227109E-02-0.209046E+000.467072E+01-0.413893E-02-0.168195E+000.261523E+01

-0.574753E-02-0.157890E+000.809014E-01

-0.731227E-02-0.152996E+000.119975E+01-0.876865E-02-0.138102E+000.130902E+01

-0.101067E-01-0.131582E+000.143611E+00

然后,再建一個文本文件用來存放三個方向的地震加速度時間點(diǎn),如time.txt,在這個數(shù)據(jù)文件里僅一列數(shù)據(jù),對應(yīng)于加速度數(shù)據(jù)文件里每一行的時間點(diǎn),這里給出數(shù)據(jù)文件中前幾行數(shù)據(jù)：0.100000E-01

0.200000E-010.300000E-010.400000E-01

0.500000E-01

0.600000E-01

編寫如下的命令流文件,并命名為acce.inp

*dim,ACCEXYZ,TABLE,2000,3!01行

*vread,ACCEXYZ<1,1>,accexyz,txt,,JIK,3,2000!02行<3e16.6>!03行

*vread,ACCEXYZ<1,0>,time,txt!04行

<e16.6>!05行

ACCEXYZ<0,1>=1!06行

ACCEXYZ<0,2>=2!07行,同上

ACCEXYZ<0,3>=3!08行,同上

finish

/SOLU

ANTYPE,trans

btime=0.01!定義計算起始時間

etime=15.00!定義計算結(jié)束時間

dtime=0.01!定義計算時間步長

*DO,itime,btime,etime,dtime

time,itime

AUTOTS,0

NSUBST,1,,,1

KBC,1

acel,ACCEXYZ<itime,1>,ACCEXYZ<itime,2>,ACCEXYZ<itime,3>!施加三個方向的地震加速度

SOLVE

*ENDDO

最后,在命令窗口里輸入/input,acce,inp即可對結(jié)構(gòu)進(jìn)行地震動力分析。

說明和討論：

1、命令流中各行說明：

01行：定義2000行,3列的數(shù)組,<行數(shù)根據(jù)數(shù)據(jù)文件里加速度點(diǎn)數(shù)來定>ACCEXYZ用來存放三個方向的加速度值。

02行：從數(shù)據(jù)文件accexyz.txt里讀加速度值到數(shù)組ACCEXYZ,2000為行數(shù),可根據(jù)情況修改。其中的JIK,3,2000非常重要,它決定著將從加速度數(shù)據(jù)文件中數(shù)據(jù)輸入到ACCEXYZ數(shù)組時的存放格式,這里用到按JIK方式變化,即讀進(jìn)來的數(shù)據(jù)依次放入ACCEXYZ<I,J>〔J從1到3,I從1到2000,K默認(rèn)從1到1。根據(jù)這行命令下面<3e16.6>的格式,每次從數(shù)據(jù)文件里讀一行三個數(shù)據(jù),分別放入ACCEXYZ<1,1>、ACCEXYZ<1,2>、ACCEXYZ<1,3>,接著再讀下一行的三個數(shù)據(jù),分別放入ACCEXYZ<2,1>、ACCEXYZ<2,2>、ACCEXYZ<2,3>,依次類推。當(dāng)然依據(jù)數(shù)據(jù)文件的格式,也可以選用IJK,IKJ,JIK,JKI,KIJ,KJI等格式,其中IJK為默認(rèn)。03行：讀數(shù)據(jù)的格式,每行三個數(shù)值。由于數(shù)據(jù)文件中的數(shù)據(jù)是用科學(xué)記數(shù)年法表示的,因此,這里也用相應(yīng)FORTRAN的科學(xué)記數(shù)法的格式。如果數(shù)據(jù)文件里的數(shù)值是如"0.28761.23332.9938"這樣的格式,此行的數(shù)據(jù)格式也就修改為"<3f10.4>"這樣的格式。

04行：從數(shù)據(jù)文件time.txt里讀時間值到數(shù)組ACCEXYZ的第零列。

06行：將數(shù)組ACCEXYZ的第零行賦值,如果不對行插值的話也可以不賦值

對于地震波的輸入,可以把荷載記錄做成文件,利用apdl的讀取功能讀入倒數(shù)據(jù)庫中。下面的例子是自己編的一個小文件。修改一下可以更簡潔。有用到的朋友自己作一下把。

fini

/config,nres,1000

*dim,aceX,TABLE,3000,1

*dim,aceY,TABLE,3000,1

*dim,aceZ,TABLE,3000,1

*creat,ff

*vread,aceX<1,1>,acex,txt,,1

<e16.6>

*vread,aceX<1,0>,ACETT,,,1

<e17.6>

ACEX<0,1>=1

*end

/input,ff

*creat,ff

*vread,aceY<1,1>,acey,txt,,1

<e16.6>.*vread,aceY<1,0>,ACETT,,,1

<e17.6>

ACEY<0,1>=1

*end

/input,ff

*creat,ff

*vread,aceZ<1,1>,acez,txt,,1

<e16.6>

*vread,aceZ<1,0>,ACETT,,,1

<e17.6>

ACEZ<0,1>=1

*end

/input,ff

！地震波時程記錄分成了3個文件,每個文件是一列。分別記錄x,y,z方向的加速度。acett是時間記錄。

這樣就可以把加速度記錄讀取倒ansys數(shù)據(jù)庫中作為數(shù)組。

也可以把加速度記錄做成一個文件,這樣程序就簡單多了。大家可以試看看修改一下。

下面是計算部分語句：

/SOLU

ANTYPE,trans

！求解其自己選了

TM_START=0.01

TM_END=15.00

TM_INCR=0.01

*DO,TM,TM_START,TM_END,TM_INCR

TIME,tm

alpha,

BETAD,

ACEL,acex<tm>,acey<tm>,acez<tm>

SOLVE

*ENDDO

fini

!借助于ANSYS動力分析基本功能和APDL語言以及有關(guān)的地震資料編寫的SEISMIC程序

NT=100

DT=0.02

*DIM,ac,,NT

/INPUT,FF,txt

/SOLU

NSUBST,1,,,1

OUTRES,ALL,1

ANTYPE,TRANS

*DO,I,1.NT

ACEL,0.01*ac<I>*1.3,0,0!ac<I>的值由地震資料獲得

TIME,I*DT

SOLVE

*ENDDO

地震分析算例<ANSYS>土木工程中除了常見的靜力分析以外,動力分析,特別是結(jié)構(gòu)在地震荷載作用下的受力分析,也是土木工程中經(jīng)常遇到的問題。結(jié)構(gòu)的地震分析根據(jù)現(xiàn)行抗震規(guī)要求,一般分為以下兩類：基于結(jié)構(gòu)自振特性的地震反應(yīng)譜分析和基于特定地震波的地震時程分析。

本算例將以一個4質(zhì)點(diǎn)的彈簧－質(zhì)點(diǎn)體系來說明如何使用有限元軟件進(jìn)行地震分析。更復(fù)雜結(jié)構(gòu)的分析其基本過程也與之類似。關(guān)鍵知識點(diǎn)：

<a>模態(tài)分析

<b>譜分析

<c>地震反應(yīng)譜輸入

<d>地震時程輸入

<e>時程動力分析

<1>在ANSYS窗口頂部靜態(tài)菜單,進(jìn)入Parameters菜單,選擇ScalarParameters選項,在輸入窗口中填入DAMPRATIO=0.02,即所有振型的阻尼比為2%

<2>ANSYS主菜單Preprocessor->Elementtype->Add/Edit/Delete,添加Beam188單元

<3>在ElementTypes窗口中,選擇Beam188單元,選擇Options,進(jìn)入Beam188的選項窗口,將第7個和第8個選項,Stress/Strain<SectPoints>K7,Stress/Strain<SectNods>K8,從None改為MaxandMinOnly。即要求Beam188單元輸出積分點(diǎn)和節(jié)點(diǎn)上的最大、最小應(yīng)力和應(yīng)變

<4>在ElementTypes窗口中,繼續(xù)添加Mass21集中質(zhì)量單元

<5>下面輸入材料參數(shù),進(jìn)入ANSYS主菜單Preprocessor->MaterialProps->MaterialModels菜單,在MaterialModelNumber1中添加Structural->Linear->Elastic->Isotropic屬性,輸入材料的彈性模量EX和泊松比PRXY分別為210E9和0.3。

<6>繼續(xù)給MaterialModelNumber1添加Density屬性,輸入密度為7800。

<7>繼續(xù)給MaterialModelNumber1添加Damping屬性,采用參數(shù)化建模,輸入阻尼類型為Constant,數(shù)值為DAMPRATIO

<8>接著建立梁單元的幾何屬性,和上一個例子一樣,采用Sections建模,進(jìn)入ANSYS主菜單Preprocessor->Sections->Beam->CommonSections,選擇Sub-Type為工字型,截面尺寸W1=0.2,W2=0.2,W3=0.5,t1=0.01,t2=0.01,t3=0.008

<9>通過實(shí)參數(shù)輸入集中質(zhì)量單元的質(zhì)量和轉(zhuǎn)動慣量,在ANSYS主菜單中選擇Preprocessor->RealConstants->Add/Edit/Delete菜單,在RealConstants窗口中選擇Add,在ElementtypeforRealConstants選擇Mass21,在RealConstantforNumber1窗口中輸入1.6E2,1.6E2,如圖所示。即該質(zhì)量單元在X和Y方向的質(zhì)量都為160,由于本例子模型為平面問題,所以不必考慮Z方向的質(zhì)量,同樣也不考慮單元的轉(zhuǎn)動慣量。

<10>繼續(xù)添加第二類集中質(zhì)量,過程和上面一樣,但是輸入的質(zhì)量數(shù)值為1.2E2,1.2E2

<11>完成以上工作就完成了模型的基本數(shù)據(jù)準(zhǔn)備,下面開始建立物理模型。

<12>在ANSYS主菜單中選擇Preprocessor->Modeling->Create->Keypoints->InActiveCS,依次輸入關(guān)鍵點(diǎn)編號和坐標(biāo)：需要說明的是,關(guān)鍵點(diǎn)6為后面建立梁單元所需的截面方向控制點(diǎn),在上一個例子中已經(jīng)做過介紹。

<13>完成關(guān)鍵點(diǎn)輸入后下面建立直線模型。在ANSYS主菜單中選擇Preprocessor->Modeling->Create->Lines->Lines->StraightLine,依次連接關(guān)鍵點(diǎn)1～5。

<14>下面給建立完的幾何模型賦予材料屬性,在ANSYS菜單中選擇Meshing->MeshAttributes->PickedLines,選中所有的直線,進(jìn)入LineAttributes窗口,選擇相關(guān)選項材料屬性,實(shí)參數(shù),單元類型和截面類型都為1,點(diǎn)擊OK后輸入關(guān)鍵點(diǎn)6作為截面方向控制點(diǎn)。

<15>在ANSYS主菜單中選擇Preprocessor->Meshing->SizeCntrls->ManualSize->Lines->PickedLines,在ElementSizesonPickedLines窗口中設(shè)定NDIVNo.ofelementdivisions為3,即將每條直線分為3段

<16>在ANSYS主菜單中選擇Preprocessor->Meshing->Mesh->Lines,選擇所有的直線,完成直線的網(wǎng)格劃分。

<17>為了便于后面操作,將網(wǎng)格劃分后的單元和節(jié)點(diǎn)編號進(jìn)行適當(dāng)?shù)那謇怼Ｔ贏NSYS主菜單中選擇Preprocessor->NumberingCtrls->MergeItems,在MergeCoincidentorEquivalentlyDefinedItems窗口中選擇All,清理所有重復(fù)的元素。同樣選擇NumberingCtrls->CompressNumber菜單,在CompressNumber中選擇All,對節(jié)點(diǎn)和單元進(jìn)行重新編號。

<18>下面建立集中質(zhì)量單元,采取直接輸入單元的方法建立。在ANSYS主菜單中選擇Preprocessor->Modeling->Create->Elements->ElemAttributes,在ElementAttributes中設(shè)定單元類型編號為2MASS21,材料編號任意,實(shí)參數(shù)編號為1。如圖所示

<19>在ANSYS主菜單中選擇Preprocessor->Modeling->Create->Elements->AutoNumbered->ThruNodes,選擇節(jié)點(diǎn)2,建立第一個集中質(zhì)量。

<20>再次進(jìn)入第18步ElementAttributes窗口,設(shè)定實(shí)參數(shù)<Realconstantsetnumber>為2。

<21>重復(fù)19步,選擇節(jié)點(diǎn)8,14,20,建立其他的三個集中質(zhì)量單元。

<22>到此完成所有建模工作,下面開始進(jìn)行結(jié)構(gòu)分析

<23>進(jìn)入ANSYS主菜單中Solution功能模塊,選擇Solution->DefineLoads->Apply->Structural->Displacement->OnNodes,選擇節(jié)點(diǎn)1,設(shè)定約束所有的自由度。

<24>首先做一次靜力分析,選擇ANSYS主菜單Solution->AnalysisType->NewAnalysis,設(shè)定分析類型為Static

<25>選擇ANSYS主菜單Solution->Solve->CurrentLS選項,進(jìn)行一次靜力分析

<26>分析完后,下面進(jìn)行模態(tài)分析,在ANSYS主菜單中選擇Solution->AnalysisType->NewAnalysis,選擇分析類型為Model

<27>在ANSYS主菜單中選擇Solution->AnalysisType->AnalysisOptions,輸入模態(tài)分析方法為子空間法<Subspace>,求解8階模態(tài),同時需要作模態(tài)擴(kuò)展,擴(kuò)展的模態(tài)為8階,并計算單元應(yīng)力和應(yīng)變,輸入窗口如圖

<28>再次選擇ANSYS主菜單Solution->Solve->CurrentLS,計算當(dāng)前問題

<29>這時,如果需要看結(jié)果,可以進(jìn)入后處理模塊,即ANSYS主菜單GeneralPostproc,可以看到計算的各階頻率和振型。

<30>完成結(jié)構(gòu)自振分析后,下面就可以進(jìn)行反應(yīng)譜分析

<31>進(jìn)入ANSYS主菜單Solution->AnalysisType->NewAnalysis,選擇分析類型為Spectrum。

<32>在ANSYS主菜單Solution->AnalysisType->AnalysisOptions中,選擇譜分析的類型為單點(diǎn)輸入<Sing-ptresp>

〔33>在ANSYS主菜單中選擇Preprocessor->Loads->LoadStepOpts->Time/Frequenc->Damping,輸入所有的阻尼為DAMPRATIO,如圖所示〔

34>下面需要定義地震的反應(yīng)譜。我國規(guī)給定的是基于加速度的反應(yīng)譜。在ANSYS主菜單中選擇Solution->LoadStepOpts->Spectrum->SinglePoint->Settings,設(shè)定反應(yīng)譜類型為地震加速度,放大系數(shù)為1,輸入方向為X方向〔1,0,0。

<35>接下來開始輸入地震反應(yīng)譜。這里輸入的反應(yīng)譜按7度多遇地震,取地震影響系數(shù)為0.08,第一組,III類場地,卓越周期Tg=0.45s。值得注意的是,我國規(guī)給的反應(yīng)譜橫坐標(biāo)是周期,ANSYS定義的反應(yīng)譜橫坐標(biāo)是頻率,應(yīng)該注意上述區(qū)別。選擇ANSYS主菜單Solution->LoadStepOpts->Spectrum->SinglePoint->FreqTable,輸入頻率反應(yīng)譜Freq1~Freq12為0.167,0.25,0.333,0.44444,0.5,0.667,1,1.25,1.667,2.222,10,100000,如圖

<36>選擇ANSYS主菜單Solution->LoadStepOpts->Spectrum->SinglePoint->SpectrValues,輸入對應(yīng)的反應(yīng)譜數(shù)值依次如下：

<37>最后選擇ANSYS主菜單Solution->LoadStepOpts->Spectrum->SinglePoint->ModeCombine,設(shè)定振型組合方式為SRSS法,如圖所示

<38>選擇ANSYS主菜單Solution->Solve->CurrentLS,計算反應(yīng)譜結(jié)果

<39>進(jìn)入ANSYS主菜單后處理模塊GeneralPostproc,在ANSYS窗口頂部菜單選擇File->Readinputfrom,選擇文件后綴名為*.mcom的文件

<40>進(jìn)入ANSYS主菜單GeneralPostproc->PlotResults->DeformedShapes,選擇繪制變形后形狀和結(jié)構(gòu)形狀,得到地震反應(yīng)譜分析的結(jié)構(gòu)變形如圖

<41>最后我們來進(jìn)行地震時程分析,進(jìn)行地震時程分析以前,首先要有一個地震時程記錄,本例子給定的地震時程記錄總長20秒,記錄點(diǎn)間隔0.02s,共有1001個記錄點(diǎn)。該地震記錄存放在RECORD.TXT文件中。<42>首先建立兩個變量,在ANSYS窗口頂部菜單選擇Parameters->Scalarparameters,在窗口中輸入NT=1001,即總共1001個記錄點(diǎn),DT=0.02,即記錄點(diǎn)間隔0.02s

<43>然后還需要定義一個數(shù)組來存放地震記錄,選擇ANSYS窗口頂部菜單Parameters->Arrayparameters->Define/Edit,在彈出的ArrayParameters窗口中點(diǎn)擊Add按鈕,在AddNewArrayParameter窗口中輸入數(shù)組的名稱為AC,行數(shù)為NT個,如圖所示。

<44>下面從數(shù)據(jù)文件中讀入地震時程記錄。選擇ANSYS窗口頂部菜單Parameters->Arrayparameters->Readfromfile,在彈出窗口中依次輸入將數(shù)據(jù)讀入的數(shù)組名稱AC,可以用Browse選擇目標(biāo)數(shù)據(jù)文件,最后要給出數(shù)據(jù)讀入的格式〔F8.3,可以參考Fortran相應(yīng)的文件輸入輸出要求。

<45>完成數(shù)據(jù)讀入后就可以進(jìn)行地震反應(yīng)分析。由于本次地震分析一共要進(jìn)行1001步,超過了ANSYS默認(rèn)的最大1000步限制,因此首先需要進(jìn)行調(diào)整。在ANSYS的命令輸入窗口,輸入以下命令：Finish/CONFIG,NRES,20000

<46>進(jìn)入ANSYS主菜單Solution->AnalysisType->NewAnalysis,指定分析類型為瞬態(tài)分析Transient,Solutionmethod可以選擇Full。

<47>進(jìn)入ANSYS主菜單Solution->AnalysisType->Sol'nControls,在SolutionControls窗口里面選擇Transient頁面,設(shè)定瑞雷阻尼的數(shù)值。由于阻尼的機(jī)理十分復(fù)雜,因此不同問題可能各不相同,這里僅介紹一種阻尼的取法。輸入質(zhì)量阻尼系數(shù)為2*DAMPRATIO*FREQ1*2*3.1415926,剛度阻尼系數(shù)為2*DAMPRATIO/<FREQ1*2*3.1415926>,FREQ1為結(jié)構(gòu)的第一階自振頻率。

<48>下面將用一組循環(huán)控制語句來定義作用在結(jié)構(gòu)上的加速度并進(jìn)行計算。讀者可以將下面這部分語句用任意編輯器〔比如windows的記事本輸入好以后復(fù)制ANSYS的輸入欄中。*DO,I,1,1001!對變量I循環(huán)1001ACEL,AC<I>,0,0!對結(jié)構(gòu)施加X方向加速度TIME,I*0.02!計算時間步長為0.02sOUTRES,ALL,ALL!輸出所有結(jié)果SOLVE!求解*ENDDO!循環(huán)結(jié)束

<49>下面介紹適用ANSYS的時程后處理器TimeHistPostpro來處理計算結(jié)果,在ANSYS主菜單上選擇TimeHistPostpro,這時會彈出時程變量窗口TimeHistoryVariables。點(diǎn)擊工具欄上第一個綠色加號按鈕,添加時程變量。在AddTime-HistoryVariable窗口中選擇要添加的變量為節(jié)點(diǎn)的X方向位移。

<50>選擇第20號節(jié)點(diǎn),回到時程變量窗口TimeHistoryVariables,點(diǎn)擊工具欄上第三個按鈕,即繪制出節(jié)點(diǎn)的位移時程曲線如圖*SET,NT,500!假如有500個地震波數(shù)據(jù)點(diǎn)

*SET,DT,0.02!每個地震波數(shù)據(jù)點(diǎn)之間的時間間隔是0.02*dim,AC1,,NT,1!定義X方向加速波的數(shù)組

!*dim,AC2,,NT,1!定義Y、Z方向加速波的數(shù)組如果有的話

!*dim,AC3,,NT,1*vread,AC1<1>,elcentro_EW,txt,,JIK,1,NT!讀入X方向地震波文件<elcentro_EW.txt>到數(shù)組

<F16.2>

!*vread,AC2<1>,elcentro_NS,txt,,JIK,1,NT!讀入Y、Z方向地震波文件、如果有的話

!<F16.2>

!*vread,AC3<1>,elcentro_EW,txt,,JIK,1,NT

!<F16.2>*DO,I,1,NT!加載求解

ACEL,AC1<I>,AC2<I>,AC3<I>

TIME,I*DT

SOLVE

!************讀入地震波數(shù)據(jù)************

*DIM,ELC,ARRAY,1,150,0,,,!定義數(shù)組

*CREATE,ansuitmp!讀入數(shù)據(jù)

*VREAD,ELC<1,1>,'ELC','txt','',150,,,,,,

<e9.3,e11.3>

*END

/INPUT,ansuitmp

!****************求解****************

FINI

allsel

/SOLU

!eqslv,pcg,1e-4

ANTYPE,4!指定分析類型為瞬態(tài)動力學(xué)分析

TRNOPT,full!瞬態(tài)動力學(xué)分析采用模態(tài)疊加法

!**************地震載荷**************

*DO,T,1,150,1!循環(huán)讀入地震數(shù)據(jù)并求解

TIME,0.001*T!設(shè)置時間步

KBC,0!指定載荷為遞增載荷

NSUB,1!設(shè)定子步數(shù)為1

DAMPRATIO=0.05

pi=3.1415!結(jié)構(gòu)的阻尼比

FREQ1=1.0721*2*pi!角頻率變圓頻率

freq2=1.0927*2*pi

!通過模態(tài)求解得出頻率

ALPHAD,2*DAMPRATIO*FREQ1*freq2/<freq1+freq2>!質(zhì)量阻尼

BETAD,2*DAMPRATIO/<FREQ1+freq2>!剛度阻尼

ACEL,ELC<1,T>/347,,!設(shè)定x,y方向加速度

ALLSEL!選中所有元素

SOLVE!求解

*ENDDO

SAVE!保存

FINI

建立acelx、acely兩個地震波文本〔txt,文本中第一列為時間,第二列為加速度,如果兩列之間為兩空格,則加速度為負(fù),一個空格則加速度為正,190行,19秒,間隔0.1記錄一個值麻煩你看看錯在哪,

*dim,tjx,array,190,2

*dim,tjy,array,190,2

*create,ansuitmp

*vread,tjx<1,1>,ACELX,txt,,jik,2,190

<E9.3,E11.3>

*end

/input,ansuitmp

*create,ansuitmp

*vread,tjy<1,1>,ACELY,txt,,jik,2,190

<e9.3,e11.3>

*end

/input,ansuitmp

!對于地震波的輸入,可以把荷載記錄做成文件,利用apdl的讀取功能讀入倒數(shù)據(jù)庫中。下面的例子是自己編的一個小文件。修改一下可以更簡潔。有用到的朋友自己作一下把。

fini

/config,nres,1000

*dim,aceX,TABLE,3000,1

*dim,aceY,TABLE,3000,1

*dim,aceZ,TABLE,3000,1

*creat,ff

*vread,aceX<1,1>,acex,txt,,1

<e16.6>

*vread,aceX<1,0>,ACETT,,,1

<e17.6>

ACEX<0,1>=1

*end

/input,ff

*creat,ff

*vread,aceY<1,1>,acey,txt,,1

<e16.6>

*vread,aceY<1,0>,ACETT,,,1

<e17.6>

ACEY<0,1>=1

*end

/input,ff

*creat,ff

*vread,aceZ<1,1>,acez,txt,,1

<e16.6>

*vread,aceZ<1,0>,ACETT,,,1

<e17.6>

ACEZ<0,1>=1

*end

/input,ff

！地震波時程記錄分成了3個文件,每個文件是一列。分別記錄x,y,z方向的加速度。acett是時間記錄。

這樣就可以把加速度記錄讀取倒ansys數(shù)據(jù)庫中作為數(shù)組。

也可以把加速度記錄做成一個文件,這樣程序就簡單多了。大家可以試看看修改一下。

下面是計算部分語句：

/SOLU

ANTYPE,trans

！求解其自己選了

TM_START=0.01

TM_END=15.00

TM_INCR=0.01

*DO,TM,TM_START,TM_END,TM_INCR

TIME,tm

alpha,

BETAD,

ACEL,acex<tm>,acey<tm>,acez<tm>

SOLVE

*ENDDO