Вывод переменной в файл txt

[Nео 0]

Здравствуйте!

Подскажите пожалуйста, как вывести переменные в файл. 

А конкретнее, есть текстовый файл с командами adpl, там переменные геометрии, ну например размеры трещины а и б.

Запускаю на счет. В результате есть значения напряжений и деформаций.

Нужно чтобы в конце счета в текстовый файл result.txt записывалось:

a=5 мм  b=15 мм Узел 1  Sig= N

где N - полученные напряжения в узле №1

После этого, нужно чтоб программа изменяла эти переменные 

а=10 мм  b = 20 мм

запускала на счет 

и по результатам счета выводила в текстовый файл,  в конец текстового файл (аppend)

a=10 мм  b=20 мм Узел 1  Sig= N1

 

Прошу пример листинга самой простой задачи - допустим 1 переменная... длина стержня... Дальше я сам разберусь под потребности своей задачи.

 

 

 

 

 

 

 

 

[Graf Kim 217]

Можно сделать по-умному. Но я так не люблю. Гораздо проще пользоваться COM'ом как print'ом в языках программирования.
 

/OUTPUT,result,txt !Перенаправляем вывод в текстовый файл. Файл будет перезаписан.

/COM, a = %a%, b = %b%, Sig = %nd_str% !Предполагается, что в переменных a, b и nd_str заранее записаны нужные значения.

/OUTPUT !Перенаправляем вывод обратно в окно вывода, или файл по-умолчанию.

 

Чтобы делать что-то многократно - пользуйтесь циклами.

*DO,i,1,10,1

    a=5*i

    b=10+5*i

    ...

    /OUTPUT,result,txt,,APPEND !Перенаправляем вывод в файл. Файл НЕ будет перезаписан.

    /GOPR

    /COM, a = %a%, b = %b%, Sig = %nd_str%

    /NOPR

    /OUTPUT

    ...

*ENDDO

    

[Борман 2 391]

*create,MODEL

	/prep7

	a=ARG1
	b=ARG2
	
	rectng,0,a,0,b
	
	et,1,42
	
	esize,a/10
	amesh,1
	
	mp,ex,1,10
	mp,nuxy,1,0.3
	
	d,1,all,0
	d,2,all,0
	f,3,fx,1
	
	/solu
	solve

	finish		
	/post1

	res=ux(4)	

	*cfopen,result.txt,,,append
	*vwrite,a,b,res
a=%f,	b=%f,	RES=%f
	*cfclose
	
	/prep7
	aclear,all
	adele,all,,,1

*end


*do,i,1,5
	*do,j,1,5
		*use,MODEL,i,j
	*enddo
*enddo

/sys,start notepad result.txt

Это все в одном файлике.

[Nео 0]

Спасибо!!!) именно это нужно было... 

Загвостка у меня была вот с этими  кодами

a=%f,	b=%f,	RES=%f

Это, насколько я понял, коды фортрана...  Не смог найти по ним информацию..

 

 

 

[Nео 0]

Разобрался с листингом, все вроде понятно и просто.

Появилось один вопрос:

В каком формате присваивать переменной res значения эквивалентных и прочих напряжений, усилий и т.д.

res=SX(1) или res=SEQV(1)

пишет ошибка, не хватает параметров. 

Вот в примере, перемещения (U) по Х,Y,Z, для узла 4- все правильно выводит:

res=ux(4)	

 

[Борман 2 391]

42 минуты назад, Nео сказал:

В каком формате присваивать переменной res значения эквивалентных и прочих напряжений, усилий и т.д.

 

Читайте справку к команде *GET, все можно выдернуть.

[Nео 0]

Спасибо. Сделал.  

*get,res,NODE,1,S,X

Из справочника программиста APDL:

"Например,  *GET,A,ELEM,5,CENT,X  возвращает 
х-координату центра тяжести элемента номер 5 и хранит результат как параметр A."

 

 

[Nео 0]

