OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
polynomial51.F File Reference
#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine polyun51 (c0, c1, c2, c3, c4, c5, gg, volume, dvol, vold, rho, masa, rhoa0, dd, mu, mup1, pold, pext, p, pm, q, rho0e, einta, viscmax, xl, ssp, qa, qb)
subroutine poly51 (c01, c11, c21, c31, c41, c51, gg1, v10, v1, v1old, mu1, mup1, eint1, pext, p1, pm1, p1i, rho1, rho10, mas1, ssp1, dvdp1, dpdv1, e_inf, flag)

Function/Subroutine Documentation

◆ poly51()

subroutine poly51 ( intent(in) c01,
intent(in) c11,
intent(in) c21,
intent(in) c31,
intent(in) c41,
intent(in) c51,
intent(in) gg1,
intent(in) v10,
intent(in) v1,
intent(in) v1old,
intent(inout) mu1,
intent(inout) mup1,
intent(inout) eint1,
intent(in) pext,
intent(inout) p1,
intent(in) pm1,
intent(in) p1i,
intent(inout) rho1,
intent(in) rho10,
intent(in) mas1,
intent(inout) ssp1,
intent(inout) dvdp1,
intent(inout) dpdv1,
intent(in) e_inf,
integer flag )

Definition at line 96 of file polynomial51.F.

101C-----------------------------------------------
102C I m p l i c i t T y p e s
103C-----------------------------------------------
104#include "implicit_f.inc"
105C-----------------------------------------------
106C I N P U T O U T P U T A r g u m e n t s
107C-----------------------------------------------
108 my_real,INTENT(IN) :: c01,c11,c21,c31,c41,c51,gg1,
109 . v10,v1,v1old,
110 . pext,pm1,p1i,
111 . rho10,mas1,e_inf
112 my_real,INTENT(INOUT) :: dvdp1,p1,ssp1,dpdv1,eint1,mu1,mup1,rho1
113 INTEGER FLAG
114C-----------------------------------------------
115C L o c a l V a r i a b l e s
116C-----------------------------------------------
117 my_real aa,bb,dvdp1i
118C------------------------
119 dvdp1i = dvdp1
120 IF (flag == 1) rho1 = mas1/v1
121 mup1 = rho1/rho10
122 mu1 = mup1 - one
123 aa = (c41 + c51*mu1)/v10
124 bb = half*(v1-v1old)
125 IF (flag == 1) eint1 = eint1 - (p1i+pext+pext)*bb
126 IF (flag == 1) THEN
127 p1 = ( c01 + c11*mu1 + max(mu1,zero)*(c21*mu1 + c31*mu1*mu1) + aa*eint1 ) / (one+aa*bb)
128 ELSE
129 p1 = ( c01 + c11*mu1 + max(mu1,zero)*(c21*mu1 + c31*mu1*mu1) + aa*eint1 )
130 ENDIF
131 p1 = max(p1,pm1)
132 IF (flag == 1) eint1 = eint1 - p1*bb
133 IF (flag == 1) eint1 = max(eint1, e_inf*v10)
134
135 dpdv1 = (-(c11 + (c21 + c21 + three*c31*mu1)*max(zero,mu1) + c51*eint1/v10 )*mup1*mup1 -(c41 + c51*mu1)*(p1+pext) ) / v10
136 ssp1 = (-dpdv1*v1/mup1 + two_third*gg1)/rho10
137 ssp1 = sqrt(abs(ssp1))
138
139 IF(abs(dpdv1)<em20)THEN
140 dvdp1 = zero
141 ELSE
142 dvdp1 = one/dpdv1
143 ENDIF
144
145 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
max
#define max(a, b)
Definition macros.h:21

◆ polyun51()

subroutine polyun51 ( intent(in) c0,
intent(in) c1,
intent(in) c2,
intent(in) c3,
intent(in) c4,
intent(in) c5,
intent(in) gg,
intent(in) volume,
intent(inout) dvol,
intent(in) vold,
intent(inout) rho,
intent(in) masa,
intent(in) rhoa0,
intent(in) dd,
intent(inout) mu,
mup1,
intent(in) pold,
intent(in) pext,
intent(inout) p,
intent(in) pm,
intent(inout) q,
intent(in) rho0e,
intent(inout) einta,
intent(inout) viscmax,
intent(in) xl,
intent(inout) ssp,
intent(in) qa,
intent(in) qb )

Definition at line 28 of file polynomial51.F.

35C-----------------------------------------------
36C I m p l i c i t T y p e s
37C-----------------------------------------------
38#include "implicit_f.inc"
39C-----------------------------------------------
40C I N P U T O U T P U T A r g u m e n t s
41C-----------------------------------------------
42 my_real,INTENT(IN) :: c0,c1,c2,c3,c5,c4,gg,volume,vold,masa ,rhoa0,dd, pold,pext,pm
43 my_real,INTENT(IN) :: rho0e,xl, qa,qb
44 my_real,INTENT(INOUT) :: einta
45 my_real,INTENT(INOUT) :: p, dvol, rho, mu, viscmax, q, ssp
46C-----------------------------------------------
47C L o c a l V a r i a b l e s
48C-----------------------------------------------
49 my_real :: aa,bb,dpdv,dpdmu,v0,qal,qbl,mup1
50C-----------------------------------------------
51C L o c a l V a r i a b l e s
52C-----------------------------------------------
53 v0 = masa / rhoa0
54 dvol = volume - vold
55 rho = masa/volume
56 mup1 = rho/rhoa0
57 mu = mup1 - one
58
59 dpdv = (-(c1 + (c2 + c2 + three*c3*mu)*max(zero,mu) + c5*rho0e )*mup1*mup1 -(c4 + c5*mu)*(pold+pext) ) / v0
60 dpdmu = -dpdv*volume/mup1
61 dpdmu = abs(dpdmu)
62 ssp = sqrt((dpdmu + two_third*gg)/rhoa0)
63 qal = qa*xl
64 qal = qal*qal
65 qbl = qb*xl
66 viscmax = rho*(qal*max(zero,dd) + qbl*ssp)
67 q = viscmax*max(zero,dd)
68
69 aa = (c4 + c5*mu)/v0
70 bb = half*(volume-vold)
71
72 !first torder
73 p = ( c0 + c1*mu+ max(mu,zero)*(c2*mu + c3*mu*mu) + aa*einta )
74 p = max(p,pm)
75 !second order
76 !P = ( PEXT+C0 + C1*MU+ MAX(MU,ZERO)*(C2*MU + C3*MU*MU) + AA*EINTA ) / (ONE+AA*BB)
77 !P = MAX(P,PM)-PEXT
78 !EINTA = EINTA - BB*(P+PEXT+Q) !2nd order integration (semi-implicit)
79
80
81
82
83 dpdv = (-(c1 + (c2 + c2 + three*c3*mu)*max(zero,mu) + c5*rho0e )*mup1*mup1 -(c4 + c5*mu)*(p+pext) ) / v0
84 dpdmu = -dpdv*volume/mup1
85 dpdmu = abs(dpdmu)
86 ssp = sqrt((dpdmu + two_third*gg)/rhoa0)
87
88 RETURN