Здравствуйте! Упростил задачу - одна переменная и один цикл. И выдает ошибку... Что тут блин не так, все уже разжевано вроде... ((((

 

 

/BATCH  
/input,menust,tmp,''

WPSTYLE,,,,,,,,0

*create,MODEL1

/PREP7  

a=0.1

BLC4,-0.35,-0.1,a,0.1

ET,1,PLANE182   

KEYOPT,1,1,0
KEYOPT,1,3,3
KEYOPT,1,6,0

MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2.1e11 
MPDATA,PRXY,1,, 
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDE,EX,1   
MPDE,PRXY,1 
MPDATA,EX,1,,2.1E+011   
MPDATA,PRXY,1,,0.3  
sect,1,shell,,  
secdata, 0.010,1,0.0,3  
secoffset,MID   
seccontrol,,,, , , ,

FINISH  
/SOL
FINISH  
/PREP7  
/UI,MESH,OFF
MSHAPE,1,2D 
MSHKEY,0
!*  
CM,_Y,AREA  
ASEL, , , ,       1 
CM,_Y1,AREA 
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FINISH  
/SOL
FLST,2,1,4,ORDE,1   
FITEM,2,1   
!*  
/GO 
DL,P51X, ,ALL,  
FLST,2,1,4,ORDE,1   
FITEM,2,3   
/GO 
!*  
SFL,P51X,PRES,2e6,  
FLST,2,1,4,ORDE,1   
FITEM,2,3   
/GO 
!*  
SFL,P51X,PRES,-4000000, 
  

R,1,0.01,   
  
FINISH  
/SOL

STATUS,SOLU

SOLVE   
FINISH  
/POST1  

SET,FIRST   
  
 /EFACET,1   
PLNSOL, S,EQV, 0,1.0

*get,aa,NODE,1,S,eqv

res=aa

*cfopen,'c:\1\result.txt',,,
*vwrite,a,res
a=%f,res=%f
*cfclose
    
/prep7
aclear,all
adele,all,,,1

*end

*do,i,1,10
        *use,MODEL1,i

*enddo

 

[Борман 2 391]

39 минут назад, Nео сказал:

/prep7
aclear,all
adele,all,,,1

*end

*do,i,1,10
        *use,MODEL1,i

*enddo

 

 

...

*end

 

*do,i,1,10
        *use,MODEL1,i

        /prep7
        aclear,all
        adele,all,,,1

*enddo

Изменено 13 октября 2020 пользователем Борман

[Борман 2 391]

А не....

[Nео 0]

Здравствуйте! По команде *Get...

(Из справочника программиста APDL: "Например,  *GET,A,ELEM,5,CENT,X  возвращает 
х-координату центра тяжести элемента номер 5 и хранит результат как параметр A.")

 

Это все хорошо, но нужно указывать номер узла или элемента и его свойства (координаты, напряжения и т.д..).

А есть ли простой способ присвоить переменной максимальные значения напряжений(деформаций) из всех существующих?

Для примера: в Ансисе выбираешь Plot Results/Nodal Solu/X-Component of stress и выдает в данном случае напряжения по Х в модели. И на экране - Ансис уже определяет максимальные значения - указанны в верхнем левом углу переменные DMX, SMN, SMX. 

Можно ли вытащить командой *GET переменную SMX? Нужно сделать запрос, в каком элементе SMX?

 

[houston707 3]

После отрисовки *get,A,plnsol,0,MAX

[Nео 0]

Спасибо все верно! Я так и раньше пробовал, но почему-то в этом случае не записываются переменные в файл - вместо цифр- звездочки... 

Причем в Scalar Parameters по итогу счета по листингу все верно - данные последнего цикла:

A=1.1.; B=3522965.75; I=10; RES= 4466102.5

А в файл не выводит..... только звездочки вместо цифр....

 

Вот листинг:

 

/BATCH  
/input,menust,tmp,''

WPSTYLE,,,,,,,,0

a=0.1

 

*create,MODEL1

/PREP7  

BLC4,-0.35,-0.1,a,0.1

ET,1,PLANE182   

KEYOPT,1,1,0
KEYOPT,1,3,3
KEYOPT,1,6,0

MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2.1e11 
MPDATA,PRXY,1,, 
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDE,EX,1   
MPDE,PRXY,1 
MPDATA,EX,1,,2.1E+011   
MPDATA,PRXY,1,,0.3  
sect,1,shell,,  
secdata, 0.010,1,0.0,3  
secoffset,MID   
seccontrol,,,, , , ,

FINISH  
/SOL
FINISH  
/PREP7  
/UI,MESH,OFF
MSHAPE,1,2D 
MSHKEY,0
!*  
CM,_Y,AREA  
ASEL, , , ,       1 
CM,_Y1,AREA 
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FINISH  
/SOL
FLST,2,1,4,ORDE,1   
FITEM,2,1   
!*  
/GO 
DL,P51X, ,ALL,  
FLST,2,1,4,ORDE,1   
FITEM,2,3   
/GO 
SFL,P51X,PRES,2e6,  
FLST,2,1,4,ORDE,1   
FITEM,2,3   
/GO 
SFL,P51X,PRES,-4000000, 
 R,1,0.01,   
 FINISH  
/SOL

STATUS,SOLU

SOLVE   
FINISH  

/POST1  
SET,NEXT
 /EFACET,1   

 

plnsol, S,EQV, 0,1.0

*get,res,plnsol,0,MAX
*get,b,plnsol,0,MIN

 

*cfopen,'c:\111\result2.txt',,,append
*vwrite,a,res,b

a=%f,    res=%f,    b=%f

*cfclose

    
*end

*do,i,1,10
        *use,MODEL1,i
        a=a+0.1
        /prep7
        aclear,all
        adele,all,,,1

*enddo

/sys,start notepad C:\111\result2.txt

[Борман 2 391]

3 минуты назад, Nео сказал:

А в файл не выводит..... только звездочки вместо цифр....

 

Это вроде когда места под цифры не хватает. Надо формат сменить как то.

[houston707 3]

В строчке a=%f,    res=%f,    b=%f поменять %f на %e

[Nео 0]

Спасибо, получилось) 

Вообще, сложновато конечно искать ошибки... 

Скачал и ознакамливаюсь:

1.Автоматизированное моделирование и расчёт конструкций в ANSYS: од-
номерные  модели:  учеб.  пособие  /  О.А. Саченков,  А.А. Саченков, 
П.В. Большаков, О.В. Герасимов. – Казань: Казан. ун-т, 2019. – 140 с. 

2. Справочник программиста APDL  для v.10

Такое чувство что очень сжатое там изложение...

[Nео 0]

Здравствуйте!

Все получилось...

Но делая параметрическую модель, с изменением параметра в цикле столкнулся с тем что вторая итерация не хочет считать. Такое чувство что переменные неправильно удаляются... 

Такое чувство, что я неправильно удаляю геометрию в конце 1 итерации цикла... 

Нужно чтоб геометрия удалялась и параметры удалялись, кроме параметра L. А у меня. на 2ой итерации цикла какие то неправильные параметры для построения геометрии и ошибка в построении.

(Причем тестовая задача - прямоугольник - считалась нормально на всех итерациях).

 

/prep7

L=0.01

* create,MODEL1

              a=(2*L)/3 – параметры зависящие от L

              ПРОЧЕЕ

              FINISH  /POST1 

              plnsol, S,EQV, 0,1.0

             *get,Sigma_EKV_Max,plnsol,0,MAX

             *cfopen,'c:\111\result2.txt',,,append

            *vwrite,L,L1,…………….

            L=%f,………………

            *cfclose

*end

 

 

*do,ii,1,3

                               *use,MODEL1

                                L=L+0.01

                               /prep7

                               aclear,all

                               adele,all,,,1

*enddo

 

 

 

Предположение:

т.е. L =0.02 в конце 1го цикла, а потом команды aclear и adele делают L=0? Как сохранить L?

 

[houston707 3]

Сейчас параметры не удаляются, для их удаления нужна команда /clear, но перед этим нужно сохранить параметр командой parsav и после очищения загрузить командой parres. Можно сохранить параметры перед вызовом model1, потом в model1 очистить и загрузить, потом в теле цикла очистить, загрузить и увеличить l,  а потом новая итерация.

[Graf Kim 217]

@Nео может быть дадите "прочее"? Без него всё отработало, L меняется.

[Nео 0]

Здравствуйте!

Да вроде наладил цикл, переменная L меняется, в файл записывается...

Теперь другая проблема - результаты пишет в файл одинаковые.... Когда я без цикла считаю при разных L, результаты, например  макс.напряжений (по Мизесу) - разные.

А когда запускаю скрипт с циклом - то в файле результаты одинаковые. Там в конце цикла выдается  табличка - "Слишком много ошибок, выйти или нет?"

Если нажать "да", то на 1м цикле счет прерывается, если "нет", то дальше считает, при другой L.

Но результаты в файле - одинаковые, даже если несколько циклов прокручиывает.

Как вариант, забыл где-то поставить "prep7"?

Вот весь листинг. Это как пример прочностной задачи, для дальнейшего мне важно прокручивать задачу именно в цикле (возможно и в нескоьких), с множеством переменных. А тут в простой задаче не могу разобраться.

 

/prep7
    L=0.01

*create,MODEL1

    h=0.1
    W=0.2
    T=0.2
    a=(2*L)/3
    P=20e6
    E=2.1e11
    v=0.3
    Ro_st=7850
    W1=W/2
    T1=T/2
    h1=h/2
    L1=L/2
    a1=a/2
    dl1=L1/6
    dl2=L1/6
    n=1
ET,1,PLANE183   
ET,2,SOLID186   

MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2.06e11 
MPDATA,PRXY,1,,0.3   
MPDATA,DENS,1,,Ro_st

TB,PLAS,1,1,7,MISO
TBTEMP,0
TBPT,,0.0,5E+6
TBPT,,4.854E-005,10E+6
TBPT,,1.456E-004,30E+6
TBPT,,2.427E-004,50E+6
TBPT,,3.398E-004,70E+6
TBPT,,4.369E-004,90E+6
TBPT,,5.34E-004,110E+6
TBPT,,6.311E-004,130E+6
TBPT,,7.282E-004,150E+6
TBPT,,8.252E-004,170E+6
TBPT,,9.223E-004,190E+6
TBPT,,1.019E-003,210E+6
TBPT,,1.604E-003,230E+6
TBPT,,3.214E-003,250E+6
TBPT,,5.167E-003,270E+6
TBPT,,7.35E-003,290E+6
TBPT,,9.712E-003,310E+6
TBPT,,0.012,330E+6
TBPT,,0.015,350E+6
TBPT,,0.018,370E+6
TBPT,,0.020,390E+6
TBPT,,0.022,400E+6

FLST,2,2,8  
FITEM,2,0,0,0   
FITEM,2,(L1/2),0,0  
CIRCLE,P51X, , , ,90,n,   

FLST,2,1,4,ORDE,1   
FITEM,2,1   
LSSCALE,P51X, , ,((L1/2)/(a1/2)),1,1, ,1,1

FLST,2,2,8  
FITEM,2,0,0,0   
FITEM,2,((L1/2)-(dl1/2)),0,0  
CIRCLE,P51X, , , ,90,n,   

FLST,2,1,4,ORDE,1   
FITEM,2,2   
LSSCALE,P51X, , ,((L1/2)-(dl1/2))/((a1/2)-(dl1/2)),1,1, ,1,1

FLST,2,2,8  
FITEM,2,0,0,0   
FITEM,2,((L1/2)+(dl1/2)),0,0  
CIRCLE,P51X, , , ,90,n,   

FLST,2,1,4,ORDE,1   
FITEM,2,3   
LSSCALE,P51X, , ,((L1/2)+(dl1/2))/((a1/2)+(dl1/2)),1,1, ,1,1

K,7,0,0,0,  
K,8,0,h1,0,
K,9,W1,h1,0,
K,10,W1,0,0,
K,11,0,0,T1,
K,12,0,h1,T1,
K,13,W1,h1,T1,
K,14,W1,0,T1,

LSTR,       7,       3  
LSTR,       3,       1  
LSTR,       1,       5  
LSTR,       5,      10  
LSTR,      10,       9  
LSTR,       9,       8  
LSTR,       8,       6  
LSTR,       6,       2  
LSTR,       2,       4  
LSTR,       4,       7  
LSTR,      10,      14  
LSTR,      14,      13  
LSTR,      13,       9  
LSTR,      14,      11  
LSTR,      11,      12  
LSTR,      12,      13  
LSTR,      12,       8  
LSTR,      11,       7  

FLST,3,13,4,ORDE,2  
FITEM,3,1   
FITEM,3,-13 
LGEN,2,P51X, , , , ,(dl1/2), ,1

LSTR,      19,       5  
LSTR,      15,       1  
LSTR,      17,       3  
LSTR,      20,       6  
LSTR,      16,       2  
LSTR,      18,       4  
LSTR,      21,       7  

FLST,2,4,4  
FITEM,2,41  
FITEM,2,25  
FITEM,2,4   
FITEM,2,37  
AL,P51X 

FLST,2,4,4  
FITEM,2,40  
FITEM,2,34  
FITEM,2,13  
FITEM,2,41  
AL,P51X 

FLST,2,4,4  
FITEM,2,40  
FITEM,2,23  
FITEM,2,2   
FITEM,2,37  
AL,P51X 

FLST,2,3,4  
FITEM,2,23  
FITEM,2,34  
FITEM,2,25  
AL,P51X 

FLST,2,3,4  
FITEM,2,2   
FITEM,2,13  
FITEM,2,4   
AL,P51X 

FLST,2,5,5,ORDE,2   
FITEM,2,1   
FITEM,2,-5  
VA,P51X 

FLST,2,4,4  
FITEM,2,12  
FITEM,2,2   
FITEM,2,1   
FITEM,2,5   
AL,P51X 

FLST,2,4,4  
FITEM,2,33  
FITEM,2,26  
FITEM,2,23  
FITEM,2,22  
AL,P51X 

FLST,2,4,4  
FITEM,2,40  
FITEM,2,39  
FITEM,2,33  
FITEM,2,12  
AL,P51X 

FLST,2,4,4  
FITEM,2,26  
FITEM,2,5   
FITEM,2,37  
FITEM,2,36  
AL,P51X 

FLST,2,4,4  
FITEM,2,39  
FITEM,2,36  
FITEM,2,22  
FITEM,2,1   
AL,P51X 

FLST,2,6,5,ORDE,3   
FITEM,2,3   
FITEM,2,6   
FITEM,2,-10 
VA,P51X 

FLST,2,4,4  
FITEM,2,11  
FITEM,2,6   
FITEM,2,1   
FITEM,2,3   
AL,P51X 

FLST,2,4,4  
FITEM,2,27  
FITEM,2,32  
FITEM,2,22  
FITEM,2,24  
AL,P51X 

FLST,2,4,4  
FITEM,2,35  
FITEM,2,3   
FITEM,2,24  
FITEM,2,38  
AL,P51X 

FLST,2,4,4  
FITEM,2,38  
FITEM,2,39  
FITEM,2,11  
FITEM,2,32  
AL,P51X 

FLST,2,4,4  
FITEM,2,6   
FITEM,2,36  
FITEM,2,35  
FITEM,2,27  
AL,P51X 

FLST,2,6,5,ORDE,2   
FITEM,2,10  
FITEM,2,-15 
VA,P51X 

FLST,2,4,4,ORDE,4   
FITEM,2,14  
FITEM,2,16  
FITEM,2,20  
FITEM,2,-21 
LDELE,P51X, , ,1

LSTR,       9,      23  
LSTR,      23,      13  
LSTR,      10,      22  
LSTR,      22,      14  
LSTR,       8,      24  
LSTR,      24,      12  
LSTR,      11,      21  

FLST,2,10,4 
FITEM,2,30  
FITEM,2,29  
FITEM,2,28  
FITEM,2,31  
FITEM,2,32  
FITEM,2,33  
FITEM,2,34  
FITEM,2,25  
FITEM,2,26  
FITEM,2,27  
AL,P51X 

FLST,2,4,4  
FITEM,2,19  
FITEM,2,15  
FITEM,2,17  
FITEM,2,18  
AL,P51X 

FLST,2,4,4  
FITEM,2,15  
FITEM,2,16  
FITEM,2,21  
FITEM,2,29  
AL,P51X 

FLST,2,7,4  
FITEM,2,18  
FITEM,2,43  
FITEM,2,44  
FITEM,2,31  
FITEM,2,32  
FITEM,2,33  
FITEM,2,34  
AL,P51X 
  
FLST,2,4,4  
FITEM,2,19  
FITEM,2,16  
FITEM,2,43  
FITEM,2,30  
AL,P51X 

FLST,2,7,4  
FITEM,2,17  
FITEM,2,21  
FITEM,2,44  
FITEM,2,28  
FITEM,2,27  
FITEM,2,26  
FITEM,2,25  
AL,P51X 

FLST,2,4,4  
FITEM,2,28  
FITEM,2,35  
FITEM,2,7   
FITEM,2,20  
AL,P51X 

FLST,2,4,4  
FITEM,2,20  
FITEM,2,8   
FITEM,2,14  
FITEM,2,29  
AL,P51X 

FLST,2,4,4  
FITEM,2,14  
FITEM,2,9   
FITEM,2,30  
FITEM,2,42  
AL,P51X 

FLST,2,4,4  
FITEM,2,42  
FITEM,2,10  
FITEM,2,31  
FITEM,2,38  
AL,P51X 

!FLST,2,4,4  
!FITEM,2,38  
!FITEM,2,35  
!FITEM,2,3   
!FITEM,2,22  
!AL,P51X 

FLST,2,5,4  
FITEM,2,31  
FITEM,2,30  
FITEM,2,28  
FITEM,2,29  
FITEM,2,24  
AL,P51X 

FLST,2,5,4  
FITEM,2,7   
FITEM,2,8   
FITEM,2,9   
FITEM,2,10  
FITEM,2,3   
AL,P51X 

FLST,2,7,5,ORDE,3   
FITEM,2,13  
FITEM,2,22  
FITEM,2,-27 
VA,P51X 

FLST,2,6,5,ORDE,2   
FITEM,2,16  
FITEM,2,-21 
VA,P51X 

FLST,2,5,6,ORDE,2   
FITEM,2,1   
FITEM,2,-5  
VGLUE,P51X  

SMRT,6  
SMRT,1  
MSHAPE,0,2D 
MSHKEY,0

FLST,5,4,5,ORDE,4   
FITEM,5,5   
FITEM,5,-6  
FITEM,5,11  
FITEM,5,27  
CM,_Y,AREA  
ASEL, , , ,P51X 
CM,_Y1,AREA 
CHKMSH,'AREA'   
CMSEL,S,_Y  
  
AMESH,_Y1   
 
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  

FLST,5,2,5,ORDE,2   
FITEM,5,6   
FITEM,5,11  
CM,_Y,AREA  
ASEL, , , ,P51X 
CM,_Y1,AREA 
CMSEL,S,_Y  
CMDELE,_Y   
 
AREFINE,_Y1, , ,5,0,1,1 
CMDELE,_Y1  
 
FLST,5,1,5,ORDE,1   
FITEM,5,5   
CM,_Y,AREA  
ASEL, , , ,P51X 
CM,_Y1,AREA 
CMSEL,S,_Y  
CMDELE,_Y   
 
AREFINE,_Y1, , ,1,0,1,1 
CMDELE,_Y1  

FLST,5,1,5,ORDE,1   
FITEM,5,27  
CM,_Y,AREA  
ASEL, , , ,P51X 
CM,_Y1,AREA 
CMSEL,S,_Y  
CMDELE,_Y   

AREFINE,_Y1, , ,1,0,1,1 
CMDELE,_Y1  

/UI,MESH,OFF

GPLOT   

FLST,5,4,6,ORDE,2   
FITEM,5,1   
FITEM,5,-4  
CM,_Y,VOLU  
VSEL, , , ,P51X 
CM,_Y1,VOLU 
CHKMSH,'VOLU'   
CMSEL,S,_Y  
VSWEEP,_Y1  
  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
MSHAPE,1,3D 
CM,_Y,VOLU  
VSEL, , , ,       6 
CM,_Y1,VOLU 
CHKMSH,'VOLU'   
CMSEL,S,_Y  
VMESH,_Y1   
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  

FINISH  
/SOL

FLST,2,1,5,ORDE,1   
FITEM,2,17  
/GO 
 
SFA,P51X,1,PRES,P 

FLST,2,2,5,ORDE,2   
FITEM,2,11  
FITEM,2,27  
FLST,2,2,5,ORDE,2   
FITEM,2,11  
FITEM,2,27  
DA,P51X,SYMM

FLST,2,5,5,ORDE,5   
FITEM,2,2   
FITEM,2,8   
FITEM,2,14  
FITEM,2,19  
FITEM,2,25  
DA,P51X,SYMM

ACEL,0,0,-9.81, 

SOLVE   
FINISH  

/POST1  
SET,FIRST

plnsol, S,EQV, 0,1.0
*get,Sigma_EKV_Max,plnsol,0,MAX

plnsol, S,INT, 0,1.0
*get,Sigma_INT_max,plnsol,0,MAX

PLNSOL, EPTO,EQV, 0,1.0 
*get,DEF_Total_EKV_Max,plnsol,0,MAX

PLNSOL, EPTO,INT, 0,1.0 
*get,DEF_Total_INT_Max,plnsol,0,MAX

PLNSOL, U,X, 0,1.0  
*get,Perem_X_Max,plnsol,0,MAX

PLNSOL, U,Y, 0,1.0  
*get,Perem_Y_Max,plnsol,0,MAX

*cfopen,'c:\111\result2.txt',,,append
*vwrite,L,L1,Sigma_EKV_Max,Sigma_INT_max,DEF_Total_EKV_Max,DEF_Total_INT_Max,Perem_X_Max,Perem_Y_Max

L=%f,    L1=%f,    Sigma_EKV_Max=%e,    Sigma_INT_max=%e,    DEF_Total_EKV_Max=%e,    DEF_Total_INT_Max=%e,    Perem_X_Max=%e,    Perem_Y_Max=%e

*cfclose

*end

 

!_____________________________________________________________________________

*do,ii,1,7
        *use,MODEL1
                /clear,start
        L=L+0.01

*enddo